ruby-dcl 1.8.1 → 1.8.2

data/ext/numru/math1_vialib.c

@@ -4,7 +4,6 @@
 
 #include <stdio.h>
 #include "ruby.h"
-#include "libtinyf2c.h"
 #include "narray.h"
 #include "cdcl.h"
 
@@ -106,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);
@@ -158,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);
@@ -210,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);
@@ -262,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);
@@ -314,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);
@@ -366,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);
@@ -418,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);
@@ -470,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);
@@ -522,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);
@@ -574,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);
@@ -626,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);
@@ -678,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);
@@ -725,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);
@@ -769,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);
@@ -813,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);
@@ -857,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);
@@ -901,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);
@@ -945,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);
@@ -989,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);
@@ -1033,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);
@@ -1077,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);
@@ -1121,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

@@ -4,7 +4,6 @@
 
 #include <stdio.h>
 #include "ruby.h"
-#include "libtinyf2c.h"
 #include "narray.h"
 #include "cdcl.h"
 
@@ -120,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);
@@ -187,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);
@@ -254,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);
@@ -316,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);
@@ -378,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);
@@ -440,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);
@@ -502,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);
@@ -564,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);
@@ -626,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);
@@ -688,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);
@@ -750,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);
@@ -812,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);
@@ -874,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);
@@ -936,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);
@@ -998,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

@@ -4,7 +4,6 @@
 
 #include <stdio.h>
 #include "ruby.h"
-#include "libtinyf2c.h"
 #include "narray.h"
 #include "cdcl.h"
 
@@ -106,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);
@@ -158,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);
@@ -210,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);
@@ -262,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);
@@ -314,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);
@@ -366,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);
@@ -418,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);
@@ -470,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);
@@ -522,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);
@@ -574,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);
@@ -626,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);
@@ -678,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);
@@ -725,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);
@@ -769,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);
@@ -813,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);
@@ -857,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);
@@ -901,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);
@@ -945,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);
@@ -989,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);
@@ -1033,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);
@@ -1077,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);
@@ -1121,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);

data/ext/numru/math1_vrblib.c

@@ -4,7 +4,6 @@
 
 #include <stdio.h>
 #include "ruby.h"
-#include "libtinyf2c.h"
 #include "narray.h"
 #include "cdcl.h"
 
@@ -120,7 +119,7 @@ dcl_vrfnb(obj, rx, ry, n, jx, jy, jz, rfnb)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrfnb_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz, &i_rfnb);
+    vrfnb__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz, &i_rfnb);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -187,7 +186,7 @@ dcl_vrfnb0(obj, rx, ry, n, jx, jy, jz, rfnb)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrfnb0_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz, &i_rfnb);
+    vrfnb0__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz, &i_rfnb);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -254,7 +253,7 @@ dcl_vrfnb1(obj, rx, ry, n, jx, jy, jz, rfnb)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrfnb1_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz, &i_rfnb);
+    vrfnb1__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz, &i_rfnb);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -316,7 +315,7 @@ dcl_vradd(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vradd_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vradd__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -378,7 +377,7 @@ dcl_vradd0(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vradd0_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vradd0__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -440,7 +439,7 @@ dcl_vradd1(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vradd1_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vradd1__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -502,7 +501,7 @@ dcl_vrsub(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrsub_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vrsub__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -564,7 +563,7 @@ dcl_vrsub0(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrsub0_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vrsub0__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -626,7 +625,7 @@ dcl_vrsub1(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrsub1_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vrsub1__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -688,7 +687,7 @@ dcl_vrmlt(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrmlt_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vrmlt__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -750,7 +749,7 @@ dcl_vrmlt0(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrmlt0_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vrmlt0__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -812,7 +811,7 @@ dcl_vrmlt1(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrmlt1_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vrmlt1__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -874,7 +873,7 @@ dcl_vrdiv(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrdiv_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vrdiv__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -936,7 +935,7 @@ dcl_vrdiv0(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrdiv0_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vrdiv0__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);
@@ -998,7 +997,7 @@ dcl_vrdiv1(obj, rx, ry, n, jx, jy, jz)
 
     o_rz= ALLOCA_N(real, i_jz*(i_n-1)+1);
 
-    vrdiv1_(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
+    vrdiv1__(i_rx, i_ry, o_rz, &i_n, &i_jx, &i_jy, &i_jz);
 
     {int array_shape[1] = {i_jz*(i_n-1)+1};
      rz = dcl_crealary2obj(o_rz, i_jz*(i_n-1)+1, 1, array_shape);