lbfgsb 0.1.0 → 0.4.0

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

@@ -44,9 +44,9 @@
  *                        March  2011
  */
 
+#include "lbfgsb.h"
 #include "blas.h"
 #include "linpack.h"
-#include "lbfgsb.h"
 
 static double c_b9 = 0.;
 static long c__1 = 1;
@@ -232,12 +232,8 @@ static double c_b282 = .1;
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-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)
-{
+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) {
   long i__1;
 
   static long ld, lr, lt, lz, lwa, lwn, lss, lxp, lws, lwt, lsy, lwy, lsnd;
@@ -287,11 +283,9 @@ int setulb_(long *n, long *m, double *x,
   lt = isave[14];
   lxp = isave[15];
   lwa = isave[16];
-  mainlb_(n, m, &x[1], &l[1], &u[1], &nbd[1], f, &g[1], factr, pgtol, &wa[lws],
-          &wa[lwy], &wa[lsy], &wa[lss], &wa[lwt], &wa[lwn], &wa[lsnd],
-          &wa[lz], &wa[lr], &wa[ld], &wa[lt], &wa[lxp], &wa[lwa], &iwa[1],
-          &iwa[*n + 1], &iwa[(*n << 1) + 1], task, iprint, csave, &lsave[1],
-          &isave[22], &dsave[1]);
+  mainlb_(n, m, &x[1], &l[1], &u[1], &nbd[1], f, &g[1], factr, pgtol, &wa[lws], &wa[lwy], &wa[lsy], &wa[lss], &wa[lwt],
+          &wa[lwn], &wa[lsnd], &wa[lz], &wa[lr], &wa[ld], &wa[lt], &wa[lxp], &wa[lwa], &iwa[1], &iwa[*n + 1],
+          &iwa[(*n << 1) + 1], task, iprint, csave, &lsave[1], &isave[22], &dsave[1]);
   return 0;
 }
 
@@ -471,20 +465,14 @@ int setulb_(long *n, long *m, double *x,
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-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)
-{
-  long ws_dim1, ws_offset, wy_dim1, wy_offset, sy_dim1, sy_offset,
-    ss_dim1, ss_offset, wt_dim1, wt_offset, wn_dim1, wn_offset,
-    snd_dim1, snd_offset, i__1;
+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) {
+  long ws_dim1, ws_offset, wy_dim1, wy_offset, sy_dim1, sy_offset, ss_dim1, ss_offset, wt_dim1, wt_offset, wn_dim1, wn_offset,
+      snd_dim1, snd_offset, i__1;
   double d__1, d__2;
-  FILE *itfptr;
+  FILE* itfptr;
   static long i__, k;
   static double gd, dr, rr, dtd;
   static long col;
@@ -610,9 +598,8 @@ int mainlb_(long *n, long *m, double *x,
     /* Check the input arguments for errors. */
     errclb_(n, m, factr, &l[1], &u[1], &nbd[1], task, &info, &k);
     if (strncmp(task, "ERROR", 5) == 0) {
-      prn3lb_(n, &x[1], f, task, iprint, &info, &itfile, &iter, &nfgv, &nintol,
-          &nskip, &nact, &sbgnrm, &c_b9, &nseg, word, &iback,
-          &stp, &xstep, &k, &cachyt, &sbtime, &lnscht);
+      prn3lb_(n, &x[1], f, task, iprint, &info, &itfile, &iter, &nfgv, &nintol, &nskip, &nact, &sbgnrm, &c_b9, &nseg, word,
+              &iback, &stp, &xstep, &k, &cachyt, &sbtime, &lnscht);
       return 0;
     }
     prn1lb_(n, m, &l[1], &u[1], &x[1], iprint, &itfile, &epsmch);
@@ -673,8 +660,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;
@@ -707,9 +694,9 @@ L222:
   }
   iword = -1;
 
-  if (! cnstnd && col > 0) {
+  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;
@@ -718,10 +705,9 @@ L222:
    * Compute the Generalized Cauchy Point (GCP).
    */
   timer_(&cpu1);
-  cauchy_(n, &x[1], &l[1], &u[1], &nbd[1], &g[1], &indx2[1], &iwhere[1], &t[1],
-      &d__[1], &z__[1], m, &wy[wy_offset], &ws[ws_offset], &sy[sy_offset],
-      &wt[wt_offset], &theta, &col, &head, &wa[1], &wa[(*m << 1) + 1],
-      &wa[(*m << 2) + 1], &wa[*m * 6 + 1], &nseg, iprint, &sbgnrm, &info, &epsmch);
+  cauchy_(n, &x[1], &l[1], &u[1], &nbd[1], &g[1], &indx2[1], &iwhere[1], &t[1], &d__[1], &z__[1], m, &wy[wy_offset],
+          &ws[ws_offset], &sy[sy_offset], &wt[wt_offset], &theta, &col, &head, &wa[1], &wa[(*m << 1) + 1], &wa[(*m << 2) + 1],
+          &wa[*m * 6 + 1], &nseg, iprint, &sbgnrm, &info, &epsmch);
   if (info != 0) {
     /* singular triangular system detected; refresh the lbfgs memory. */
     if (*iprint >= 1) {
@@ -762,9 +748,8 @@ L333:
   /*   where     E = [-I  0] */
   /*                 [ 0  I] */
   if (wrk) {
-    formk_(n, &nfree, &index[1], &nenter, &ileave, &indx2[1], &iupdat, &updatd,
-        &wn[wn_offset], &snd[snd_offset], m, &ws[ws_offset], &wy[wy_offset],
-        &sy[sy_offset], &theta, &col, &head, &info);
+    formk_(n, &nfree, &index[1], &nenter, &ileave, &indx2[1], &iupdat, &updatd, &wn[wn_offset], &snd[snd_offset], m,
+           &ws[ws_offset], &wy[wy_offset], &sy[sy_offset], &theta, &col, &head, &info);
   }
   if (info != 0) {
     /* nonpositive definiteness in Cholesky factorization; */
@@ -786,16 +771,14 @@ L333:
   }
   /* compute r=-Z'B(xcp-xk)-Z'g (using wa(2m+1)=W'(xcp-x) */
   /*                                            from 'cauchy'). */
-  cmprlb_(n, m, &x[1], &g[1], &ws[ws_offset], &wy[wy_offset], &sy[sy_offset],
-      &wt[wt_offset], &z__[1], &r__[1], &wa[1], &index[1], &theta, &col,
-      &head, &nfree, &cnstnd, &info);
+  cmprlb_(n, m, &x[1], &g[1], &ws[ws_offset], &wy[wy_offset], &sy[sy_offset], &wt[wt_offset], &z__[1], &r__[1], &wa[1],
+          &index[1], &theta, &col, &head, &nfree, &cnstnd, &info);
   if (info != 0) {
     goto L444;
   }
   /* jlm-jn call the direct method. */
-  subsm_(n, m, &nfree, &index[1], &l[1], &u[1], &nbd[1], &z__[1], &r__[1], &xp[1],
-      &ws[ws_offset], &wy[wy_offset], &theta, &x[1], &g[1], &col,
-      &head, &iword, &wa[1], &wn[wn_offset], iprint, &info);
+  subsm_(n, m, &nfree, &index[1], &l[1], &u[1], &nbd[1], &z__[1], &r__[1], &xp[1], &ws[ws_offset], &wy[wy_offset], &theta,
+         &x[1], &g[1], &col, &head, &iword, &wa[1], &wn[wn_offset], iprint, &info);
 L444:
   if (info != 0) {
     /* singular triangular system detected; */
@@ -828,14 +811,12 @@ L555:
   }
   timer_(&cpu1);
 L666:
-  lnsrlb_(n, &l[1], &u[1], &nbd[1], &x[1], f, &fold, &gd, &gdold, &g[1],
-      &d__[1], &r__[1], &t[1], &z__[1], &stp, &dnorm, &dtd, &xstep,
-      &stpmx, &iter, &ifun, &iback, &nfgv, &info, task, &boxed, &cnstnd,
-      csave, &isave[22], &dsave[17]);
+  lnsrlb_(n, &l[1], &u[1], &nbd[1], &x[1], f, &fold, &gd, &gdold, &g[1], &d__[1], &r__[1], &t[1], &z__[1], &stp, &dnorm, &dtd,
+          &xstep, &stpmx, &iter, &ifun, &iback, &nfgv, &info, task, &boxed, &cnstnd, 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. */
@@ -881,8 +862,7 @@ L666:
     /* Compute the infinity norm of the projected (-)gradient. */
     projgr_(n, &l[1], &u[1], &nbd[1], &x[1], &g[1], &sbgnrm);
     /* Print iteration information. */
-    prn2lb_(n, &x[1], f, &g[1], iprint, &itfile, &iter, &nfgv, &nact,
-        &sbgnrm, &nseg, word, &iword, &iback, &stp, &xstep);
+    prn2lb_(n, &x[1], f, &g[1], iprint, &itfile, &iter, &nfgv, &nact, &sbgnrm, &nseg, word, &iword, &iback, &stp, &xstep);
     goto L1000;
   }
 L777:
@@ -910,13 +890,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) {
@@ -934,9 +914,8 @@ L777:
   updatd = TRUE_;
   ++iupdat;
   /* Update matrices WS and WY and form the middle matrix in B. */
-  matupd_(n, m, &ws[ws_offset], &wy[wy_offset], &sy[sy_offset],
-      &ss[ss_offset], &d__[1], &r__[1], &itail, &iupdat, &col, &head,
-      &theta, &rr, &dr, &stp, &dtd);
+  matupd_(n, m, &ws[ws_offset], &wy[wy_offset], &sy[sy_offset], &ss[ss_offset], &d__[1], &r__[1], &itail, &iupdat, &col, &head,
+          &theta, &rr, &dr, &stp, &dtd);
   /* Form the upper half of the pds T = theta*SS + L*D^(-1)*L'; */
   /*    Store T in the upper triangular of the array wt; */
   /*    Cholesky factorize T to J*J' with */
@@ -967,9 +946,8 @@ L888:
 L999:
   timer_(&time2);
   time = time2 - time1;
-  prn3lb_(n, &x[1], f, task, iprint, &info, &itfile, &iter, &nfgv, &nintol,
-    &nskip, &nact, &sbgnrm, &time, &nseg, word, &iback, &stp, &xstep,
-    &k, &cachyt, &sbtime, &lnscht);
+  prn3lb_(n, &x[1], f, task, iprint, &info, &itfile, &iter, &nfgv, &nintol, &nskip, &nact, &sbgnrm, &time, &nseg, word, &iback,
+          &stp, &xstep, &k, &cachyt, &sbtime, &lnscht);
 L1000:
   /* Save local variables. */
   lsave[1] = prjctd;
@@ -1035,10 +1013,8 @@ L1000:
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-int active_(long *n, double *l, double *u,
-  long *nbd, double *x, long *iwhere, long *iprint,
-  long *prjctd, long *cnstnd, long *boxed)
-{
+int active_(long* n, double* l, double* u, long* nbd, double* x, long* iwhere, long* iprint, long* prjctd, long* cnstnd,
+            long* boxed) {
   long i__1;
   static long i__, nbdd;
   --iwhere;
@@ -1080,7 +1056,7 @@ int active_(long *n, double *l, double *u,
     if (nbd[i__] == 0) {
       /* this variable is always free */
       iwhere[i__] = -1;
-    /* otherwise set x(i)=mid(x(i), u(i), l(i)). */
+      /* otherwise set x(i)=mid(x(i), u(i), l(i)). */
     } else {
       *cnstnd = TRUE_;
       if (nbd[i__] == 2 && u[i__] - l[i__] <= 0.) {
@@ -1095,7 +1071,7 @@ int active_(long *n, double *l, double *u,
     if (*prjctd) {
       fprintf(stdout, " The initial X is infeasible.  Restart with its projection.\n");
     }
-    if (! (*cnstnd)) {
+    if (!(*cnstnd)) {
       fprintf(stdout, " This problem is unconstrained.\n");
     }
   }
@@ -1159,9 +1135,7 @@ int active_(long *n, double *l, double *u,
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-int bmv_(long *m, double *sy, double *wt, long
-  *col, double *v, double *p, long *info)
-{
+int bmv_(long* m, double* sy, double* wt, long* col, double* v, double* p, long* info) {
   long sy_dim1, sy_offset, wt_dim1, wt_offset, i__1, i__2;
   static long i__, k, i2;
   static double sum;
@@ -1193,7 +1167,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 +1179,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;
   }
@@ -1400,14 +1374,10 @@ int bmv_(long *m, double *sy, double *wt, long
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-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)
-{
+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) {
   long wy_dim1, wy_offset, ws_dim1, ws_offset, sy_dim1, sy_offset, wt_dim1, wt_offset, i__1, i__2;
   double d__1;
   static long i__, j;
@@ -1463,7 +1433,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 +1532,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 +1566,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 +1655,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 +1670,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 +1707,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");
@@ -1782,12 +1752,8 @@ L999:
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-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)
-{
+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) {
   long ws_dim1, ws_offset, wy_dim1, wy_offset, sy_dim1, sy_offset, wt_dim1, wt_offset, i__1, i__2;
   static long i__, j, k;
   static double a1, a2;
@@ -1812,7 +1778,7 @@ int cmprlb_(long *n, long *m, double *x,
   ws_offset = 1 + ws_dim1;
   ws -= ws_offset;
 
-  if (! (*cnstnd) && *col > 0) {
+  if (!(*cnstnd) && *col > 0) {
     i__1 = *n;
     for (i__ = 1; i__ <= i__1; ++i__) {
       r__[i__] = -g[i__];
@@ -1858,9 +1824,7 @@ int cmprlb_(long *n, long *m, double *x,
  *                       Ciyou Zhu
  *    in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-int errclb_(long *n, long *m, double *factr,
-  double *l, double *u, long *nbd, char *task, long *info, long *k)
-{
+int errclb_(long* n, long* m, double* factr, double* l, double* u, long* nbd, char* task, long* info, long* k) {
   long i__1;
   static long i__;
   --nbd;
@@ -2017,14 +1981,9 @@ int errclb_(long *n, long *m, double *factr,
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-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)
-{
-  long wn_dim1, wn_offset, wn1_dim1, wn1_offset, ws_dim1, ws_offset,
-    wy_dim1, wy_offset, sy_dim1, sy_offset, i__1, i__2, i__3;
+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) {
+  long wn_dim1, wn_offset, wn1_dim1, wn1_offset, ws_dim1, ws_offset, wy_dim1, wy_offset, sy_dim1, sy_offset, i__1, i__2, i__3;
   static long i__, k, k1, m2, is, js, iy, jy, is1, js1, col2, dend, pend;
   static long upcl;
   static double temp1, temp2, temp3, temp4;
@@ -2060,11 +2019,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 +2174,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 +2183,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 +2191,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;
@@ -2265,9 +2224,7 @@ int formk_(long *n, long *nsub, long *ind, long *nenter,
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-int formt_(long *m, double *wt, double *sy,
-  double *ss, long *col, double *theta, long *info)
-{
+int formt_(long* m, double* wt, double* sy, double* ss, long* col, double* theta, long* info) {
   long wt_dim1, wt_offset, sy_dim1, sy_offset, ss_dim1, ss_offset, i__1, i__2, i__3;
   static long i__, j, k, k1;
   static double ddum;
@@ -2303,7 +2260,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;
   }
@@ -2343,11 +2300,8 @@ int formt_(long *m, double *wt, double *sy,
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-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)
-{
+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) {
   long i__1;
   static long i__, k, iact;
 
@@ -2385,7 +2339,7 @@ int freev_(long *n, long *nfree, long *index,
     }
     if (*iprint >= 99) {
       i__1 = *n + 1 - *ileave;
-      fprintf(stdout,  " %2ld variables leave; %2ld variables enter\n", i__1, *nenter);
+      fprintf(stdout, " %2ld variables leave; %2ld variables enter\n", i__1, *nenter);
     }
   }
   *wrk = *ileave < *n + 1 || *nenter > 0 || *updatd;
@@ -2448,8 +2402,7 @@ int freev_(long *n, long *nfree, long *index,
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-int hpsolb_(long *n, double *t, long *iorder, long *iheap)
-{
+int hpsolb_(long* n, double* t, long* iorder, long* iheap) {
   long i__1;
   static long i__, j, k;
   static double out, ddum;
@@ -2466,7 +2419,7 @@ int hpsolb_(long *n, double *t, long *iorder, long *iheap)
       indxin = iorder[k];
       /* Add ddum to the heap. */
       i__ = k;
-L10:
+    L10:
       if (i__ > 1) {
         j = i__ / 2;
         if (ddum < t[j]) {
@@ -2490,7 +2443,7 @@ L10:
     ddum = t[*n];
     indxin = iorder[*n];
     /* Restore the heap */
-L30:
+  L30:
     j = i__ + i__;
     if (j <= *n - 1) {
       if (t[j + 1] < t[j]) {
@@ -2534,15 +2487,10 @@ L30:
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-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)
-{
+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) {
   long i__1;
   double d__1;
   static long i__;
@@ -2563,7 +2511,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;
@@ -2594,20 +2542,20 @@ int lnsrlb_(long *n, double *l, double *u,
       }
     }
   }
-  if (*iter == 0 && ! (*boxed)) {
+  if (*iter == 0 && !(*boxed)) {
     d__1 = 1. / *dnorm;
     *stp = d__1 <= *stpmx ? d__1 : *stpmx;
   } 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 +2574,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__) {
@@ -2658,12 +2606,8 @@ L556:
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-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)
-{
+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) {
   long ws_dim1, ws_offset, wy_dim1, wy_offset, sy_dim1, sy_offset, ss_dim1, ss_offset, i__1, i__2;
   static long j;
   static long pointr;
@@ -2693,8 +2637,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 +2648,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 +2658,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.) {
@@ -2743,12 +2687,9 @@ int matupd_(long *n, long *m, double *ws,
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-int prn1lb_(long *n, long *m, double *l,
-  double *u, double *x, long *iprint, long *itfile,
-  double *epsmch)
-{
+int prn1lb_(long* n, long* m, double* l, double* u, double* x, long* iprint, long* itfile, double* epsmch) {
   long i__1;
-  FILE *itfptr;
+  FILE* itfptr;
   static long i__;
 
   --x;
@@ -2840,14 +2781,11 @@ int prn1lb_(long *n, long *m, double *l,
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-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)
-{
+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) {
   long i__1;
   static long i__, imod;
-  FILE *itfptr;
+  FILE* itfptr;
   --g;
   --x;
 
@@ -2890,8 +2828,8 @@ int prn2lb_(long *n, double *x, double *f,
   }
   if (*iprint >= 1) {
     itfptr = fopen("iterate.dat", "a");
-    fprintf(itfptr, " %4ld %4ld %5ld %5ld  %3s %4ld  %7.1E  %7.1E %10.3E %10.3E\n",
-        *iter, *nfgv, *nseg, *nact, word, *iback, *stp, *xstep, *sbgnrm, *f);
+    fprintf(itfptr, " %4ld %4ld %5ld %5ld  %3s %4ld  %7.1E  %7.1E %10.3E %10.3E\n", *iter, *nfgv, *nseg, *nact, word, *iback,
+            *stp, *xstep, *sbgnrm, *f);
     fclose(itfptr);
   }
   return 0;
@@ -2913,15 +2851,11 @@ int prn2lb_(long *n, double *x, double *f,
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-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)
-{
+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) {
   long i__1;
-  FILE *itfptr;
+  FILE* itfptr;
   static long i__;
 
   --x;
@@ -3023,8 +2957,8 @@ L999:
     if (*iprint >= 1) {
       itfptr = fopen("iterate.dat", "a");
       if (*info == -4 || *info == -9) {
-        fprintf(itfptr, " %4ld %4ld %5ld %5ld  %3s %4ld  %7.1E  %7.1E      -          -\n",
-            *iter, *nfgv, *nseg, *nact, word, *iback, *stp, *xstep);
+        fprintf(itfptr, " %4ld %4ld %5ld %5ld  %3s %4ld  %7.1E  %7.1E      -          -\n", *iter, *nfgv, *nseg, *nact, word,
+                *iback, *stp, *xstep);
       }
       fprintf(itfptr, "\n");
       fprintf(itfptr, "%s\n", task);
@@ -3090,9 +3024,7 @@ L999:
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal.
  */
-int projgr_(long *n, double *l, double *u,
-  long *nbd, double *x, double *g, double *sbgnrm)
-{
+int projgr_(long* n, double* l, double* u, long* nbd, double* x, double* g, double* sbgnrm) {
   long i__1;
   double d__1, d__2;
   static long i__;
@@ -3299,13 +3231,9 @@ int projgr_(long *n, double *l, double *u,
  *                        Ciyou Zhu
  *     in collaboration with R.H. Byrd, P. Lu-Chen and J. Nocedal
  */
-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)
-{
+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) {
   long ws_dim1, ws_offset, wy_dim1, wy_offset, wn_dim1, wn_offset, i__1, i__2;
   double d__1, d__2;
   static long i__, j, k, m2;
@@ -3361,7 +3289,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 +3297,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;
   }
@@ -3381,18 +3309,17 @@ int subsm_(long *n, long *m, long *nsub, long *ind,
     i__2 = *nsub;
     for (i__ = 1; i__ <= i__2; ++i__) {
       k = ind[i__];
-      d__[i__] = d__[i__] + wy[k + pointr * wy_dim1] * wv[jy] / *theta
-        + ws[k + pointr * ws_dim1] * wv[js];
+      d__[i__] = d__[i__] + wy[k + pointr * wy_dim1] * wv[jy] / *theta + ws[k + pointr * ws_dim1] * wv[js];
     }
     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 +3372,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 {
@@ -3637,11 +3564,8 @@ L911:
  *     Argonne National Laboratory and University of Minnesota.
  *     Brett M. Averick, Richard G. Carter, and Jorge J. More'.
  */
-int dcsrch_(double *f, double *g, double *stp,
-  double *ftol, double *gtol, double *xtol,
-  double *stpmin, double *stpmax,
-  char *task, long *isave, double *dsave)
-{
+int dcsrch_(double* f, double* g, double* stp, double* ftol, double* gtol, double* xtol, double* stpmin, double* stpmax,
+            char* task, long* isave, double* dsave) {
 
   double d__1;
   static double fm, gm, fx, fy, gx, gy, fxm, fym, gxm, gym, stx, sty;
@@ -3917,11 +3841,8 @@ L1000:
  *     Argonne National Laboratory and University of Minnesota.
  *     Brett M. Averick and Jorge J. More'.
  */
-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)
-{
+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) {
   double d__1, d__2, d__3;
   static double p, q, r__, s, sgnd, stpc, stpf, stpq, gamma, theta;
 
@@ -3953,10 +3874,10 @@ int dcstep_(double *stx, double *fx, double *dx,
       stpf = stpc + (stpq - stpc) / 2.;
     }
     *brackt = TRUE_;
-  /* Second case: A lower function value and derivatives of opposite */
-  /* sign. The minimum is bracketed. If the cubic step is farther from */
-  /* stp than the secant step, the cubic step is taken, otherwise the */
-  /* secant step is taken. */
+    /* Second case: A lower function value and derivatives of opposite */
+    /* sign. The minimum is bracketed. If the cubic step is farther from */
+    /* stp than the secant step, the cubic step is taken, otherwise the */
+    /* secant step is taken. */
   } else if (sgnd < 0.) {
     theta = (*fx - *fp) * 3. / (*stp - *stx) + *dx + *dp;
     d__1 = fabs(theta);
@@ -3980,8 +3901,8 @@ int dcstep_(double *stx, double *fx, double *dx,
       stpf = stpq;
     }
     *brackt = TRUE_;
-  /* Third case: A lower function value, derivatives of the same sign, */
-  /* and the magnitude of the derivative decreases. */
+    /* Third case: A lower function value, derivatives of the same sign, */
+    /* and the magnitude of the derivative decreases. */
   } else if (fabs(*dp) < fabs(*dx)) {
     /* The cubic step is computed only if the cubic tends to infinity */
     /* in the direction of the step or if the minimum of the cubic */
@@ -4041,10 +3962,10 @@ int dcstep_(double *stx, double *fx, double *dx,
       stpf = *stpmax <= stpf ? *stpmax : stpf;
       stpf = *stpmin >= stpf ? *stpmin : stpf;
     }
-  /* Fourth case: A lower function value, derivatives of the same sign, */
-  /* and the magnitude of the derivative does not decrease. If the */
-  /* minimum is not bracketed, the step is either stpmin or stpmax, */
-  /* otherwise the cubic step is taken. */
+    /* Fourth case: A lower function value, derivatives of the same sign, */
+    /* and the magnitude of the derivative does not decrease. If the */
+    /* minimum is not bracketed, the step is either stpmin or stpmax, */
+    /* otherwise the cubic step is taken. */
   } else {
     if (*brackt) {
       theta = (*fp - *fy) * 3. / (*sty - *stp) + *dy + *dp;
@@ -4052,7 +3973,7 @@ int dcstep_(double *stx, double *fx, double *dx,
       d__2 = fabs(*dy);
       d__1 = d__1 >= d__2 ? d__1 : d__2;
       d__2 = fabs(*dp);
-      s = d__1 >= d__2 ? d__1: d__2;
+      s = d__1 >= d__2 ? d__1 : d__2;
       d__1 = theta / s;
       gamma = s * sqrt(d__1 * d__1 - *dy / s * (*dp / s));
       if (*stp > *sty) {
@@ -4089,8 +4010,7 @@ int dcstep_(double *stx, double *fx, double *dx,
   return 0;
 }
 
-int timer_(double *ttime)
-{
+int timer_(double* ttime) {
   *ttime = (double)clock() / CLOCKS_PER_SEC;
   return 0;
 }