ruby-dcl 1.7.1 → 1.8.2

data/ext/numru/math1_vialib.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 *);
@@ -102,7 +105,7 @@ dcl_vifna(obj, ix, n, jx, jy, ifna)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    vifna_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ifna);
+    vifna__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ifna);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -154,7 +157,7 @@ dcl_vifna0(obj, ix, n, jx, jy, ifna)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    vifna0_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ifna);
+    vifna0__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ifna);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -206,7 +209,7 @@ dcl_vifna1(obj, ix, n, jx, jy, ifna)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    vifna1_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ifna);
+    vifna1__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ifna);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -258,7 +261,7 @@ dcl_viinc(obj, ix, n, jx, jy, ii)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    viinc_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
+    viinc__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -310,7 +313,7 @@ dcl_viinc0(obj, ix, n, jx, jy, ii)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    viinc0_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
+    viinc0__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -362,7 +365,7 @@ dcl_viinc1(obj, ix, n, jx, jy, ii)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    viinc1_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
+    viinc1__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -414,7 +417,7 @@ dcl_vifct(obj, ix, n, jx, jy, ii)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    vifct_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
+    vifct__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -466,7 +469,7 @@ dcl_vifct0(obj, ix, n, jx, jy, ii)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    vifct0_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
+    vifct0__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -518,7 +521,7 @@ dcl_vifct1(obj, ix, n, jx, jy, ii)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    vifct1_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
+    vifct1__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -570,7 +573,7 @@ dcl_vicon(obj, ix, n, jx, jy, ii)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    vicon_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
+    vicon__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -622,7 +625,7 @@ dcl_vicon0(obj, ix, n, jx, jy, ii)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    vicon0_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
+    vicon0__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -674,7 +677,7 @@ dcl_vicon1(obj, ix, n, jx, jy, ii)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    vicon1_(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
+    vicon1__(i_ix, o_iy, &i_n, &i_jx, &i_jy, &i_ii);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -721,7 +724,7 @@ dcl_viset(obj, ix, n, jx, jy)
 
     o_iy= ALLOCA_N(integer, i_jy*(i_n-1)+1);
 
-    viset_(i_ix, o_iy, &i_n, &i_jx, &i_jy);
+    viset__(i_ix, o_iy, &i_n, &i_jx, &i_jy);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      iy = dcl_cintegerary2obj(o_iy, i_jy*(i_n-1)+1, 1, array_shape);
@@ -765,7 +768,7 @@ dcl_iadd(obj, ix, n, jx, ii)
     io_ix = dcl_obj2cintegerary(ix);
 
 
-    iadd_(io_ix, &i_n, &i_jx, &i_ii);
+    iadd__(io_ix, &i_n, &i_jx, &i_ii);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      ix = dcl_cintegerary2obj(io_ix, i_jx*(i_n-1)+1, 1, array_shape);
@@ -809,7 +812,7 @@ dcl_iadd0(obj, ix, n, jx, ii)
     io_ix = dcl_obj2cintegerary(ix);
 
 
-    iadd0_(io_ix, &i_n, &i_jx, &i_ii);
+    iadd0__(io_ix, &i_n, &i_jx, &i_ii);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      ix = dcl_cintegerary2obj(io_ix, i_jx*(i_n-1)+1, 1, array_shape);
@@ -853,7 +856,7 @@ dcl_iadd1(obj, ix, n, jx, ii)
     io_ix = dcl_obj2cintegerary(ix);
 
 
-    iadd1_(io_ix, &i_n, &i_jx, &i_ii);
+    iadd1__(io_ix, &i_n, &i_jx, &i_ii);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      ix = dcl_cintegerary2obj(io_ix, i_jx*(i_n-1)+1, 1, array_shape);
@@ -897,7 +900,7 @@ dcl_imlt(obj, ix, n, jx, ii)
     io_ix = dcl_obj2cintegerary(ix);
 
 
-    imlt_(io_ix, &i_n, &i_jx, &i_ii);
+    imlt__(io_ix, &i_n, &i_jx, &i_ii);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      ix = dcl_cintegerary2obj(io_ix, i_jx*(i_n-1)+1, 1, array_shape);
@@ -941,7 +944,7 @@ dcl_imlt0(obj, ix, n, jx, ii)
     io_ix = dcl_obj2cintegerary(ix);
 
 
-    imlt0_(io_ix, &i_n, &i_jx, &i_ii);
+    imlt0__(io_ix, &i_n, &i_jx, &i_ii);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      ix = dcl_cintegerary2obj(io_ix, i_jx*(i_n-1)+1, 1, array_shape);
@@ -985,7 +988,7 @@ dcl_imlt1(obj, ix, n, jx, ii)
     io_ix = dcl_obj2cintegerary(ix);
 
 
-    imlt1_(io_ix, &i_n, &i_jx, &i_ii);
+    imlt1__(io_ix, &i_n, &i_jx, &i_ii);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      ix = dcl_cintegerary2obj(io_ix, i_jx*(i_n-1)+1, 1, array_shape);
@@ -1029,7 +1032,7 @@ dcl_iset(obj, ix, n, jx, ii)
     io_ix = dcl_obj2cintegerary(ix);
 
 
-    iset_(io_ix, &i_n, &i_jx, &i_ii);
+    iset__(io_ix, &i_n, &i_jx, &i_ii);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      ix = dcl_cintegerary2obj(io_ix, i_jx*(i_n-1)+1, 1, array_shape);
@@ -1073,7 +1076,7 @@ dcl_iset0(obj, ix, n, jx, ii)
     io_ix = dcl_obj2cintegerary(ix);
 
 
-    iset0_(io_ix, &i_n, &i_jx, &i_ii);
+    iset0__(io_ix, &i_n, &i_jx, &i_ii);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      ix = dcl_cintegerary2obj(io_ix, i_jx*(i_n-1)+1, 1, array_shape);
@@ -1117,7 +1120,7 @@ dcl_iset1(obj, ix, n, jx, ii)
     io_ix = dcl_obj2cintegerary(ix);
 
 
-    iset1_(io_ix, &i_n, &i_jx, &i_ii);
+    iset1__(io_ix, &i_n, &i_jx, &i_ii);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      ix = dcl_cintegerary2obj(io_ix, i_jx*(i_n-1)+1, 1, array_shape);

data/ext/numru/math1_viblib.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 *);
@@ -116,7 +119,7 @@ dcl_vifnb(obj, ix, iy, n, jx, jy, jz, ifnb)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    vifnb_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz, &i_ifnb);
+    vifnb__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz, &i_ifnb);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -183,7 +186,7 @@ dcl_vifnb0(obj, ix, iy, n, jx, jy, jz, ifnb)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    vifnb0_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz, &i_ifnb);
+    vifnb0__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz, &i_ifnb);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -250,7 +253,7 @@ dcl_vifnb1(obj, ix, iy, n, jx, jy, jz, ifnb)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    vifnb1_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz, &i_ifnb);
+    vifnb1__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz, &i_ifnb);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -312,7 +315,7 @@ dcl_viadd(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    viadd_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    viadd__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -374,7 +377,7 @@ dcl_viadd0(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    viadd0_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    viadd0__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -436,7 +439,7 @@ dcl_viadd1(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    viadd1_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    viadd1__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -498,7 +501,7 @@ dcl_visub(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    visub_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    visub__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -560,7 +563,7 @@ dcl_visub0(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    visub0_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    visub0__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -622,7 +625,7 @@ dcl_visub1(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    visub1_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    visub1__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -684,7 +687,7 @@ dcl_vimlt(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    vimlt_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    vimlt__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -746,7 +749,7 @@ dcl_vimlt0(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    vimlt0_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    vimlt0__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -808,7 +811,7 @@ dcl_vimlt1(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    vimlt1_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    vimlt1__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -870,7 +873,7 @@ dcl_vidiv(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    vidiv_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    vidiv__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -932,7 +935,7 @@ dcl_vidiv0(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    vidiv0_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    vidiv0__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -994,7 +997,7 @@ dcl_vidiv1(obj, ix, iy, n, jx, jy, jz)
 
     o_iz= ALLOCA_N(integer, i_jz*(i_n-1)+1);
 
-    vidiv1_(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
+    vidiv1__(i_ix, i_iy, o_iz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      iz = dcl_cintegerary2obj(o_iz, i_jz*(i_n-1)+1, 1, array_shape);

data/ext/numru/math1_vralib.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 *);
@@ -102,7 +105,7 @@ dcl_vrfna(obj, rx, n, jx, jy, rfna)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrfna_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rfna);
+    vrfna__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rfna);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -154,7 +157,7 @@ dcl_vrfna0(obj, rx, n, jx, jy, rfna)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrfna0_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rfna);
+    vrfna0__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rfna);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -206,7 +209,7 @@ dcl_vrfna1(obj, rx, n, jx, jy, rfna)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrfna1_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rfna);
+    vrfna1__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rfna);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -258,7 +261,7 @@ dcl_vrinc(obj, rx, n, jx, jy, rr)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrinc_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
+    vrinc__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -310,7 +313,7 @@ dcl_vrinc0(obj, rx, n, jx, jy, rr)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrinc0_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
+    vrinc0__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -362,7 +365,7 @@ dcl_vrinc1(obj, rx, n, jx, jy, rr)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrinc1_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
+    vrinc1__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -414,7 +417,7 @@ dcl_vrfct(obj, rx, n, jx, jy, rr)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrfct_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
+    vrfct__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -466,7 +469,7 @@ dcl_vrfct0(obj, rx, n, jx, jy, rr)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrfct0_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
+    vrfct0__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -518,7 +521,7 @@ dcl_vrfct1(obj, rx, n, jx, jy, rr)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrfct1_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
+    vrfct1__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -570,7 +573,7 @@ dcl_vrcon(obj, rx, n, jx, jy, rr)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrcon_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
+    vrcon__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -622,7 +625,7 @@ dcl_vrcon0(obj, rx, n, jx, jy, rr)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrcon0_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
+    vrcon0__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -674,7 +677,7 @@ dcl_vrcon1(obj, rx, n, jx, jy, rr)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrcon1_(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
+    vrcon1__(i_rx, o_ry, &i_n, &i_jx, &i_jy, &i_rr);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -721,7 +724,7 @@ dcl_vrset(obj, rx, n, jx, jy)
 
     o_ry= ALLOCA_N(real, i_jy*(i_n-1)+1);
 
-    vrset_(i_rx, o_ry, &i_n, &i_jx, &i_jy);
+    vrset__(i_rx, o_ry, &i_n, &i_jx, &i_jy);
 
     {int array_shape[1] = {i_jy*(i_n-1)+1};
      ry = dcl_crealary2obj(o_ry, i_jy*(i_n-1)+1, 1, array_shape);
@@ -765,7 +768,7 @@ dcl_radd(obj, rx, n, jx, rr)
     io_rx = dcl_obj2crealary(rx);
 
 
-    radd_(io_rx, &i_n, &i_jx, &i_rr);
+    radd__(io_rx, &i_n, &i_jx, &i_rr);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      rx = dcl_crealary2obj(io_rx, i_jx*(i_n-1)+1, 1, array_shape);
@@ -809,7 +812,7 @@ dcl_radd0(obj, rx, n, jx, rr)
     io_rx = dcl_obj2crealary(rx);
 
 
-    radd0_(io_rx, &i_n, &i_jx, &i_rr);
+    radd0__(io_rx, &i_n, &i_jx, &i_rr);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      rx = dcl_crealary2obj(io_rx, i_jx*(i_n-1)+1, 1, array_shape);
@@ -853,7 +856,7 @@ dcl_radd1(obj, rx, n, jx, rr)
     io_rx = dcl_obj2crealary(rx);
 
 
-    radd1_(io_rx, &i_n, &i_jx, &i_rr);
+    radd1__(io_rx, &i_n, &i_jx, &i_rr);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      rx = dcl_crealary2obj(io_rx, i_jx*(i_n-1)+1, 1, array_shape);
@@ -897,7 +900,7 @@ dcl_rmlt(obj, rx, n, jx, rr)
     io_rx = dcl_obj2crealary(rx);
 
 
-    rmlt_(io_rx, &i_n, &i_jx, &i_rr);
+    rmlt__(io_rx, &i_n, &i_jx, &i_rr);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      rx = dcl_crealary2obj(io_rx, i_jx*(i_n-1)+1, 1, array_shape);
@@ -941,7 +944,7 @@ dcl_rmlt0(obj, rx, n, jx, rr)
     io_rx = dcl_obj2crealary(rx);
 
 
-    rmlt0_(io_rx, &i_n, &i_jx, &i_rr);
+    rmlt0__(io_rx, &i_n, &i_jx, &i_rr);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      rx = dcl_crealary2obj(io_rx, i_jx*(i_n-1)+1, 1, array_shape);
@@ -985,7 +988,7 @@ dcl_rmlt1(obj, rx, n, jx, rr)
     io_rx = dcl_obj2crealary(rx);
 
 
-    rmlt1_(io_rx, &i_n, &i_jx, &i_rr);
+    rmlt1__(io_rx, &i_n, &i_jx, &i_rr);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      rx = dcl_crealary2obj(io_rx, i_jx*(i_n-1)+1, 1, array_shape);
@@ -1029,7 +1032,7 @@ dcl_rset(obj, rx, n, jx, rr)
     io_rx = dcl_obj2crealary(rx);
 
 
-    rset_(io_rx, &i_n, &i_jx, &i_rr);
+    rset__(io_rx, &i_n, &i_jx, &i_rr);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      rx = dcl_crealary2obj(io_rx, i_jx*(i_n-1)+1, 1, array_shape);
@@ -1073,7 +1076,7 @@ dcl_rset0(obj, rx, n, jx, rr)
     io_rx = dcl_obj2crealary(rx);
 
 
-    rset0_(io_rx, &i_n, &i_jx, &i_rr);
+    rset0__(io_rx, &i_n, &i_jx, &i_rr);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      rx = dcl_crealary2obj(io_rx, i_jx*(i_n-1)+1, 1, array_shape);
@@ -1117,7 +1120,7 @@ dcl_rset1(obj, rx, n, jx, rr)
     io_rx = dcl_obj2crealary(rx);
 
 
-    rset1_(io_rx, &i_n, &i_jx, &i_rr);
+    rset1__(io_rx, &i_n, &i_jx, &i_rr);
 
     {int array_shape[1] = {i_jx*(i_n-1)+1};
      rx = dcl_crealary2obj(io_rx, i_jx*(i_n-1)+1, 1, array_shape);