ruby-dcl 1.7.1 → 1.8.2

data/ext/numru/grph2_udpack.c

@@ -1,10 +1,9 @@
 /*
- * $Id: p_header,v 1.3 2015/03/18 03:27:42 horinout Exp $
+ * $Id$
  */
 
 #include <stdio.h>
 #include "ruby.h"
-#include "libtinyf2c.h"
 #include "narray.h"
 #include "cdcl.h"
 
@@ -15,6 +14,10 @@
 
 #define DFLT_SIZE 32
 
+#ifndef NARRAY_BIGMEM
+typedef int    na_shape_t;
+#endif
+
 extern char    *dcl_obj2ccharary(VALUE, int, int);
 extern integer *dcl_obj2cintegerary(VALUE);
 extern real    *dcl_obj2crealary(VALUE);
@@ -22,10 +25,10 @@ extern complex *dcl_obj2ccomplexary(VALUE);
 extern logical *dcl_obj2clogicalary(VALUE);
 
 extern VALUE dcl_ccharary2obj(char *, int, int);
-extern VALUE dcl_cintegerary2obj(integer *, int, int, int *);
-extern VALUE dcl_crealary2obj(real *, int, int, int *);
-extern VALUE dcl_ccomplexary2obj(complex *, int, char *);
-extern VALUE dcl_clogicalary2obj(logical *, int, int, int *);
+extern VALUE dcl_cintegerary2obj(integer *, na_shape_t, int, na_shape_t *);
+extern VALUE dcl_crealary2obj(real *, na_shape_t, int, na_shape_t *);
+/* extern VALUE dcl_ccomplexary2obj(complex *, int, char *); */
+extern VALUE dcl_clogicalary2obj(logical *, na_shape_t, int, na_shape_t *);
 
 extern void dcl_freeccharary(char *);
 extern void dcl_freecintegerary(integer *);
@@ -94,7 +97,7 @@ dcl_udcntr(obj, z, mx, nx, ny)
     i_z = dcl_obj2crealary(z);
 
 
-    udcntr_(i_z, &i_mx, &i_nx, &i_ny);
+    udcntr__(i_z, &i_mx, &i_nx, &i_ny);
 
     dcl_freecrealary(i_z);
 
@@ -141,7 +144,7 @@ dcl_udcntz(obj, z, mx, nx, ny, nbr2)
 
     w_ibr= ALLOCA_N(integer, (i_nbr2));
 
-    udcntz_(i_z, &i_mx, &i_nx, &i_ny, w_ibr, &i_nbr2);
+    udcntz__(i_z, &i_mx, &i_nx, &i_ny, w_ibr, &i_nbr2);
 
     dcl_freecrealary(i_z);
 
@@ -172,7 +175,7 @@ dcl_udgcla(obj, xmin, xmax, dx)
     i_dx = (real)NUM2DBL(dx);
 
 
-    udgcla_(&i_xmin, &i_xmax, &i_dx);
+    udgcla__(&i_xmin, &i_xmax, &i_dx);
 
     return Qnil;
 
@@ -215,7 +218,7 @@ dcl_udgclb(obj, z, mx, nx, ny, dx)
     i_z = dcl_obj2crealary(z);
 
 
-    udgclb_(i_z, &i_mx, &i_nx, &i_ny, &i_dx);
+    udgclb__(i_z, &i_mx, &i_nx, &i_ny, &i_dx);
 
     dcl_freecrealary(i_z);
 
@@ -227,7 +230,7 @@ static VALUE
 dcl_udiclv(obj)
     VALUE obj;
 {
-    udiclv_();
+    udiclv__();
 
     return Qnil;
 
@@ -266,7 +269,7 @@ dcl_udsclv(obj, zlev, indx, ityp, clv, hl)
     i_hl = (real)NUM2DBL(hl);
 
 
-    udsclv_(&i_zlev, &i_indx, &i_ityp, i_clv, &i_hl, (ftnlen)strlen(i_clv));
+    udsclv__(&i_zlev, &i_indx, &i_ityp, i_clv, &i_hl, (ftnlen)strlen(i_clv));
 
     return Qnil;
 
@@ -297,7 +300,7 @@ dcl_udqclv(obj, nlev)
     o_clv= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_clv, '\0', DFLT_SIZE+1);
 
-    udqclv_(&o_zlev, &o_indx, &o_ityp, o_clv, &o_hl, &i_nlev, (ftnlen)DFLT_SIZE);
+    udqclv__(&o_zlev, &o_indx, &o_ityp, o_clv, &o_hl, &i_nlev, (ftnlen)DFLT_SIZE);
 
     zlev = rb_float_new((double)o_zlev);
     indx = INT2NUM(o_indx);
@@ -317,7 +320,7 @@ dcl_udqcln(obj)
     integer o_nlev;
     VALUE nlev;
 
-    udqcln_(&o_nlev);
+    udqcln__(&o_nlev);
 
     nlev = INT2NUM(o_nlev);
 
@@ -339,7 +342,7 @@ dcl_uddclv(obj, zlev)
     i_zlev = (real)NUM2DBL(zlev);
 
 
-    uddclv_(&i_zlev);
+    uddclv__(&i_zlev);
 
     return Qnil;
 
@@ -360,7 +363,7 @@ dcl_rudlev(obj, nlev)
     i_nlev = NUM2INT(nlev);
 
 
-    o_rtn_val = rudlev_(&i_nlev);
+    o_rtn_val = rudlev__(&i_nlev);
 
     rtn_val = rb_float_new((double)o_rtn_val);
 
@@ -385,7 +388,7 @@ dcl_udiclr(obj, n)
 
     o_ix= ALLOCA_N(integer, (i_n));
 
-    udiclr_(o_ix, &i_n);
+    udiclr__(o_ix, &i_n);
 
     {int array_shape[1] = {i_n};
      ix = dcl_cintegerary2obj(o_ix, (i_n), 1, array_shape);
@@ -413,7 +416,7 @@ dcl_udlabl(obj, val)
     o_cval= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cval, '\0', DFLT_SIZE+1);
 
-    udlabl_(&i_val, o_cval, (ftnlen)DFLT_SIZE);
+    udlabl__(&i_val, o_cval, (ftnlen)DFLT_SIZE);
 
     cval = rb_str_new2(o_cval);
 
@@ -435,7 +438,7 @@ dcl_udsfmt(obj, cfmt)
     i_cfmt = StringValuePtr(cfmt);
 
 
-    udsfmt_(i_cfmt, (ftnlen)strlen(i_cfmt));
+    udsfmt__(i_cfmt, (ftnlen)strlen(i_cfmt));
 
     return Qnil;
 
@@ -451,7 +454,7 @@ dcl_udqfmt(obj)
     o_cfmt= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cfmt, '\0', DFLT_SIZE+1);
 
-    udqfmt_(o_cfmt, (ftnlen)DFLT_SIZE);
+    udqfmt__(o_cfmt, (ftnlen)DFLT_SIZE);
 
     cfmt = rb_str_new2(o_cfmt);
 
@@ -467,7 +470,7 @@ dcl_udpqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    udpqnp_(&o_ncp);
+    udpqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -491,7 +494,7 @@ dcl_udpqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udpqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    udpqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -517,7 +520,7 @@ dcl_udpqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    udpqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    udpqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -543,7 +546,7 @@ dcl_udpqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    udpqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    udpqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -567,7 +570,7 @@ dcl_udpqit(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    udpqit_(&i_idx, &o_itp);
+    udpqit__(&i_idx, &o_itp);
 
     itp = INT2NUM(o_itp);
 
@@ -591,7 +594,7 @@ dcl_udpqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    udpqvl_(&i_idx, &o_ipara);
+    udpqvl__(&i_idx, &o_ipara);
 
     ipara = INT2NUM(o_ipara);
 
@@ -618,7 +621,7 @@ dcl_udpsvl(obj, idx, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    udpsvl_(&i_idx, &i_ipara);
+    udpsvl__(&i_idx, &i_ipara);
 
     return Qnil;
 
@@ -639,7 +642,7 @@ dcl_udpqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udpqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    udpqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);
 
@@ -663,7 +666,7 @@ dcl_udiget(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udiget_(i_cp, &o_ipara, (ftnlen)strlen(i_cp));
+    udiget__(i_cp, &o_ipara, (ftnlen)strlen(i_cp));
 
     ipara = INT2NUM(o_ipara);
 
@@ -690,7 +693,7 @@ dcl_udiset(obj, cp, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    udiset_(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
+    udiset__(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -714,7 +717,7 @@ dcl_udistx(obj, cp, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    udistx_(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
+    udistx__(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -727,7 +730,7 @@ dcl_udiqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    udiqnp_(&o_ncp);
+    udiqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -751,7 +754,7 @@ dcl_udiqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udiqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    udiqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -777,7 +780,7 @@ dcl_udiqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    udiqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    udiqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -803,7 +806,7 @@ dcl_udiqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    udiqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    udiqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -827,7 +830,7 @@ dcl_udiqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    udiqvl_(&i_idx, &o_ipara);
+    udiqvl__(&i_idx, &o_ipara);
 
     ipara = INT2NUM(o_ipara);
 
@@ -854,7 +857,7 @@ dcl_udisvl(obj, idx, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    udisvl_(&i_idx, &i_ipara);
+    udisvl__(&i_idx, &i_ipara);
 
     return Qnil;
 
@@ -875,7 +878,7 @@ dcl_udiqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udiqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    udiqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);
 
@@ -899,7 +902,7 @@ dcl_udlget(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udlget_(i_cp, &o_lpara, (ftnlen)strlen(i_cp));
+    udlget__(i_cp, &o_lpara, (ftnlen)strlen(i_cp));
 
     lpara = (o_lpara == FALSE_) ? Qfalse : Qtrue;
 
@@ -923,7 +926,7 @@ dcl_udlset(obj, cp, lpara)
     i_lpara = ((lpara == Qnil)||(lpara == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    udlset_(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
+    udlset__(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -944,7 +947,7 @@ dcl_udlstx(obj, cp, lpara)
     i_lpara = ((lpara == Qnil)||(lpara == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    udlstx_(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
+    udlstx__(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -957,7 +960,7 @@ dcl_udlqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    udlqnp_(&o_ncp);
+    udlqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -981,7 +984,7 @@ dcl_udlqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udlqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    udlqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -1007,7 +1010,7 @@ dcl_udlqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    udlqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    udlqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1033,7 +1036,7 @@ dcl_udlqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    udlqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    udlqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1057,7 +1060,7 @@ dcl_udlqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    udlqvl_(&i_idx, &o_lpara);
+    udlqvl__(&i_idx, &o_lpara);
 
     lpara = (o_lpara == FALSE_) ? Qfalse : Qtrue;
 
@@ -1081,7 +1084,7 @@ dcl_udlsvl(obj, idx, lpara)
     i_lpara = ((lpara == Qnil)||(lpara == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    udlsvl_(&i_idx, &i_lpara);
+    udlsvl__(&i_idx, &i_lpara);
 
     return Qnil;
 
@@ -1102,7 +1105,7 @@ dcl_udlqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udlqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    udlqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);
 
@@ -1126,7 +1129,7 @@ dcl_udrget(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udrget_(i_cp, &o_rpara, (ftnlen)strlen(i_cp));
+    udrget__(i_cp, &o_rpara, (ftnlen)strlen(i_cp));
 
     rpara = rb_float_new((double)o_rpara);
 
@@ -1153,7 +1156,7 @@ dcl_udrset(obj, cp, rpara)
     i_rpara = (real)NUM2DBL(rpara);
 
 
-    udrset_(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
+    udrset__(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -1177,7 +1180,7 @@ dcl_udrstx(obj, cp, rpara)
     i_rpara = (real)NUM2DBL(rpara);
 
 
-    udrstx_(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
+    udrstx__(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -1190,7 +1193,7 @@ dcl_udrqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    udrqnp_(&o_ncp);
+    udrqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -1214,7 +1217,7 @@ dcl_udrqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udrqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    udrqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -1240,7 +1243,7 @@ dcl_udrqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    udrqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    udrqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1266,7 +1269,7 @@ dcl_udrqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    udrqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    udrqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1290,7 +1293,7 @@ dcl_udrqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    udrqvl_(&i_idx, &o_rpara);
+    udrqvl__(&i_idx, &o_rpara);
 
     rpara = rb_float_new((double)o_rpara);
 
@@ -1317,7 +1320,7 @@ dcl_udrsvl(obj, idx, rpara)
     i_rpara = (real)NUM2DBL(rpara);
 
 
-    udrsvl_(&i_idx, &i_rpara);
+    udrsvl__(&i_idx, &i_rpara);
 
     return Qnil;
 
@@ -1338,7 +1341,7 @@ dcl_udrqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    udrqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    udrqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);