numo-gsl-alt 0.1.3

data/ext/numo/gsl/histogram/tmpl/c_DFloat_f_2d_get.c

@@ -0,0 +1,54 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   n, i;
+    char    *p1, *p2, *p3;
+    ssize_t  s1, s2, s3;
+    ssize_t  x1, x2;
+    size_t   y1, y2;
+    double   y3;
+    <%=struct%> *w = (<%=struct%>*)(lp->opt_ptr);
+
+    INIT_COUNTER(lp, n);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+    INIT_PTR(lp, 2, p3, s3);
+
+    for (i=0; i<n; i++) {
+        GET_DATA_STRIDE(p1,s1,ssize_t,x1);
+        GET_DATA_STRIDE(p2,s2,ssize_t,x2);
+        if (x1 < 0) {
+            y1 = x1 + w->nx;
+        } else {
+            y1 = x1;
+        }
+        if (x2 < 0) {
+            y2 = x2 + w->ny;
+        } else {
+            y2 = x2;
+        }
+        y3 = <%=func_name%>(<%=get(:postpose) ? "y1,y2,w" : "w,y1,y2"%>);
+        SET_DATA_STRIDE(p3,s3,double,y3);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args[1][1]%>, <%=args[2][1]%>)
+  @param  [Int32/64]  <%=args[1][1]%>
+  @param  [Int32/64]  <%=args[2][1]%>
+  @return [<%=class_name%>]  self
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(2)%>(VALUE self, VALUE v1, VALUE v2)
+{
+    <%=struct%> *w;
+    ndfunc_arg_in_t ain[2] = {{cSSZ,0},{cSSZ,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 2,1, ain,aout};
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    return na_ndloop3(&ndf, w, 2, v1, v2);
+}

data/ext/numo/gsl/histogram/tmpl/c_DFloat_f_field.c

@@ -0,0 +1,35 @@
+<% nd = (/histogram2d/ =~ struct) ? 2 : 1 %>
+/*
+  @overload <%=name%>()
+  @return [DFloat]    narray of <%=name%> field in <%=struct%>.
+
+  returns <%=name%> field in <%=struct%> struct.
+  <%= description %>
+*/
+static VALUE
+<%=c_func(0)%>(VALUE self)
+{
+    <%=struct%> *w;
+    double *d;
+    VALUE v;
+    size_t n, i, shape[<%=nd%>];
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    <%if /([xy]?)range/ =~ name %>
+    shape[0] = n = w->n<%=$1%>+1;
+    <%elsif nd==1%>
+    shape[0] = n = w->n;
+    <%else%>
+    n = w->nx * w->ny;
+    shape[0] = w->nx;
+    shape[1] = w->ny;
+    <%end%>
+
+    v = rb_narray_new(cDF, <%=nd%>, shape);
+    d = (double*)na_get_pointer_for_write(v);
+    for (i=0; i<n; i++) {
+        d[i] = w-><%=name%>[i];
+    }
+    return v;
+}

data/ext/numo/gsl/histogram/tmpl/c_DFloat_f_get.c

@@ -0,0 +1,47 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   n, i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    ssize_t  x1;
+    size_t   y1;
+    double   y2;
+    <%=struct%> *w = (<%=struct%>*)(lp->opt_ptr);
+
+    INIT_COUNTER(lp, n);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    for (i=0; i<n; i++) {
+        GET_DATA_STRIDE(p1,s1,ssize_t,x1);
+        if (x1 < 0) {
+            y1 = x1 + w->n;
+        } else {
+            y1 = x1;
+        }
+        y2 = <%=func_name%>(<%=get(:postpose) ? "y1,w" : "w,y1"%>);
+        SET_DATA_STRIDE(p2,s2,double,y2);
+    }
+}
+#undef FUNC
+
+/*
+  @overload <%=name%>(<%=args[1][1]%>)
+  @param  [Int32/64]  <%=args[1][1]%>
+  @return [<%=class_name%>]  self
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE self, VALUE v1)
+{
+    <%=struct%> *w;
+    ndfunc_arg_in_t ain[1] = {{cSSZ,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    return na_ndloop3(&ndf, w, 1, v1);
+}

data/ext/numo/gsl/histogram/tmpl/c_DFloat_x2_f_get_range.c

@@ -0,0 +1,56 @@
+<%
+n = case name
+    when /_xrange/; "nx"
+    when /_yrange/; "ny"
+    when /_range/;  "n"
+    end
+%>
+
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   n, i;
+    char    *p1, *p2, *p3;
+    ssize_t  s1, s2, s3;
+    ssize_t  x1;
+    size_t   y1;
+    double   y2, y3;
+    <%=struct%> *w = (<%=struct%>*)(lp->opt_ptr);
+
+    INIT_COUNTER(lp, n);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+    INIT_PTR(lp, 2, p3, s3);
+
+    for (i=0; i<n; i++) {
+        GET_DATA_STRIDE(p1,s1,ssize_t,x1);
+        if (x1 < 0) {
+            y1 = x1 + w-><%=n%>;
+        } else {
+            y1 = x1;
+        }
+        <%=func_name%>(<%=get(:postpose) ? "y1,&y2,&y3,w" : "w,y1,&y2,&y3"%>);
+        SET_DATA_STRIDE(p2,s2,double,y2);
+        SET_DATA_STRIDE(p3,s3,double,y3);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args[1][1]%>)
+  @param  [Int32/64]  <%=args[1][1]%>  Int32/64 NArray
+  @return [Array]     array of [[DFloat] <%=args[2][1]%>, [DFloat] <%=args[3][1]%>]
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE self, VALUE v1)
+{
+    <%=struct%> *w;
+    ndfunc_arg_in_t ain[1] = {{cSSZ,0}};
+    ndfunc_arg_out_t aout[2] = {{cDF,0},{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,2, ain,aout};
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    return na_ndloop3(&ndf, w, 1, v1);
+}

data/ext/numo/gsl/histogram/tmpl/c_bool_f_other.c

@@ -0,0 +1,19 @@
+/*
+  @overload <%=name%>(<%=args[1][1]%>)
+  @param  [<%=class_name%>] <%=args[1][1]%>
+  @return [Bool]
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE self, VALUE v1)
+{
+    int stat;
+    <%=struct%> *w, *w1;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+    TypedData_Get_Struct(v1, <%=struct%>, &<%=data_type_var%>, w1);
+
+    stat = <%=func_name%>(w, w1);
+    return (stat) ? Qtrue: Qfalse;
+}

data/ext/numo/gsl/histogram/tmpl/c_bool_f_pdf_init.c

@@ -0,0 +1,24 @@
+<%
+h_struct = struct.sub(/_pdf/,"")
+h_data_type_var = data_type_var.sub(/_pdf/,"")
+%>
+/*
+  @overload <%=name%>(<%=args[1][1]%>)
+  @param  [<%=class_name%>] <%=args[1][1]%>
+  @return [Bool]
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE self, VALUE v1)
+{
+    int stat;
+    <%=struct%> *w;
+    <%=h_struct%> *w1;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+    TypedData_Get_Struct(v1, <%=h_struct%>, &<%=h_data_type_var%>, w1);
+
+    stat = <%=func_name%>(w, w1);
+    return (stat) ? Qtrue: Qfalse;
+}

data/ext/numo/gsl/histogram/tmpl/c_self_f_2d_set_ranges.c

@@ -0,0 +1,26 @@
+/*
+  @overload <%=name%>(<%=args[1][1]%>,<%=args[3][1]%>)
+  @param  [DFloat]   <%=args[1][1]%>
+  @param  [DFloat]   <%=args[3][1]%>
+  @return [<%=class_name%>]  self
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(2)%>(VALUE self, VALUE v1, VALUE v2)
+{
+    <%=struct%> *w;
+    double *p1, *p2;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    v1 = cast_1d_contiguous(v1, cDF);
+    p1 = (double*)na_get_pointer_for_read(v1);
+    v2 = cast_1d_contiguous(v2, cDF);
+    p2 = (double*)na_get_pointer_for_read(v2);
+
+    <%=func_name%>(w, p1, RNARRAY_SIZE(v1), p2, RNARRAY_SIZE(v2));
+    RB_GC_GUARD(v1);
+    RB_GC_GUARD(v2);
+    return self;
+}

data/ext/numo/gsl/histogram/tmpl/c_self_f_set_ranges.c

@@ -0,0 +1,21 @@
+/*
+  @overload <%=name%>(<%=args[1][1]%>)
+  @param  [DFloat]   <%=args[1][1]%>
+  @return [<%=class_name%>]  self
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE self, VALUE v1)
+{
+    <%=struct%> *w;
+    double *p1;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    v1 = cast_1d_contiguous(v1, cDF);
+    p1 = (double*)na_get_pointer_for_read(v1);
+    <%=func_name%>(w, p1, RNARRAY_SIZE(v1));
+    RB_GC_GUARD(v1);
+    return self;
+}

data/ext/numo/gsl/histogram/tmpl/c_sizet_f_field.c

@@ -0,0 +1,15 @@
+/*
+  @overload <%=name%>
+  @return [Integer]
+
+  returns <%=name%> field in <%=struct%> struct.
+*/
+static VALUE
+<%=c_func(0)%>(VALUE self)
+{
+    <%=struct%> *w;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    return SIZET2NUM(w-><%=name%>);
+}

data/ext/numo/gsl/interp/depend

@@ -0,0 +1,17 @@
+COGEN=ruby erbpp_interp.rb
+GENSRC=gsl_interp.c
+GENDEPS=erbpp_interp.rb ../gen/*.rb ../gen/tmpl/*.c tmpl/*.c
+
+src : $(GENSRC)
+
+$(GENSRC) : $(GENDEPS)
+	$(COGEN) > $@
+
+doc : $(GENSRC)
+	yard doc $(GENSRC)
+
+clean: cleansrc cleandoc
+cleansrc:
+	-$(Q)$(RM) $(GENSRC)
+cleandoc:
+	-$(Q)$(RM_RF) doc .yardoc

data/ext/numo/gsl/interp/erbpp_interp.rb

@@ -0,0 +1,60 @@
+require_relative "../gen/erbpp_gsl"
+require_relative "parse_interp"
+
+gsl_list = ErbppGsl.read_func_pattern(
+ [/^gsl_interp_accel_(\w+)$/,    interp_accel_list=[]],
+ [/^gsl_interp2d_(\w+)$/,    interp2d_list=[]],
+ [/^gsl_interp_(\w+)$/,    interp_list=[]],
+ [/^gsl_spline2d_(\w+)$/,    spline2d_list=[]],
+ [/^gsl_spline_(\w+)$/,    spline_list=[]],
+)
+
+desc = {}
+gsl_list.each do |h|
+  desc[h[:func_name]] = h[:desc]
+end
+gsl_list.each do |h|
+  if /^gsl_spline/ =~ h[:func_name] && h[:desc].empty?
+    h[:desc] = desc[h[:func_name].sub(/_spline/,"_interp")]
+  end
+end
+
+DefLib.new do
+  set line_number: true
+  set erb_dir: %w[tmpl ../gen/tmpl]
+  set erb_suffix: ".c"
+  set ns_var: "mG"
+
+  ErbPP.new(self,"cast_1d_contiguous")
+  ErbPP.new(self,"cast_2d_contiguous")
+
+  name = "Interp"
+  set file_name: "gsl_#{name}.c"
+  set include_files: %w[gsl/gsl_interp.h gsl/gsl_interp2d.h
+                        gsl/gsl_spline.h gsl/gsl_spline2d.h]
+  set lib_name: name.downcase
+
+  [
+   #["Interp",      "interp",       interp_list],
+   ["InterpAccel", "interp_accel", interp_accel_list],
+   #["Interp2D",    "interp2d",     interp2d_list],
+   ["Spline",      "spline",       spline_list],
+   ["Spline2D",    "spline2d",     spline2d_list],
+  ].each do |name,base,list|
+
+    DefInterp.new(self) do
+      set name: base
+      set class_name: name
+      set class_var: "c"+name
+      set full_class_name: "Numo::GSL::"+name
+      set struct: "gsl_"+base
+
+      undef_alloc_func
+      init = nil
+      list.each do |h|
+        check_func(h)
+      end
+    end
+
+  end
+end.run

data/ext/numo/gsl/interp/extconf.rb

@@ -0,0 +1,9 @@
+require_relative '../extconf_gsl.rb'
+
+# source file to compile
+srcs = %w(
+ gsl_interp
+)
+$objs = srcs.collect{|i| i+".o"}
+
+create_makefile('numo/gsl/interp')