numo-gsl-alt 0.1.3

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

@@ -0,0 +1,38 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2, *p3;
+    ssize_t  s1, s2, s3;
+    double   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,double,x);
+        GET_DATA_STRIDE(p2,s2,double,y);
+        z = <%=func_name%>(x,y);
+        SET_DATA_STRIDE(p3,s3,double,z);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args[0][1]%>, <%=args[1][1]%>)
+  @param  [DFloat]   <%=args[0][1]%>
+  @param  [DFloat]   <%=args[1][1]%>
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(2)%>(VALUE mod, VALUE v0, VALUE v1)
+{
+    ndfunc_arg_in_t ain[2] = {{cDF,0},{cDF,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 2,1, ain,aout};
+
+    return na_ndloop(&ndf, 2, v0, v1);
+}

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

@@ -0,0 +1,41 @@
+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, y;
+
+    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);
+        y = <%=func_name%>(x1,x2,x3);
+        SET_DATA_STRIDE(p4,s4,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args[0][1]%>, <%=args[1][1]%>, <%=args[2][1]%>)
+  @param  [DFloat]   <%=args[0][1]%>
+  @param  [DFloat]   <%=args[1][1]%>
+  @param  [DFloat]   <%=args[2][1]%>
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(3)%>(VALUE mod, VALUE v0, VALUE v1, VALUE v2)
+{
+    ndfunc_arg_in_t ain[3] = {{cDF,0},{cDF,0},{cDF,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 3,1, ain,aout};
+
+    return na_ndloop(&ndf, 3, v0, v1, v2);
+}

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

@@ -0,0 +1,36 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    int      x;
+    double   y;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    for (; i--; ) {
+        GET_DATA_STRIDE(p1,s1,int,x);
+        y = <%=func_name%>(x);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args[0][1]%>)
+  @param  [Int]  <%=args[0][1]%>
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE mod, VALUE v0)
+{
+    ndfunc_arg_in_t ain[1] = {{cInt,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+
+    return na_ndloop(&ndf, 1, v0);
+}

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

@@ -0,0 +1,36 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    unsigned int x;
+    double   y;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,unsigned int,x);
+        y = <%=func_name%>(x);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [UInt]     <%=args[0][1]%>
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE mod, VALUE v0)
+{
+    ndfunc_arg_in_t ain[1] = {{cUI,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+
+    return na_ndloop(&ndf, 1, v0);
+}

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

@@ -0,0 +1,36 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    unsigned int x;
+    double   y;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,unsigned int,x);
+        y = <%=func_name%>(x);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [UInt32]   <%=args[0][1]%>
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(1)%>(VALUE mod, VALUE v0)
+{
+    ndfunc_arg_in_t ain[1] = {{numo_cUInt32,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+
+    return na_ndloop(&ndf, 1, v0);
+}

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

@@ -0,0 +1,43 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    unsigned int x;
+    double   y;
+    double   c1;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    c1 = *(double*)(lp->opt_ptr);
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,unsigned int,x);
+        y = <%=func_name%>(x, c1);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [UInt32]   <%=args[0][1]%>
+  @param  [Float]    <%=args[1][1]%> parameter
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(2)%>(VALUE mod, VALUE v0, VALUE v1)
+{
+    ndfunc_arg_in_t ain[1] = {{numo_cUInt32,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+    double c1;
+
+    c1 = NUM2DBL(v1);
+
+    return na_ndloop3(&ndf, &c1, 1, v0);
+}

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

@@ -0,0 +1,49 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    unsigned int x;
+    double   y;
+    struct opt_d_u *opt;
+    double c1;
+    unsigned int c2;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    opt = (struct opt_d_u*)(lp->opt_ptr);
+    c1 = opt->d;
+    c2 = opt->u;
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,double,x);
+        y = <%=func_name%>(x, c1, c2);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [UInt32]   <%=args[0][1]%>
+  @param  [Integer]  <%=args[1][1]%> parameter
+  @param  [Float]    <%=args[2][1]%> parameter
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(3)%>(VALUE mod, VALUE v0, VALUE v1, VALUE v2)
+{
+    ndfunc_arg_in_t ain[1] = {{numo_cUInt32,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+    struct opt_d_u opt;
+
+    opt.d = NUM2DBL(v1);
+    opt.u = NUM2UINT(v2);
+
+    return na_ndloop3(&ndf, &opt, 1, v0);
+}

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

@@ -0,0 +1,48 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    unsigned int x;
+    double   y;
+    double  *opt;
+    double   c1, c2;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    opt = (double*)(lp->opt_ptr);
+    c1 = opt[0];
+    c2 = opt[1];
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,double,x);
+        y = <%=func_name%>(x, c1, c2);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [UInt32]   <%=args[0][1]%>
+  @param  [Float]    <%=args[1][1]%> parameter
+  @param  [Float]    <%=args[2][1]%> parameter
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(3)%>(VALUE mod , VALUE v0, VALUE v1, VALUE v2)
+{
+    ndfunc_arg_in_t ain[1] = {{numo_cUInt32,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+    double opt[2];
+
+    opt[0] = NUM2DBL(v1);
+    opt[1] = NUM2DBL(v2);
+
+    return na_ndloop3(&ndf, opt, 1, v0);
+}

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

@@ -0,0 +1,51 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    unsigned int x;
+    double   y;
+    unsigned int *opt;
+    unsigned int c1, c2, c3;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    opt = (unsigned int*)(lp->opt_ptr);
+    c1 = opt[0];
+    c2 = opt[1];
+    c3 = opt[2];
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,double,x);
+        y = <%=func_name%>(x, c1, c2, c3);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [UInt32]   <%=args[0][1]%>
+  @param  [Integer]  <%=args[1][1]%> parameter
+  @param  [Integer]  <%=args[2][1]%> parameter
+  @param  [Integer]  <%=args[3][1]%> parameter
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(4)%>(VALUE mod, VALUE v0, VALUE v1, VALUE v2, VALUE v3)
+{
+    ndfunc_arg_in_t ain[1] = {{numo_cUInt32,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+    unsigned int opt[3];
+
+    opt[0] = NUM2UINT(v1);
+    opt[1] = NUM2UINT(v2);
+    opt[2] = NUM2UINT(v3);
+
+    return na_ndloop3(&ndf, opt, 1, v0);
+}

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

@@ -0,0 +1,43 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    unsigned int x;
+    double   y;
+    double   c1;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    c1 = *(double*)(lp->opt_ptr);
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,unsigned int,x);
+        y = <%=func_name%>(x, c1);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [UInt]     <%=args[0][1]%>
+  @param  [Float]    <%=args[1][1]%> parameter
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(2)%>(VALUE mod, VALUE v0, VALUE v1)
+{
+    ndfunc_arg_in_t ain[1] = {{cUI,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+    double c1;
+
+    c1 = NUM2DBL(v1);
+
+    return na_ndloop3(&ndf, &c1, 1, v0);
+}

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

@@ -0,0 +1,50 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    unsigned int x;
+    double   y;
+    char *opt;
+    double c1;
+    unsigned int c2;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    opt = (char*)(lp->opt_ptr);
+    c1 = *(double*)opt;
+    c2 = *(unsigned int*)(opt+sizeof(double));
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,double,x);
+        y = <%=func_name%>(x, c1, c2);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [UInt]     <%=args[0][1]%>
+  @param  [Integer]  <%=args[1][1]%> parameter
+  @param  [Float]    <%=args[2][1]%> parameter
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(3)%>(VALUE mod, VALUE v0, VALUE v1, VALUE v2)
+{
+    ndfunc_arg_in_t ain[1] = {{cUI,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+    char *opt;
+    opt = ALLOCA_N(char,sizeof(double)+sizeof(unsigned int));
+
+    *(double*)opt = NUM2DBL(v1);
+    *(unsigned int*)(opt+sizeof(double)) = NUM2UINT(v2);
+
+    return na_ndloop3(&ndf, opt, 1, v0);
+}

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

@@ -0,0 +1,48 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    unsigned int x;
+    double   y;
+    double  *opt;
+    double   c1, c2;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    opt = (double*)(lp->opt_ptr);
+    c1 = opt[0];
+    c2 = opt[1];
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,double,x);
+        y = <%=func_name%>(x, c1, c2);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [UInt]     <%=args[0][1]%>
+  @param  [Float]    <%=args[1][1]%> parameter
+  @param  [Float]    <%=args[2][1]%> parameter
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(3)%>(VALUE mod , VALUE v0, VALUE v1, VALUE v2)
+{
+    ndfunc_arg_in_t ain[1] = {{cUI,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+    double opt[2];
+
+    opt[0] = NUM2DBL(v1);
+    opt[1] = NUM2DBL(v2);
+
+    return na_ndloop3(&ndf, opt, 1, v0);
+}

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

@@ -0,0 +1,51 @@
+static void
+iter_<%=c_func%>(na_loop_t *const lp)
+{
+    size_t   i;
+    char    *p1, *p2;
+    ssize_t  s1, s2;
+    unsigned int x;
+    double   y;
+    unsigned int *opt;
+    unsigned int c1, c2, c3;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    opt = (unsigned int*)(lp->opt_ptr);
+    c1 = opt[0];
+    c2 = opt[1];
+    c3 = opt[2];
+
+    for (; i--;) {
+        GET_DATA_STRIDE(p1,s1,double,x);
+        y = <%=func_name%>(x, c1, c2, c3);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [UInt]     <%=args[0][1]%>
+  @param  [Integer]  <%=args[1][1]%> parameter
+  @param  [Integer]  <%=args[2][1]%> parameter
+  @param  [Integer]  <%=args[3][1]%> parameter
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(4)%>(VALUE mod, VALUE v0, VALUE v1, VALUE v2, VALUE v3)
+{
+    ndfunc_arg_in_t ain[1] = {{cUI,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+    unsigned int opt[3];
+
+    opt[0] = NUM2UINT(v1);
+    opt[1] = NUM2UINT(v2);
+    opt[2] = NUM2UINT(v3);
+
+    return na_ndloop3(&ndf, opt, 1, v0);
+}

data/ext/numo/gsl/gen/tmpl/m_DFloat_f_int_DFloat.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;
+    int      c0;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+    c0 = *(int*)(lp->opt_ptr);
+
+    for (; i--; ) {
+        GET_DATA_STRIDE(p1,s1,double,x);
+        y = <%=func_name%>(c0, x);
+        SET_DATA_STRIDE(p2,s2,double,y);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args[0][1]%>, <%=args[1][1]%>)
+  @param  [Integer]   <%=args[0][1]%>
+  @param  [DFloat]    <%=args[1][1]%>
+  @return [DFloat]    result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(2)%>(VALUE mod, VALUE v0, VALUE v1)
+{
+    ndfunc_arg_in_t ain[1] = {{cDF,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 1,1, ain,aout};
+    int c0;
+
+    c0 = NUM2INT(v0);
+
+    return na_ndloop3(&ndf, &c0, 1, v1);
+}

data/ext/numo/gsl/gen/tmpl/m_DFloat_f_int_DFloat_x2.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   x, y, z;
+    int     *opts;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR(lp, 0, p1, s1);
+    INIT_PTR(lp, 1, p2, s2);
+    INIT_PTR(lp, 2, p3, s3);
+    opts = (int*)(lp->opt_ptr);
+
+    for (; i--; ) {
+        GET_DATA_STRIDE(p1,s1,double,x);
+        GET_DATA_STRIDE(p2,s2,double,y);
+        z = <%=func_name%>(opts[0],x,y);
+        SET_DATA_STRIDE(p3,s3,double,z);
+    }
+}
+
+/*
+  @overload <%=name%>(<%=args.map{|a| a[1]}.join(",")%>)
+  @param  [Integer]  <%=args[0][1]%> parameter
+  @param  [DFloat]   <%=args[1][1]%>
+  @param  [DFloat]   <%=args[2][1]%>
+  @return [DFloat]   result
+
+  <%= description %>
+*/
+static VALUE
+<%=c_func(3)%>(VALUE mod, VALUE v0, VALUE v1, VALUE v2)
+{
+    ndfunc_arg_in_t ain[2] = {{cDF,0},{cDF,0}};
+    ndfunc_arg_out_t aout[1] = {{cDF,0}};
+    ndfunc_t ndf = {iter_<%=c_func%>, STRIDE_LOOP|NDF_EXTRACT, 2,1, ain,aout};
+    int opts[1];
+
+    opts[0] = NUM2INT(v0);
+
+    return na_ndloop3(&ndf, opts, 2, v1, v2);
+}