ruby-dcl 1.7.1 → 1.8.2

data/ext/numru/grph2_umpack.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 *);
@@ -75,7 +78,7 @@ dcl_umpmap(obj, cdsn)
     i_cdsn = StringValuePtr(cdsn);
 
 
-    umpmap_(i_cdsn, (ftnlen)strlen(i_cdsn));
+    umpmap__(i_cdsn, (ftnlen)strlen(i_cdsn));
 
     return Qnil;
 
@@ -98,7 +101,7 @@ dcl_umqfnm(obj, cpara)
     o_cfname= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cfname, '\0', DFLT_SIZE+1);
 
-    umqfnm_(i_cpara, o_cfname, (ftnlen)strlen(i_cpara), (ftnlen)DFLT_SIZE);
+    umqfnm__(i_cpara, o_cfname, (ftnlen)strlen(i_cpara), (ftnlen)DFLT_SIZE);
 
     cfname = rb_str_new2(o_cfname);
 
@@ -111,7 +114,7 @@ static VALUE
 dcl_umpglb(obj)
     VALUE obj;
 {
-    umpglb_();
+    umpglb__();
 
     return Qnil;
 
@@ -121,7 +124,7 @@ static VALUE
 dcl_umpgrd(obj)
     VALUE obj;
 {
-    umpgrd_();
+    umpgrd__();
 
     return Qnil;
 
@@ -131,7 +134,7 @@ static VALUE
 dcl_umplim(obj)
     VALUE obj;
 {
-    umplim_();
+    umplim__();
 
     return Qnil;
 
@@ -141,7 +144,7 @@ static VALUE
 dcl_uminit(obj)
     VALUE obj;
 {
-    uminit_();
+    uminit__();
 
     return Qnil;
 
@@ -170,7 +173,7 @@ dcl_umscnt(obj, xcnt, ycnt, rot)
     i_rot = (real)NUM2DBL(rot);
 
 
-    umscnt_(&i_xcnt, &i_ycnt, &i_rot);
+    umscnt__(&i_xcnt, &i_ycnt, &i_rot);
 
     return Qnil;
 
@@ -187,7 +190,7 @@ dcl_umqcnt(obj)
     VALUE ycnt;
     VALUE rot;
 
-    umqcnt_(&o_xcnt, &o_ycnt, &o_rot);
+    umqcnt__(&o_xcnt, &o_ycnt, &o_rot);
 
     xcnt = rb_float_new((double)o_xcnt);
     ycnt = rb_float_new((double)o_ycnt);
@@ -221,7 +224,7 @@ dcl_umscwd(obj, xcntr, ycntr, r)
     i_r = (real)NUM2DBL(r);
 
 
-    umscwd_(&i_xcntr, &i_ycntr, &i_r);
+    umscwd__(&i_xcntr, &i_ycntr, &i_r);
 
     return Qnil;
 
@@ -238,7 +241,7 @@ dcl_umqcwd(obj)
     VALUE ycntr;
     VALUE r;
 
-    umqcwd_(&o_xcntr, &o_ycntr, &o_r);
+    umqcwd__(&o_xcntr, &o_ycntr, &o_r);
 
     xcntr = rb_float_new((double)o_xcntr);
     ycntr = rb_float_new((double)o_ycntr);
@@ -280,7 +283,7 @@ dcl_umspnt(obj, n, ux, uy)
     i_uy = dcl_obj2crealary(uy);
 
 
-    umspnt_(&i_n, i_ux, i_uy);
+    umspnt__(&i_n, i_ux, i_uy);
 
     dcl_freecrealary(i_ux);
     dcl_freecrealary(i_uy);
@@ -306,7 +309,7 @@ dcl_umqpnt(obj, n)
     i_n = NUM2INT(n);
 
 
-    umqpnt_(&i_n, &o_uxa, &o_uya);
+    umqpnt__(&i_n, &o_uxa, &o_uya);
 
     uxa = rb_float_new((double)o_uxa);
     uya = rb_float_new((double)o_uya);
@@ -323,7 +326,7 @@ dcl_umqptn(obj)
     integer o_n;
     VALUE n;
 
-    umqptn_(&o_n);
+    umqptn__(&o_n);
 
     n = INT2NUM(o_n);
 
@@ -336,7 +339,7 @@ static VALUE
 dcl_umrpnt(obj)
     VALUE obj;
 {
-    umrpnt_();
+    umrpnt__();
 
     return Qnil;
 
@@ -346,7 +349,7 @@ static VALUE
 dcl_umpfit(obj)
     VALUE obj;
 {
-    umpfit_();
+    umpfit__();
 
     return Qnil;
 
@@ -359,7 +362,7 @@ dcl_umpqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    umpqnp_(&o_ncp);
+    umpqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -383,7 +386,7 @@ dcl_umpqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umpqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    umpqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -409,7 +412,7 @@ dcl_umpqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    umpqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    umpqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -435,7 +438,7 @@ dcl_umpqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    umpqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    umpqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -459,7 +462,7 @@ dcl_umpqit(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    umpqit_(&i_idx, &o_itp);
+    umpqit__(&i_idx, &o_itp);
 
     itp = INT2NUM(o_itp);
 
@@ -483,7 +486,7 @@ dcl_umpqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    umpqvl_(&i_idx, &o_ipara);
+    umpqvl__(&i_idx, &o_ipara);
 
     ipara = INT2NUM(o_ipara);
 
@@ -510,7 +513,7 @@ dcl_umpsvl(obj, idx, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    umpsvl_(&i_idx, &i_ipara);
+    umpsvl__(&i_idx, &i_ipara);
 
     return Qnil;
 
@@ -531,7 +534,7 @@ dcl_umpqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umpqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    umpqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);
 
@@ -544,7 +547,7 @@ static VALUE
 dcl_umspct(obj)
     VALUE obj;
 {
-    umspct_();
+    umspct__();
 
     return Qnil;
 
@@ -554,7 +557,7 @@ static VALUE
 dcl_umspcw(obj)
     VALUE obj;
 {
-    umspcw_();
+    umspcw__();
 
     return Qnil;
 
@@ -564,7 +567,7 @@ static VALUE
 dcl_umspwd(obj)
     VALUE obj;
 {
-    umspwd_();
+    umspwd__();
 
     return Qnil;
 
@@ -574,7 +577,7 @@ static VALUE
 dcl_umsppt(obj)
     VALUE obj;
 {
-    umsppt_();
+    umsppt__();
 
     return Qnil;
 
@@ -584,7 +587,7 @@ static VALUE
 dcl_umspdf(obj)
     VALUE obj;
 {
-    umspdf_();
+    umspdf__();
 
     return Qnil;
 
@@ -609,7 +612,7 @@ dcl_umbndc(obj, xmin)
     io_xmin = (real)NUM2DBL(xmin);
 
 
-    umbndc_(&io_xmin, &o_xmax, &o_ymin, &o_ymax);
+    umbndc__(&io_xmin, &o_xmax, &o_ymin, &o_ymax);
 
     xmin = rb_float_new((double)io_xmin);
     xmax = rb_float_new((double)o_xmax);
@@ -640,7 +643,7 @@ dcl_umbndp(obj, vxmin)
     io_vxmin = (real)NUM2DBL(vxmin);
 
 
-    umbndp_(&io_vxmin, &o_vxmax, &o_vymin, &o_vymax);
+    umbndp__(&io_vxmin, &o_vxmax, &o_vymin, &o_vymax);
 
     vxmin = rb_float_new((double)io_vxmin);
     vxmax = rb_float_new((double)o_vxmax);
@@ -678,7 +681,7 @@ dcl_umbndr(obj, func, ftr)
     i_ftr = (real)NUM2DBL(ftr);
 
 
-    umbndr_(&i_func, &i_ftr, &o_xmin, &o_xmax, &o_ymin, &o_ymax);
+    umbndr__(&i_func, &i_ftr, &o_xmin, &o_xmax, &o_ymin, &o_ymax);
 
     xmin = rb_float_new((double)o_xmin);
     xmax = rb_float_new((double)o_xmax);
@@ -703,7 +706,7 @@ dcl_umqtxy(obj)
     VALUE tyminz;
     VALUE tymaxz;
 
-    umqtxy_(&o_txminz, &o_txmaxz, &o_tyminz, &o_tymaxz);
+    umqtxy__(&o_txminz, &o_txmaxz, &o_tyminz, &o_tymaxz);
 
     txminz = rb_float_new((double)o_txminz);
     txmaxz = rb_float_new((double)o_txmaxz);
@@ -719,7 +722,7 @@ static VALUE
 dcl_umstvz(obj)
     VALUE obj;
 {
-    umstvz_();
+    umstvz__();
 
     return Qnil;
 
@@ -729,7 +732,7 @@ static VALUE
 dcl_umsgrd(obj)
     VALUE obj;
 {
-    umsgrd_();
+    umsgrd__();
 
     return Qnil;
 
@@ -739,7 +742,7 @@ static VALUE
 dcl_umscom(obj)
     VALUE obj;
 {
-    umscom_();
+    umscom__();
 
     return Qnil;
 
@@ -760,7 +763,7 @@ dcl_umiget(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umiget_(i_cp, &o_ipara, (ftnlen)strlen(i_cp));
+    umiget__(i_cp, &o_ipara, (ftnlen)strlen(i_cp));
 
     ipara = INT2NUM(o_ipara);
 
@@ -787,7 +790,7 @@ dcl_umiset(obj, cp, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    umiset_(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
+    umiset__(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -811,7 +814,7 @@ dcl_umistx(obj, cp, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    umistx_(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
+    umistx__(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -824,7 +827,7 @@ dcl_umiqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    umiqnp_(&o_ncp);
+    umiqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -848,7 +851,7 @@ dcl_umiqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umiqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    umiqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -874,7 +877,7 @@ dcl_umiqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    umiqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    umiqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -900,7 +903,7 @@ dcl_umiqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    umiqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    umiqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -924,7 +927,7 @@ dcl_umiqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    umiqvl_(&i_idx, &o_ipara);
+    umiqvl__(&i_idx, &o_ipara);
 
     ipara = INT2NUM(o_ipara);
 
@@ -951,7 +954,7 @@ dcl_umisvl(obj, idx, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    umisvl_(&i_idx, &i_ipara);
+    umisvl__(&i_idx, &i_ipara);
 
     return Qnil;
 
@@ -972,7 +975,7 @@ dcl_umiqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umiqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    umiqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);
 
@@ -996,7 +999,7 @@ dcl_umlget(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umlget_(i_cp, &o_lpara, (ftnlen)strlen(i_cp));
+    umlget__(i_cp, &o_lpara, (ftnlen)strlen(i_cp));
 
     lpara = (o_lpara == FALSE_) ? Qfalse : Qtrue;
 
@@ -1020,7 +1023,7 @@ dcl_umlset(obj, cp, lpara)
     i_lpara = ((lpara == Qnil)||(lpara == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    umlset_(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
+    umlset__(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -1041,7 +1044,7 @@ dcl_umlstx(obj, cp, lpara)
     i_lpara = ((lpara == Qnil)||(lpara == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    umlstx_(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
+    umlstx__(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -1054,7 +1057,7 @@ dcl_umlqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    umlqnp_(&o_ncp);
+    umlqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -1078,7 +1081,7 @@ dcl_umlqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umlqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    umlqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -1104,7 +1107,7 @@ dcl_umlqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    umlqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    umlqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1130,7 +1133,7 @@ dcl_umlqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    umlqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    umlqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1154,7 +1157,7 @@ dcl_umlqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    umlqvl_(&i_idx, &o_lpara);
+    umlqvl__(&i_idx, &o_lpara);
 
     lpara = (o_lpara == FALSE_) ? Qfalse : Qtrue;
 
@@ -1178,7 +1181,7 @@ dcl_umlsvl(obj, idx, lpara)
     i_lpara = ((lpara == Qnil)||(lpara == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    umlsvl_(&i_idx, &i_lpara);
+    umlsvl__(&i_idx, &i_lpara);
 
     return Qnil;
 
@@ -1199,7 +1202,7 @@ dcl_umlqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umlqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    umlqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);
 
@@ -1223,7 +1226,7 @@ dcl_umrget(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umrget_(i_cp, &o_rpara, (ftnlen)strlen(i_cp));
+    umrget__(i_cp, &o_rpara, (ftnlen)strlen(i_cp));
 
     rpara = rb_float_new((double)o_rpara);
 
@@ -1250,7 +1253,7 @@ dcl_umrset(obj, cp, rpara)
     i_rpara = (real)NUM2DBL(rpara);
 
 
-    umrset_(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
+    umrset__(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -1274,7 +1277,7 @@ dcl_umrstx(obj, cp, rpara)
     i_rpara = (real)NUM2DBL(rpara);
 
 
-    umrstx_(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
+    umrstx__(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -1287,7 +1290,7 @@ dcl_umrqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    umrqnp_(&o_ncp);
+    umrqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -1311,7 +1314,7 @@ dcl_umrqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umrqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    umrqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -1337,7 +1340,7 @@ dcl_umrqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    umrqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    umrqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1363,7 +1366,7 @@ dcl_umrqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    umrqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    umrqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1387,7 +1390,7 @@ dcl_umrqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    umrqvl_(&i_idx, &o_rpara);
+    umrqvl__(&i_idx, &o_rpara);
 
     rpara = rb_float_new((double)o_rpara);
 
@@ -1414,7 +1417,7 @@ dcl_umrsvl(obj, idx, rpara)
     i_rpara = (real)NUM2DBL(rpara);
 
 
-    umrsvl_(&i_idx, &i_rpara);
+    umrsvl__(&i_idx, &i_rpara);
 
     return Qnil;
 
@@ -1435,7 +1438,7 @@ dcl_umrqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    umrqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    umrqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);
 
@@ -1459,7 +1462,7 @@ dcl_umfmap(obj, cdsn)
     i_cdsn = StringValuePtr(cdsn);
 
 
-    umfmap_(i_cdsn, (ftnlen)strlen(i_cdsn));
+    umfmap__(i_cdsn, (ftnlen)strlen(i_cdsn));
 
     return Qnil;