ext/lbfgsb/src/linpack.c in lbfgsb-0.2.0 vs ext/lbfgsb/src/linpack.c in lbfgsb-0.3.0

- old
+ new

@@ -42,11 +42,11 @@ * of order k is not positive definite. * * 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; static double s, t; static long jm1; @@ -64,11 +64,11 @@ 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 + 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: @@ -131,11 +131,11 @@ * the first zero diagonal element of t. * * 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__; static double temp; @@ -177,11 +177,11 @@ } i__1 = *n; 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: goto L140; @@ -193,11 +193,11 @@ } i__1 = *n; 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: goto L140; @@ -209,11 +209,11 @@ } i__1 = *n; 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: goto L140; @@ -224,10 +224,10 @@ goto L130; } 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: L140: L150: