ruby-dcl 1.7.1 → 1.8.2

data/ext/numru/grph1_scpack.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_scsobj(obj, xobj3, yobj3, zobj3)
     i_zobj3 = (real)NUM2DBL(zobj3);
 
 
-    scsobj_(&i_xobj3, &i_yobj3, &i_zobj3);
+    scsobj__(&i_xobj3, &i_yobj3, &i_zobj3);
 
     return Qnil;
 
@@ -102,7 +105,7 @@ dcl_scqobj(obj)
     VALUE yobj3;
     VALUE zobj3;
 
-    scqobj_(&o_xobj3, &o_yobj3, &o_zobj3);
+    scqobj__(&o_xobj3, &o_yobj3, &o_zobj3);
 
     xobj3 = rb_float_new((double)o_xobj3);
     yobj3 = rb_float_new((double)o_yobj3);
@@ -136,7 +139,7 @@ dcl_scseye(obj, xeye3, yeye3, zeye3)
     i_zeye3 = (real)NUM2DBL(zeye3);
 
 
-    scseye_(&i_xeye3, &i_yeye3, &i_zeye3);
+    scseye__(&i_xeye3, &i_yeye3, &i_zeye3);
 
     return Qnil;
 
@@ -153,7 +156,7 @@ dcl_scqeye(obj)
     VALUE yeye3;
     VALUE zeye3;
 
-    scqeye_(&o_xeye3, &o_yeye3, &o_zeye3);
+    scqeye__(&o_xeye3, &o_yeye3, &o_zeye3);
 
     xeye3 = rb_float_new((double)o_xeye3);
     yeye3 = rb_float_new((double)o_yeye3);
@@ -187,7 +190,7 @@ dcl_scspln(obj, ixax, iyax, sect)
     i_sect = (real)NUM2DBL(sect);
 
 
-    scspln_(&i_ixax, &i_iyax, &i_sect);
+    scspln__(&i_ixax, &i_iyax, &i_sect);
 
     return Qnil;
 
@@ -204,7 +207,7 @@ dcl_scqpln(obj)
     VALUE iyax;
     VALUE sect;
 
-    scqpln_(&o_ixax, &o_iyax, &o_sect);
+    scqpln__(&o_ixax, &o_iyax, &o_sect);
 
     ixax = INT2NUM(o_ixax);
     iyax = INT2NUM(o_iyax);
@@ -219,7 +222,7 @@ static VALUE
 dcl_scsprj(obj)
     VALUE obj;
 {
-    scsprj_();
+    scsprj__();
 
     return Qnil;
 
@@ -263,7 +266,7 @@ dcl_scsvpt(obj, vxmin, vxmax, vymin, vymax, vzmin, vzmax)
     i_vzmax = (real)NUM2DBL(vzmax);
 
 
-    scsvpt_(&i_vxmin, &i_vxmax, &i_vymin, &i_vymax, &i_vzmin, &i_vzmax);
+    scsvpt__(&i_vxmin, &i_vxmax, &i_vymin, &i_vymax, &i_vzmin, &i_vzmax);
 
     return Qnil;
 
@@ -286,7 +289,7 @@ dcl_scqvpt(obj)
     VALUE vzmin;
     VALUE vzmax;
 
-    scqvpt_(&o_vxmin, &o_vxmax, &o_vymin, &o_vymax, &o_vzmin, &o_vzmax);
+    scqvpt__(&o_vxmin, &o_vxmax, &o_vymin, &o_vymax, &o_vzmin, &o_vzmax);
 
     vxmin = rb_float_new((double)o_vxmin);
     vxmax = rb_float_new((double)o_vxmax);
@@ -338,7 +341,7 @@ dcl_scswnd(obj, uxmin, uxmax, uymin, uymax, uzmin, uzmax)
     i_uzmax = (real)NUM2DBL(uzmax);
 
 
-    scswnd_(&i_uxmin, &i_uxmax, &i_uymin, &i_uymax, &i_uzmin, &i_uzmax);
+    scswnd__(&i_uxmin, &i_uxmax, &i_uymin, &i_uymax, &i_uzmin, &i_uzmax);
 
     return Qnil;
 
@@ -361,7 +364,7 @@ dcl_scqwnd(obj)
     VALUE uzmin;
     VALUE uzmax;
 
-    scqwnd_(&o_uxmin, &o_uxmax, &o_uymin, &o_uymax, &o_uzmin, &o_uzmax);
+    scqwnd__(&o_uxmin, &o_uxmax, &o_uymin, &o_uymax, &o_uzmin, &o_uzmax);
 
     uxmin = rb_float_new((double)o_uxmin);
     uxmax = rb_float_new((double)o_uxmax);
@@ -388,7 +391,7 @@ dcl_scslog(obj, lxlog3, lylog3, lzlog3)
     i_lzlog3 = ((lzlog3 == Qnil)||(lzlog3 == Qfalse)) ? FALSE_ : TRUE_;
 
 
-    scslog_(&i_lxlog3, &i_lylog3, &i_lzlog3);
+    scslog__(&i_lxlog3, &i_lylog3, &i_lzlog3);
 
     return Qnil;
 
@@ -405,7 +408,7 @@ dcl_scqlog(obj)
     VALUE lylog3;
     VALUE lzlog3;
 
-    scqlog_(&o_lxlog3, &o_lylog3, &o_lzlog3);
+    scqlog__(&o_lxlog3, &o_lylog3, &o_lzlog3);
 
     lxlog3 = (o_lxlog3 == FALSE_) ? Qfalse : Qtrue;
     lylog3 = (o_lylog3 == FALSE_) ? Qfalse : Qtrue;
@@ -444,7 +447,7 @@ dcl_scsorg(obj, simfac, vxorg3, vyorg3, vzorg3)
     i_vzorg3 = (real)NUM2DBL(vzorg3);
 
 
-    scsorg_(&i_simfac, &i_vxorg3, &i_vyorg3, &i_vzorg3);
+    scsorg__(&i_simfac, &i_vxorg3, &i_vyorg3, &i_vzorg3);
 
     return Qnil;
 
@@ -463,7 +466,7 @@ dcl_scqorg(obj)
     VALUE vyorg3;
     VALUE vzorg3;
 
-    scqorg_(&o_simfac, &o_vxorg3, &o_vyorg3, &o_vzorg3);
+    scqorg__(&o_simfac, &o_vxorg3, &o_vyorg3, &o_vzorg3);
 
     simfac = rb_float_new((double)o_simfac);
     vxorg3 = rb_float_new((double)o_vxorg3);
@@ -488,7 +491,7 @@ dcl_scstrn(obj, itr3)
     i_itr3 = NUM2INT(itr3);
 
 
-    scstrn_(&i_itr3);
+    scstrn__(&i_itr3);
 
     return Qnil;
 
@@ -501,7 +504,7 @@ dcl_scqtrn(obj)
     integer o_itr3;
     VALUE itr3;
 
-    scqtrn_(&o_itr3);
+    scqtrn__(&o_itr3);
 
     itr3 = INT2NUM(o_itr3);
 
@@ -514,7 +517,7 @@ static VALUE
 dcl_scstrf(obj)
     VALUE obj;
 {
-    scstrf_();
+    scstrf__();
 
     return Qnil;
 
@@ -560,7 +563,7 @@ dcl_scplu(obj, n, upx, upy, upz)
     i_upz = dcl_obj2crealary(upz);
 
 
-    scplu_(&i_n, i_upx, i_upy, i_upz);
+    scplu__(&i_n, i_upx, i_upy, i_upz);
 
     dcl_freecrealary(i_upx);
     dcl_freecrealary(i_upy);
@@ -610,7 +613,7 @@ dcl_scplv(obj, n, vpx, vpy, vpz)
     i_vpz = dcl_obj2crealary(vpz);
 
 
-    scplv_(&i_n, i_vpx, i_vpy, i_vpz);
+    scplv__(&i_n, i_vpx, i_vpy, i_vpz);
 
     dcl_freecrealary(i_vpx);
     dcl_freecrealary(i_vpy);
@@ -633,7 +636,7 @@ dcl_scspli(obj, index)
     i_index = NUM2INT(index);
 
 
-    scspli_(&i_index);
+    scspli__(&i_index);
 
     return Qnil;
 
@@ -646,7 +649,7 @@ dcl_scqpli(obj)
     integer o_index;
     VALUE index;
 
-    scqpli_(&o_index);
+    scqpli__(&o_index);
 
     index = INT2NUM(o_index);
 
@@ -700,7 +703,7 @@ dcl_scplzu(obj, n, upx, upy, upz, index)
     i_upz = dcl_obj2crealary(upz);
 
 
-    scplzu_(&i_n, i_upx, i_upy, i_upz, &i_index);
+    scplzu__(&i_n, i_upx, i_upy, i_upz, &i_index);
 
     dcl_freecrealary(i_upx);
     dcl_freecrealary(i_upy);
@@ -755,7 +758,7 @@ dcl_scplzv(obj, n, vpx, vpy, vpz, index)
     i_vpz = dcl_obj2crealary(vpz);
 
 
-    scplzv_(&i_n, i_vpx, i_vpy, i_vpz, &i_index);
+    scplzv__(&i_n, i_vpx, i_vpy, i_vpz, &i_index);
 
     dcl_freecrealary(i_vpx);
     dcl_freecrealary(i_vpy);
@@ -805,7 +808,7 @@ dcl_scpmu(obj, n, upx, upy, upz)
     i_upz = dcl_obj2crealary(upz);
 
 
-    scpmu_(&i_n, i_upx, i_upy, i_upz);
+    scpmu__(&i_n, i_upx, i_upy, i_upz);
 
     dcl_freecrealary(i_upx);
     dcl_freecrealary(i_upy);
@@ -855,7 +858,7 @@ dcl_scpmv(obj, n, vpx, vpy, vpz)
     i_vpz = dcl_obj2crealary(vpz);
 
 
-    scpmv_(&i_n, i_vpx, i_vpy, i_vpz);
+    scpmv__(&i_n, i_vpx, i_vpy, i_vpz);
 
     dcl_freecrealary(i_vpx);
     dcl_freecrealary(i_vpy);
@@ -878,7 +881,7 @@ dcl_scspmt(obj, itype)
     i_itype = NUM2INT(itype);
 
 
-    scspmt_(&i_itype);
+    scspmt__(&i_itype);
 
     return Qnil;
 
@@ -891,7 +894,7 @@ dcl_scqpmt(obj)
     integer o_itype;
     VALUE itype;
 
-    scqpmt_(&o_itype);
+    scqpmt__(&o_itype);
 
     itype = INT2NUM(o_itype);
 
@@ -913,7 +916,7 @@ dcl_scspmi(obj, index)
     i_index = NUM2INT(index);
 
 
-    scspmi_(&i_index);
+    scspmi__(&i_index);
 
     return Qnil;
 
@@ -926,7 +929,7 @@ dcl_scqpmi(obj)
     integer o_index;
     VALUE index;
 
-    scqpmi_(&o_index);
+    scqpmi__(&o_index);
 
     index = INT2NUM(o_index);
 
@@ -948,7 +951,7 @@ dcl_scspms(obj, rsize)
     i_rsize = (real)NUM2DBL(rsize);
 
 
-    scspms_(&i_rsize);
+    scspms__(&i_rsize);
 
     return Qnil;
 
@@ -961,7 +964,7 @@ dcl_scqpms(obj)
     real o_rsize;
     VALUE rsize;
 
-    scqpms_(&o_rsize);
+    scqpms__(&o_rsize);
 
     rsize = rb_float_new((double)o_rsize);
 
@@ -1025,7 +1028,7 @@ dcl_scpmzu(obj, n, upx, upy, upz, itype, index, rsize)
     i_upz = dcl_obj2crealary(upz);
 
 
-    scpmzu_(&i_n, i_upx, i_upy, i_upz, &i_itype, &i_index, &i_rsize);
+    scpmzu__(&i_n, i_upx, i_upy, i_upz, &i_itype, &i_index, &i_rsize);
 
     dcl_freecrealary(i_upx);
     dcl_freecrealary(i_upy);
@@ -1090,7 +1093,7 @@ dcl_scpmzv(obj, n, vpx, vpy, vpz, itype, index, rsize)
     i_vpz = dcl_obj2crealary(vpz);
 
 
-    scpmzv_(&i_n, i_vpx, i_vpy, i_vpz, &i_itype, &i_index, &i_rsize);
+    scpmzv__(&i_n, i_vpx, i_vpy, i_vpz, &i_itype, &i_index, &i_rsize);
 
     dcl_freecrealary(i_vpx);
     dcl_freecrealary(i_vpy);
@@ -1135,7 +1138,7 @@ dcl_sctnu(obj, upx, upy, upz)
     i_upz = dcl_obj2crealary(upz);
 
 
-    sctnu_(i_upx, i_upy, i_upz);
+    sctnu__(i_upx, i_upy, i_upz);
 
     dcl_freecrealary(i_upx);
     dcl_freecrealary(i_upy);
@@ -1180,7 +1183,7 @@ dcl_sctnv(obj, vpx, vpy, vpz)
     i_vpz = dcl_obj2crealary(vpz);
 
 
-    sctnv_(i_vpx, i_vpy, i_vpz);
+    sctnv__(i_vpx, i_vpy, i_vpz);
 
     dcl_freecrealary(i_vpx);
     dcl_freecrealary(i_vpy);
@@ -1208,7 +1211,7 @@ dcl_scstnp(obj, itpat1, itpat2)
     i_itpat2 = NUM2INT(itpat2);
 
 
-    scstnp_(&i_itpat1, &i_itpat2);
+    scstnp__(&i_itpat1, &i_itpat2);
 
     return Qnil;
 
@@ -1223,7 +1226,7 @@ dcl_scqtnp(obj)
     VALUE itpat1;
     VALUE itpat2;
 
-    scqtnp_(&o_itpat1, &o_itpat2);
+    scqtnp__(&o_itpat1, &o_itpat2);
 
     itpat1 = INT2NUM(o_itpat1);
     itpat2 = INT2NUM(o_itpat2);
@@ -1278,7 +1281,7 @@ dcl_sctnzu(obj, upx, upy, upz, itpat1, itpat2)
     i_upz = dcl_obj2crealary(upz);
 
 
-    sctnzu_(i_upx, i_upy, i_upz, &i_itpat1, &i_itpat2);
+    sctnzu__(i_upx, i_upy, i_upz, &i_itpat1, &i_itpat2);
 
     dcl_freecrealary(i_upx);
     dcl_freecrealary(i_upy);
@@ -1333,7 +1336,7 @@ dcl_sctnzv(obj, vpx, vpy, vpz, itpat1, itpat2)
     i_vpz = dcl_obj2crealary(vpz);
 
 
-    sctnzv_(i_vpx, i_vpy, i_vpz, &i_itpat1, &i_itpat2);
+    sctnzv__(i_vpx, i_vpy, i_vpz, &i_itpat1, &i_itpat2);
 
     dcl_freecrealary(i_vpx);
     dcl_freecrealary(i_vpy);