numo-gsl 0.1.0

data/ext/numo/gsl/gen/tmpl/c_other_f_void.c

@@ -0,0 +1,19 @@
+/*
+  @overload <%=name%>()
+  @return [<%=full_class_name%>]
+
+ <%= desc %>
+ */
+static VALUE
+<%=c_func(0)%>(VALUE self)
+{
+    <%=struct%> *x, *w;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, x);
+
+    w = <%=func_name%>(x);
+    if (!w) {
+        rb_raise(rb_eNoMemError,"fail to allocate struct");
+    }
+    return TypedData_Wrap_Struct(<%=class_var%>, &<%=data_type_var%>, (void*)w);
+}

data/ext/numo/gsl/gen/tmpl/c_self_f_DFloat.c

@@ -0,0 +1,46 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1;
+    ssize_t  s1;
+    size_t  *idx1;
+    double   x;
+    <%=struct%> *w = (<%=struct%>*)(lp->opt_ptr);
+    <% c_args = get(:postpose) ? "x,w" : "w,x" %>
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR_IDX(lp, 0, p1, s1, idx1);
+
+    if (idx1) {
+        for (; i--;) {
+            GET_DATA_INDEX(p1,idx1,double,x);
+            <%=func_name%>(<%=c_args%>);
+        }
+    } else {
+        for (; i--;) {
+            GET_DATA_STRIDE(p1,s1,double,x);
+            <%=func_name%>(<%=c_args%>);
+        }
+    }
+}
+
+/*
+  @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;
+    ndfunc_arg_in_t ain[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, FULL_LOOP, 1,0, ain,0};
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    na_ndloop3(&ndf, w, 1, v1);
+    return self;
+}

data/ext/numo/gsl/gen/tmpl/c_self_f_DFloat_x2.c

@@ -0,0 +1,41 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    double   x, y;
+    <%=struct%> *w = (<%=struct%>*)(lp->opt_ptr);
+    <% c_args = get(:postpose) ? "x,y,w" : "w,x,y" %>
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,double,x);
+        GET_DATA_STRIDE(p2,s2,double,y);
+        <%=func_name%>(<%=c_args%>);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args[1..2].map{|a| a[1]}.join(",")%>)
+  @param  [DFloat]   <%=args[1][1]%>
+  @param  [DFloat]   <%=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] = {{cDF,0},{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP, 2,0, ain,0};
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    na_ndloop3(&ndf, w, 2, v1, v2);
+    return self;
+}

data/ext/numo/gsl/gen/tmpl/c_self_f_DFloat_x3.c

@@ -0,0 +1,44 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2, *p3;
+    ssize_t  s1, s2, s3;
+    double   x1, x2, x3;
+    <%=struct%> *w = (<%=struct%>*)(lp->opt_ptr);
+    <% c_args = get(:postpose) ? "x1,x2,x3,w" : "w,x1,x2,x3" %>
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+    INIT_PTR(lp, 2, p3, s3);
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,double,x1);
+        GET_DATA_STRIDE(p2,s2,double,x2);
+        GET_DATA_STRIDE(p3,s3,double,x3);
+        <%=func_name%>(<%=c_args%>);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args[1..3].map{|a| a[1]}.join(",")%>)
+  @param  [DFloat]   <%=args[1][1]%>
+  @param  [DFloat]   <%=args[2][1]%>
+  @param  [DFloat]   <%=args[3][1]%>
+  @return [<%=class_name%>]  self
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(3)%>(VALUE self, VALUE v1, VALUE v2, VALUE v3)
+{
+    <%=struct%> *w;
+    ndfunc_arg_in_t ain[3] = {{cDF,0},{cDF,0},{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP, 3,0, ain,0};
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    na_ndloop3(&ndf, w, 3, v1, v2, v3);
+    return self;
+}

data/ext/numo/gsl/gen/tmpl/c_self_f_DFloat_x4.c

@@ -0,0 +1,47 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2, *p3, *p4;
+    ssize_t  s1, s2, s3, s4;
+    double   x1, x2, x3, x4;
+    <%=struct%> *w = (<%=struct%>*)(lp->opt_ptr);
+    <% c_args = get(:postpose) ? "x1,x2,x3,x4,w" : "w,x1,x2,x3,x4" %>
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+    INIT_PTR(lp, 2, p3, s3);
+    INIT_PTR(lp, 3, p4, s4);
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,double,x1);
+        GET_DATA_STRIDE(p2,s2,double,x2);
+        GET_DATA_STRIDE(p3,s3,double,x3);
+        GET_DATA_STRIDE(p4,s4,double,x4);
+        <%=func_name%>(<%=c_args%>);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args[1..4].map{|a| a[1]}.join(",")%>)
+  @param  [DFloat]   <%=args[1][1]%>
+  @param  [DFloat]   <%=args[2][1]%>
+  @param  [DFloat]   <%=args[3][1]%>
+  @param  [DFloat]   <%=args[4][1]%>
+  @return [<%=class_name%>]  self
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(4)%>(VALUE self, VALUE v1, VALUE v2, VALUE v3, VALUE v4)
+{
+    <%=struct%> *w;
+    ndfunc_arg_in_t ain[4] = {{cDF,0},{cDF,0},{cDF,0},{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP, 4,0, ain,0};
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    na_ndloop3(&ndf, w, 4, v1, v2, v3, v4);
+    return self;
+}

data/ext/numo/gsl/gen/tmpl/c_self_f_SZ_x2_DFloat.c

@@ -0,0 +1,45 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2, *p3;
+    ssize_t  s1, s2, s3;
+    size_t   x, y;
+    double   z;
+    <%=struct%> *w = (<%=struct%>*)(lp->opt_ptr);
+    <% c_args = get(:postpose) ? "x,y,z,w" : "w,x,y,z" %>
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+    INIT_PTR(lp, 2, p3, s3);
+
+    for (; i--; ) {
+        GET_DATA_STRIDE(p1,s1,size_t,x);
+        GET_DATA_STRIDE(p2,s2,size_t,y);
+        GET_DATA_STRIDE(p3,s3,double,z);
+        <%=func_name%>(<%=c_args%>);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args[1][1]%>,<%=args[2][1]%>,<%=args[2][1]%>)
+  @param  [UInt]      <%=args[1][1]%>
+  @param  [UInt]      <%=args[2][1]%>
+  @param  [DFloat]    <%=args[3][1]%>
+  @return [DFloat]    self
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(3)%>(VALUE self, VALUE v1, VALUE v2, VALUE v3)
+{
+    <%=struct%> *w;
+    ndfunc_arg_in_t ain[3] = {{cSZ,0},{cSZ,0},{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 3,0, ain,0};
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    na_ndloop3(&ndf, w, 3, v1, v2, v3);
+    return self;
+}

data/ext/numo/gsl/gen/tmpl/c_self_f_double.c

@@ -0,0 +1,17 @@
+/*
+  @overload <%=name%>(<%=args[1][1]%>)
+  @param  [Float]     <%=args[1][1]%>
+  @return [<%=full_class_name%>]   self
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE self, VALUE v1)
+{
+    <%=struct%> *w;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    <%=func_name%>(w, NUM2DBL(v1));
+    return self;
+}

data/ext/numo/gsl/gen/tmpl/c_self_f_other.c

@@ -0,0 +1,17 @@
+/*
+  @overload <%=name%>(<%=args[1][1]%>)
+  @param  [<%=full_class_name%>]  <%=args[1][1]%> other
+  @return [<%=full_class_name%>]  self
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE self, VALUE other)
+{
+    <%=struct%> *w, *x;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+    TypedData_Get_Struct(other, <%=struct%>, &<%=data_type_var%>, x);
+    <%=func_name%>(w,x);
+    return self;
+}

data/ext/numo/gsl/gen/tmpl/c_self_f_sizet.c

@@ -0,0 +1,21 @@
+/*
+  @overload <%=name%>(<%=args[1][1]%>)
+  @param  [Integer]   <%=args[1][1]%>
+  @return [<%=full_class_name%>]   self
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE self, VALUE v1)
+{
+    <%=struct%> *w;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    <% if get(:postpose) %>
+    <%=func_name%>(NUM2SIZET(v1),w);
+    <% else %>
+    <%=func_name%>(w, NUM2SIZET(v1));
+    <% end %>
+    return self;
+}

data/ext/numo/gsl/gen/tmpl/c_self_f_void.c

@@ -0,0 +1,14 @@
+/*
+  @overload <%=name%>
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(0)%>(VALUE self)
+{
+    <%=struct%> *w;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+    <%=func_name%>(w);
+    return self;
+}

data/ext/numo/gsl/gen/tmpl/c_sizet_f_void.c

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

data/ext/numo/gsl/gen/tmpl/c_sizet_x2_f_void.c

@@ -0,0 +1,17 @@
+/*
+  @overload <%=name%>()
+  @return [Array]   array of [<%=args[0][1]%>,<%=args[1][1]%>]
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(0)%>(VALUE self)
+{
+    size_t i, j;
+    <%=struct%> *w;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    <%=func_name%>(w, &i, &j);
+    return rb_assoc_new(SIZET2NUM(i), SIZET2NUM(j));
+}

data/ext/numo/gsl/gen/tmpl/c_str_f_void.c

@@ -0,0 +1,15 @@
+/*
+  @overload <%=name%>
+  @return [String]
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(0)%>(VALUE self)
+{
+    <%=struct%> *w;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    return rb_str_new_cstr(<%=func_name%>(w));
+}

data/ext/numo/gsl/gen/tmpl/c_super_new_void.c

@@ -0,0 +1,18 @@
+<%
+set name: "new"
+set singleton: true
+%>
+/*
+  :nodoc:
+ */
+static VALUE
+<%=c_func(1)%>(VALUE self, const <%=type_struct%> *t)
+{
+    <%=struct%> *w;
+
+    w = <%=func_name%>(t);
+    if (!w) {
+        rb_raise(rb_eNoMemError,"fail to allocate struct");
+    }
+    return TypedData_Wrap_Struct(<%=class_var%>, &<%=data_type_var%>, (void*)w);
+}

data/ext/numo/gsl/gen/tmpl/c_uint_f_void.c

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

data/ext/numo/gsl/gen/tmpl/c_ulong_f_void.c

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

data/ext/numo/gsl/gen/tmpl/c_void_f_double.c

@@ -0,0 +1,17 @@
+/*
+  @overload <%=name%>(<%=args[1][1]%>)
+  @param  [Float]  <%=args[1][1]%>
+  @return [Qnil]
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE self, VALUE v1)
+{
+    <%=struct%> *w;
+
+    TypedData_Get_Struct(self, <%=struct%>, &<%=data_type_var%>, w);
+
+    <%=func_name%>(w, NUM2DBL(v1));
+    return Qnil;
+}