ruby-dcl 1.7.1 → 1.8.2

data/ext/numru/grph2_uepack.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_uetone(obj, z, mx, nx, ny)
     i_z = dcl_obj2crealary(z);
 
 
-    uetone_(i_z, &i_mx, &i_nx, &i_ny);
+    uetone__(i_z, &i_mx, &i_nx, &i_ny);
 
     dcl_freecrealary(i_z);
 
@@ -134,7 +137,7 @@ dcl_uetonf(obj, z, mx, nx, ny)
     i_z = dcl_obj2crealary(z);
 
 
-    uetonf_(i_z, &i_mx, &i_nx, &i_ny);
+    uetonf__(i_z, &i_mx, &i_nx, &i_ny);
 
     dcl_freecrealary(i_z);
 
@@ -176,7 +179,7 @@ dcl_uetonb(obj, z, mx, nx, ny)
     i_z = dcl_obj2crealary(z);
 
 
-    uetonb_(i_z, &i_mx, &i_nx, &i_ny);
+    uetonb__(i_z, &i_mx, &i_nx, &i_ny);
 
     dcl_freecrealary(i_z);
 
@@ -218,7 +221,7 @@ dcl_uetonc(obj, z, mx, nx, ny)
     i_z = dcl_obj2crealary(z);
 
 
-    uetonc_(i_z, &i_mx, &i_nx, &i_ny);
+    uetonc__(i_z, &i_mx, &i_nx, &i_ny);
 
     dcl_freecrealary(i_z);
 
@@ -265,7 +268,7 @@ dcl_uezchk(obj, z, mx, nx, ny, cname)
     i_z = dcl_obj2crealary(z);
 
 
-    uezchk_(i_z, &i_mx, &i_nx, &i_ny, i_cname, &o_istat, (ftnlen)strlen(i_cname));
+    uezchk__(i_z, &i_mx, &i_nx, &i_ny, i_cname, &o_istat, (ftnlen)strlen(i_cname));
 
     istat = INT2NUM(o_istat);
 
@@ -298,7 +301,7 @@ dcl_uegtla(obj, xmin, xmax, dx)
     i_dx = (real)NUM2DBL(dx);
 
 
-    uegtla_(&i_xmin, &i_xmax, &i_dx);
+    uegtla__(&i_xmin, &i_xmax, &i_dx);
 
     return Qnil;
 
@@ -341,7 +344,7 @@ dcl_uegtlb(obj, z, mx, nx, ny, dx)
     i_z = dcl_obj2crealary(z);
 
 
-    uegtlb_(i_z, &i_mx, &i_nx, &i_ny, &i_dx);
+    uegtlb__(i_z, &i_mx, &i_nx, &i_ny, &i_dx);
 
     dcl_freecrealary(i_z);
 
@@ -353,7 +356,7 @@ static VALUE
 dcl_ueitlv(obj)
     VALUE obj;
 {
-    ueitlv_();
+    ueitlv__();
 
     return Qnil;
 
@@ -382,7 +385,7 @@ dcl_uestlv(obj, tlev1, tlev2, ipat)
     i_ipat = NUM2INT(ipat);
 
 
-    uestlv_(&i_tlev1, &i_tlev2, &i_ipat);
+    uestlv__(&i_tlev1, &i_tlev2, &i_ipat);
 
     return Qnil;
 
@@ -407,7 +410,7 @@ dcl_ueqtlv(obj, iton)
     i_iton = NUM2INT(iton);
 
 
-    ueqtlv_(&o_tlev1, &o_tlev2, &o_ipat, &i_iton);
+    ueqtlv__(&o_tlev1, &o_tlev2, &o_ipat, &i_iton);
 
     tlev1 = rb_float_new((double)o_tlev1);
     tlev2 = rb_float_new((double)o_tlev2);
@@ -425,7 +428,7 @@ dcl_ueqntl(obj)
     integer o_nton;
     VALUE nton;
 
-    ueqntl_(&o_nton);
+    ueqntl__(&o_nton);
 
     nton = INT2NUM(o_nton);
 
@@ -465,7 +468,7 @@ dcl_uestln(obj, tlevn, ipatn, nton)
     i_ipatn = dcl_obj2cintegerary(ipatn);
 
 
-    uestln_(i_tlevn, i_ipatn, &i_nton);
+    uestln__(i_tlevn, i_ipatn, &i_nton);
 
     dcl_freecrealary(i_tlevn);
     dcl_freecintegerary(i_ipatn);
@@ -489,7 +492,7 @@ dcl_iueton(obj, zlev)
     i_zlev = (real)NUM2DBL(zlev);
 
 
-    o_rtn_val = iueton_(&i_zlev);
+    o_rtn_val = iueton__(&i_zlev);
 
     rtn_val = INT2NUM(o_rtn_val);
 
@@ -505,7 +508,7 @@ dcl_uepqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    uepqnp_(&o_ncp);
+    uepqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -529,7 +532,7 @@ dcl_uepqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    uepqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    uepqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -555,7 +558,7 @@ dcl_uepqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    uepqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    uepqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -581,7 +584,7 @@ dcl_uepqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    uepqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    uepqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -605,7 +608,7 @@ dcl_uepqit(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    uepqit_(&i_idx, &o_itp);
+    uepqit__(&i_idx, &o_itp);
 
     itp = INT2NUM(o_itp);
 
@@ -629,7 +632,7 @@ dcl_uepqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    uepqvl_(&i_idx, &o_ipara);
+    uepqvl__(&i_idx, &o_ipara);
 
     ipara = INT2NUM(o_ipara);
 
@@ -656,7 +659,7 @@ dcl_uepsvl(obj, idx, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    uepsvl_(&i_idx, &i_ipara);
+    uepsvl__(&i_idx, &i_ipara);
 
     return Qnil;
 
@@ -677,7 +680,7 @@ dcl_uepqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    uepqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    uepqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);
 
@@ -701,7 +704,7 @@ dcl_ueiget(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    ueiget_(i_cp, &o_ipara, (ftnlen)strlen(i_cp));
+    ueiget__(i_cp, &o_ipara, (ftnlen)strlen(i_cp));
 
     ipara = INT2NUM(o_ipara);
 
@@ -728,7 +731,7 @@ dcl_ueiset(obj, cp, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    ueiset_(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
+    ueiset__(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -752,7 +755,7 @@ dcl_ueistx(obj, cp, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    ueistx_(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
+    ueistx__(i_cp, &i_ipara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -765,7 +768,7 @@ dcl_ueiqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    ueiqnp_(&o_ncp);
+    ueiqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -789,7 +792,7 @@ dcl_ueiqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    ueiqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    ueiqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -815,7 +818,7 @@ dcl_ueiqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    ueiqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    ueiqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -841,7 +844,7 @@ dcl_ueiqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    ueiqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    ueiqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -865,7 +868,7 @@ dcl_ueiqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    ueiqvl_(&i_idx, &o_ipara);
+    ueiqvl__(&i_idx, &o_ipara);
 
     ipara = INT2NUM(o_ipara);
 
@@ -892,7 +895,7 @@ dcl_ueisvl(obj, idx, ipara)
     i_ipara = NUM2INT(ipara);
 
 
-    ueisvl_(&i_idx, &i_ipara);
+    ueisvl__(&i_idx, &i_ipara);
 
     return Qnil;
 
@@ -913,7 +916,7 @@ dcl_ueiqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    ueiqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    ueiqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);
 
@@ -937,7 +940,7 @@ dcl_uelget(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    uelget_(i_cp, &o_lpara, (ftnlen)strlen(i_cp));
+    uelget__(i_cp, &o_lpara, (ftnlen)strlen(i_cp));
 
     lpara = (o_lpara == FALSE_) ? Qfalse : Qtrue;
 
@@ -961,7 +964,7 @@ dcl_uelset(obj, cp, lpara)
     i_lpara = ((lpara == Qnil)||(lpara == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    uelset_(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
+    uelset__(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -982,7 +985,7 @@ dcl_uelstx(obj, cp, lpara)
     i_lpara = ((lpara == Qnil)||(lpara == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    uelstx_(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
+    uelstx__(i_cp, &i_lpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -995,7 +998,7 @@ dcl_uelqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    uelqnp_(&o_ncp);
+    uelqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -1019,7 +1022,7 @@ dcl_uelqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    uelqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    uelqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -1045,7 +1048,7 @@ dcl_uelqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    uelqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    uelqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1071,7 +1074,7 @@ dcl_uelqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    uelqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    uelqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1095,7 +1098,7 @@ dcl_uelqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    uelqvl_(&i_idx, &o_lpara);
+    uelqvl__(&i_idx, &o_lpara);
 
     lpara = (o_lpara == FALSE_) ? Qfalse : Qtrue;
 
@@ -1119,7 +1122,7 @@ dcl_uelsvl(obj, idx, lpara)
     i_lpara = ((lpara == Qnil)||(lpara == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    uelsvl_(&i_idx, &i_lpara);
+    uelsvl__(&i_idx, &i_lpara);
 
     return Qnil;
 
@@ -1140,7 +1143,7 @@ dcl_uelqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    uelqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    uelqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);
 
@@ -1164,7 +1167,7 @@ dcl_uerget(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    uerget_(i_cp, &o_rpara, (ftnlen)strlen(i_cp));
+    uerget__(i_cp, &o_rpara, (ftnlen)strlen(i_cp));
 
     rpara = rb_float_new((double)o_rpara);
 
@@ -1191,7 +1194,7 @@ dcl_uerset(obj, cp, rpara)
     i_rpara = (real)NUM2DBL(rpara);
 
 
-    uerset_(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
+    uerset__(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -1215,7 +1218,7 @@ dcl_uerstx(obj, cp, rpara)
     i_rpara = (real)NUM2DBL(rpara);
 
 
-    uerstx_(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
+    uerstx__(i_cp, &i_rpara, (ftnlen)strlen(i_cp));
 
     return Qnil;
 
@@ -1228,7 +1231,7 @@ dcl_uerqnp(obj)
     integer o_ncp;
     VALUE ncp;
 
-    uerqnp_(&o_ncp);
+    uerqnp__(&o_ncp);
 
     ncp = INT2NUM(o_ncp);
 
@@ -1252,7 +1255,7 @@ dcl_uerqid(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    uerqid_(i_cp, &o_idx, (ftnlen)strlen(i_cp));
+    uerqid__(i_cp, &o_idx, (ftnlen)strlen(i_cp));
 
     idx = INT2NUM(o_idx);
 
@@ -1278,7 +1281,7 @@ dcl_uerqcp(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    uerqcp_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    uerqcp__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1304,7 +1307,7 @@ dcl_uerqcl(obj, idx)
     o_cp= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cp, '\0', DFLT_SIZE+1);
 
-    uerqcl_(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
+    uerqcl__(&i_idx, o_cp, (ftnlen)DFLT_SIZE);
 
     cp = rb_str_new2(o_cp);
 
@@ -1328,7 +1331,7 @@ dcl_uerqvl(obj, idx)
     i_idx = NUM2INT(idx);
 
 
-    uerqvl_(&i_idx, &o_rpara);
+    uerqvl__(&i_idx, &o_rpara);
 
     rpara = rb_float_new((double)o_rpara);
 
@@ -1355,7 +1358,7 @@ dcl_uersvl(obj, idx, rpara)
     i_rpara = (real)NUM2DBL(rpara);
 
 
-    uersvl_(&i_idx, &i_rpara);
+    uersvl__(&i_idx, &i_rpara);
 
     return Qnil;
 
@@ -1376,7 +1379,7 @@ dcl_uerqin(obj, cp)
     i_cp = StringValuePtr(cp);
 
 
-    uerqin_(i_cp, &o_in, (ftnlen)strlen(i_cp));
+    uerqin__(i_cp, &o_in, (ftnlen)strlen(i_cp));
 
     in = INT2NUM(o_in);