lbfgsb 0.3.1 → 0.5.0

data/ext/lbfgsb/src/lbfgsb.h

@@ -1,122 +1,83 @@
 #ifndef LBFGSB_RB_LBFGSB_H_
 #define LBFGSB_RB_LBFGSB_H_
 
-#include <stdio.h>
-#include <stdlib.h>
 #include <float.h>
 #include <math.h>
+#include <stdio.h>
+#include <stdlib.h>
 #include <string.h>
 #include <time.h>
 
+#include "common.h"
+
 #define TRUE_ (1)
 #define FALSE_ (0)
 
-extern int setulb_(long *n, long *m, double *x,
-  double *l, double *u, long *nbd, double *f, double *g,
-  double *factr, double *pgtol, double *wa, long *iwa,
-  char *task, long *iprint, char *csave, long *lsave,
-  long *isave, double *dsave);
-
-extern int mainlb_(long *n, long *m, double *x,
-  double *l, double *u, long *nbd, double *f, double *g,
-  double *factr, double *pgtol, double *ws, double *wy,
-  double *sy, double *ss, double *wt, double *wn,
-  double *snd, double *z__, double *r__, double *d__,
-  double *t, double *xp, double *wa, long *index,
-  long *iwhere, long *indx2, char *task, long *iprint,
-  char *csave, long *lsave, long *isave, double *dsave);
-
-extern int active_(long *n, double *l, double *u,
-  long *nbd, double *x, long *iwhere, long *iprint,
-  long *prjctd, long *cnstnd, long *boxed);
-
-extern int bmv_(long *m, double *sy, double *wt, long
-  *col, double *v, double *p, long *info);
-
-extern int cauchy_(long *n, double *x, double *l,
-  double *u, long *nbd, double *g, long *iorder, long *iwhere,
-  double *t, double *d__, double *xcp, long *m,
-  double *wy, double *ws, double *sy, double *wt,
-  double *theta, long *col, long *head, double *p,
-  double *c__, double *wbp, double *v, long *nseg,
-  long *iprint, double *sbgnrm, long *info, double *epsmch);
-
-extern int cmprlb_(long *n, long *m, double *x,
-  double *g, double *ws, double *wy, double *sy,
-  double *wt, double *z__, double *r__, double *wa,
-  long *index, double *theta, long *col, long *head,
-  long *nfree, long *cnstnd, long *info);
-
-extern int errclb_(long *n, long *m, double *factr,
-  double *l, double *u, long *nbd, char *task, long *info, long *k);
-
-extern int formk_(long *n, long *nsub, long *ind, long *nenter,
-  long *ileave, long *indx2, long *iupdat, long *updatd,
-  double *wn, double *wn1, long *m, double *ws,
-  double *wy, double *sy, double *theta, long *col,
-  long *head, long *info);
-
-extern int formt_(long *m, double *wt, double *sy,
-  double *ss, long *col, double *theta, long *info);
-
-extern int freev_(long *n, long *nfree, long *index,
-  long *nenter, long *ileave, long *indx2, long *iwhere,
-  long *wrk, long *updatd, long *cnstnd, long *iprint,
-  long *iter);
-
-extern int hpsolb_(long *n, double *t, long *iorder, long *iheap);
-
-extern int lnsrlb_(long *n, double *l, double *u,
-  long *nbd, double *x, double *f, double *fold,
-  double *gd, double *gdold, double *g, double *d__,
-  double *r__, double *t, double *z__, double *stp,
-  double *dnorm, double *dtd, double *xstep, double *stpmx,
-  long *iter, long *ifun, long *iback, long *nfgv,
-  long *info, char *task, long *boxed, long *cnstnd,
-  char *csave, long *isave, double *dsave);
-
-extern int matupd_(long *n, long *m, double *ws,
-  double *wy, double *sy, double *ss, double *d__,
-  double *r__, long *itail, long *iupdat, long *col,
-  long *head, double *theta, double *rr, double *dr,
-  double *stp, double *dtd);
-
-extern int prn1lb_(long *n, long *m, double *l,
-  double *u, double *x, long *iprint, long *itfile,
-  double *epsmch);
-
-extern int prn2lb_(long *n, double *x, double *f,
-  double *g, long *iprint, long *itfile, long *iter,
-  long *nfgv, long *nact, double *sbgnrm, long *nseg, char*word,
-  long *iword, long *iback, double *stp, double *xstep);
-
-extern int prn3lb_(long *n, double *x, double *f, char *task,
-  long *iprint, long *info, long *itfile, long *iter,
-  long *nfgv, long *nintol, long *nskip, long *nact,
-  double *sbgnrm, double *time, long *nseg, char *word,
-  long *iback, double *stp, double *xstep, long *k,
-  double *cachyt, double *sbtime, double *lnscht);
-
-extern int projgr_(long *n, double *l, double *u,
-  long *nbd, double *x, double *g, double *sbgnrm);
-
-extern int subsm_(long *n, long *m, long *nsub, long *ind,
-  double *l, double *u, long *nbd, double *x,
-  double *d__, double *xp, double *ws, double *wy,
-  double *theta, double *xx, double *gg, long *col,
-  long *head, long *iword, double *wv, double *wn,
-  long *iprint, long *info);
-
-extern int dcsrch_(double *f, double *g, double *stp,
-  double *ftol, double *gtol, double *xtol,
-  double *stpmin, double *stpmax,
-  char *task, long *isave, double *dsave);
-
-extern int dcstep_(double *stx, double *fx, double *dx,
-  double *sty, double *fy, double *dy, double *stp,
-  double *fp, double *dp, long *brackt, double *stpmin,
-  double *stpmax);
-
-extern int timer_(double *ttime);
+extern void setulb_(F77_int* n, F77_int* m, double* x, double* l, double* u, F77_int* nbd, double* f, double* g, double* factr,
+                   double* pgtol, double* wa, F77_int* iwa, char* task, F77_int* iprint, char* csave, F77_int* lsave, F77_int* isave,
+                   double* dsave);
+
+extern void mainlb_(F77_int* n, F77_int* m, double* x, double* l, double* u, F77_int* nbd, double* f, double* g, double* factr,
+                   double* pgtol, double* ws, double* wy, double* sy, double* ss, double* wt, double* wn, double* snd,
+                   double* z__, double* r__, double* d__, double* t, double* xp, double* wa, F77_int* index, F77_int* iwhere,
+                   F77_int* indx2, char* task, F77_int* iprint, char* csave, F77_int* lsave, F77_int* isave, double* dsave);
+
+extern void active_(F77_int* n, double* l, double* u, F77_int* nbd, double* x, F77_int* iwhere, F77_int* iprint, F77_int* prjctd, F77_int* cnstnd,
+                   F77_int* boxed);
+
+extern void bmv_(F77_int* m, double* sy, double* wt, F77_int* col, double* v, double* p, F77_int* info);
+
+extern void cauchy_(F77_int* n, double* x, double* l, double* u, F77_int* nbd, double* g, F77_int* iorder, F77_int* iwhere, double* t,
+                   double* d__, double* xcp, F77_int* m, double* wy, double* ws, double* sy, double* wt, double* theta, F77_int* col,
+                   F77_int* head, double* p, double* c__, double* wbp, double* v, F77_int* nseg, F77_int* iprint, double* sbgnrm,
+                   F77_int* info, double* epsmch);
+
+extern void cmprlb_(F77_int* n, F77_int* m, double* x, double* g, double* ws, double* wy, double* sy, double* wt, double* z__,
+                   double* r__, double* wa, F77_int* index, double* theta, F77_int* col, F77_int* head, F77_int* nfree, F77_int* cnstnd,
+                   F77_int* info);
+
+extern void errclb_(F77_int* n, F77_int* m, double* factr, double* l, double* u, F77_int* nbd, char* task, F77_int* info, F77_int* k);
+
+extern void formk_(F77_int* n, F77_int* nsub, F77_int* ind, F77_int* nenter, F77_int* ileave, F77_int* indx2, F77_int* iupdat, F77_int* updatd,
+                  double* wn, double* wn1, F77_int* m, double* ws, double* wy, double* sy, double* theta, F77_int* col, F77_int* head,
+                  F77_int* info);
+
+extern void formt_(F77_int* m, double* wt, double* sy, double* ss, F77_int* col, double* theta, F77_int* info);
+
+extern void freev_(F77_int* n, F77_int* nfree, F77_int* index, F77_int* nenter, F77_int* ileave, F77_int* indx2, F77_int* iwhere, F77_int* wrk,
+                  F77_int* updatd, F77_int* cnstnd, F77_int* iprint, F77_int* iter);
+
+extern void hpsolb_(F77_int* n, double* t, F77_int* iorder, F77_int* iheap);
+
+extern void lnsrlb_(F77_int* n, double* l, double* u, F77_int* nbd, double* x, double* f, double* fold, double* gd, double* gdold,
+                   double* g, double* d__, double* r__, double* t, double* z__, double* stp, double* dnorm, double* dtd,
+                   double* xstep, double* stpmx, F77_int* iter, F77_int* ifun, F77_int* iback, F77_int* nfgv, F77_int* info, char* task,
+                   F77_int* boxed, F77_int* cnstnd, char* csave, F77_int* isave, double* dsave);
+
+extern void matupd_(F77_int* n, F77_int* m, double* ws, double* wy, double* sy, double* ss, double* d__, double* r__, F77_int* itail,
+                   F77_int* iupdat, F77_int* col, F77_int* head, double* theta, double* rr, double* dr, double* stp, double* dtd);
+
+extern void prn1lb_(F77_int* n, F77_int* m, double* l, double* u, double* x, F77_int* iprint, F77_int* itfile, double* epsmch);
+
+extern void prn2lb_(F77_int* n, double* x, double* f, double* g, F77_int* iprint, F77_int* itfile, F77_int* iter, F77_int* nfgv, F77_int* nact,
+                   double* sbgnrm, F77_int* nseg, char* word, F77_int* iword, F77_int* iback, double* stp, double* xstep);
+
+extern void prn3lb_(F77_int* n, double* x, double* f, char* task, F77_int* iprint, F77_int* info, F77_int* itfile, F77_int* iter, F77_int* nfgv,
+                   F77_int* nintol, F77_int* nskip, F77_int* nact, double* sbgnrm, double* time, F77_int* nseg, char* word, F77_int* iback,
+                   double* stp, double* xstep, F77_int* k, double* cachyt, double* sbtime, double* lnscht);
+
+extern void projgr_(F77_int* n, double* l, double* u, F77_int* nbd, double* x, double* g, double* sbgnrm);
+
+extern void subsm_(F77_int* n, F77_int* m, F77_int* nsub, F77_int* ind, double* l, double* u, F77_int* nbd, double* x, double* d__, double* xp,
+                  double* ws, double* wy, double* theta, double* xx, double* gg, F77_int* col, F77_int* head, F77_int* iword, double* wv,
+                  double* wn, F77_int* iprint, F77_int* info);
+
+extern void dcsrch_(double* f, double* g, double* stp, double* ftol, double* gtol, double* xtol, double* stpmin, double* stpmax,
+                   char* task, F77_int* isave, double* dsave);
+
+extern void dcstep_(double* stx, double* fx, double* dx, double* sty, double* fy, double* dy, double* stp, double* fp,
+                   double* dp, F77_int* brackt, double* stpmin, double* stpmax);
+
+extern void timer_(double* ttime);
 
 #endif /* LBFGSB_RB_LBFGSB_H_ */

data/ext/lbfgsb/src/linpack.c

@@ -3,11 +3,10 @@
  * or “3-clause license”)
  *  Please read attached file License.txt
  */
-
 #include "blas.h"
 #include "linpack.h"
 
-static long c__1 = 1;
+static F77_int c__1 = 1;
 
 /**
  *     dpofa factors a double precision symmetric positive definite
@@ -23,10 +22,10 @@ static long c__1 = 1;
  *                the symmetric matrix to be factored.  only the
  *                diagonal and upper triangle are used.
  *
- *        lda     long
+ *        lda     integer
  *                the leading dimension of the array  a .
  *
- *        n       long
+ *        n       integer
  *                the order of the matrix  a .
  *
  *     on return
@@ -36,7 +35,7 @@ static long c__1 = 1;
  *                the strict lower triangle is unaltered.
  *                if  info .ne. 0 , the factorization is not complete.
  *
- *        info    long
+ *        info    integer
  *                = 0  for normal return.
  *                = k  signals an error condition.  the leading minor
  *                     of order  k  is not positive definite.
@@ -44,12 +43,11 @@ static long c__1 = 1;
  *     linpack.  this version dated 08/14/78 .
  *     cleve moler, university of new mexico, argonne national lab.
  */
-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;
+void dpofa_(double* a, F77_int* lda, F77_int* n, F77_int* info) {
+  F77_int a_dim1, a_offset, i__1, i__2, i__3;
+  static F77_int j, k;
   static double s, t;
-  static long jm1;
+  static F77_int jm1;
 
   a_dim1 = *lda;
   a_offset = 1 + a_dim1;
@@ -57,31 +55,31 @@ int lbfgsb_rb_dpofa_(double *a, long *lda, long *n, long *info)
 
   i__1 = *n;
   for (j = 1; j <= i__1; ++j) {
-	  *info = j;
-	  s = 0.;
-	  jm1 = j - 1;
-	  if (jm1 < 1) {
-	    goto L20;
-	  }
-	  i__2 = jm1;
-	  for (k = 1; k <= i__2; ++k) {
-	    i__3 = k - 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;
-	  }
-L20:
-	  s = a[j + j * a_dim1] - s;
+    *info = j;
+    s = 0.;
+    jm1 = j - 1;
+    if (jm1 < 1) {
+      goto L20;
+    }
+    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 + k * a_dim1];
+      a[k + j * a_dim1] = t;
+      s += t * t;
+    }
+  L20:
+    s = a[j + j * a_dim1] - s;
     /* ......exit */
-	  if (s <= 0.) {
-	    goto L40;
-	  }
-	  a[j + j * a_dim1] = sqrt(s);
+    if (s <= 0.) {
+      goto L40;
+    }
+    a[j + j * a_dim1] = sqrt(s);
   }
   *info = 0;
 L40:
-  return 0;
+  return;
 }
 
 /**
@@ -102,16 +100,16 @@ L40:
  *                   the corresponding elements of the array can be
  *                   used to store other information.
  *
- *         ldt       long
+ *         ldt       integer
  *                   ldt is the leading dimension of the array t.
  *
- *         n         long
+ *         n         integer
  *                   n is the order of the system.
  *
  *         b         double precision(n).
  *                   b contains the right hand side of the system.
  *
- *         job       long
+ *         job       integer
  *                   job specifies what kind of system is to be solved.
  *                   if job is
  *
@@ -125,7 +123,7 @@ L40:
  *         b         b contains the solution, if info .eq. 0.
  *                   otherwise b is unaltered.
  *
- *         info      long
+ *         info      integer
  *                   info contains zero if the system is nonsingular.
  *                   otherwise info contains the index of
  *                   the first zero diagonal element of t.
@@ -133,10 +131,9 @@ L40:
  *     linpack. this version dated 08/14/78 .
  *     g. w. stewart, university of maryland, argonne national lab.
  */
-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__;
+void dtrsl_(double* t, F77_int* ldt, F77_int* n, double* b, F77_int* job, F77_int* info) {
+  F77_int t_dim1, t_offset, i__1, i__2;
+  static F77_int j, jj, case__;
   static double temp;
 
   /* check for zero diagonal elements. */
@@ -148,39 +145,43 @@ int lbfgsb_rb_dtrsl_(double *t, long *ldt, long *n, double *b, long *job, long *
   i__1 = *n;
   for (*info = 1; *info <= i__1; ++(*info)) {
     /* ......exit */
-	  if (t[*info + *info * t_dim1] == 0.) {
-	    goto L150;
-	  }
+    if (t[*info + *info * t_dim1] == 0.) {
+      goto L150;
+    }
   }
   *info = 0;
 
   /* determine the task and go to it. */
   case__ = 1;
   if (*job % 10 != 0) {
-	  case__ = 2;
+    case__ = 2;
   }
   if (*job % 100 / 10 != 0) {
-	  case__ += 2;
+    case__ += 2;
   }
   switch (case__) {
-	  case 1: goto L20;
-	  case 2: goto L50;
-	  case 3: goto L80;
-	  case 4: goto L110;
+  case 1:
+    goto L20;
+  case 2:
+    goto L50;
+  case 3:
+    goto L80;
+  case 4:
+    goto L110;
   }
 
   /* solve t*x=b for t lower triangular */
 L20:
   b[1] /= t[t_dim1 + 1];
   if (*n < 2) {
-	  goto L40;
+    goto L40;
   }
   i__1 = *n;
   for (j = 2; j <= i__1; ++j) {
-	  temp = -b[j - 1];
-	  i__2 = *n - j + 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];
+    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);
+    b[j] /= t[j + j * t_dim1];
   }
 L40:
   goto L140;
@@ -189,14 +190,14 @@ L40:
 L50:
   b[*n] /= t[*n + *n * t_dim1];
   if (*n < 2) {
-	  goto L70;
+    goto L70;
   }
   i__1 = *n;
   for (jj = 2; jj <= i__1; ++jj) {
-	  j = *n - jj + 1;
-	  temp = -b[j + 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];
+    j = *n - jj + 1;
+    temp = -b[j + 1];
+    daxpy_(&j, &temp, &t[(j + 1) * t_dim1 + 1], &c__1, &b[1], &c__1);
+    b[j] /= t[j + j * t_dim1];
   }
 L70:
   goto L140;
@@ -205,14 +206,14 @@ L70:
 L80:
   b[*n] /= t[*n + *n * t_dim1];
   if (*n < 2) {
-	  goto L100;
+    goto L100;
   }
   i__1 = *n;
   for (jj = 2; jj <= i__1; ++jj) {
-	  j = *n - jj + 1;
-	  i__2 = jj - 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];
+    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] /= t[j + j * t_dim1];
   }
 L100:
   goto L140;
@@ -221,16 +222,16 @@ L100:
 L110:
   b[1] /= t[t_dim1 + 1];
   if (*n < 2) {
-	  goto L130;
+    goto L130;
   }
   i__1 = *n;
   for (j = 2; j <= i__1; ++j) {
-	  i__2 = j - 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];
+    i__2 = j - 1;
+    b[j] -= ddot_(&i__2, &t[j * t_dim1 + 1], &c__1, &b[1], &c__1);
+    b[j] /= t[j + j * t_dim1];
   }
 L130:
 L140:
 L150:
-  return 0;
+  return;
 }

data/ext/lbfgsb/src/linpack.h

@@ -1,9 +1,9 @@
 #ifndef LBFGSB_RB_LINPACK_H_
 #define LBFGSB_RB_LINPACK_H_
 
-#include <math.h>
+#include "common.h"
 
-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);
+extern void dpofa_(double* a, F77_int* lda, F77_int* n, F77_int* info);
+extern void dtrsl_(double* t, F77_int* ldt, F77_int* n, double* b, F77_int* job, F77_int* 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.3.1'
+  VERSION = '0.5.0'
 end

data/lib/lbfgsb.rb

@@ -43,7 +43,7 @@ module Lbfgsb
     n_elements = x_init.size
     l = Numo::DFloat.zeros(n_elements)
     u = Numo::DFloat.zeros(n_elements)
-    nbd = Numo::Int64.zeros(n_elements)
+    nbd = SZ_F77_INTEGER == 64 ? Numo::Int64.zeros(n_elements) : Numo::Int32.zeros(n_elements)
 
     unless bounds.nil?
       n_elements.times do |n|

data/sig/lbfgsb.rbs

@@ -0,0 +1,31 @@
+module Lbfgsb
+  type result = {
+    task: String,
+    x: Numo::DFloat,
+    fnc: Float,
+    jcb: Numo::DFloat,
+    n_iter: Integer,
+    n_fev: Integer,
+    n_jev: Integer,
+    success: bool
+  }
+
+  VERSION: String
+
+  DBL_EPSILON: Float
+
+  def self?.minimize: (fnc: Method | Proc fnc, x_init: Numo::DFloat x_init, jcb: Method | Proc | bool jcb,
+    ?args: untyped args, ?bounds: Numo::DFloat? bounds, ?factr: Float factr, ?pgtol: Float pgtol,
+    ?maxcor: Integer maxcor, ?maxiter: Integer maxiter,
+    ?verbose: Integer? verbose) -> result
+
+  private
+
+  def self?.fnc: (Method | Proc fnc, Numo::DFloat x, untyped args) -> Float
+
+  def self?.jcb: (Method | Proc jcb, Numo::DFloat x, untyped args) -> Numo::DFloat
+
+  def self?.min_l_bfgs_b: (Method | Proc fnc, Numo::DFloat x, Method | Proc | bool jcb, untyped args,
+    Numo::DFloat l, Numo::DFloat u, Numo::DFloat nbd, Integer maxcor, Float ftol, Float gtol, Integer maxiter,
+    Integer? disp) -> result
+end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: lbfgsb
 version: !ruby/object:Gem::Version
-  version: 0.3.1
+  version: 0.5.0
 platform: ruby
 authors:
 - yoshoku
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2021-01-16 00:00:00.000000000 Z
+date: 2022-02-13 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: numo-narray
@@ -24,7 +24,9 @@ dependencies:
     - - ">="
       - !ruby/object:Gem::Version
         version: 0.9.1
-description: Lbfgsb.rb is a Ruby binding for L-BFGS-B
+description: |
+  Lbfgsb.rb is a Ruby binding for L-BFGS-B that is a limited-memory algorithm for
+  solving large nonlinear optimization problems subject to simple bounds on the variables.
 email:
 - yoshoku@outlook.com
 executables: []
@@ -32,28 +34,23 @@ extensions:
 - ext/lbfgsb/extconf.rb
 extra_rdoc_files: []
 files:
-- ".github/workflows/build.yml"
-- ".gitignore"
-- ".rspec"
-- ".yardopts"
 - CHANGELOG.md
-- Gemfile
 - LICENSE.txt
 - README.md
-- Rakefile
 - ext/lbfgsb/extconf.rb
 - ext/lbfgsb/lbfgsbext.c
 - ext/lbfgsb/lbfgsbext.h
 - ext/lbfgsb/src/License.txt
 - ext/lbfgsb/src/blas.c
 - ext/lbfgsb/src/blas.h
+- ext/lbfgsb/src/common.h
 - ext/lbfgsb/src/lbfgsb.c
 - ext/lbfgsb/src/lbfgsb.h
 - ext/lbfgsb/src/linpack.c
 - ext/lbfgsb/src/linpack.h
-- lbfgsb.gemspec
 - lib/lbfgsb.rb
 - lib/lbfgsb/version.rb
+- sig/lbfgsb.rbs
 homepage: https://github.com/yoshoku/lbfgsb.rb
 licenses:
 - BSD-3-Clause
@@ -77,8 +74,8 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.1.4
+rubygems_version: 3.3.7
 signing_key:
 specification_version: 4
-summary: Lbfgsb.rb is a Ruby binding for L-BFGS-B
+summary: Lbfgsb.rb is a Ruby binding for L-BFGS-B.
 test_files: []

data/.github/workflows/build.yml

@@ -1,21 +0,0 @@
-name: build
-
-on: [push, pull_request]
-
-jobs:
-  build:
-    runs-on: ubuntu-latest
-    strategy:
-      matrix:
-        ruby: [ '2.5', '2.6', '2.7' ]
-    steps:
-      - uses: actions/checkout@v2
-      - name: Set up Ruby ${{ matrix.ruby }}
-        uses: actions/setup-ruby@v1
-        with:
-          ruby-version: ${{ matrix.ruby }}
-      - name: Build and test with Rake
-        run: |
-          gem install bundler
-          bundle install --jobs 4 --retry 3
-          bundle exec rake

data/.gitignore

@@ -1,17 +0,0 @@
-/.bundle/
-/.yardoc
-/_yardoc/
-/coverage/
-/doc/
-/pkg/
-/spec/reports/
-/tmp/
-*.bundle
-*.so
-*.o
-*.a
-mkmf.log
-iterate.dat
-
-# rspec failure tracking
-.rspec_status

data/.rspec

@@ -1,3 +0,0 @@
---format documentation
---color
---require spec_helper

data/.yardopts

@@ -1 +0,0 @@
---exclude ext/lbfgsb/src/

data/Gemfile

@@ -1,8 +0,0 @@
-source "https://rubygems.org"
-
-# Specify your gem's dependencies in lbfgsb.gemspec
-gemspec
-
-gem "rake", "~> 12.0"
-gem "rake-compiler"
-gem "rspec", "~> 3.0"