ruby-dcl 1.7.1 → 1.8.2

data/ext/numru/grph1_slpack.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 *);
@@ -85,7 +88,7 @@ dcl_slinit(obj, wxmax, wymax, fact)
     i_fact = (real)NUM2DBL(fact);
 
 
-    slinit_(&i_wxmax, &i_wymax, &i_fact);
+    slinit__(&i_wxmax, &i_wymax, &i_fact);
 
     return Qnil;
 
@@ -104,7 +107,7 @@ dcl_slsize(obj, cszez)
     i_cszez = StringValuePtr(cszez);
 
 
-    slsize_(i_cszez, (ftnlen)strlen(i_cszez));
+    slsize__(i_cszez, (ftnlen)strlen(i_cszez));
 
     return Qnil;
 
@@ -128,7 +131,7 @@ dcl_slform(obj, dxa, dya)
     i_dya = (real)NUM2DBL(dya);
 
 
-    slform_(&i_dxa, &i_dya);
+    slform__(&i_dxa, &i_dya);
 
     return Qnil;
 
@@ -157,7 +160,7 @@ dcl_sldiv(obj, cform, ix, iy)
     i_iy = NUM2INT(iy);
 
 
-    sldiv_(i_cform, &i_ix, &i_iy, (ftnlen)strlen(i_cform));
+    sldiv__(i_cform, &i_ix, &i_iy, (ftnlen)strlen(i_cform));
 
     return Qnil;
 
@@ -191,7 +194,7 @@ dcl_slmgn(obj, xl, xr, yb, yt)
     i_yt = (real)NUM2DBL(yt);
 
 
-    slmgn_(&i_xl, &i_xr, &i_yb, &i_yt);
+    slmgn__(&i_xl, &i_xr, &i_yb, &i_yt);
 
     return Qnil;
 
@@ -215,7 +218,7 @@ dcl_slrat(obj, rx, ry)
     i_ry = (real)NUM2DBL(ry);
 
 
-    slrat_(&i_rx, &i_ry);
+    slrat__(&i_rx, &i_ry);
 
     return Qnil;
 
@@ -259,7 +262,7 @@ dcl_slsttl(obj, cttl, cside, px, py, ht, nt)
     i_nt = NUM2INT(nt);
 
 
-    slsttl_(i_cttl, i_cside, &i_px, &i_py, &i_ht, &i_nt, (ftnlen)strlen(i_cttl), (ftnlen)strlen(i_cside));
+    slsttl__(i_cttl, i_cside, &i_px, &i_py, &i_ht, &i_nt, (ftnlen)strlen(i_cttl), (ftnlen)strlen(i_cside));
 
     return Qnil;
 
@@ -278,7 +281,7 @@ dcl_sldttl(obj, nt)
     i_nt = NUM2INT(nt);
 
 
-    sldttl_(&i_nt);
+    sldttl__(&i_nt);
 
     return Qnil;
 
@@ -297,7 +300,7 @@ dcl_slpvpr(obj, index)
     i_index = NUM2INT(index);
 
 
-    slpvpr_(&i_index);
+    slpvpr__(&i_index);
 
     return Qnil;
 
@@ -316,7 +319,7 @@ dcl_slpwwr(obj, index)
     i_index = NUM2INT(index);
 
 
-    slpwwr_(&i_index);
+    slpwwr__(&i_index);
 
     return Qnil;
 
@@ -335,7 +338,7 @@ dcl_slpwvr(obj, index)
     i_index = NUM2INT(index);
 
 
-    slpwvr_(&i_index);
+    slpwvr__(&i_index);
 
     return Qnil;
 
@@ -359,7 +362,7 @@ dcl_slpvpc(obj, index, rc)
     i_rc = (real)NUM2DBL(rc);
 
 
-    slpvpc_(&i_index, &i_rc);
+    slpvpc__(&i_index, &i_rc);
 
     return Qnil;
 
@@ -383,7 +386,7 @@ dcl_slpwwc(obj, index, rc)
     i_rc = (real)NUM2DBL(rc);
 
 
-    slpwwc_(&i_index, &i_rc);
+    slpwwc__(&i_index, &i_rc);
 
     return Qnil;
 
@@ -407,7 +410,7 @@ dcl_slpwvc(obj, index, rc)
     i_rc = (real)NUM2DBL(rc);
 
 
-    slpwvc_(&i_index, &i_rc);
+    slpwvc__(&i_index, &i_rc);
 
     return Qnil;
 

data/ext/numru/grph1_stpack.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 *);
@@ -69,7 +72,7 @@ dcl_stqtrf(obj)
     logical o_lmapa;
     VALUE lmapa;
 
-    stqtrf_(&o_lmapa);
+    stqtrf__(&o_lmapa);
 
     lmapa = (o_lmapa == FALSE_) ? Qfalse : Qtrue;
 
@@ -87,7 +90,7 @@ dcl_ststrf(obj, lmapa)
     i_lmapa = ((lmapa == Qnil)||(lmapa == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    ststrf_(&i_lmapa);
+    ststrf__(&i_lmapa);
 
     return Qnil;
 
@@ -115,7 +118,7 @@ dcl_stftrn(obj, ux, uy)
     i_uy = (real)NUM2DBL(uy);
 
 
-    stftrn_(&i_ux, &i_uy, &o_vx, &o_vy);
+    stftrn__(&i_ux, &i_uy, &o_vx, &o_vy);
 
     vx = rb_float_new((double)o_vx);
     vy = rb_float_new((double)o_vy);
@@ -147,7 +150,7 @@ dcl_stitrn(obj, vx, vy)
     i_vy = (real)NUM2DBL(vy);
 
 
-    stitrn_(&i_vx, &i_vy, &o_ux, &o_uy);
+    stitrn__(&i_vx, &i_vy, &o_ux, &o_uy);
 
     ux = rb_float_new((double)o_ux);
     uy = rb_float_new((double)o_uy);
@@ -172,7 +175,7 @@ dcl_ststri(obj, itr)
     i_itr = NUM2INT(itr);
 
 
-    ststri_(&i_itr);
+    ststri__(&i_itr);
 
     return Qnil;
 
@@ -210,7 +213,7 @@ dcl_ststrp(obj, cxa, cya, vxoff, vyoff)
     i_vyoff = (real)NUM2DBL(vyoff);
 
 
-    ststrp_(&i_cxa, &i_cya, &i_vxoff, &i_vyoff);
+    ststrp__(&i_cxa, &i_cya, &i_vxoff, &i_vyoff);
 
     return Qnil;
 
@@ -251,7 +254,7 @@ dcl_ststrn(obj, itr, cxa, cya, vxoff, vyoff)
     i_vyoff = (real)NUM2DBL(vyoff);
 
 
-    ststrn_(&i_itr, &i_cxa, &i_cya, &i_vxoff, &i_vyoff);
+    ststrn__(&i_itr, &i_cxa, &i_cya, &i_vxoff, &i_vyoff);
 
     return Qnil;
 
@@ -279,7 +282,7 @@ dcl_stfrot(obj, ux, uy)
     i_uy = (real)NUM2DBL(uy);
 
 
-    stfrot_(&i_ux, &i_uy, &o_tx, &o_ty);
+    stfrot__(&i_ux, &i_uy, &o_tx, &o_ty);
 
     tx = rb_float_new((double)o_tx);
     ty = rb_float_new((double)o_ty);
@@ -311,7 +314,7 @@ dcl_stirot(obj, tx, ty)
     i_ty = (real)NUM2DBL(ty);
 
 
-    stirot_(&i_tx, &i_ty, &o_ux, &o_uy);
+    stirot__(&i_tx, &i_ty, &o_ux, &o_uy);
 
     ux = rb_float_new((double)o_ux);
     uy = rb_float_new((double)o_uy);
@@ -344,7 +347,7 @@ dcl_stsrot(obj, theta, phi, psi)
     i_psi = (real)NUM2DBL(psi);
 
 
-    stsrot_(&i_theta, &i_phi, &i_psi);
+    stsrot__(&i_theta, &i_phi, &i_psi);
 
     return Qnil;
 
@@ -372,7 +375,7 @@ dcl_stfrad(obj, x, y)
     i_y = (real)NUM2DBL(y);
 
 
-    stfrad_(&i_x, &i_y, &o_rx, &o_ry);
+    stfrad__(&i_x, &i_y, &o_rx, &o_ry);
 
     rx = rb_float_new((double)o_rx);
     ry = rb_float_new((double)o_ry);
@@ -404,7 +407,7 @@ dcl_stirad(obj, rx, ry)
     i_ry = (real)NUM2DBL(ry);
 
 
-    stirad_(&i_rx, &i_ry, &o_x, &o_y);
+    stirad__(&i_rx, &i_ry, &o_x, &o_y);
 
     x = rb_float_new((double)o_x);
     y = rb_float_new((double)o_y);
@@ -425,7 +428,7 @@ dcl_stsrad(obj, lxdeg, lydeg)
     i_lydeg = ((lydeg == Qnil)||(lydeg == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    stsrad_(&i_lxdeg, &i_lydeg);
+    stsrad__(&i_lxdeg, &i_lydeg);
 
     return Qnil;
 
@@ -455,7 +458,7 @@ dcl_stfusr(obj, xlon, ylat)
     i_ylat = (real)NUM2DBL(ylat);
 
 
-    stfusr_(&i_xlon, &i_ylat, &o_xx, &o_yy);
+    stfusr__(&i_xlon, &i_ylat, &o_xx, &o_yy);
 
     xx = rb_float_new((double)o_xx);
     yy = rb_float_new((double)o_yy);
@@ -489,7 +492,7 @@ dcl_stfusr(obj, ux, uy)
     i_uy = (real)NUM2DBL(uy);
 
 
-    stfusr_(&i_ux, &i_uy, &o_xx, &o_yy);
+    stfusr__(&i_ux, &i_uy, &o_xx, &o_yy);
 
     xx = rb_float_new((double)o_xx);
     yy = rb_float_new((double)o_yy);
@@ -525,7 +528,7 @@ dcl_stiusr(obj, xx, yy)
     i_yy = (real)NUM2DBL(yy);
 
 
-    stiusr_(&i_xx, &i_yy, &o_xlon, &o_ylat);
+    stiusr__(&i_xx, &i_yy, &o_xlon, &o_ylat);
 
     xlon = rb_float_new((double)o_xlon);
     ylat = rb_float_new((double)o_ylat);
@@ -559,7 +562,7 @@ dcl_stiusr(obj, xx, yy)
     i_yy = (real)NUM2DBL(yy);
 
 
-    stiusr_(&i_xx, &i_yy, &o_ux, &o_uy);
+    stiusr__(&i_xx, &i_yy, &o_ux, &o_uy);
 
     ux = rb_float_new((double)o_ux);
     uy = rb_float_new((double)o_uy);
@@ -575,7 +578,7 @@ static VALUE
 dcl_stsusr(obj)
     VALUE obj;
 {
-    stsusr_();
+    stsusr__();
 
     return Qnil;
 
@@ -596,7 +599,7 @@ dcl_stnusr(obj, kproj)
     i_kproj = NUM2INT(kproj);
 
 
-    stnusr_(&i_kproj);
+    stnusr__(&i_kproj);
 
     return Qnil;
 
@@ -623,7 +626,7 @@ dcl_stcusr(obj, i)
     o_cttl= ALLOCA_N(char, (DFLT_SIZE+1));
     memset(o_cttl, '\0', DFLT_SIZE+1);
 
-    stcusr_(&i_i, o_cttl, (ftnlen)DFLT_SIZE);
+    stcusr__(&i_i, o_cttl, (ftnlen)DFLT_SIZE);
 
     cttl = rb_str_new2(o_cttl);
 
@@ -656,7 +659,7 @@ dcl_stfwtr(obj, rx, ry)
     i_ry = (real)NUM2DBL(ry);
 
 
-    stfwtr_(&i_rx, &i_ry, &o_wx, &o_wy);
+    stfwtr__(&i_rx, &i_ry, &o_wx, &o_wy);
 
     wx = rb_float_new((double)o_wx);
     wy = rb_float_new((double)o_wy);
@@ -688,7 +691,7 @@ dcl_stiwtr(obj, wx, wy)
     i_wy = (real)NUM2DBL(wy);
 
 
-    stiwtr_(&i_wx, &i_wy, &o_rx, &o_ry);
+    stiwtr__(&i_wx, &i_wy, &o_rx, &o_ry);
 
     rx = rb_float_new((double)o_rx);
     ry = rb_float_new((double)o_ry);
@@ -751,7 +754,7 @@ dcl_stswtr(obj, rxmin, rxmax, rymin, rymax, wxmin, wxmax, wymin, wymax, iwtrf)
     i_iwtrf = NUM2INT(iwtrf);
 
 
-    stswtr_(&i_rxmin, &i_rxmax, &i_rymin, &i_rymax, &i_wxmin, &i_wxmax, &i_wymin, &i_wymax, &i_iwtrf);
+    stswtr__(&i_rxmin, &i_rxmax, &i_rymin, &i_rymax, &i_wxmin, &i_wxmax, &i_wymin, &i_wymax, &i_iwtrf);
 
     return Qnil;
 
@@ -780,7 +783,7 @@ dcl_stqwtr(obj)
     VALUE wymax;
     VALUE iwtrf;
 
-    stqwtr_(&o_rxmin, &o_rxmax, &o_rymin, &o_rymax, &o_wxmin, &o_wxmax, &o_wymin, &o_wymax, &o_iwtrf);
+    stqwtr__(&o_rxmin, &o_rxmax, &o_rymin, &o_rymax, &o_wxmin, &o_wxmax, &o_wymin, &o_wymax, &o_iwtrf);
 
     rxmin = rb_float_new((double)o_rxmin);
     rxmax = rb_float_new((double)o_rxmax);
@@ -825,7 +828,7 @@ dcl_stswrc(obj, wsxmn, wsxmx, wsymn, wsymx)
     i_wsymx = (real)NUM2DBL(wsymx);
 
 
-    stswrc_(&i_wsxmn, &i_wsxmx, &i_wsymn, &i_wsymx);
+    stswrc__(&i_wsxmn, &i_wsxmx, &i_wsymn, &i_wsymx);
 
     return Qnil;
 
@@ -844,7 +847,7 @@ dcl_stqwrc(obj)
     VALUE wsymn;
     VALUE wsymx;
 
-    stqwrc_(&o_wsxmn, &o_wsxmx, &o_wsymn, &o_wsymx);
+    stqwrc__(&o_wsxmn, &o_wsxmx, &o_wsymn, &o_wsymx);
 
     wsxmn = rb_float_new((double)o_wsxmn);
     wsxmx = rb_float_new((double)o_wsxmx);
@@ -883,7 +886,7 @@ dcl_stfpr3(obj, x, y, z)
     i_z = (real)NUM2DBL(z);
 
 
-    stfpr3_(&i_x, &i_y, &i_z, &o_rx, &o_ry);
+    stfpr3__(&i_x, &i_y, &i_z, &o_rx, &o_ry);
 
     rx = rb_float_new((double)o_rx);
     ry = rb_float_new((double)o_ry);
@@ -951,7 +954,7 @@ dcl_stspr3(obj, xfc, yfc, zfc, theta, phi, psi, fac, zview, rxoff, ryoff)
     i_ryoff = (real)NUM2DBL(ryoff);
 
 
-    stspr3_(&i_xfc, &i_yfc, &i_zfc, &i_theta, &i_phi, &i_psi, &i_fac, &i_zview, &i_rxoff, &i_ryoff);
+    stspr3__(&i_xfc, &i_yfc, &i_zfc, &i_theta, &i_phi, &i_psi, &i_fac, &i_zview, &i_rxoff, &i_ryoff);
 
     return Qnil;
 
@@ -980,7 +983,7 @@ dcl_stspr2(obj, ix, iy, sect)
     i_sect = (real)NUM2DBL(sect);
 
 
-    stspr2_(&i_ix, &i_iy, &i_sect);
+    stspr2__(&i_ix, &i_iy, &i_sect);
 
     return Qnil;
 
@@ -990,7 +993,7 @@ static VALUE
 dcl_stepr2(obj)
     VALUE obj;
 {
-    stepr2_();
+    stepr2__();
 
     return Qnil;
 
@@ -1000,7 +1003,7 @@ static VALUE
 dcl_strpr2(obj)
     VALUE obj;
 {
-    strpr2_();
+    strpr2__();
 
     return Qnil;
 
@@ -1035,7 +1038,7 @@ dcl_stftr3(obj, ux, uy, uz)
     i_uz = (real)NUM2DBL(uz);
 
 
-    stftr3_(&i_ux, &i_uy, &i_uz, &o_vx, &o_vy, &o_vz);
+    stftr3__(&i_ux, &i_uy, &i_uz, &o_vx, &o_vy, &o_vz);
 
     vx = rb_float_new((double)o_vx);
     vy = rb_float_new((double)o_vy);
@@ -1089,7 +1092,7 @@ dcl_ststr3(obj, itr, cxa, cya, cza, vxoff, vyoff, vzoff)
     i_vzoff = (real)NUM2DBL(vzoff);
 
 
-    ststr3_(&i_itr, &i_cxa, &i_cya, &i_cza, &i_vxoff, &i_vyoff, &i_vzoff);
+    ststr3__(&i_itr, &i_cxa, &i_cya, &i_cza, &i_vxoff, &i_vyoff, &i_vzoff);
 
     return Qnil;
 
@@ -1108,7 +1111,7 @@ dcl_stsrd3(obj, lxrd, lyrd, lzrd)
     i_lzrd = ((lzrd == Qnil)||(lzrd == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    stsrd3_(&i_lxrd, &i_lyrd, &i_lzrd);
+    stsrd3__(&i_lxrd, &i_lyrd, &i_lzrd);
 
     return Qnil;
 
@@ -1127,7 +1130,7 @@ dcl_stslg3(obj, lxlg, lylg, lzlg)
     i_lzlg = ((lzlg == Qnil)||(lzlg == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    stslg3_(&i_lxlg, &i_lylg, &i_lzlg);
+    stslg3__(&i_lxlg, &i_lylg, &i_lzlg);
 
     return Qnil;