lbfgsb 0.2.0 → 0.3.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 272abcfabc883d993bd0b9eae08b4e0d590524b4a5aea38bc7b1f8c92fb0a5b5
-  data.tar.gz: 315ffe16183e0a08d236a9563cc6daa29db0931310b8d0cd2ccbd10810b2353e
+  metadata.gz: 38b3fbdad3f3976d24e542cdb86a2a8c2c4fdadf85948526c51e513dc2165343
+  data.tar.gz: 3a42da751483bb5b8d3940208b47da081388ba39c44c9cfc9b3eec0f41b6ec67
 SHA512:
-  metadata.gz: 73b05d103b55055894fec49d587c78a798b94d7cc3333324a5e1331e1dd8ce30bb0ee0280f190e111cc04e3d31c3c16e06fce51c91fe004a7f59bcc05aa8f1c5
-  data.tar.gz: 0ebcc71411a6fd7b9339d7c9a335277ba69eb74cb83f301de3c651bc243d1daef396b84e9904a680511d7f5f07ac1ff4af469d344f68b2ad2b41ff003b2c5534
+  metadata.gz: 00ca5c403a6871799f94540e91fb8ee8cbdf921ab81510408c6c75465f5db7ea9b911048b78c56933bac7ada91a6264b5f704f6f805ba126bc25c572e7327020
+  data.tar.gz: 36e5bfa35fc9bf58bfb5cb008c5c2efa4460a9fc76119d48d7dd9f92f0e4d61b851281e8a752a4dbb89b41bbbb2e8abdcc186f55366d3f115b2c861dad2617f8

data/CHANGELOG.md

@@ -1,3 +1,7 @@
+## 0.3.0
+- Add library name prefix to functions derived from blas and linpack
+to avoid errors due to collision with the blas functions called by numo-linalg.
+
 ## 0.2.0
 - Add an option to minimize method for calculating both function value and gradient vector using fnc.
 If true is given to jcb, fnc is assumed to return the function value and gardient vector as [f, g] array.

data/ext/lbfgsb/lbfgsbext.c

@@ -109,6 +109,7 @@ VALUE lbfgsb_min_l_bfgs_b(VALUE self,
   wa = ALLOC_N(double, (2 * m + 5) * n + 12 * m * m + 12 * m);
   iwa = ALLOC_N(long, 3 * n);
 
+  g_val = Qnil;
   f = 0.0;
   for (i = 0; i < n; g[i++] = 0.0);
   strcpy(task, "START");

data/ext/lbfgsb/src/blas.c

@@ -5,37 +5,7 @@
  */
 #include "blas.h"
 
-double dnrm2_(long *n, double *x, long *incx)
-{
-  long i__1, i__2;
-  double ret_val, d__1, d__2, d__3;
-  static long i__;
-  static double scale;
-
-  --x;
-
-  ret_val = 0.;
-  scale = 0.;
-  i__1 = *n;
-  i__2 = *incx;
-  for (i__ = 1; i__2 < 0 ? i__ >= i__1 : i__ <= i__1; i__ += i__2) {
-    d__2 = scale, d__3 = (d__1 = x[i__], fabs(d__1));
-    scale = d__2 >= d__3 ? d__2 : d__3;
-  }
-  if (scale == 0.) {
-    return ret_val;
-  }
-  i__2 = *n;
-  i__1 = *incx;
-  for (i__ = 1; i__1 < 0 ? i__ >= i__2 : i__ <= i__2; i__ += i__1) {
-    d__1 = x[i__] / scale;
-    ret_val += d__1 * d__1;
-  }
-  ret_val = scale * sqrt(ret_val);
-  return ret_val;
-}
-
-int daxpy_(long *n, double *da, double *dx, long *incx, double *dy, long *incy)
+int lbfgsb_rb_daxpy_(long *n, double *da, double *dx, long *incx, double *dy, long *incy)
 {
   long i__1;
   static long i__, m, ix, iy, mp1;
@@ -100,7 +70,7 @@ L40:
   return 0;
 }
 
-int dcopy_(long *n, double *dx, long *incx, double *dy, long *incy)
+int lbfgsb_rb_dcopy_(long *n, double *dx, long *incx, double *dy, long *incy)
 {
   long i__1;
   static long i__, m, ix, iy, mp1;
@@ -165,7 +135,7 @@ L40:
   return 0;
 }
 
-double ddot_(long *n, double *dx, long *incx, double *dy, long *incy)
+double lbfgsb_rb_ddot_(long *n, double *dx, long *incx, double *dy, long *incy)
 {
   long i__1;
   double ret_val;
@@ -232,7 +202,7 @@ L60:
   return ret_val;
 }
 
-int dscal_(long *n, double *da, double *dx, long *incx)
+int lbfgsb_rb_dscal_(long *n, double *da, double *dx, long *incx)
 {
   long i__1, i__2;
   static long i__, m, mp1, nincx;

data/ext/lbfgsb/src/blas.h

@@ -3,10 +3,9 @@
 
 #include <math.h>
 
-extern double dnrm2_(long *n, double *x, long *incx);
-extern int daxpy_(long *n, double *da, double *dx, long *incx, double *dy, long *incy);
-extern int dcopy_(long *n, double *dx, long *incx, double *dy, long *incy);
-extern double ddot_(long *n, double *dx, long *incx, double *dy, long *incy);
-extern int dscal_(long *n, double *da, double *dx, long *incx);
+extern int lbfgsb_rb_daxpy_(long *n, double *da, double *dx, long *incx, double *dy, long *incy);
+extern int lbfgsb_rb_dcopy_(long *n, double *dx, long *incx, double *dy, long *incy);
+extern double lbfgsb_rb_ddot_(long *n, double *dx, long *incx, double *dy, long *incy);
+extern int lbfgsb_rb_dscal_(long *n, double *da, double *dx, long *incx);
 
 #endif /* LBFGSB_RB_BLAS_H_ */

data/ext/lbfgsb/src/lbfgsb.c

@@ -673,8 +673,8 @@ int mainlb_(long *n, long *m, double *x,
     if (strncmp(task, "STOP", 4) == 0) {
       if (strncmp(task + 6, "CPU", 3) == 0) {
         /* restore the previous iterate. */
-        dcopy_(n, &t[1], &c__1, &x[1], &c__1);
-        dcopy_(n, &r__[1], &c__1, &g[1], &c__1);
+        lbfgsb_rb_dcopy_(n, &t[1], &c__1, &x[1], &c__1);
+        lbfgsb_rb_dcopy_(n, &r__[1], &c__1, &g[1], &c__1);
         *f = fold;
       }
       goto L999;
@@ -709,7 +709,7 @@ L222:
 
   if (! cnstnd && col > 0) {
     /* skip the search for GCP. */
-    dcopy_(n, &x[1], &c__1, &z__[1], &c__1);
+    lbfgsb_rb_dcopy_(n, &x[1], &c__1, &z__[1], &c__1);
     wrk = updatd;
     nseg = 0;
     goto L333;
@@ -834,8 +834,8 @@ L666:
       csave, &isave[22], &dsave[17]);
   if (info != 0 || iback >= 20) {
     /* restore the previous iterate. */
-    dcopy_(n, &t[1], &c__1, &x[1], &c__1);
-    dcopy_(n, &r__[1], &c__1, &g[1], &c__1);
+    lbfgsb_rb_dcopy_(n, &t[1], &c__1, &x[1], &c__1);
+    lbfgsb_rb_dcopy_(n, &r__[1], &c__1, &g[1], &c__1);
     *f = fold;
     if (col == 0) {
       /* abnormal termination. */
@@ -910,13 +910,13 @@ L777:
   for (i__ = 1; i__ <= i__1; ++i__) {
     r__[i__] = g[i__] - r__[i__];
   }
-  rr = ddot_(n, &r__[1], &c__1, &r__[1], &c__1);
+  rr = lbfgsb_rb_ddot_(n, &r__[1], &c__1, &r__[1], &c__1);
   if (stp == 1.) {
     dr = gd - gdold;
     ddum = -gdold;
   } else {
     dr = (gd - gdold) * stp;
-    dscal_(n, &stp, &d__[1], &c__1);
+    lbfgsb_rb_dscal_(n, &stp, &d__[1], &c__1);
     ddum = -gdold * stp;
   }
   if (dr <= epsmch * ddum) {
@@ -1193,7 +1193,7 @@ int bmv_(long *m, double *sy, double *wt, long
     p[i2] = v[i2] + sum;
   }
   /* Solve the triangular system */
-  dtrsl_(&wt[wt_offset], m, col, &p[*col + 1], &c__11, info);
+  lbfgsb_rb_dtrsl_(&wt[wt_offset], m, col, &p[*col + 1], &c__11, info);
   if (*info != 0) {
     return 0;
   }
@@ -1205,7 +1205,7 @@ int bmv_(long *m, double *sy, double *wt, long
   /* PART II: solve [ -D^(1/2)   D^(-1/2)*L'  ] [ p1 ] = [ p1 ] */
   /*                [  0         J'           ] [ p2 ]   [ p2 ]. */
   /*   solve J^Tp2=p2. */
-  dtrsl_(&wt[wt_offset], m, col, &p[*col + 1], &c__1, info);
+  lbfgsb_rb_dtrsl_(&wt[wt_offset], m, col, &p[*col + 1], &c__1, info);
   if (*info != 0) {
     return 0;
   }
@@ -1463,7 +1463,7 @@ int cauchy_(long *n, double *x, double *l,
     if (*iprint >= 0) {
       fprintf(stdout, " Subgnorm = 0.  GCP = X.\n");
     }
-    dcopy_(n, &x[1], &c__1, &xcp[1], &c__1);
+    lbfgsb_rb_dcopy_(n, &x[1], &c__1, &xcp[1], &c__1);
     return 0;
   }
   bnded = TRUE_;
@@ -1562,10 +1562,10 @@ int cauchy_(long *n, double *x, double *l,
   /*   The smallest of the nbreak breakpoints is in t(ibkmin)=bkmin. */
   if (*theta != 1.) {
     /* complete the initialization of p for theta not= one. */
-    dscal_(col, theta, &p[*col + 1], &c__1);
+    lbfgsb_rb_dscal_(col, theta, &p[*col + 1], &c__1);
   }
   /* Initialize GCP xcp = x. */
-  dcopy_(n, &x[1], &c__1, &xcp[1], &c__1);
+  lbfgsb_rb_dcopy_(n, &x[1], &c__1, &xcp[1], &c__1);
   if (nbreak == 0 && nfree == *n + 1) {
     /* is a zero vector, return with the initial xcp as GCP. */
     if (*iprint > 100) {
@@ -1596,7 +1596,7 @@ int cauchy_(long *n, double *x, double *l,
     if (*info != 0) {
       return 0;
     }
-    f2 -= ddot_(&col2, &v[1], &c__1, &p[1], &c__1);
+    f2 -= lbfgsb_rb_ddot_(&col2, &v[1], &c__1, &p[1], &c__1);
   }
   dtm = -f1 / f2;
   tsum = 0.;
@@ -1685,7 +1685,7 @@ L777:
   f2 -= *theta * dibp2;
   if (*col > 0) {
     /* update c = c + dt*p. */
-    daxpy_(&col2, &dt, &p[1], &c__1, &c__[1], &c__1);
+    lbfgsb_rb_daxpy_(&col2, &dt, &p[1], &c__1, &c__[1], &c__1);
     /* choose wbp, */
     /* the row of W corresponding to the breakpoint encountered. */
     pointr = *head;
@@ -1700,12 +1700,12 @@ L777:
     if (*info != 0) {
       return 0;
     }
-    wmc = ddot_(&col2, &c__[1], &c__1, &v[1], &c__1);
-    wmp = ddot_(&col2, &p[1], &c__1, &v[1], &c__1);
-    wmw = ddot_(&col2, &wbp[1], &c__1, &v[1], &c__1);
+    wmc = lbfgsb_rb_ddot_(&col2, &c__[1], &c__1, &v[1], &c__1);
+    wmp = lbfgsb_rb_ddot_(&col2, &p[1], &c__1, &v[1], &c__1);
+    wmw = lbfgsb_rb_ddot_(&col2, &wbp[1], &c__1, &v[1], &c__1);
     /* update p = p - dibp*wbp. */
     d__1 = -dibp;
-    daxpy_(&col2, &d__1, &wbp[1], &c__1, &p[1], &c__1);
+    lbfgsb_rb_daxpy_(&col2, &d__1, &wbp[1], &c__1, &p[1], &c__1);
     /* complete updating f1 and f2 while col > 0. */
     f1 += dibp * wmc;
     f2 = f2 + dibp * 2. * wmp - dibp2 * wmw;
@@ -1737,12 +1737,12 @@ L888:
   tsum += dtm;
   /* Move free variables (i.e., the ones w/o breakpoints) and */
   /*   the variables whose breakpoints haven't been reached. */
-  daxpy_(n, &tsum, &d__[1], &c__1, &xcp[1], &c__1);
+  lbfgsb_rb_daxpy_(n, &tsum, &d__[1], &c__1, &xcp[1], &c__1);
 L999:
   /* Update c = c + dtm*p = W'(x^c - x) */
   /*   which will be used in computing r = Z'(B(x^c - x) + g). */
   if (*col > 0) {
-    daxpy_(&col2, &dtm, &p[1], &c__1, &c__[1], &c__1);
+    lbfgsb_rb_daxpy_(&col2, &dtm, &p[1], &c__1, &c__[1], &c__1);
   }
   if (*iprint > 100) {
     fprintf(stdout, "Cauchy X =  \n");
@@ -2060,11 +2060,11 @@ int formk_(long *n, long *nsub, long *ind, long *nenter,
       for (jy = 1; jy <= i__1; ++jy) {
         js = *m + jy;
         i__2 = *m - jy;
-        dcopy_(&i__2, &wn1[jy + 1 + (jy + 1) * wn1_dim1], &c__1, &wn1[jy + jy * wn1_dim1], &c__1);
+        lbfgsb_rb_dcopy_(&i__2, &wn1[jy + 1 + (jy + 1) * wn1_dim1], &c__1, &wn1[jy + jy * wn1_dim1], &c__1);
         i__2 = *m - jy;
-        dcopy_(&i__2, &wn1[js + 1 + (js + 1) * wn1_dim1], &c__1, &wn1[js + js * wn1_dim1], &c__1);
+        lbfgsb_rb_dcopy_(&i__2, &wn1[js + 1 + (js + 1) * wn1_dim1], &c__1, &wn1[js + js * wn1_dim1], &c__1);
         i__2 = *m - 1;
-        dcopy_(&i__2, &wn1[*m + 2 + (jy + 1) * wn1_dim1], &c__1, &wn1[*m + 1 + jy * wn1_dim1], &c__1);
+        lbfgsb_rb_dcopy_(&i__2, &wn1[*m + 2 + (jy + 1) * wn1_dim1], &c__1, &wn1[*m + 1 + jy * wn1_dim1], &c__1);
       }
     }
     /* put new rows in blocks (1,1), (2,1) and (2,2). */
@@ -2215,7 +2215,7 @@ int formk_(long *n, long *nsub, long *ind, long *nenter,
   /*                                [(-L_a +R_z)L'^-1   S'AA'S*theta  ] */
   /*    first Cholesky factor (1,1) block of wn to get LL' */
   /*                      with L' stored in the upper triangle of wn. */
-  dpofa_(&wn[wn_offset], &m2, col, info);
+  lbfgsb_rb_dpofa_(&wn[wn_offset], &m2, col, info);
   if (*info != 0) {
     *info = -1;
     return 0;
@@ -2224,7 +2224,7 @@ int formk_(long *n, long *nsub, long *ind, long *nenter,
   col2 = *col << 1;
   i__1 = col2;
   for (js = *col + 1; js <= i__1; ++js) {
-    dtrsl_(&wn[wn_offset], &m2, col, &wn[js * wn_dim1 + 1], &c__11, info);
+    lbfgsb_rb_dtrsl_(&wn[wn_offset], &m2, col, &wn[js * wn_dim1 + 1], &c__11, info);
   }
   /* Form S'AA'S*theta + (L^-1(-L_a'+R_z'))'L^-1(-L_a'+R_z') in the */
   /*    upper triangle of (2,2) block of wn. */
@@ -2232,11 +2232,11 @@ int formk_(long *n, long *nsub, long *ind, long *nenter,
   for (is = *col + 1; is <= i__1; ++is) {
     i__2 = col2;
     for (js = is; js <= i__2; ++js) {
-      wn[is + js * wn_dim1] += ddot_(col, &wn[is * wn_dim1 + 1], &c__1, &wn[js * wn_dim1 + 1], &c__1);
+      wn[is + js * wn_dim1] += lbfgsb_rb_ddot_(col, &wn[is * wn_dim1 + 1], &c__1, &wn[js * wn_dim1 + 1], &c__1);
     }
   }
   /* Cholesky factorization of (2,2) block of wn. */
-  dpofa_(&wn[*col + 1 + (*col + 1) * wn_dim1], &m2, col, info);
+  lbfgsb_rb_dpofa_(&wn[*col + 1 + (*col + 1) * wn_dim1], &m2, col, info);
   if (*info != 0) {
     *info = -2;
     return 0;
@@ -2303,7 +2303,7 @@ int formt_(long *m, double *wt, double *sy,
   }
   /* Cholesky factorize T to J*J' with */
   /*    J' stored in the upper triangle of wt. */
-  dpofa_(&wt[wt_offset], m, col, info);
+  lbfgsb_rb_dpofa_(&wt[wt_offset], m, col, info);
   if (*info != 0) {
     *info = -3;
   }
@@ -2563,7 +2563,7 @@ int lnsrlb_(long *n, double *l, double *u,
   if (strncmp(task, "FG_LN", 5) == 0) {
     goto L556;
   }
-  *dtd = ddot_(n, &d__[1], &c__1, &d__[1], &c__1);
+  *dtd = lbfgsb_rb_ddot_(n, &d__[1], &c__1, &d__[1], &c__1);
   *dnorm = sqrt(*dtd);
   /* Determine the maximum step length. */
   *stpmx = 1e10;
@@ -2600,14 +2600,14 @@ int lnsrlb_(long *n, double *l, double *u,
   } else {
     *stp = 1.;
   }
-  dcopy_(n, &x[1], &c__1, &t[1], &c__1);
-  dcopy_(n, &g[1], &c__1, &r__[1], &c__1);
+  lbfgsb_rb_dcopy_(n, &x[1], &c__1, &t[1], &c__1);
+  lbfgsb_rb_dcopy_(n, &g[1], &c__1, &r__[1], &c__1);
   *fold = *f;
   *ifun = 0;
   *iback = 0;
   strcpy(csave, "START");
 L556:
-  *gd = ddot_(n, &g[1], &c__1, &d__[1], &c__1);
+  *gd = lbfgsb_rb_ddot_(n, &g[1], &c__1, &d__[1], &c__1);
   if (*ifun == 0) {
     *gdold = *gd;
     if (*gd >= 0.) {
@@ -2626,7 +2626,7 @@ L556:
     ++(*nfgv);
     *iback = *ifun - 1;
     if (*stp == 1.) {
-      dcopy_(n, &z__[1], &c__1, &x[1], &c__1);
+      lbfgsb_rb_dcopy_(n, &z__[1], &c__1, &x[1], &c__1);
     } else {
       i__1 = *n;
       for (i__ = 1; i__ <= i__1; ++i__) {
@@ -2693,8 +2693,8 @@ int matupd_(long *n, long *m, double *ws,
     *head = *head % *m + 1;
   }
   /* Update matrices WS and WY. */
-  dcopy_(n, &d__[1], &c__1, &ws[*itail * ws_dim1 + 1], &c__1);
-  dcopy_(n, &r__[1], &c__1, &wy[*itail * wy_dim1 + 1], &c__1);
+  lbfgsb_rb_dcopy_(n, &d__[1], &c__1, &ws[*itail * ws_dim1 + 1], &c__1);
+  lbfgsb_rb_dcopy_(n, &r__[1], &c__1, &wy[*itail * wy_dim1 + 1], &c__1);
   /* Set theta=yy/ys. */
   *theta = *rr / *dr;
   /* Form the middle matrix in B. */
@@ -2704,9 +2704,9 @@ int matupd_(long *n, long *m, double *ws,
     /* move old information */
     i__1 = *col - 1;
     for (j = 1; j <= i__1; ++j) {
-      dcopy_(&j, &ss[(j + 1) * ss_dim1 + 2], &c__1, &ss[j * ss_dim1 + 1], &c__1);
+      lbfgsb_rb_dcopy_(&j, &ss[(j + 1) * ss_dim1 + 2], &c__1, &ss[j * ss_dim1 + 1], &c__1);
       i__2 = *col - j;
-      dcopy_(&i__2, &sy[j + 1 + (j + 1) * sy_dim1], &c__1, &sy[j + j * sy_dim1], &c__1);
+      lbfgsb_rb_dcopy_(&i__2, &sy[j + 1 + (j + 1) * sy_dim1], &c__1, &sy[j + j * sy_dim1], &c__1);
     }
   }
   /* add new information: the last row of SY */
@@ -2714,8 +2714,8 @@ int matupd_(long *n, long *m, double *ws,
   pointr = *head;
   i__1 = *col - 1;
   for (j = 1; j <= i__1; ++j) {
-    sy[*col + j * sy_dim1] = ddot_(n, &d__[1], &c__1, &wy[pointr * wy_dim1 + 1], &c__1);
-    ss[j + *col * ss_dim1] = ddot_(n, &ws[pointr * ws_dim1 + 1], &c__1, &d__[1], &c__1);
+    sy[*col + j * sy_dim1] = lbfgsb_rb_ddot_(n, &d__[1], &c__1, &wy[pointr * wy_dim1 + 1], &c__1);
+    ss[j + *col * ss_dim1] = lbfgsb_rb_ddot_(n, &ws[pointr * ws_dim1 + 1], &c__1, &d__[1], &c__1);
     pointr = pointr % *m + 1;
   }
   if (*stp == 1.) {
@@ -3361,7 +3361,7 @@ int subsm_(long *n, long *m, long *nsub, long *ind,
   /* Compute wv:=K^(-1)wv. */
   m2 = *m << 1;
   col2 = *col << 1;
-  dtrsl_(&wn[wn_offset], &m2, &col2, &wv[1], &c__11, info);
+  lbfgsb_rb_dtrsl_(&wn[wn_offset], &m2, &col2, &wv[1], &c__11, info);
   if (*info != 0) {
     return 0;
   }
@@ -3369,7 +3369,7 @@ int subsm_(long *n, long *m, long *nsub, long *ind,
   for (i__ = 1; i__ <= i__1; ++i__) {
     wv[i__] = -wv[i__];
   }
-  dtrsl_(&wn[wn_offset], &m2, &col2, &wv[1], &c__1, info);
+  lbfgsb_rb_dtrsl_(&wn[wn_offset], &m2, &col2, &wv[1], &c__1, info);
   if (*info != 0) {
     return 0;
   }
@@ -3387,12 +3387,12 @@ int subsm_(long *n, long *m, long *nsub, long *ind,
     pointr = pointr % *m + 1;
   }
   d__1 = 1. / *theta;
-  dscal_(nsub, &d__1, &d__[1], &c__1);
+  lbfgsb_rb_dscal_(nsub, &d__1, &d__[1], &c__1);
 
   /* ----------------------------------------------------------------- */
   /* Let us try the projection, d is the Newton direction */
   *iword = 0;
-  dcopy_(n, &x[1], &c__1, &xp[1], &c__1);
+  lbfgsb_rb_dcopy_(n, &x[1], &c__1, &xp[1], &c__1);
 
   i__1 = *nsub;
   for (i__ = 1; i__ <= i__1; ++i__) {
@@ -3445,7 +3445,7 @@ int subsm_(long *n, long *m, long *nsub, long *ind,
     dd_p__ += (x[i__] - xx[i__]) * gg[i__];
   }
   if (dd_p__ > 0.) {
-    dcopy_(n, &xp[1], &c__1, &x[1], &c__1);
+    lbfgsb_rb_dcopy_(n, &xp[1], &c__1, &x[1], &c__1);
     fprintf(stderr, "  Positive dir derivative in projection\n");
     fprintf(stderr, "  Using the backtracking step\n");
   } else {

data/ext/lbfgsb/src/linpack.c

@@ -44,7 +44,7 @@ static long c__1 = 1;
  *     linpack.  this version dated 08/14/78 .
  *     cleve moler, university of new mexico, argonne national lab.
  */
-int dpofa_(double *a, long *lda, long *n, long *info)
+int lbfgsb_rb_dpofa_(double *a, long *lda, long *n, long *info)
 {
   long a_dim1, a_offset, i__1, i__2, i__3;
   static long j, k;
@@ -66,7 +66,7 @@ int dpofa_(double *a, long *lda, long *n, long *info)
 	  i__2 = jm1;
 	  for (k = 1; k <= i__2; ++k) {
 	    i__3 = k - 1;
-	    t = a[k + j * a_dim1] - ddot_(&i__3, &a[k * a_dim1 + 1], &c__1, &a[j * a_dim1 + 1], &c__1);
+	    t = a[k + j * a_dim1] - lbfgsb_rb_ddot_(&i__3, &a[k * a_dim1 + 1], &c__1, &a[j * a_dim1 + 1], &c__1);
 	    t /= a[k + k * a_dim1];
 	    a[k + j * a_dim1] = t;
 	    s += t * t;
@@ -133,7 +133,7 @@ L40:
  *     linpack. this version dated 08/14/78 .
  *     g. w. stewart, university of maryland, argonne national lab.
  */
-int dtrsl_(double *t, long *ldt, long *n, double *b, long *job, long *info)
+int lbfgsb_rb_dtrsl_(double *t, long *ldt, long *n, double *b, long *job, long *info)
 {
   long t_dim1, t_offset, i__1, i__2;
   static long j, jj, case__;
@@ -179,7 +179,7 @@ L20:
   for (j = 2; j <= i__1; ++j) {
 	  temp = -b[j - 1];
 	  i__2 = *n - j + 1;
-	  daxpy_(&i__2, &temp, &t[j + (j - 1) * t_dim1], &c__1, &b[j], &c__1);
+	  lbfgsb_rb_daxpy_(&i__2, &temp, &t[j + (j - 1) * t_dim1], &c__1, &b[j], &c__1);
 	  b[j] /= t[j + j * t_dim1];
   }
 L40:
@@ -195,7 +195,7 @@ L50:
   for (jj = 2; jj <= i__1; ++jj) {
 	  j = *n - jj + 1;
 	  temp = -b[j + 1];
-	  daxpy_(&j, &temp, &t[(j + 1) * t_dim1 + 1], &c__1, &b[1], &c__1);
+	  lbfgsb_rb_daxpy_(&j, &temp, &t[(j + 1) * t_dim1 + 1], &c__1, &b[1], &c__1);
 	  b[j] /= t[j + j * t_dim1];
   }
 L70:
@@ -211,7 +211,7 @@ L80:
   for (jj = 2; jj <= i__1; ++jj) {
 	  j = *n - jj + 1;
 	  i__2 = jj - 1;
-	  b[j] -= ddot_(&i__2, &t[j + 1 + j * t_dim1], &c__1, &b[j + 1], &c__1);
+	  b[j] -= lbfgsb_rb_ddot_(&i__2, &t[j + 1 + j * t_dim1], &c__1, &b[j + 1], &c__1);
 	  b[j] /= t[j + j * t_dim1];
   }
 L100:
@@ -226,7 +226,7 @@ L110:
   i__1 = *n;
   for (j = 2; j <= i__1; ++j) {
 	  i__2 = j - 1;
-	  b[j] -= ddot_(&i__2, &t[j * t_dim1 + 1], &c__1, &b[1], &c__1);
+	  b[j] -= lbfgsb_rb_ddot_(&i__2, &t[j * t_dim1 + 1], &c__1, &b[1], &c__1);
 	  b[j] /= t[j + j * t_dim1];
   }
 L130:

data/ext/lbfgsb/src/linpack.h

@@ -3,7 +3,7 @@
 
 #include <math.h>
 
-extern int dpofa_(double *a, long *lda, long *n, long *info);
-extern int dtrsl_(double *t, long *ldt, long *n, double *b, long *job, long *info);
+extern int lbfgsb_rb_dpofa_(double *a, long *lda, long *n, long *info);
+extern int lbfgsb_rb_dtrsl_(double *t, long *ldt, long *n, double *b, long *job, long *info);
 
 #endif /* LBFGSB_RB_LINPACK_H_ */

data/lib/lbfgsb/version.rb

@@ -3,5 +3,5 @@
 # Lbfgsb.rb is a Ruby binding for L-BFGS-B with Numo::NArray.
 module Lbfgsb
   # The version of Lbfgsb.rb you are using.
-  VERSION = '0.2.0'
+  VERSION = '0.3.0'
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: lbfgsb
 version: !ruby/object:Gem::Version
-  version: 0.2.0
+  version: 0.3.0
 platform: ruby
 authors:
 - yoshoku
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2020-11-14 00:00:00.000000000 Z
+date: 2020-11-22 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: numo-narray
@@ -77,7 +77,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.1.2
+rubygems_version: 3.1.4
 signing_key:
 specification_version: 4
 summary: Lbfgsb.rb is a Ruby binding for L-BFGS-B