File fix-CBLAS-warn.diff of Package superlu
Index: SuperLU_3.0/CBLAS/isamax.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/isamax.c
+++ SuperLU_3.0/CBLAS/isamax.c
@@ -10,9 +10,9 @@ integer isamax_(integer *n, real *sx, in
{
/* System generated locals */
- integer ret_val, i__1;
+ integer ret_val; // , i__1;
real r__1;
/* Local variables */
static real smax;
@@ -47,9 +47,9 @@ integer isamax_(integer *n, real *sx, in
ix = 1;
smax = dabs(SX(1));
ix += *incx;
- i__1 = *n;
+ //i__1 = *n;
for (i = 2; i <= *n; ++i) {
if ((r__1 = SX(ix), dabs(r__1)) <= smax) {
goto L5;
}
@@ -64,9 +64,9 @@ L5:
/* code for increment equal to 1 */
L20:
smax = dabs(SX(1));
- i__1 = *n;
+ //i__1 = *n;
for (i = 2; i <= *n; ++i) {
if ((r__1 = SX(i), dabs(r__1)) <= smax) {
goto L30;
}
Index: SuperLU_3.0/CBLAS/sasum.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/sasum.c
+++ SuperLU_3.0/CBLAS/sasum.c
@@ -10,9 +10,9 @@ real sasum_(integer *n, real *sx, intege
{
/* System generated locals */
- integer i__1, i__2;
+ //integer i__1, i__2;
real ret_val, r__1, r__2, r__3, r__4, r__5, r__6;
/* Local variables */
static integer i, m, nincx;
@@ -44,10 +44,10 @@ real sasum_(integer *n, real *sx, intege
/* code for increment not equal to 1 */
nincx = *n * *incx;
- i__1 = nincx;
- i__2 = *incx;
+ //i__1 = nincx;
+ //i__2 = *incx;
for (i = 1; *incx < 0 ? i >= nincx : i <= nincx; i += *incx) {
stemp += (r__1 = SX(i), dabs(r__1));
/* L10: */
}
@@ -63,9 +63,9 @@ L20:
m = *n % 6;
if (m == 0) {
goto L40;
}
- i__2 = m;
+ //i__2 = m;
for (i = 1; i <= m; ++i) {
stemp += (r__1 = SX(i), dabs(r__1));
/* L30: */
}
@@ -73,9 +73,9 @@ L20:
goto L60;
}
L40:
mp1 = m + 1;
- i__2 = *n;
+ //i__2 = *n;
for (i = mp1; i <= *n; i += 6) {
stemp = stemp + (r__1 = SX(i), dabs(r__1)) + (r__2 = SX(i + 1), dabs(
r__2)) + (r__3 = SX(i + 2), dabs(r__3)) + (r__4 = SX(i + 3),
dabs(r__4)) + (r__5 = SX(i + 4), dabs(r__5)) + (r__6 = SX(i +
Index: SuperLU_3.0/CBLAS/saxpy.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/saxpy.c
+++ SuperLU_3.0/CBLAS/saxpy.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1;
+ //integer i__1;
/* Local variables */
static integer i, m, ix, iy, mp1;
@@ -51,9 +51,9 @@
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
SY(iy) += *sa * SX(ix);
ix += *incx;
iy += *incy;
@@ -70,9 +70,9 @@ L20:
m = *n % 4;
if (m == 0) {
goto L40;
}
- i__1 = m;
+ //i__1 = m;
for (i = 1; i <= m; ++i) {
SY(i) += *sa * SX(i);
/* L30: */
}
@@ -80,9 +80,9 @@ L20:
return 0;
}
L40:
mp1 = m + 1;
- i__1 = *n;
+ //i__1 = *n;
for (i = mp1; i <= *n; i += 4) {
SY(i) += *sa * SX(i);
SY(i + 1) += *sa * SX(i + 1);
SY(i + 2) += *sa * SX(i + 2);
Index: SuperLU_3.0/CBLAS/scopy.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/scopy.c
+++ SuperLU_3.0/CBLAS/scopy.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1;
+ //integer i__1;
/* Local variables */
static integer i, m, ix, iy, mp1;
@@ -48,9 +48,9 @@
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
SY(iy) = SX(ix);
ix += *incx;
iy += *incy;
@@ -67,9 +67,9 @@ L20:
m = *n % 7;
if (m == 0) {
goto L40;
}
- i__1 = m;
+ //i__1 = m;
for (i = 1; i <= m; ++i) {
SY(i) = SX(i);
/* L30: */
}
@@ -77,9 +77,9 @@ L20:
return 0;
}
L40:
mp1 = m + 1;
- i__1 = *n;
+ //i__1 = *n;
for (i = mp1; i <= *n; i += 7) {
SY(i) = SX(i);
SY(i + 1) = SX(i + 1);
SY(i + 2) = SX(i + 2);
Index: SuperLU_3.0/CBLAS/sdot.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/sdot.c
+++ SuperLU_3.0/CBLAS/sdot.c
@@ -10,9 +10,9 @@ real sdot_(integer *n, real *sx, integer
{
/* System generated locals */
- integer i__1;
+ //integer i__1;
real ret_val;
/* Local variables */
static integer i, m;
@@ -52,9 +52,9 @@ real sdot_(integer *n, real *sx, integer
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
stemp += SX(ix) * SY(iy);
ix += *incx;
iy += *incy;
@@ -72,9 +72,9 @@ L20:
m = *n % 5;
if (m == 0) {
goto L40;
}
- i__1 = m;
+ //i__1 = m;
for (i = 1; i <= m; ++i) {
stemp += SX(i) * SY(i);
/* L30: */
}
@@ -82,9 +82,9 @@ L20:
goto L60;
}
L40:
mp1 = m + 1;
- i__1 = *n;
+ //i__1 = *n;
for (i = mp1; i <= *n; i += 5) {
stemp = stemp + SX(i) * SY(i) + SX(i + 1) * SY(i + 1) + SX(i + 2) *
SY(i + 2) + SX(i + 3) * SY(i + 3) + SX(i + 4) * SY(i + 4);
/* L50: */
Index: SuperLU_3.0/CBLAS/sgemv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/sgemv.c
+++ SuperLU_3.0/CBLAS/sgemv.c
@@ -12,9 +12,9 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2;
+ //integer a_dim1, a_offset, i__1, i__2;
/* Local variables */
static integer info;
static real temp;
@@ -150,9 +150,9 @@
}
/* Quick return if possible. */
- if (*m == 0 || *n == 0 || *alpha == 0.f && *beta == 1.f) {
+ if (*m == 0 || *n == 0 || (*alpha == 0.f && *beta == 1.f)) {
return 0;
}
/* Set LENX and LENY, the lengths of the vectors x and y, and set
@@ -184,31 +184,31 @@
if (*beta != 1.f) {
if (*incy == 1) {
if (*beta == 0.f) {
- i__1 = leny;
+ //i__1 = leny;
for (i = 1; i <= leny; ++i) {
Y(i) = 0.f;
/* L10: */
}
} else {
- i__1 = leny;
+ //i__1 = leny;
for (i = 1; i <= leny; ++i) {
Y(i) = *beta * Y(i);
/* L20: */
}
}
} else {
iy = ky;
if (*beta == 0.f) {
- i__1 = leny;
+ //i__1 = leny;
for (i = 1; i <= leny; ++i) {
Y(iy) = 0.f;
iy += *incy;
/* L30: */
}
} else {
- i__1 = leny;
+ //i__1 = leny;
for (i = 1; i <= leny; ++i) {
Y(iy) = *beta * Y(iy);
iy += *incy;
/* L40: */
@@ -224,13 +224,13 @@
/* Form y := alpha*A*x + y. */
jx = kx;
if (*incy == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(jx) != 0.f) {
temp = *alpha * X(jx);
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
Y(i) += temp * A(i,j);
/* L50: */
}
@@ -238,14 +238,14 @@
jx += *incx;
/* L60: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(jx) != 0.f) {
temp = *alpha * X(jx);
iy = ky;
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
Y(iy) += temp * A(i,j);
iy += *incy;
/* L70: */
@@ -260,12 +260,12 @@
/* Form y := alpha*A'*x + y. */
jy = ky;
if (*incx == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp = 0.f;
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
temp += A(i,j) * X(i);
/* L90: */
}
@@ -273,13 +273,13 @@
jy += *incy;
/* L100: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp = 0.f;
ix = kx;
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
temp += A(i,j) * X(ix);
ix += *incx;
/* L110: */
Index: SuperLU_3.0/CBLAS/snrm2.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/snrm2.c
+++ SuperLU_3.0/CBLAS/snrm2.c
@@ -10,9 +10,9 @@ real snrm2_(integer *n, real *x, integer
{
/* System generated locals */
- integer i__1, i__2;
+ //integer i__1, i__2;
real ret_val, r__1;
/* Builtin functions */
double sqrt(doublereal);
@@ -52,10 +52,10 @@ real snrm2_(integer *n, real *x, integer
auxiliary routine:
CALL SLASSQ( N, X, INCX, SCALE, SSQ ) */
- i__1 = (*n - 1) * *incx + 1;
- i__2 = *incx;
+ //i__1 = (*n - 1) * *incx + 1;
+ //i__2 = *incx;
for (ix = 1; *incx < 0 ? ix >= (*n-1)**incx+1 : ix <= (*n-1)**incx+1; ix += *incx) {
if (X(ix) != 0.f) {
absxi = (r__1 = X(ix), dabs(r__1));
if (scale < absxi) {
Index: SuperLU_3.0/CBLAS/srot.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/srot.c
+++ SuperLU_3.0/CBLAS/srot.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1;
+ //integer i__1;
/* Local variables */
static integer i;
static real stemp;
@@ -49,9 +49,9 @@
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
stemp = *c * SX(ix) + *s * SY(iy);
SY(iy) = *c * SY(iy) - *s * SX(ix);
SX(ix) = stemp;
@@ -63,9 +63,9 @@
/* code for both increments equal to 1 */
L20:
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
stemp = *c * SX(i) + *s * SY(i);
SY(i) = *c * SY(i) - *s * SX(i);
SX(i) = stemp;
Index: SuperLU_3.0/CBLAS/sscal.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/sscal.c
+++ SuperLU_3.0/CBLAS/sscal.c
@@ -10,9 +10,9 @@
{
/* System generated locals */
- integer i__1, i__2;
+ //integer i__1, i__2;
/* Local variables */
static integer i, m, nincx, mp1;
@@ -39,10 +39,10 @@
/* code for increment not equal to 1 */
nincx = *n * *incx;
- i__1 = nincx;
- i__2 = *incx;
+ //i__1 = nincx;
+ //i__2 = *incx;
for (i = 1; *incx < 0 ? i >= nincx : i <= nincx; i += *incx) {
SX(i) = *sa * SX(i);
/* L10: */
}
@@ -57,9 +57,9 @@ L20:
m = *n % 5;
if (m == 0) {
goto L40;
}
- i__2 = m;
+ //i__2 = m;
for (i = 1; i <= m; ++i) {
SX(i) = *sa * SX(i);
/* L30: */
}
@@ -67,9 +67,9 @@ L20:
return 0;
}
L40:
mp1 = m + 1;
- i__2 = *n;
+ //i__2 = *n;
for (i = mp1; i <= *n; i += 5) {
SX(i) = *sa * SX(i);
SX(i + 1) = *sa * SX(i + 1);
SX(i + 2) = *sa * SX(i + 2);
Index: SuperLU_3.0/CBLAS/ssymv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/ssymv.c
+++ SuperLU_3.0/CBLAS/ssymv.c
@@ -12,9 +12,9 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2;
+ //integer a_dim1, a_offset, i__1, i__2;
/* Local variables */
static integer info;
static real temp1, temp2;
@@ -145,9 +145,9 @@
}
/* Quick return if possible. */
- if (*n == 0 || *alpha == 0.f && *beta == 1.f) {
+ if (*n == 0 || (*alpha == 0.f && *beta == 1.f)) {
return 0;
}
/* Set up the start points in X and Y. */
@@ -171,31 +171,31 @@
if (*beta != 1.f) {
if (*incy == 1) {
if (*beta == 0.f) {
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
Y(i) = 0.f;
/* L10: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
Y(i) = *beta * Y(i);
/* L20: */
}
}
} else {
iy = ky;
if (*beta == 0.f) {
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
Y(iy) = 0.f;
iy += *incy;
/* L30: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
Y(iy) = *beta * Y(iy);
iy += *incy;
/* L40: */
@@ -210,13 +210,13 @@
/* Form y when A is stored in upper triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp1 = *alpha * X(j);
temp2 = 0.f;
- i__2 = j - 1;
+ //i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
Y(i) += temp1 * A(i,j);
temp2 += A(i,j) * X(i);
/* L50: */
@@ -226,15 +226,15 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp1 = *alpha * X(jx);
temp2 = 0.f;
ix = kx;
iy = ky;
- i__2 = j - 1;
+ //i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
Y(iy) += temp1 * A(i,j);
temp2 += A(i,j) * X(ix);
ix += *incx;
@@ -251,14 +251,14 @@
/* Form y when A is stored in lower triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp1 = *alpha * X(j);
temp2 = 0.f;
Y(j) += temp1 * A(j,j);
- i__2 = *n;
+ //i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
Y(i) += temp1 * A(i,j);
temp2 += A(i,j) * X(i);
/* L90: */
@@ -268,16 +268,16 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp1 = *alpha * X(jx);
temp2 = 0.f;
Y(jy) += temp1 * A(j,j);
ix = jx;
iy = jy;
- i__2 = *n;
+ //i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
iy += *incy;
Y(iy) += temp1 * A(i,j);
Index: SuperLU_3.0/CBLAS/dasum.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dasum.c
+++ SuperLU_3.0/CBLAS/dasum.c
@@ -10,9 +10,9 @@ doublereal dasum_(integer *n, doublereal
{
/* System generated locals */
- integer i__1, i__2;
+ //integer i__1, i__2;
doublereal ret_val, d__1, d__2, d__3, d__4, d__5, d__6;
/* Local variables */
static integer i, m;
@@ -43,10 +43,10 @@ doublereal dasum_(integer *n, doublereal
/* code for increment not equal to 1 */
nincx = *n * *incx;
- i__1 = nincx;
- i__2 = *incx;
+ //i__1 = nincx;
+ //i__2 = *incx;
for (i = 1; *incx < 0 ? i >= nincx : i <= nincx; i += *incx) {
dtemp += (d__1 = DX(i), abs(d__1));
/* L10: */
}
@@ -62,9 +62,9 @@ L20:
m = *n % 6;
if (m == 0) {
goto L40;
}
- i__2 = m;
+ //i__2 = m;
for (i = 1; i <= m; ++i) {
dtemp += (d__1 = DX(i), abs(d__1));
/* L30: */
}
@@ -72,9 +72,9 @@ L20:
goto L60;
}
L40:
mp1 = m + 1;
- i__2 = *n;
+ //i__2 = *n;
for (i = mp1; i <= *n; i += 6) {
dtemp = dtemp + (d__1 = DX(i), abs(d__1)) + (d__2 = DX(i + 1), abs(
d__2)) + (d__3 = DX(i + 2), abs(d__3)) + (d__4 = DX(i + 3),
abs(d__4)) + (d__5 = DX(i + 4), abs(d__5)) + (d__6 = DX(i + 5)
Index: SuperLU_3.0/CBLAS/daxpy.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/daxpy.c
+++ SuperLU_3.0/CBLAS/daxpy.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1;
+ //integer i__1;
/* Local variables */
static integer i, m, ix, iy, mp1;
@@ -51,9 +51,9 @@
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
DY(iy) += *da * DX(ix);
ix += *incx;
iy += *incy;
@@ -70,9 +70,9 @@ L20:
m = *n % 4;
if (m == 0) {
goto L40;
}
- i__1 = m;
+ //i__1 = m;
for (i = 1; i <= m; ++i) {
DY(i) += *da * DX(i);
/* L30: */
}
@@ -80,9 +80,9 @@ L20:
return 0;
}
L40:
mp1 = m + 1;
- i__1 = *n;
+ //i__1 = *n;
for (i = mp1; i <= *n; i += 4) {
DY(i) += *da * DX(i);
DY(i + 1) += *da * DX(i + 1);
DY(i + 2) += *da * DX(i + 2);
Index: SuperLU_3.0/CBLAS/dcopy.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dcopy.c
+++ SuperLU_3.0/CBLAS/dcopy.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1;
+ //integer i__1;
/* Local variables */
static integer i, m, ix, iy, mp1;
@@ -48,9 +48,9 @@
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
DY(iy) = DX(ix);
ix += *incx;
iy += *incy;
@@ -67,9 +67,9 @@ L20:
m = *n % 7;
if (m == 0) {
goto L40;
}
- i__1 = m;
+ //i__1 = m;
for (i = 1; i <= m; ++i) {
DY(i) = DX(i);
/* L30: */
}
@@ -77,9 +77,9 @@ L20:
return 0;
}
L40:
mp1 = m + 1;
- i__1 = *n;
+ //i__1 = *n;
for (i = mp1; i <= *n; i += 7) {
DY(i) = DX(i);
DY(i + 1) = DX(i + 1);
DY(i + 2) = DX(i + 2);
Index: SuperLU_3.0/CBLAS/ddot.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/ddot.c
+++ SuperLU_3.0/CBLAS/ddot.c
@@ -11,9 +11,9 @@ doublereal ddot_(integer *n, doublereal
{
/* System generated locals */
- integer i__1;
+ //integer i__1;
doublereal ret_val;
/* Local variables */
static integer i, m;
@@ -53,9 +53,9 @@ doublereal ddot_(integer *n, doublereal
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
dtemp += DX(ix) * DY(iy);
ix += *incx;
iy += *incy;
@@ -73,9 +73,9 @@ L20:
m = *n % 5;
if (m == 0) {
goto L40;
}
- i__1 = m;
+ //i__1 = m;
for (i = 1; i <= m; ++i) {
dtemp += DX(i) * DY(i);
/* L30: */
}
@@ -83,9 +83,9 @@ L20:
goto L60;
}
L40:
mp1 = m + 1;
- i__1 = *n;
+ //i__1 = *n;
for (i = mp1; i <= *n; i += 5) {
dtemp = dtemp + DX(i) * DY(i) + DX(i + 1) * DY(i + 1) + DX(i + 2) *
DY(i + 2) + DX(i + 3) * DY(i + 3) + DX(i + 4) * DY(i + 4);
/* L50: */
Index: SuperLU_3.0/CBLAS/idamax.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/idamax.c
+++ SuperLU_3.0/CBLAS/idamax.c
@@ -10,9 +10,9 @@ integer idamax_(integer *n, doublereal *
{
/* System generated locals */
- integer ret_val, i__1;
+ integer ret_val; // , i__1;
doublereal d__1;
/* Local variables */
static doublereal dmax__;
@@ -47,9 +47,9 @@ integer idamax_(integer *n, doublereal *
ix = 1;
dmax__ = abs(DX(1));
ix += *incx;
- i__1 = *n;
+ //i__1 = *n;
for (i = 2; i <= *n; ++i) {
if ((d__1 = DX(ix), abs(d__1)) <= dmax__) {
goto L5;
}
@@ -64,9 +64,9 @@ L5:
/* code for increment equal to 1 */
L20:
dmax__ = abs(DX(1));
- i__1 = *n;
+ //i__1 = *n;
for (i = 2; i <= *n; ++i) {
if ((d__1 = DX(i), abs(d__1)) <= dmax__) {
goto L30;
}
Index: SuperLU_3.0/CBLAS/sger.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/sger.c
+++ SuperLU_3.0/CBLAS/sger.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2;
+ //integer a_dim1, a_offset, i__1, i__2;
/* Local variables */
static integer info;
static real temp;
@@ -136,13 +136,13 @@
} else {
jy = 1 - (*n - 1) * *incy;
}
if (*incx == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (Y(jy) != 0.f) {
temp = *alpha * Y(jy);
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
A(i,j) += X(i) * temp;
/* L10: */
}
@@ -155,14 +155,14 @@
kx = 1;
} else {
kx = 1 - (*m - 1) * *incx;
}
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (Y(jy) != 0.f) {
temp = *alpha * Y(jy);
ix = kx;
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
A(i,j) += X(ix) * temp;
ix += *incx;
/* L30: */
Index: SuperLU_3.0/CBLAS/ssyr2.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/ssyr2.c
+++ SuperLU_3.0/CBLAS/ssyr2.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2;
+ //integer a_dim1, a_offset, i__1, i__2;
/* Local variables */
static integer info;
static real temp1, temp2;
@@ -175,14 +175,14 @@
/* Form A when A is stored in the upper triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(j) != 0.f || Y(j) != 0.f) {
temp1 = *alpha * Y(j);
temp2 = *alpha * X(j);
- i__2 = j;
+ //i__2 = j;
for (i = 1; i <= j; ++i) {
A(i,j) = A(i,j) + X(i) * temp1
+ Y(i) * temp2;
/* L10: */
@@ -190,16 +190,16 @@
}
/* L20: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(jx) != 0.f || Y(jy) != 0.f) {
temp1 = *alpha * Y(jy);
temp2 = *alpha * X(jx);
ix = kx;
iy = ky;
- i__2 = j;
+ //i__2 = j;
for (i = 1; i <= j; ++i) {
A(i,j) = A(i,j) + X(ix) * temp1
+ Y(iy) * temp2;
ix += *incx;
@@ -216,14 +216,14 @@
/* Form A when A is stored in the lower triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(j) != 0.f || Y(j) != 0.f) {
temp1 = *alpha * Y(j);
temp2 = *alpha * X(j);
- i__2 = *n;
+ //i__2 = *n;
for (i = j; i <= *n; ++i) {
A(i,j) = A(i,j) + X(i) * temp1
+ Y(i) * temp2;
/* L50: */
@@ -231,16 +231,16 @@
}
/* L60: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(jx) != 0.f || Y(jy) != 0.f) {
temp1 = *alpha * Y(jy);
temp2 = *alpha * X(jx);
ix = jx;
iy = jy;
- i__2 = *n;
+ //i__2 = *n;
for (i = j; i <= *n; ++i) {
A(i,j) = A(i,j) + X(ix) * temp1
+ Y(iy) * temp2;
ix += *incx;
Index: SuperLU_3.0/CBLAS/strsv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/strsv.c
+++ SuperLU_3.0/CBLAS/strsv.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2;
+ //integer a_dim1, a_offset, i__1, i__2;
/* Local variables */
static integer info;
static real temp;
@@ -214,16 +214,16 @@
}
}
} else {
if (*incx == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(j) != 0.f) {
if (nounit) {
X(j) /= A(j,j);
}
temp = X(j);
- i__2 = *n;
+ //i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
X(i) -= temp * A(i,j);
/* L50: */
}
@@ -231,17 +231,17 @@
/* L60: */
}
} else {
jx = kx;
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(jx) != 0.f) {
if (nounit) {
X(jx) /= A(j,j);
}
temp = X(jx);
ix = jx;
- i__2 = *n;
+ //i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
X(ix) -= temp * A(i,j);
/* L70: */
@@ -257,12 +257,12 @@
/* Form x := inv( A' )*x. */
if (lsame_(uplo, "U")) {
if (*incx == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp = X(j);
- i__2 = j - 1;
+ //i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
temp -= A(i,j) * X(i);
/* L90: */
}
@@ -273,13 +273,13 @@
/* L100: */
}
} else {
jx = kx;
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp = X(jx);
ix = kx;
- i__2 = j - 1;
+ //i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
temp -= A(i,j) * X(ix);
ix += *incx;
/* L110: */
@@ -295,9 +295,9 @@
} else {
if (*incx == 1) {
for (j = *n; j >= 1; --j) {
temp = X(j);
- i__1 = j + 1;
+ //i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
temp -= A(i,j) * X(i);
/* L130: */
}
@@ -312,9 +312,9 @@
jx = kx;
for (j = *n; j >= 1; --j) {
temp = X(jx);
ix = kx;
- i__1 = j + 1;
+ //i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
temp -= A(i,j) * X(ix);
ix -= *incx;
/* L150: */
Index: SuperLU_3.0/CBLAS/dgemv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dgemv.c
+++ SuperLU_3.0/CBLAS/dgemv.c
@@ -12,9 +12,9 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2;
+ //integer a_dim1, a_offset, i__1, i__2;
/* Local variables */
static integer info;
static doublereal temp;
@@ -150,9 +150,9 @@
}
/* Quick return if possible. */
- if (*m == 0 || *n == 0 || *alpha == 0. && *beta == 1.) {
+ if (*m == 0 || *n == 0 || (*alpha == 0. && *beta == 1.)) {
return 0;
}
/* Set LENX and LENY, the lengths of the vectors x and y, and set
@@ -184,31 +184,31 @@
if (*beta != 1.) {
if (*incy == 1) {
if (*beta == 0.) {
- i__1 = leny;
+ //i__1 = leny;
for (i = 1; i <= leny; ++i) {
Y(i) = 0.;
/* L10: */
}
} else {
- i__1 = leny;
+ //i__1 = leny;
for (i = 1; i <= leny; ++i) {
Y(i) = *beta * Y(i);
/* L20: */
}
}
} else {
iy = ky;
if (*beta == 0.) {
- i__1 = leny;
+ //i__1 = leny;
for (i = 1; i <= leny; ++i) {
Y(iy) = 0.;
iy += *incy;
/* L30: */
}
} else {
- i__1 = leny;
+ //i__1 = leny;
for (i = 1; i <= leny; ++i) {
Y(iy) = *beta * Y(iy);
iy += *incy;
/* L40: */
@@ -224,13 +224,13 @@
/* Form y := alpha*A*x + y. */
jx = kx;
if (*incy == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(jx) != 0.) {
temp = *alpha * X(jx);
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
Y(i) += temp * A(i,j);
/* L50: */
}
@@ -238,14 +238,14 @@
jx += *incx;
/* L60: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(jx) != 0.) {
temp = *alpha * X(jx);
iy = ky;
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
Y(iy) += temp * A(i,j);
iy += *incy;
/* L70: */
@@ -260,12 +260,12 @@
/* Form y := alpha*A'*x + y. */
jy = ky;
if (*incx == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp = 0.;
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
temp += A(i,j) * X(i);
/* L90: */
}
@@ -273,13 +273,13 @@
jy += *incy;
/* L100: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp = 0.;
ix = kx;
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
temp += A(i,j) * X(ix);
ix += *incx;
/* L110: */
Index: SuperLU_3.0/CBLAS/dnrm2.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dnrm2.c
+++ SuperLU_3.0/CBLAS/dnrm2.c
@@ -10,9 +10,9 @@ doublereal dnrm2_(integer *n, doublereal
{
/* System generated locals */
- integer i__1, i__2;
+ //integer i__1, i__2;
doublereal ret_val, d__1;
/* Builtin functions */
double sqrt(doublereal);
@@ -52,10 +52,10 @@ doublereal dnrm2_(integer *n, doublereal
auxiliary routine:
CALL DLASSQ( N, X, INCX, SCALE, SSQ ) */
- i__1 = (*n - 1) * *incx + 1;
- i__2 = *incx;
+ //i__1 = (*n - 1) * *incx + 1;
+ //i__2 = *incx;
for (ix = 1; *incx < 0 ? ix >= (*n-1)**incx+1 : ix <= (*n-1)**incx+1; ix += *incx) {
if (X(ix) != 0.) {
absxi = (d__1 = X(ix), abs(d__1));
if (scale < absxi) {
Index: SuperLU_3.0/CBLAS/drot.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/drot.c
+++ SuperLU_3.0/CBLAS/drot.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1;
+ //integer i__1;
/* Local variables */
static integer i;
static doublereal dtemp;
@@ -49,9 +49,9 @@
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
dtemp = *c * DX(ix) + *s * DY(iy);
DY(iy) = *c * DY(iy) - *s * DX(ix);
DX(ix) = dtemp;
@@ -63,9 +63,9 @@
/* code for both increments equal to 1 */
L20:
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
dtemp = *c * DX(i) + *s * DY(i);
DY(i) = *c * DY(i) - *s * DX(i);
DX(i) = dtemp;
Index: SuperLU_3.0/CBLAS/dscal.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dscal.c
+++ SuperLU_3.0/CBLAS/dscal.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1, i__2;
+ //integer i__1, i__2;
/* Local variables */
static integer i, m, nincx, mp1;
@@ -40,10 +40,10 @@
/* code for increment not equal to 1 */
nincx = *n * *incx;
- i__1 = nincx;
- i__2 = *incx;
+ //i__1 = nincx;
+ //i__2 = *incx;
for (i = 1; *incx < 0 ? i >= nincx : i <= nincx; i += *incx) {
DX(i) = *da * DX(i);
/* L10: */
}
@@ -58,9 +58,9 @@ L20:
m = *n % 5;
if (m == 0) {
goto L40;
}
- i__2 = m;
+ //i__2 = m;
for (i = 1; i <= m; ++i) {
DX(i) = *da * DX(i);
/* L30: */
}
@@ -68,9 +68,9 @@ L20:
return 0;
}
L40:
mp1 = m + 1;
- i__2 = *n;
+ //i__2 = *n;
for (i = mp1; i <= *n; i += 5) {
DX(i) = *da * DX(i);
DX(i + 1) = *da * DX(i + 1);
DX(i + 2) = *da * DX(i + 2);
Index: SuperLU_3.0/CBLAS/dsymv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dsymv.c
+++ SuperLU_3.0/CBLAS/dsymv.c
@@ -12,9 +12,9 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2;
+ //integer a_dim1, a_offset, i__1, i__2;
/* Local variables */
static integer info;
static doublereal temp1, temp2;
@@ -145,9 +145,9 @@
}
/* Quick return if possible. */
- if (*n == 0 || *alpha == 0. && *beta == 1.) {
+ if (*n == 0 || (*alpha == 0. && *beta == 1.)) {
return 0;
}
/* Set up the start points in X and Y. */
@@ -171,31 +171,31 @@
if (*beta != 1.) {
if (*incy == 1) {
if (*beta == 0.) {
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
Y(i) = 0.;
/* L10: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
Y(i) = *beta * Y(i);
/* L20: */
}
}
} else {
iy = ky;
if (*beta == 0.) {
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
Y(iy) = 0.;
iy += *incy;
/* L30: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
Y(iy) = *beta * Y(iy);
iy += *incy;
/* L40: */
@@ -210,13 +210,13 @@
/* Form y when A is stored in upper triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp1 = *alpha * X(j);
temp2 = 0.;
- i__2 = j - 1;
+ //i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
Y(i) += temp1 * A(i,j);
temp2 += A(i,j) * X(i);
/* L50: */
@@ -226,15 +226,15 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp1 = *alpha * X(jx);
temp2 = 0.;
ix = kx;
iy = ky;
- i__2 = j - 1;
+ //i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
Y(iy) += temp1 * A(i,j);
temp2 += A(i,j) * X(ix);
ix += *incx;
@@ -251,14 +251,14 @@
/* Form y when A is stored in lower triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp1 = *alpha * X(j);
temp2 = 0.;
Y(j) += temp1 * A(j,j);
- i__2 = *n;
+ //i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
Y(i) += temp1 * A(i,j);
temp2 += A(i,j) * X(i);
/* L90: */
@@ -268,16 +268,16 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp1 = *alpha * X(jx);
temp2 = 0.;
Y(jy) += temp1 * A(j,j);
ix = jx;
iy = jy;
- i__2 = *n;
+ //i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
iy += *incy;
Y(iy) += temp1 * A(i,j);
Index: SuperLU_3.0/CBLAS/caxpy.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/caxpy.c
+++ SuperLU_3.0/CBLAS/caxpy.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1, i__2, i__3, i__4;
+ //integer i__1, i__2, i__3, i__4;
real r__1, r__2;
complex q__1, q__2;
/* Builtin functions */
@@ -55,13 +55,13 @@
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
- i__3 = iy;
- i__4 = ix;
+ //i__2 = iy;
+ //i__3 = iy;
+ //i__4 = ix;
q__2.r = ca->r * CX(ix).r - ca->i * CX(ix).i, q__2.i = ca->r * CX(
ix).i + ca->i * CX(ix).r;
q__1.r = CY(iy).r + q__2.r, q__1.i = CY(iy).i + q__2.i;
CY(iy).r = q__1.r, CY(iy).i = q__1.i;
@@ -73,13 +73,13 @@
/* code for both increments equal to 1 */
L20:
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
- i__3 = i;
- i__4 = i;
+ //i__2 = i;
+ //i__3 = i;
+ //i__4 = i;
q__2.r = ca->r * CX(i).r - ca->i * CX(i).i, q__2.i = ca->r * CX(
i).i + ca->i * CX(i).r;
q__1.r = CY(i).r + q__2.r, q__1.i = CY(i).i + q__2.i;
CY(i).r = q__1.r, CY(i).i = q__1.i;
Index: SuperLU_3.0/CBLAS/ccopy.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/ccopy.c
+++ SuperLU_3.0/CBLAS/ccopy.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1, i__2, i__3;
+ //integer i__1, i__2, i__3;
/* Local variables */
static integer i, ix, iy;
@@ -47,12 +47,12 @@
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
- i__3 = ix;
+ //i__2 = iy;
+ //i__3 = ix;
CY(iy).r = CX(ix).r, CY(iy).i = CX(ix).i;
ix += *incx;
iy += *incy;
/* L10: */
@@ -61,12 +61,12 @@
/* code for both increments equal to 1 */
L20:
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
- i__3 = i;
+ //i__2 = i;
+ //i__3 = i;
CY(i).r = CX(i).r, CY(i).i = CX(i).i;
/* L30: */
}
return 0;
Index: SuperLU_3.0/CBLAS/dger.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dger.c
+++ SuperLU_3.0/CBLAS/dger.c
@@ -12,9 +12,9 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2;
+ //integer a_dim1, a_offset, i__1, i__2;
/* Local variables */
static integer info;
static doublereal temp;
@@ -137,13 +137,13 @@
} else {
jy = 1 - (*n - 1) * *incy;
}
if (*incx == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (Y(jy) != 0.) {
temp = *alpha * Y(jy);
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
A(i,j) += X(i) * temp;
/* L10: */
}
@@ -156,14 +156,14 @@
kx = 1;
} else {
kx = 1 - (*m - 1) * *incx;
}
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (Y(jy) != 0.) {
temp = *alpha * Y(jy);
ix = kx;
- i__2 = *m;
+ //i__2 = *m;
for (i = 1; i <= *m; ++i) {
A(i,j) += X(ix) * temp;
ix += *incx;
/* L30: */
Index: SuperLU_3.0/CBLAS/dsyr2.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dsyr2.c
+++ SuperLU_3.0/CBLAS/dsyr2.c
@@ -12,9 +12,9 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2;
+ //integer a_dim1, a_offset, i__1, i__2;
/* Local variables */
static integer info;
static doublereal temp1, temp2;
@@ -176,14 +176,14 @@
/* Form A when A is stored in the upper triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(j) != 0. || Y(j) != 0.) {
temp1 = *alpha * Y(j);
temp2 = *alpha * X(j);
- i__2 = j;
+ //i__2 = j;
for (i = 1; i <= j; ++i) {
A(i,j) = A(i,j) + X(i) * temp1
+ Y(i) * temp2;
/* L10: */
@@ -191,16 +191,16 @@
}
/* L20: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(jx) != 0. || Y(jy) != 0.) {
temp1 = *alpha * Y(jy);
temp2 = *alpha * X(jx);
ix = kx;
iy = ky;
- i__2 = j;
+ //i__2 = j;
for (i = 1; i <= j; ++i) {
A(i,j) = A(i,j) + X(ix) * temp1
+ Y(iy) * temp2;
ix += *incx;
@@ -217,14 +217,14 @@
/* Form A when A is stored in the lower triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(j) != 0. || Y(j) != 0.) {
temp1 = *alpha * Y(j);
temp2 = *alpha * X(j);
- i__2 = *n;
+ //i__2 = *n;
for (i = j; i <= *n; ++i) {
A(i,j) = A(i,j) + X(i) * temp1
+ Y(i) * temp2;
/* L50: */
@@ -232,16 +232,16 @@
}
/* L60: */
}
} else {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(jx) != 0. || Y(jy) != 0.) {
temp1 = *alpha * Y(jy);
temp2 = *alpha * X(jx);
ix = jx;
iy = jy;
- i__2 = *n;
+ //i__2 = *n;
for (i = j; i <= *n; ++i) {
A(i,j) = A(i,j) + X(ix) * temp1
+ Y(iy) * temp2;
ix += *incx;
Index: SuperLU_3.0/CBLAS/dtrsv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dtrsv.c
+++ SuperLU_3.0/CBLAS/dtrsv.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2;
+ //integer a_dim1, a_offset, i__1, i__2;
/* Local variables */
static integer info;
static doublereal temp;
@@ -214,16 +214,16 @@
}
}
} else {
if (*incx == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(j) != 0.) {
if (nounit) {
X(j) /= A(j,j);
}
temp = X(j);
- i__2 = *n;
+ //i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
X(i) -= temp * A(i,j);
/* L50: */
}
@@ -231,17 +231,17 @@
/* L60: */
}
} else {
jx = kx;
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
if (X(jx) != 0.) {
if (nounit) {
X(jx) /= A(j,j);
}
temp = X(jx);
ix = jx;
- i__2 = *n;
+ //i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
X(ix) -= temp * A(i,j);
/* L70: */
@@ -257,12 +257,12 @@
/* Form x := inv( A' )*x. */
if (lsame_(uplo, "U")) {
if (*incx == 1) {
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp = X(j);
- i__2 = j - 1;
+ //i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
temp -= A(i,j) * X(i);
/* L90: */
}
@@ -273,13 +273,13 @@
/* L100: */
}
} else {
jx = kx;
- i__1 = *n;
+ //i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp = X(jx);
ix = kx;
- i__2 = j - 1;
+ //i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
temp -= A(i,j) * X(ix);
ix += *incx;
/* L110: */
@@ -295,9 +295,9 @@
} else {
if (*incx == 1) {
for (j = *n; j >= 1; --j) {
temp = X(j);
- i__1 = j + 1;
+ //i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
temp -= A(i,j) * X(i);
/* L130: */
}
@@ -312,9 +312,9 @@
jx = kx;
for (j = *n; j >= 1; --j) {
temp = X(jx);
ix = kx;
- i__1 = j + 1;
+ //i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
temp -= A(i,j) * X(ix);
ix -= *incx;
/* L150: */
Index: SuperLU_3.0/CBLAS/icamax.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/icamax.c
+++ SuperLU_3.0/CBLAS/icamax.c
@@ -2,9 +2,9 @@
integer icamax_(integer *n, complex *cx, integer *incx)
{
/* System generated locals */
- integer ret_val, i__1, i__2;
+ integer ret_val; //, i__1, i__2;
real r__1, r__2;
/* Builtin functions */
double r_imag(complex *);
/* Local variables */
@@ -32,17 +32,17 @@ integer icamax_(integer *n, complex *cx,
/* code for increment not equal to 1 */
ix = 1;
smax = (r__1 = CX(1).r, dabs(r__1)) + (r__2 = r_imag(&CX(1)), dabs(r__2));
ix += *incx;
- i__1 = *n;
+ //i__1 = *n;
for (i = 2; i <= *n; ++i) {
- i__2 = ix;
+ //i__2 = ix;
if ((r__1 = CX(ix).r, dabs(r__1)) + (r__2 = r_imag(&CX(ix)), dabs(
r__2)) <= smax) {
goto L5;
}
ret_val = i;
- i__2 = ix;
+ //i__2 = ix;
smax = (r__1 = CX(ix).r, dabs(r__1)) + (r__2 = r_imag(&CX(ix)),
dabs(r__2));
L5:
ix += *incx;
@@ -51,17 +51,17 @@ L5:
return ret_val;
/* code for increment equal to 1 */
L20:
smax = (r__1 = CX(1).r, dabs(r__1)) + (r__2 = r_imag(&CX(1)), dabs(r__2));
- i__1 = *n;
+ //i__1 = *n;
for (i = 2; i <= *n; ++i) {
- i__2 = i;
+ //i__2 = i;
if ((r__1 = CX(i).r, dabs(r__1)) + (r__2 = r_imag(&CX(i)), dabs(
r__2)) <= smax) {
goto L30;
}
ret_val = i;
- i__2 = i;
+ //i__2 = i;
smax = (r__1 = CX(i).r, dabs(r__1)) + (r__2 = r_imag(&CX(i)), dabs(
r__2));
L30:
;
Index: SuperLU_3.0/CBLAS/scasum.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/scasum.c
+++ SuperLU_3.0/CBLAS/scasum.c
@@ -10,9 +10,9 @@ real scasum_(integer *n, complex *cx, in
{
/* System generated locals */
- integer i__1, i__2, i__3;
+ //integer i__1, i__2, i__3;
real ret_val, r__1, r__2;
/* Builtin functions */
double r_imag(complex *);
@@ -46,12 +46,12 @@ real scasum_(integer *n, complex *cx, in
/* code for increment not equal to 1 */
nincx = *n * *incx;
- i__1 = nincx;
- i__2 = *incx;
+ //i__1 = nincx;
+ //i__2 = *incx;
for (i = 1; *incx < 0 ? i >= nincx : i <= nincx; i += *incx) {
- i__3 = i;
+ //i__3 = i;
stemp = stemp + (r__1 = CX(i).r, dabs(r__1)) + (r__2 = r_imag(&CX(
i)), dabs(r__2));
/* L10: */
}
@@ -60,11 +60,11 @@ real scasum_(integer *n, complex *cx, in
/* code for increment equal to 1 */
L20:
- i__2 = *n;
+ //i__2 = *n;
for (i = 1; i <= *n; ++i) {
- i__1 = i;
+ //i__1 = i;
stemp = stemp + (r__1 = CX(i).r, dabs(r__1)) + (r__2 = r_imag(&CX(
i)), dabs(r__2));
/* L30: */
}
Index: SuperLU_3.0/CBLAS/cgemv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/cgemv.c
+++ SuperLU_3.0/CBLAS/cgemv.c
@@ -156,10 +156,10 @@
}
/* Quick return if possible. */
- if (*m == 0 || *n == 0 || alpha->r == 0.f && alpha->i == 0.f && (beta->r
- == 1.f && beta->i == 0.f)) {
+ if (*m == 0 || *n == 0 || (alpha->r == 0.f && alpha->i == 0.f && (beta->r
+ == 1.f && beta->i == 0.f))) {
return 0;
}
noconj = lsame_(trans, "T");
Index: SuperLU_3.0/CBLAS/cgerc.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/cgerc.c
+++ SuperLU_3.0/CBLAS/cgerc.c
@@ -126,9 +126,9 @@
}
/* Quick return if possible. */
- if (*m == 0 || *n == 0 || alpha->r == 0.f && alpha->i == 0.f) {
+ if (*m == 0 || *n == 0 || (alpha->r == 0.f && alpha->i == 0.f)) {
return 0;
}
/* Start the operations. In this version the elements of A are
Index: SuperLU_3.0/CBLAS/chemv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/chemv.c
+++ SuperLU_3.0/CBLAS/chemv.c
@@ -152,10 +152,10 @@
}
/* Quick return if possible. */
- if (*n == 0 || alpha->r == 0.f && alpha->i == 0.f && (beta->r == 1.f &&
- beta->i == 0.f)) {
+ if (*n == 0 || (alpha->r == 0.f && alpha->i == 0.f && (beta->r == 1.f &&
+ beta->i == 0.f))) {
return 0;
}
/* Set up the start points in X and Y. */
Index: SuperLU_3.0/CBLAS/cher2.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/cher2.c
+++ SuperLU_3.0/CBLAS/cher2.c
@@ -151,9 +151,9 @@
}
/* Quick return if possible. */
- if (*n == 0 || alpha->r == 0.f && alpha->i == 0.f) {
+ if (*n == 0 || (alpha->r == 0.f && alpha->i == 0.f)) {
return 0;
}
/* Set up the start points in X and Y if the increments are not both
Index: SuperLU_3.0/CBLAS/cscal.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/cscal.c
+++ SuperLU_3.0/CBLAS/cscal.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1, i__2, i__3, i__4;
+ //integer i__1, i__2, i__3, i__4;
complex q__1;
/* Local variables */
static integer i, nincx;
@@ -40,13 +40,13 @@
/* code for increment not equal to 1 */
nincx = *n * *incx;
- i__1 = nincx;
- i__2 = *incx;
+ //i__1 = nincx;
+ //i__2 = *incx;
for (i = 1; *incx < 0 ? i >= nincx : i <= nincx; i += *incx) {
- i__3 = i;
- i__4 = i;
+ //i__3 = i;
+ //i__4 = i;
q__1.r = ca->r * CX(i).r - ca->i * CX(i).i, q__1.i = ca->r * CX(
i).i + ca->i * CX(i).r;
CX(i).r = q__1.r, CX(i).i = q__1.i;
/* L10: */
@@ -55,12 +55,12 @@
/* code for increment equal to 1 */
L20:
- i__2 = *n;
+ //i__2 = *n;
for (i = 1; i <= *n; ++i) {
- i__1 = i;
- i__3 = i;
+ //i__1 = i;
+ //i__3 = i;
q__1.r = ca->r * CX(i).r - ca->i * CX(i).i, q__1.i = ca->r * CX(
i).i + ca->i * CX(i).r;
CX(i).r = q__1.r, CX(i).i = q__1.i;
/* L30: */
Index: SuperLU_3.0/CBLAS/izamax.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/izamax.c
+++ SuperLU_3.0/CBLAS/izamax.c
@@ -10,9 +10,9 @@ integer izamax_(integer *n, doublecomple
{
/* System generated locals */
- integer ret_val, i__1;
+ integer ret_val; //, i__1;
/* Local variables */
static doublereal smax;
static integer i;
@@ -48,9 +48,9 @@ integer izamax_(integer *n, doublecomple
ix = 1;
smax = dcabs1_(&ZX(1));
ix += *incx;
- i__1 = *n;
+ //i__1 = *n;
for (i = 2; i <= *n; ++i) {
if (dcabs1_(&ZX(ix)) <= smax) {
goto L5;
}
@@ -65,9 +65,9 @@ L5:
/* code for increment equal to 1 */
L20:
smax = dcabs1_(&ZX(1));
- i__1 = *n;
+ //i__1 = *n;
for (i = 2; i <= *n; ++i) {
if (dcabs1_(&ZX(i)) <= smax) {
goto L30;
}
Index: SuperLU_3.0/CBLAS/scnrm2.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/scnrm2.c
+++ SuperLU_3.0/CBLAS/scnrm2.c
@@ -10,9 +10,9 @@ real scnrm2_(integer *n, complex *x, int
{
/* System generated locals */
- integer i__1, i__2, i__3;
+ //integer i__1, i__2, i__3;
real ret_val, r__1;
/* Builtin functions */
double r_imag(complex *), sqrt(doublereal);
@@ -50,14 +50,14 @@ real scnrm2_(integer *n, complex *x, int
auxiliary routine:
CALL CLASSQ( N, X, INCX, SCALE, SSQ ) */
- i__1 = (*n - 1) * *incx + 1;
- i__2 = *incx;
+ //i__1 = (*n - 1) * *incx + 1;
+ //i__2 = *incx;
for (ix = 1; *incx < 0 ? ix >= (*n-1)**incx+1 : ix <= (*n-1)**incx+1; ix += *incx) {
- i__3 = ix;
+ //i__3 = ix;
if (X(ix).r != 0.f) {
- i__3 = ix;
+ //i__3 = ix;
temp = (r__1 = X(ix).r, dabs(r__1));
if (scale < temp) {
/* Computing 2nd power */
r__1 = scale / temp;
Index: SuperLU_3.0/CBLAS/dzasum.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dzasum.c
+++ SuperLU_3.0/CBLAS/dzasum.c
@@ -10,9 +10,9 @@ doublereal dzasum_(integer *n, doublecom
{
/* System generated locals */
- integer i__1;
+ //integer i__1;
doublereal ret_val;
/* Local variables */
static integer i;
@@ -44,9 +44,9 @@ doublereal dzasum_(integer *n, doublecom
/* code for increment not equal to 1 */
ix = 1;
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
stemp += dcabs1_(&ZX(ix));
ix += *incx;
/* L10: */
@@ -56,9 +56,9 @@ doublereal dzasum_(integer *n, doublecom
/* code for increment equal to 1 */
L20:
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
stemp += dcabs1_(&ZX(i));
/* L30: */
}
Index: SuperLU_3.0/CBLAS/dznrm2.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/dznrm2.c
+++ SuperLU_3.0/CBLAS/dznrm2.c
@@ -10,9 +10,9 @@ doublereal dznrm2_(integer *n, doublecom
{
/* System generated locals */
- integer i__1, i__2, i__3;
+ //integer i__1, i__2, i__3;
doublereal ret_val, d__1;
/* Builtin functions */
double d_imag(doublecomplex *), sqrt(doublereal);
@@ -50,14 +50,14 @@ doublereal dznrm2_(integer *n, doublecom
auxiliary routine:
CALL ZLASSQ( N, X, INCX, SCALE, SSQ ) */
- i__1 = (*n - 1) * *incx + 1;
- i__2 = *incx;
+ //i__1 = (*n - 1) * *incx + 1;
+ //i__2 = *incx;
for (ix = 1; *incx < 0 ? ix >= (*n-1)**incx+1 : ix <= (*n-1)**incx+1; ix += *incx) {
- i__3 = ix;
+ //i__3 = ix;
if (X(ix).r != 0.) {
- i__3 = ix;
+ //i__3 = ix;
temp = (d__1 = X(ix).r, abs(d__1));
if (scale < temp) {
/* Computing 2nd power */
d__1 = scale / temp;
Index: SuperLU_3.0/CBLAS/zaxpy.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zaxpy.c
+++ SuperLU_3.0/CBLAS/zaxpy.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1, i__2, i__3, i__4;
+ //integer i__1, i__2, i__3, i__4;
doublecomplex z__1, z__2;
/* Local variables */
static integer i;
@@ -52,13 +52,13 @@
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
- i__3 = iy;
- i__4 = ix;
+ //i__2 = iy;
+ //i__3 = iy;
+ //i__4 = ix;
z__2.r = za->r * ZX(ix).r - za->i * ZX(ix).i, z__2.i = za->r * ZX(
ix).i + za->i * ZX(ix).r;
z__1.r = ZY(iy).r + z__2.r, z__1.i = ZY(iy).i + z__2.i;
ZY(iy).r = z__1.r, ZY(iy).i = z__1.i;
@@ -70,13 +70,13 @@
/* code for both increments equal to 1 */
L20:
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
- i__3 = i;
- i__4 = i;
+ //i__2 = i;
+ //i__3 = i;
+ //i__4 = i;
z__2.r = za->r * ZX(i).r - za->i * ZX(i).i, z__2.i = za->r * ZX(
i).i + za->i * ZX(i).r;
z__1.r = ZY(i).r + z__2.r, z__1.i = ZY(i).i + z__2.i;
ZY(i).r = z__1.r, ZY(i).i = z__1.i;
Index: SuperLU_3.0/CBLAS/zcopy.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zcopy.c
+++ SuperLU_3.0/CBLAS/zcopy.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1, i__2, i__3;
+ //integer i__1, i__2, i__3;
/* Local variables */
static integer i, ix, iy;
@@ -47,12 +47,12 @@
}
if (*incy < 0) {
iy = (-(*n) + 1) * *incy + 1;
}
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
- i__3 = ix;
+ //i__2 = iy;
+ //i__3 = ix;
ZY(iy).r = ZX(ix).r, ZY(iy).i = ZX(ix).i;
ix += *incx;
iy += *incy;
/* L10: */
@@ -61,12 +61,12 @@
/* code for both increments equal to 1 */
L20:
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
- i__3 = i;
+ //i__2 = i;
+ //i__3 = i;
ZY(i).r = ZX(i).r, ZY(i).i = ZX(i).i;
/* L30: */
}
return 0;
Index: SuperLU_3.0/CBLAS/zgemv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zgemv.c
+++ SuperLU_3.0/CBLAS/zgemv.c
@@ -157,10 +157,10 @@
}
/* Quick return if possible. */
- if (*m == 0 || *n == 0 || alpha->r == 0. && alpha->i == 0. && (beta->r ==
- 1. && beta->i == 0.)) {
+ if (*m == 0 || *n == 0 || (alpha->r == 0. && alpha->i == 0. && (beta->r ==
+ 1. && beta->i == 0.))) {
return 0;
}
noconj = lsame_(trans, "T");
Index: SuperLU_3.0/CBLAS/zgerc.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zgerc.c
+++ SuperLU_3.0/CBLAS/zgerc.c
@@ -127,9 +127,9 @@
}
/* Quick return if possible. */
- if (*m == 0 || *n == 0 || alpha->r == 0. && alpha->i == 0.) {
+ if (*m == 0 || *n == 0 || (alpha->r == 0. && alpha->i == 0.)) {
return 0;
}
/* Start the operations. In this version the elements of A are
Index: SuperLU_3.0/CBLAS/zhemv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zhemv.c
+++ SuperLU_3.0/CBLAS/zhemv.c
@@ -152,10 +152,10 @@
}
/* Quick return if possible. */
- if (*n == 0 || alpha->r == 0. && alpha->i == 0. && (beta->r == 1. &&
- beta->i == 0.)) {
+ if (*n == 0 || (alpha->r == 0. && alpha->i == 0. && (beta->r == 1. &&
+ beta->i == 0.))) {
return 0;
}
/* Set up the start points in X and Y. */
Index: SuperLU_3.0/CBLAS/zher2.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zher2.c
+++ SuperLU_3.0/CBLAS/zher2.c
@@ -152,9 +152,9 @@
}
/* Quick return if possible. */
- if (*n == 0 || alpha->r == 0. && alpha->i == 0.) {
+ if (*n == 0 || (alpha->r == 0. && alpha->i == 0.)) {
return 0;
}
/* Set up the start points in X and Y if the increments are not both
Index: SuperLU_3.0/CBLAS/zscal.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zscal.c
+++ SuperLU_3.0/CBLAS/zscal.c
@@ -11,9 +11,9 @@
{
/* System generated locals */
- integer i__1, i__2, i__3;
+ //integer i__1, i__2, i__3;
doublecomplex z__1;
/* Local variables */
static integer i, ix;
@@ -40,12 +40,12 @@
/* code for increment not equal to 1 */
ix = 1;
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = ix;
- i__3 = ix;
+ //i__2 = ix;
+ //i__3 = ix;
z__1.r = za->r * ZX(ix).r - za->i * ZX(ix).i, z__1.i = za->r * ZX(
ix).i + za->i * ZX(ix).r;
ZX(ix).r = z__1.r, ZX(ix).i = z__1.i;
ix += *incx;
@@ -55,12 +55,12 @@
/* code for increment equal to 1 */
L20:
- i__1 = *n;
+ //i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
- i__3 = i;
+ //i__2 = i;
+ //i__3 = i;
z__1.r = za->r * ZX(i).r - za->i * ZX(i).i, z__1.i = za->r * ZX(
i).i + za->i * ZX(i).r;
ZX(i).r = z__1.r, ZX(i).i = z__1.i;
/* L30: */
