pool/superlu
SHA256
2
0
Fork 1
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity
4a58f37a6b
superlu/superlu-initialize.diff

2669 lines
76 KiB
Diff
Raw Blame History

Index: SuperLU_3.0/CBLAS/cgemv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/cgemv.c
+++ SuperLU_3.0/CBLAS/cgemv.c
@@ -12,9 +12,8 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5;
complex q__1, q__2, q__3;
/* Builtin functions */
void r_cnjg(complex *, complex *);
@@ -194,19 +193,14 @@
if (beta->r != 1.f || beta->i != 0.f) {
if (*incy == 1) {
if (beta->r == 0.f && beta->i == 0.f) {
- i__1 = leny;
for (i = 1; i <= leny; ++i) {
- i__2 = i;
Y(i).r = 0.f, Y(i).i = 0.f;
/* L10: */
}
} else {
- i__1 = leny;
for (i = 1; i <= leny; ++i) {
- i__2 = i;
- i__3 = i;
q__1.r = beta->r * Y(i).r - beta->i * Y(i).i,
q__1.i = beta->r * Y(i).i + beta->i * Y(i)
.r;
Y(i).r = q__1.r, Y(i).i = q__1.i;
@@ -215,20 +209,15 @@
}
} else {
iy = ky;
if (beta->r == 0.f && beta->i == 0.f) {
- i__1 = leny;
for (i = 1; i <= leny; ++i) {
- i__2 = iy;
Y(iy).r = 0.f, Y(iy).i = 0.f;
iy += *incy;
/* L30: */
}
} else {
- i__1 = leny;
for (i = 1; i <= leny; ++i) {
- i__2 = iy;
- i__3 = iy;
q__1.r = beta->r * Y(iy).r - beta->i * Y(iy).i,
q__1.i = beta->r * Y(iy).i + beta->i * Y(iy)
.r;
Y(iy).r = q__1.r, Y(iy).i = q__1.i;
@@ -246,22 +235,15 @@
/* Form y := alpha*A*x + y. */
jx = kx;
if (*incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
if (X(jx).r != 0.f || X(jx).i != 0.f) {
- i__2 = jx;
q__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i,
q__1.i = alpha->r * X(jx).i + alpha->i * X(jx)
.r;
temp.r = q__1.r, temp.i = q__1.i;
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
q__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
q__2.i = temp.r * A(i,j).i + temp.i * A(i,j)
.r;
q__1.r = Y(i).r + q__2.r, q__1.i = Y(i).i +
@@ -273,23 +255,16 @@
jx += *incx;
/* L60: */
}
} else {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
if (X(jx).r != 0.f || X(jx).i != 0.f) {
- i__2 = jx;
q__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i,
q__1.i = alpha->r * X(jx).i + alpha->i * X(jx)
.r;
temp.r = q__1.r, temp.i = q__1.i;
iy = ky;
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = iy;
- i__4 = iy;
- i__5 = i + j * a_dim1;
q__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
q__2.i = temp.r * A(i,j).i + temp.i * A(i,j)
.r;
q__1.r = Y(iy).r + q__2.r, q__1.i = Y(iy).i +
@@ -309,55 +284,43 @@
*/
jy = ky;
if (*incx == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp.r = 0.f, temp.i = 0.f;
if (noconj) {
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i;
q__2.r = A(i,j).r * X(i).r - A(i,j).i * X(i)
.i, q__2.i = A(i,j).r * X(i).i + A(i,j)
.i * X(i).r;
q__1.r = temp.r + q__2.r, q__1.i = temp.i + q__2.i;
temp.r = q__1.r, temp.i = q__1.i;
/* L90: */
}
} else {
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
r_cnjg(&q__3, &A(i,j));
- i__3 = i;
q__2.r = q__3.r * X(i).r - q__3.i * X(i).i,
q__2.i = q__3.r * X(i).i + q__3.i * X(i)
.r;
q__1.r = temp.r + q__2.r, q__1.i = temp.i + q__2.i;
temp.r = q__1.r, temp.i = q__1.i;
/* L100: */
}
}
- i__2 = jy;
- i__3 = jy;
q__2.r = alpha->r * temp.r - alpha->i * temp.i, q__2.i =
alpha->r * temp.i + alpha->i * temp.r;
q__1.r = Y(jy).r + q__2.r, q__1.i = Y(jy).i + q__2.i;
Y(jy).r = q__1.r, Y(jy).i = q__1.i;
jy += *incy;
/* L110: */
}
} else {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp.r = 0.f, temp.i = 0.f;
ix = kx;
if (noconj) {
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = ix;
q__2.r = A(i,j).r * X(ix).r - A(i,j).i * X(ix)
.i, q__2.i = A(i,j).r * X(ix).i + A(i,j)
.i * X(ix).r;
q__1.r = temp.r + q__2.r, q__1.i = temp.i + q__2.i;
@@ -365,12 +328,10 @@
ix += *incx;
/* L120: */
}
} else {
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
r_cnjg(&q__3, &A(i,j));
- i__3 = ix;
q__2.r = q__3.r * X(ix).r - q__3.i * X(ix).i,
q__2.i = q__3.r * X(ix).i + q__3.i * X(ix)
.r;
q__1.r = temp.r + q__2.r, q__1.i = temp.i + q__2.i;
@@ -378,10 +339,8 @@
ix += *incx;
/* L130: */
}
}
- i__2 = jy;
- i__3 = jy;
q__2.r = alpha->r * temp.r - alpha->i * temp.i, q__2.i =
alpha->r * temp.i + alpha->i * temp.r;
q__1.r = Y(jy).r + q__2.r, q__1.i = Y(jy).i + q__2.i;
Y(jy).r = q__1.r, Y(jy).i = q__1.i;
Index: SuperLU_3.0/CBLAS/cgerc.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/cgerc.c
+++ SuperLU_3.0/CBLAS/cgerc.c
@@ -11,9 +11,8 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5;
complex q__1, q__2;
/* Builtin functions */
void r_cnjg(complex *, complex *);
@@ -140,21 +139,15 @@
} else {
jy = 1 - (*n - 1) * *incy;
}
if (*incx == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jy;
if (Y(jy).r != 0.f || Y(jy).i != 0.f) {
r_cnjg(&q__2, &Y(jy));
q__1.r = alpha->r * q__2.r - alpha->i * q__2.i, q__1.i =
alpha->r * q__2.i + alpha->i * q__2.r;
temp.r = q__1.r, temp.i = q__1.i;
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = i;
q__2.r = X(i).r * temp.r - X(i).i * temp.i, q__2.i =
X(i).r * temp.i + X(i).i * temp.r;
q__1.r = A(i,j).r + q__2.r, q__1.i = A(i,j).i + q__2.i;
A(i,j).r = q__1.r, A(i,j).i = q__1.i;
@@ -169,22 +162,16 @@
kx = 1;
} else {
kx = 1 - (*m - 1) * *incx;
}
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jy;
if (Y(jy).r != 0.f || Y(jy).i != 0.f) {
r_cnjg(&q__2, &Y(jy));
q__1.r = alpha->r * q__2.r - alpha->i * q__2.i, q__1.i =
alpha->r * q__2.i + alpha->i * q__2.r;
temp.r = q__1.r, temp.i = q__1.i;
ix = kx;
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = ix;
q__2.r = X(ix).r * temp.r - X(ix).i * temp.i, q__2.i =
X(ix).r * temp.i + X(ix).i * temp.r;
q__1.r = A(i,j).r + q__2.r, q__1.i = A(i,j).i + q__2.i;
A(i,j).r = q__1.r, A(i,j).i = q__1.i;
Index: SuperLU_3.0/CBLAS/chemv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/chemv.c
+++ SuperLU_3.0/CBLAS/chemv.c
@@ -12,9 +12,8 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5;
doublereal d__1;
complex q__1, q__2, q__3, q__4;
/* Builtin functions */
@@ -180,19 +179,14 @@
if (beta->r != 1.f || beta->i != 0.f) {
if (*incy == 1) {
if (beta->r == 0.f && beta->i == 0.f) {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
Y(i).r = 0.f, Y(i).i = 0.f;
/* L10: */
}
} else {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
- i__3 = i;
q__1.r = beta->r * Y(i).r - beta->i * Y(i).i,
q__1.i = beta->r * Y(i).i + beta->i * Y(i)
.r;
Y(i).r = q__1.r, Y(i).i = q__1.i;
@@ -201,20 +195,15 @@
}
} else {
iy = ky;
if (beta->r == 0.f && beta->i == 0.f) {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
Y(iy).r = 0.f, Y(iy).i = 0.f;
iy += *incy;
/* L30: */
}
} else {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
- i__3 = iy;
q__1.r = beta->r * Y(iy).r - beta->i * Y(iy).i,
q__1.i = beta->r * Y(iy).i + beta->i * Y(iy)
.r;
Y(iy).r = q__1.r, Y(iy).i = q__1.i;
@@ -231,36 +220,26 @@
/* Form y when A is stored in upper triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
q__1.r = alpha->r * X(j).r - alpha->i * X(j).i, q__1.i =
alpha->r * X(j).i + alpha->i * X(j).r;
temp1.r = q__1.r, temp1.i = q__1.i;
temp2.r = 0.f, temp2.i = 0.f;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
q__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
q__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
q__1.r = Y(i).r + q__2.r, q__1.i = Y(i).i + q__2.i;
Y(i).r = q__1.r, Y(i).i = q__1.i;
r_cnjg(&q__3, &A(i,j));
- i__3 = i;
q__2.r = q__3.r * X(i).r - q__3.i * X(i).i, q__2.i =
q__3.r * X(i).i + q__3.i * X(i).r;
q__1.r = temp2.r + q__2.r, q__1.i = temp2.i + q__2.i;
temp2.r = q__1.r, temp2.i = q__1.i;
/* L50: */
}
- i__2 = j;
- i__3 = j;
- i__4 = j + j * a_dim1;
d__1 = A(j,j).r;
q__3.r = d__1 * temp1.r, q__3.i = d__1 * temp1.i;
q__2.r = Y(j).r + q__3.r, q__2.i = Y(j).i + q__3.i;
q__4.r = alpha->r * temp2.r - alpha->i * temp2.i, q__4.i =
@@ -271,40 +250,30 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
q__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i, q__1.i =
alpha->r * X(jx).i + alpha->i * X(jx).r;
temp1.r = q__1.r, temp1.i = q__1.i;
temp2.r = 0.f, temp2.i = 0.f;
ix = kx;
iy = ky;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = iy;
- i__4 = iy;
- i__5 = i + j * a_dim1;
q__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
q__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
q__1.r = Y(iy).r + q__2.r, q__1.i = Y(iy).i + q__2.i;
Y(iy).r = q__1.r, Y(iy).i = q__1.i;
r_cnjg(&q__3, &A(i,j));
- i__3 = ix;
q__2.r = q__3.r * X(ix).r - q__3.i * X(ix).i, q__2.i =
q__3.r * X(ix).i + q__3.i * X(ix).r;
q__1.r = temp2.r + q__2.r, q__1.i = temp2.i + q__2.i;
temp2.r = q__1.r, temp2.i = q__1.i;
ix += *incx;
iy += *incy;
/* L70: */
}
- i__2 = jy;
- i__3 = jy;
- i__4 = j + j * a_dim1;
d__1 = A(j,j).r;
q__3.r = d__1 * temp1.r, q__3.i = d__1 * temp1.i;
q__2.r = Y(jy).r + q__3.r, q__2.i = Y(jy).i + q__3.i;
q__4.r = alpha->r * temp2.r - alpha->i * temp2.i, q__4.i =
@@ -320,42 +289,30 @@
/* Form y when A is stored in lower triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
q__1.r = alpha->r * X(j).r - alpha->i * X(j).i, q__1.i =
alpha->r * X(j).i + alpha->i * X(j).r;
temp1.r = q__1.r, temp1.i = q__1.i;
temp2.r = 0.f, temp2.i = 0.f;
- i__2 = j;
- i__3 = j;
- i__4 = j + j * a_dim1;
d__1 = A(j,j).r;
q__2.r = d__1 * temp1.r, q__2.i = d__1 * temp1.i;
q__1.r = Y(j).r + q__2.r, q__1.i = Y(j).i + q__2.i;
Y(j).r = q__1.r, Y(j).i = q__1.i;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
q__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
q__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
q__1.r = Y(i).r + q__2.r, q__1.i = Y(i).i + q__2.i;
Y(i).r = q__1.r, Y(i).i = q__1.i;
r_cnjg(&q__3, &A(i,j));
- i__3 = i;
q__2.r = q__3.r * X(i).r - q__3.i * X(i).i, q__2.i =
q__3.r * X(i).i + q__3.i * X(i).r;
q__1.r = temp2.r + q__2.r, q__1.i = temp2.i + q__2.i;
temp2.r = q__1.r, temp2.i = q__1.i;
/* L90: */
}
- i__2 = j;
- i__3 = j;
q__2.r = alpha->r * temp2.r - alpha->i * temp2.i, q__2.i =
alpha->r * temp2.i + alpha->i * temp2.r;
q__1.r = Y(j).r + q__2.r, q__1.i = Y(j).i + q__2.i;
Y(j).r = q__1.r, Y(j).i = q__1.i;
@@ -363,46 +320,34 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
q__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i, q__1.i =
alpha->r * X(jx).i + alpha->i * X(jx).r;
temp1.r = q__1.r, temp1.i = q__1.i;
temp2.r = 0.f, temp2.i = 0.f;
- i__2 = jy;
- i__3 = jy;
- i__4 = j + j * a_dim1;
d__1 = A(j,j).r;
q__2.r = d__1 * temp1.r, q__2.i = d__1 * temp1.i;
q__1.r = Y(jy).r + q__2.r, q__1.i = Y(jy).i + q__2.i;
Y(jy).r = q__1.r, Y(jy).i = q__1.i;
ix = jx;
iy = jy;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
iy += *incy;
- i__3 = iy;
- i__4 = iy;
- i__5 = i + j * a_dim1;
q__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
q__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
q__1.r = Y(iy).r + q__2.r, q__1.i = Y(iy).i + q__2.i;
Y(iy).r = q__1.r, Y(iy).i = q__1.i;
r_cnjg(&q__3, &A(i,j));
- i__3 = ix;
q__2.r = q__3.r * X(ix).r - q__3.i * X(ix).i, q__2.i =
q__3.r * X(ix).i + q__3.i * X(ix).r;
q__1.r = temp2.r + q__2.r, q__1.i = temp2.i + q__2.i;
temp2.r = q__1.r, temp2.i = q__1.i;
/* L110: */
}
- i__2 = jy;
- i__3 = jy;
q__2.r = alpha->r * temp2.r - alpha->i * temp2.i, q__2.i =
alpha->r * temp2.i + alpha->i * temp2.r;
q__1.r = Y(jy).r + q__2.r, q__1.i = Y(jy).i + q__2.i;
Y(jy).r = q__1.r, Y(jy).i = q__1.i;
Index: SuperLU_3.0/CBLAS/cher2.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/cher2.c
+++ SuperLU_3.0/CBLAS/cher2.c
@@ -11,9 +11,8 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5, i__6;
doublereal d__1;
complex q__1, q__2, q__3, q__4;
/* Builtin functions */
@@ -184,93 +183,69 @@
/* Form A when A is stored in the upper triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
- i__3 = j;
if (X(j).r != 0.f || X(j).i != 0.f || (Y(j).r != 0.f
|| Y(j).i != 0.f)) {
r_cnjg(&q__2, &Y(j));
q__1.r = alpha->r * q__2.r - alpha->i * q__2.i, q__1.i =
alpha->r * q__2.i + alpha->i * q__2.r;
temp1.r = q__1.r, temp1.i = q__1.i;
- i__2 = j;
q__2.r = alpha->r * X(j).r - alpha->i * X(j).i,
q__2.i = alpha->r * X(j).i + alpha->i * X(j)
.r;
r_cnjg(&q__1, &q__2);
temp2.r = q__1.r, temp2.i = q__1.i;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = i;
q__3.r = X(i).r * temp1.r - X(i).i * temp1.i,
q__3.i = X(i).r * temp1.i + X(i).i *
temp1.r;
q__2.r = A(i,j).r + q__3.r, q__2.i = A(i,j).i +
q__3.i;
- i__6 = i;
q__4.r = Y(i).r * temp2.r - Y(i).i * temp2.i,
q__4.i = Y(i).r * temp2.i + Y(i).i *
temp2.r;
q__1.r = q__2.r + q__4.r, q__1.i = q__2.i + q__4.i;
A(i,j).r = q__1.r, A(i,j).i = q__1.i;
/* L10: */
}
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
- i__4 = j;
q__2.r = X(j).r * temp1.r - X(j).i * temp1.i,
q__2.i = X(j).r * temp1.i + X(j).i *
temp1.r;
- i__5 = j;
q__3.r = Y(j).r * temp2.r - Y(j).i * temp2.i,
q__3.i = Y(j).r * temp2.i + Y(j).i *
temp2.r;
q__1.r = q__2.r + q__3.r, q__1.i = q__2.i + q__3.i;
d__1 = A(j,j).r + q__1.r;
A(j,j).r = d__1, A(j,j).i = 0.f;
} else {
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
d__1 = A(j,j).r;
A(j,j).r = d__1, A(j,j).i = 0.f;
}
/* L20: */
}
} else {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
- i__3 = jy;
if (X(jx).r != 0.f || X(jx).i != 0.f || (Y(jy).r != 0.f
|| Y(jy).i != 0.f)) {
r_cnjg(&q__2, &Y(jy));
q__1.r = alpha->r * q__2.r - alpha->i * q__2.i, q__1.i =
alpha->r * q__2.i + alpha->i * q__2.r;
temp1.r = q__1.r, temp1.i = q__1.i;
- i__2 = jx;
q__2.r = alpha->r * X(jx).r - alpha->i * X(jx).i,
q__2.i = alpha->r * X(jx).i + alpha->i * X(jx)
.r;
r_cnjg(&q__1, &q__2);
temp2.r = q__1.r, temp2.i = q__1.i;
ix = kx;
iy = ky;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = ix;
q__3.r = X(ix).r * temp1.r - X(ix).i * temp1.i,
q__3.i = X(ix).r * temp1.i + X(ix).i *
temp1.r;
q__2.r = A(i,j).r + q__3.r, q__2.i = A(i,j).i +
q__3.i;
- i__6 = iy;
q__4.r = Y(iy).r * temp2.r - Y(iy).i * temp2.i,
q__4.i = Y(iy).r * temp2.i + Y(iy).i *
temp2.r;
q__1.r = q__2.r + q__4.r, q__1.i = q__2.i + q__4.i;
@@ -278,24 +253,18 @@
ix += *incx;
iy += *incy;
/* L30: */
}
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
- i__4 = jx;
q__2.r = X(jx).r * temp1.r - X(jx).i * temp1.i,
q__2.i = X(jx).r * temp1.i + X(jx).i *
temp1.r;
- i__5 = jy;
q__3.r = Y(jy).r * temp2.r - Y(jy).i * temp2.i,
q__3.i = Y(jy).r * temp2.i + Y(jy).i *
temp2.r;
q__1.r = q__2.r + q__3.r, q__1.i = q__2.i + q__3.i;
d__1 = A(j,j).r + q__1.r;
A(j,j).r = d__1, A(j,j).i = 0.f;
} else {
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
d__1 = A(j,j).r;
A(j,j).r = d__1, A(j,j).i = 0.f;
}
jx += *incx;
@@ -307,118 +276,88 @@
/* Form A when A is stored in the lower triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
- i__3 = j;
if (X(j).r != 0.f || X(j).i != 0.f || (Y(j).r != 0.f
|| Y(j).i != 0.f)) {
r_cnjg(&q__2, &Y(j));
q__1.r = alpha->r * q__2.r - alpha->i * q__2.i, q__1.i =
alpha->r * q__2.i + alpha->i * q__2.r;
temp1.r = q__1.r, temp1.i = q__1.i;
- i__2 = j;
q__2.r = alpha->r * X(j).r - alpha->i * X(j).i,
q__2.i = alpha->r * X(j).i + alpha->i * X(j)
.r;
r_cnjg(&q__1, &q__2);
temp2.r = q__1.r, temp2.i = q__1.i;
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
- i__4 = j;
q__2.r = X(j).r * temp1.r - X(j).i * temp1.i,
q__2.i = X(j).r * temp1.i + X(j).i *
temp1.r;
- i__5 = j;
q__3.r = Y(j).r * temp2.r - Y(j).i * temp2.i,
q__3.i = Y(j).r * temp2.i + Y(j).i *
temp2.r;
q__1.r = q__2.r + q__3.r, q__1.i = q__2.i + q__3.i;
d__1 = A(j,j).r + q__1.r;
A(j,j).r = d__1, A(j,j).i = 0.f;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = i;
q__3.r = X(i).r * temp1.r - X(i).i * temp1.i,
q__3.i = X(i).r * temp1.i + X(i).i *
temp1.r;
q__2.r = A(i,j).r + q__3.r, q__2.i = A(i,j).i +
q__3.i;
- i__6 = i;
q__4.r = Y(i).r * temp2.r - Y(i).i * temp2.i,
q__4.i = Y(i).r * temp2.i + Y(i).i *
temp2.r;
q__1.r = q__2.r + q__4.r, q__1.i = q__2.i + q__4.i;
A(i,j).r = q__1.r, A(i,j).i = q__1.i;
/* L50: */
}
} else {
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
d__1 = A(j,j).r;
A(j,j).r = d__1, A(j,j).i = 0.f;
}
/* L60: */
}
} else {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
- i__3 = jy;
if (X(jx).r != 0.f || X(jx).i != 0.f || (Y(jy).r != 0.f
|| Y(jy).i != 0.f)) {
r_cnjg(&q__2, &Y(jy));
q__1.r = alpha->r * q__2.r - alpha->i * q__2.i, q__1.i =
alpha->r * q__2.i + alpha->i * q__2.r;
temp1.r = q__1.r, temp1.i = q__1.i;
- i__2 = jx;
q__2.r = alpha->r * X(jx).r - alpha->i * X(jx).i,
q__2.i = alpha->r * X(jx).i + alpha->i * X(jx)
.r;
r_cnjg(&q__1, &q__2);
temp2.r = q__1.r, temp2.i = q__1.i;
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
- i__4 = jx;
q__2.r = X(jx).r * temp1.r - X(jx).i * temp1.i,
q__2.i = X(jx).r * temp1.i + X(jx).i *
temp1.r;
- i__5 = jy;
q__3.r = Y(jy).r * temp2.r - Y(jy).i * temp2.i,
q__3.i = Y(jy).r * temp2.i + Y(jy).i *
temp2.r;
q__1.r = q__2.r + q__3.r, q__1.i = q__2.i + q__3.i;
d__1 = A(j,j).r + q__1.r;
A(j,j).r = d__1, A(j,j).i = 0.f;
ix = jx;
iy = jy;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
iy += *incy;
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = ix;
q__3.r = X(ix).r * temp1.r - X(ix).i * temp1.i,
q__3.i = X(ix).r * temp1.i + X(ix).i *
temp1.r;
q__2.r = A(i,j).r + q__3.r, q__2.i = A(i,j).i +
q__3.i;
- i__6 = iy;
q__4.r = Y(iy).r * temp2.r - Y(iy).i * temp2.i,
q__4.i = Y(iy).r * temp2.i + Y(iy).i *
temp2.r;
q__1.r = q__2.r + q__4.r, q__1.i = q__2.i + q__4.i;
A(i,j).r = q__1.r, A(i,j).i = q__1.i;
/* L70: */
}
} else {
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
d__1 = A(j,j).r;
A(j,j).r = d__1, A(j,j).i = 0.f;
}
jx += *incx;
Index: SuperLU_3.0/CBLAS/ctrsv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/ctrsv.c
+++ SuperLU_3.0/CBLAS/ctrsv.c
@@ -11,9 +11,8 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5;
complex q__1, q__2, q__3;
/* Builtin functions */
void c_div(complex *, complex *, complex *), r_cnjg(complex *, complex *);
@@ -186,21 +185,15 @@
if (lsame_(uplo, "U")) {
if (*incx == 1) {
for (j = *n; j >= 1; --j) {
- i__1 = j;
if (X(j).r != 0.f || X(j).i != 0.f) {
if (nounit) {
- i__1 = j;
c_div(&q__1, &X(j), &A(j,j));
X(j).r = q__1.r, X(j).i = q__1.i;
}
- i__1 = j;
temp.r = X(j).r, temp.i = X(j).i;
for (i = j - 1; i >= 1; --i) {
- i__1 = i;
- i__2 = i;
- i__3 = i + j * a_dim1;
q__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
q__2.i = temp.r * A(i,j).i + temp.i * A(i,j).r;
q__1.r = X(i).r - q__2.r, q__1.i = X(i).i -
q__2.i;
@@ -212,23 +205,17 @@
}
} else {
jx = kx + (*n - 1) * *incx;
for (j = *n; j >= 1; --j) {
- i__1 = jx;
if (X(jx).r != 0.f || X(jx).i != 0.f) {
if (nounit) {
- i__1 = jx;
c_div(&q__1, &X(jx), &A(j,j));
X(jx).r = q__1.r, X(jx).i = q__1.i;
}
- i__1 = jx;
temp.r = X(jx).r, temp.i = X(jx).i;
ix = jx;
for (i = j - 1; i >= 1; --i) {
ix -= *incx;
- i__1 = ix;
- i__2 = ix;
- i__3 = i + j * a_dim1;
q__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
q__2.i = temp.r * A(i,j).i + temp.i * A(i,j).r;
q__1.r = X(ix).r - q__2.r, q__1.i = X(ix).i -
q__2.i;
@@ -241,24 +228,16 @@
}
}
} else {
if (*incx == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
if (X(j).r != 0.f || X(j).i != 0.f) {
if (nounit) {
- i__2 = j;
c_div(&q__1, &X(j), &A(j,j));
X(j).r = q__1.r, X(j).i = q__1.i;
}
- i__2 = j;
temp.r = X(j).r, temp.i = X(j).i;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
q__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
q__2.i = temp.r * A(i,j).i + temp.i * A(i,j).r;
q__1.r = X(i).r - q__2.r, q__1.i = X(i).i -
q__2.i;
@@ -269,26 +248,18 @@
/* L60: */
}
} else {
jx = kx;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
if (X(jx).r != 0.f || X(jx).i != 0.f) {
if (nounit) {
- i__2 = jx;
c_div(&q__1, &X(jx), &A(j,j));
X(jx).r = q__1.r, X(jx).i = q__1.i;
}
- i__2 = jx;
temp.r = X(jx).r, temp.i = X(jx).i;
ix = jx;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
- i__3 = ix;
- i__4 = ix;
- i__5 = i + j * a_dim1;
q__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
q__2.i = temp.r * A(i,j).i + temp.i * A(i,j).r;
q__1.r = X(ix).r - q__2.r, q__1.i = X(ix).i -
q__2.i;
@@ -306,17 +277,12 @@
/* Form x := inv( A' )*x or x := inv( conjg( A' ) )*x. */
if (lsame_(uplo, "U")) {
if (*incx == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
temp.r = X(j).r, temp.i = X(j).i;
if (noconj) {
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i;
q__2.r = A(i,j).r * X(i).r - A(i,j).i * X(
i).i, q__2.i = A(i,j).r * X(i).i +
A(i,j).i * X(i).r;
q__1.r = temp.r - q__2.r, q__1.i = temp.i -
@@ -328,12 +294,10 @@
c_div(&q__1, &temp, &A(j,j));
temp.r = q__1.r, temp.i = q__1.i;
}
} else {
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
r_cnjg(&q__3, &A(i,j));
- i__3 = i;
q__2.r = q__3.r * X(i).r - q__3.i * X(i).i,
q__2.i = q__3.r * X(i).i + q__3.i * X(
i).r;
q__1.r = temp.r - q__2.r, q__1.i = temp.i -
@@ -346,24 +310,18 @@
c_div(&q__1, &temp, &q__2);
temp.r = q__1.r, temp.i = q__1.i;
}
}
- i__2 = j;
X(j).r = temp.r, X(j).i = temp.i;
/* L110: */
}
} else {
jx = kx;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
ix = kx;
- i__2 = jx;
temp.r = X(jx).r, temp.i = X(jx).i;
if (noconj) {
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = ix;
q__2.r = A(i,j).r * X(ix).r - A(i,j).i * X(
ix).i, q__2.i = A(i,j).r * X(ix).i +
A(i,j).i * X(ix).r;
q__1.r = temp.r - q__2.r, q__1.i = temp.i -
@@ -376,12 +334,10 @@
c_div(&q__1, &temp, &A(j,j));
temp.r = q__1.r, temp.i = q__1.i;
}
} else {
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
r_cnjg(&q__3, &A(i,j));
- i__3 = ix;
q__2.r = q__3.r * X(ix).r - q__3.i * X(ix).i,
q__2.i = q__3.r * X(ix).i + q__3.i * X(
ix).r;
q__1.r = temp.r - q__2.r, q__1.i = temp.i -
@@ -395,24 +351,19 @@
c_div(&q__1, &temp, &q__2);
temp.r = q__1.r, temp.i = q__1.i;
}
}
- i__2 = jx;
X(jx).r = temp.r, X(jx).i = temp.i;
jx += *incx;
/* L140: */
}
}
} else {
if (*incx == 1) {
for (j = *n; j >= 1; --j) {
- i__1 = j;
temp.r = X(j).r, temp.i = X(j).i;
if (noconj) {
- i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
- i__2 = i + j * a_dim1;
- i__3 = i;
q__2.r = A(i,j).r * X(i).r - A(i,j).i * X(
i).i, q__2.i = A(i,j).r * X(i).i +
A(i,j).i * X(i).r;
q__1.r = temp.r - q__2.r, q__1.i = temp.i -
@@ -424,12 +375,10 @@
c_div(&q__1, &temp, &A(j,j));
temp.r = q__1.r, temp.i = q__1.i;
}
} else {
- i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
r_cnjg(&q__3, &A(i,j));
- i__2 = i;
q__2.r = q__3.r * X(i).r - q__3.i * X(i).i,
q__2.i = q__3.r * X(i).i + q__3.i * X(
i).r;
q__1.r = temp.r - q__2.r, q__1.i = temp.i -
@@ -442,24 +391,19 @@
c_div(&q__1, &temp, &q__2);
temp.r = q__1.r, temp.i = q__1.i;
}
}
- i__1 = j;
X(j).r = temp.r, X(j).i = temp.i;
/* L170: */
}
} else {
kx += (*n - 1) * *incx;
jx = kx;
for (j = *n; j >= 1; --j) {
ix = kx;
- i__1 = jx;
temp.r = X(jx).r, temp.i = X(jx).i;
if (noconj) {
- i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
- i__2 = i + j * a_dim1;
- i__3 = ix;
q__2.r = A(i,j).r * X(ix).r - A(i,j).i * X(
ix).i, q__2.i = A(i,j).r * X(ix).i +
A(i,j).i * X(ix).r;
q__1.r = temp.r - q__2.r, q__1.i = temp.i -
@@ -472,12 +416,10 @@
c_div(&q__1, &temp, &A(j,j));
temp.r = q__1.r, temp.i = q__1.i;
}
} else {
- i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
r_cnjg(&q__3, &A(i,j));
- i__2 = ix;
q__2.r = q__3.r * X(ix).r - q__3.i * X(ix).i,
q__2.i = q__3.r * X(ix).i + q__3.i * X(
ix).r;
q__1.r = temp.r - q__2.r, q__1.i = temp.i -
@@ -491,9 +433,8 @@
c_div(&q__1, &temp, &q__2);
temp.r = q__1.r, temp.i = q__1.i;
}
}
- i__1 = jx;
X(jx).r = temp.r, X(jx).i = temp.i;
jx -= *incx;
/* L200: */
}
Index: SuperLU_3.0/CBLAS/zgemv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zgemv.c
+++ SuperLU_3.0/CBLAS/zgemv.c
@@ -13,9 +13,8 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5;
doublecomplex z__1, z__2, z__3;
/* Builtin functions */
void d_cnjg(doublecomplex *, doublecomplex *);
@@ -195,19 +194,14 @@
if (beta->r != 1. || beta->i != 0.) {
if (*incy == 1) {
if (beta->r == 0. && beta->i == 0.) {
- i__1 = leny;
for (i = 1; i <= leny; ++i) {
- i__2 = i;
Y(i).r = 0., Y(i).i = 0.;
/* L10: */
}
} else {
- i__1 = leny;
for (i = 1; i <= leny; ++i) {
- i__2 = i;
- i__3 = i;
z__1.r = beta->r * Y(i).r - beta->i * Y(i).i,
z__1.i = beta->r * Y(i).i + beta->i * Y(i)
.r;
Y(i).r = z__1.r, Y(i).i = z__1.i;
@@ -216,20 +210,15 @@
}
} else {
iy = ky;
if (beta->r == 0. && beta->i == 0.) {
- i__1 = leny;
for (i = 1; i <= leny; ++i) {
- i__2 = iy;
Y(iy).r = 0., Y(iy).i = 0.;
iy += *incy;
/* L30: */
}
} else {
- i__1 = leny;
for (i = 1; i <= leny; ++i) {
- i__2 = iy;
- i__3 = iy;
z__1.r = beta->r * Y(iy).r - beta->i * Y(iy).i,
z__1.i = beta->r * Y(iy).i + beta->i * Y(iy)
.r;
Y(iy).r = z__1.r, Y(iy).i = z__1.i;
@@ -247,22 +236,15 @@
/* Form y := alpha*A*x + y. */
jx = kx;
if (*incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
if (X(jx).r != 0. || X(jx).i != 0.) {
- i__2 = jx;
z__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i,
z__1.i = alpha->r * X(jx).i + alpha->i * X(jx)
.r;
temp.r = z__1.r, temp.i = z__1.i;
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
z__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
z__2.i = temp.r * A(i,j).i + temp.i * A(i,j)
.r;
z__1.r = Y(i).r + z__2.r, z__1.i = Y(i).i +
@@ -274,23 +256,16 @@
jx += *incx;
/* L60: */
}
} else {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
if (X(jx).r != 0. || X(jx).i != 0.) {
- i__2 = jx;
z__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i,
z__1.i = alpha->r * X(jx).i + alpha->i * X(jx)
.r;
temp.r = z__1.r, temp.i = z__1.i;
iy = ky;
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = iy;
- i__4 = iy;
- i__5 = i + j * a_dim1;
z__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
z__2.i = temp.r * A(i,j).i + temp.i * A(i,j)
.r;
z__1.r = Y(iy).r + z__2.r, z__1.i = Y(iy).i +
@@ -310,55 +285,43 @@
*/
jy = ky;
if (*incx == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp.r = 0., temp.i = 0.;
if (noconj) {
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i;
z__2.r = A(i,j).r * X(i).r - A(i,j).i * X(i)
.i, z__2.i = A(i,j).r * X(i).i + A(i,j)
.i * X(i).r;
z__1.r = temp.r + z__2.r, z__1.i = temp.i + z__2.i;
temp.r = z__1.r, temp.i = z__1.i;
/* L90: */
}
} else {
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
d_cnjg(&z__3, &A(i,j));
- i__3 = i;
z__2.r = z__3.r * X(i).r - z__3.i * X(i).i,
z__2.i = z__3.r * X(i).i + z__3.i * X(i)
.r;
z__1.r = temp.r + z__2.r, z__1.i = temp.i + z__2.i;
temp.r = z__1.r, temp.i = z__1.i;
/* L100: */
}
}
- i__2 = jy;
- i__3 = jy;
z__2.r = alpha->r * temp.r - alpha->i * temp.i, z__2.i =
alpha->r * temp.i + alpha->i * temp.r;
z__1.r = Y(jy).r + z__2.r, z__1.i = Y(jy).i + z__2.i;
Y(jy).r = z__1.r, Y(jy).i = z__1.i;
jy += *incy;
/* L110: */
}
} else {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
temp.r = 0., temp.i = 0.;
ix = kx;
if (noconj) {
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = ix;
z__2.r = A(i,j).r * X(ix).r - A(i,j).i * X(ix)
.i, z__2.i = A(i,j).r * X(ix).i + A(i,j)
.i * X(ix).r;
z__1.r = temp.r + z__2.r, z__1.i = temp.i + z__2.i;
@@ -366,12 +329,10 @@
ix += *incx;
/* L120: */
}
} else {
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
d_cnjg(&z__3, &A(i,j));
- i__3 = ix;
z__2.r = z__3.r * X(ix).r - z__3.i * X(ix).i,
z__2.i = z__3.r * X(ix).i + z__3.i * X(ix)
.r;
z__1.r = temp.r + z__2.r, z__1.i = temp.i + z__2.i;
@@ -379,10 +340,8 @@
ix += *incx;
/* L130: */
}
}
- i__2 = jy;
- i__3 = jy;
z__2.r = alpha->r * temp.r - alpha->i * temp.i, z__2.i =
alpha->r * temp.i + alpha->i * temp.r;
z__1.r = Y(jy).r + z__2.r, z__1.i = Y(jy).i + z__2.i;
Y(jy).r = z__1.r, Y(jy).i = z__1.i;
Index: SuperLU_3.0/CBLAS/zgerc.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zgerc.c
+++ SuperLU_3.0/CBLAS/zgerc.c
@@ -12,9 +12,8 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5;
doublecomplex z__1, z__2;
/* Builtin functions */
void d_cnjg(doublecomplex *, doublecomplex *);
@@ -141,21 +140,15 @@
} else {
jy = 1 - (*n - 1) * *incy;
}
if (*incx == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jy;
if (Y(jy).r != 0. || Y(jy).i != 0.) {
d_cnjg(&z__2, &Y(jy));
z__1.r = alpha->r * z__2.r - alpha->i * z__2.i, z__1.i =
alpha->r * z__2.i + alpha->i * z__2.r;
temp.r = z__1.r, temp.i = z__1.i;
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = i;
z__2.r = X(i).r * temp.r - X(i).i * temp.i, z__2.i =
X(i).r * temp.i + X(i).i * temp.r;
z__1.r = A(i,j).r + z__2.r, z__1.i = A(i,j).i + z__2.i;
A(i,j).r = z__1.r, A(i,j).i = z__1.i;
@@ -170,22 +163,16 @@
kx = 1;
} else {
kx = 1 - (*m - 1) * *incx;
}
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jy;
if (Y(jy).r != 0. || Y(jy).i != 0.) {
d_cnjg(&z__2, &Y(jy));
z__1.r = alpha->r * z__2.r - alpha->i * z__2.i, z__1.i =
alpha->r * z__2.i + alpha->i * z__2.r;
temp.r = z__1.r, temp.i = z__1.i;
ix = kx;
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = ix;
z__2.r = X(ix).r * temp.r - X(ix).i * temp.i, z__2.i =
X(ix).r * temp.i + X(ix).i * temp.r;
z__1.r = A(i,j).r + z__2.r, z__1.i = A(i,j).i + z__2.i;
A(i,j).r = z__1.r, A(i,j).i = z__1.i;
Index: SuperLU_3.0/CBLAS/zhemv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zhemv.c
+++ SuperLU_3.0/CBLAS/zhemv.c
@@ -12,9 +12,8 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5;
doublereal d__1;
doublecomplex z__1, z__2, z__3, z__4;
/* Builtin functions */
@@ -180,19 +179,14 @@
if (beta->r != 1. || beta->i != 0.) {
if (*incy == 1) {
if (beta->r == 0. && beta->i == 0.) {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
Y(i).r = 0., Y(i).i = 0.;
/* L10: */
}
} else {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
- i__3 = i;
z__1.r = beta->r * Y(i).r - beta->i * Y(i).i,
z__1.i = beta->r * Y(i).i + beta->i * Y(i)
.r;
Y(i).r = z__1.r, Y(i).i = z__1.i;
@@ -201,20 +195,15 @@
}
} else {
iy = ky;
if (beta->r == 0. && beta->i == 0.) {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
Y(iy).r = 0., Y(iy).i = 0.;
iy += *incy;
/* L30: */
}
} else {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
- i__3 = iy;
z__1.r = beta->r * Y(iy).r - beta->i * Y(iy).i,
z__1.i = beta->r * Y(iy).i + beta->i * Y(iy)
.r;
Y(iy).r = z__1.r, Y(iy).i = z__1.i;
@@ -231,36 +220,26 @@
/* Form y when A is stored in upper triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
z__1.r = alpha->r * X(j).r - alpha->i * X(j).i, z__1.i =
alpha->r * X(j).i + alpha->i * X(j).r;
temp1.r = z__1.r, temp1.i = z__1.i;
temp2.r = 0., temp2.i = 0.;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
z__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
z__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
z__1.r = Y(i).r + z__2.r, z__1.i = Y(i).i + z__2.i;
Y(i).r = z__1.r, Y(i).i = z__1.i;
d_cnjg(&z__3, &A(i,j));
- i__3 = i;
z__2.r = z__3.r * X(i).r - z__3.i * X(i).i, z__2.i =
z__3.r * X(i).i + z__3.i * X(i).r;
z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i;
temp2.r = z__1.r, temp2.i = z__1.i;
/* L50: */
}
- i__2 = j;
- i__3 = j;
- i__4 = j + j * a_dim1;
d__1 = A(j,j).r;
z__3.r = d__1 * temp1.r, z__3.i = d__1 * temp1.i;
z__2.r = Y(j).r + z__3.r, z__2.i = Y(j).i + z__3.i;
z__4.r = alpha->r * temp2.r - alpha->i * temp2.i, z__4.i =
@@ -271,40 +250,30 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
z__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i, z__1.i =
alpha->r * X(jx).i + alpha->i * X(jx).r;
temp1.r = z__1.r, temp1.i = z__1.i;
temp2.r = 0., temp2.i = 0.;
ix = kx;
iy = ky;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = iy;
- i__4 = iy;
- i__5 = i + j * a_dim1;
z__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
z__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
z__1.r = Y(iy).r + z__2.r, z__1.i = Y(iy).i + z__2.i;
Y(iy).r = z__1.r, Y(iy).i = z__1.i;
d_cnjg(&z__3, &A(i,j));
- i__3 = ix;
z__2.r = z__3.r * X(ix).r - z__3.i * X(ix).i, z__2.i =
z__3.r * X(ix).i + z__3.i * X(ix).r;
z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i;
temp2.r = z__1.r, temp2.i = z__1.i;
ix += *incx;
iy += *incy;
/* L70: */
}
- i__2 = jy;
- i__3 = jy;
- i__4 = j + j * a_dim1;
d__1 = A(j,j).r;
z__3.r = d__1 * temp1.r, z__3.i = d__1 * temp1.i;
z__2.r = Y(jy).r + z__3.r, z__2.i = Y(jy).i + z__3.i;
z__4.r = alpha->r * temp2.r - alpha->i * temp2.i, z__4.i =
@@ -320,42 +289,30 @@
/* Form y when A is stored in lower triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
z__1.r = alpha->r * X(j).r - alpha->i * X(j).i, z__1.i =
alpha->r * X(j).i + alpha->i * X(j).r;
temp1.r = z__1.r, temp1.i = z__1.i;
temp2.r = 0., temp2.i = 0.;
- i__2 = j;
- i__3 = j;
- i__4 = j + j * a_dim1;
d__1 = A(j,j).r;
z__2.r = d__1 * temp1.r, z__2.i = d__1 * temp1.i;
z__1.r = Y(j).r + z__2.r, z__1.i = Y(j).i + z__2.i;
Y(j).r = z__1.r, Y(j).i = z__1.i;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
z__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
z__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
z__1.r = Y(i).r + z__2.r, z__1.i = Y(i).i + z__2.i;
Y(i).r = z__1.r, Y(i).i = z__1.i;
d_cnjg(&z__3, &A(i,j));
- i__3 = i;
z__2.r = z__3.r * X(i).r - z__3.i * X(i).i, z__2.i =
z__3.r * X(i).i + z__3.i * X(i).r;
z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i;
temp2.r = z__1.r, temp2.i = z__1.i;
/* L90: */
}
- i__2 = j;
- i__3 = j;
z__2.r = alpha->r * temp2.r - alpha->i * temp2.i, z__2.i =
alpha->r * temp2.i + alpha->i * temp2.r;
z__1.r = Y(j).r + z__2.r, z__1.i = Y(j).i + z__2.i;
Y(j).r = z__1.r, Y(j).i = z__1.i;
@@ -363,46 +320,34 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
z__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i, z__1.i =
alpha->r * X(jx).i + alpha->i * X(jx).r;
temp1.r = z__1.r, temp1.i = z__1.i;
temp2.r = 0., temp2.i = 0.;
- i__2 = jy;
- i__3 = jy;
- i__4 = j + j * a_dim1;
d__1 = A(j,j).r;
z__2.r = d__1 * temp1.r, z__2.i = d__1 * temp1.i;
z__1.r = Y(jy).r + z__2.r, z__1.i = Y(jy).i + z__2.i;
Y(jy).r = z__1.r, Y(jy).i = z__1.i;
ix = jx;
iy = jy;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
iy += *incy;
- i__3 = iy;
- i__4 = iy;
- i__5 = i + j * a_dim1;
z__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
z__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
z__1.r = Y(iy).r + z__2.r, z__1.i = Y(iy).i + z__2.i;
Y(iy).r = z__1.r, Y(iy).i = z__1.i;
d_cnjg(&z__3, &A(i,j));
- i__3 = ix;
z__2.r = z__3.r * X(ix).r - z__3.i * X(ix).i, z__2.i =
z__3.r * X(ix).i + z__3.i * X(ix).r;
z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i;
temp2.r = z__1.r, temp2.i = z__1.i;
/* L110: */
}
- i__2 = jy;
- i__3 = jy;
z__2.r = alpha->r * temp2.r - alpha->i * temp2.i, z__2.i =
alpha->r * temp2.i + alpha->i * temp2.r;
z__1.r = Y(jy).r + z__2.r, z__1.i = Y(jy).i + z__2.i;
Y(jy).r = z__1.r, Y(jy).i = z__1.i;
Index: SuperLU_3.0/CBLAS/zher2.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/zher2.c
+++ SuperLU_3.0/CBLAS/zher2.c
@@ -12,9 +12,8 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5, i__6;
doublereal d__1;
doublecomplex z__1, z__2, z__3, z__4;
/* Builtin functions */
@@ -185,93 +184,69 @@
/* Form A when A is stored in the upper triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
- i__3 = j;
if (X(j).r != 0. || X(j).i != 0. || (Y(j).r != 0. ||
Y(j).i != 0.)) {
d_cnjg(&z__2, &Y(j));
z__1.r = alpha->r * z__2.r - alpha->i * z__2.i, z__1.i =
alpha->r * z__2.i + alpha->i * z__2.r;
temp1.r = z__1.r, temp1.i = z__1.i;
- i__2 = j;
z__2.r = alpha->r * X(j).r - alpha->i * X(j).i,
z__2.i = alpha->r * X(j).i + alpha->i * X(j)
.r;
d_cnjg(&z__1, &z__2);
temp2.r = z__1.r, temp2.i = z__1.i;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = i;
z__3.r = X(i).r * temp1.r - X(i).i * temp1.i,
z__3.i = X(i).r * temp1.i + X(i).i *
temp1.r;
z__2.r = A(i,j).r + z__3.r, z__2.i = A(i,j).i +
z__3.i;
- i__6 = i;
z__4.r = Y(i).r * temp2.r - Y(i).i * temp2.i,
z__4.i = Y(i).r * temp2.i + Y(i).i *
temp2.r;
z__1.r = z__2.r + z__4.r, z__1.i = z__2.i + z__4.i;
A(i,j).r = z__1.r, A(i,j).i = z__1.i;
/* L10: */
}
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
- i__4 = j;
z__2.r = X(j).r * temp1.r - X(j).i * temp1.i,
z__2.i = X(j).r * temp1.i + X(j).i *
temp1.r;
- i__5 = j;
z__3.r = Y(j).r * temp2.r - Y(j).i * temp2.i,
z__3.i = Y(j).r * temp2.i + Y(j).i *
temp2.r;
z__1.r = z__2.r + z__3.r, z__1.i = z__2.i + z__3.i;
d__1 = A(j,j).r + z__1.r;
A(j,j).r = d__1, A(j,j).i = 0.;
} else {
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
d__1 = A(j,j).r;
A(j,j).r = d__1, A(j,j).i = 0.;
}
/* L20: */
}
} else {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
- i__3 = jy;
if (X(jx).r != 0. || X(jx).i != 0. || (Y(jy).r != 0. ||
Y(jy).i != 0.)) {
d_cnjg(&z__2, &Y(jy));
z__1.r = alpha->r * z__2.r - alpha->i * z__2.i, z__1.i =
alpha->r * z__2.i + alpha->i * z__2.r;
temp1.r = z__1.r, temp1.i = z__1.i;
- i__2 = jx;
z__2.r = alpha->r * X(jx).r - alpha->i * X(jx).i,
z__2.i = alpha->r * X(jx).i + alpha->i * X(jx)
.r;
d_cnjg(&z__1, &z__2);
temp2.r = z__1.r, temp2.i = z__1.i;
ix = kx;
iy = ky;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = ix;
z__3.r = X(ix).r * temp1.r - X(ix).i * temp1.i,
z__3.i = X(ix).r * temp1.i + X(ix).i *
temp1.r;
z__2.r = A(i,j).r + z__3.r, z__2.i = A(i,j).i +
z__3.i;
- i__6 = iy;
z__4.r = Y(iy).r * temp2.r - Y(iy).i * temp2.i,
z__4.i = Y(iy).r * temp2.i + Y(iy).i *
temp2.r;
z__1.r = z__2.r + z__4.r, z__1.i = z__2.i + z__4.i;
@@ -279,24 +254,18 @@
ix += *incx;
iy += *incy;
/* L30: */
}
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
- i__4 = jx;
z__2.r = X(jx).r * temp1.r - X(jx).i * temp1.i,
z__2.i = X(jx).r * temp1.i + X(jx).i *
temp1.r;
- i__5 = jy;
z__3.r = Y(jy).r * temp2.r - Y(jy).i * temp2.i,
z__3.i = Y(jy).r * temp2.i + Y(jy).i *
temp2.r;
z__1.r = z__2.r + z__3.r, z__1.i = z__2.i + z__3.i;
d__1 = A(j,j).r + z__1.r;
A(j,j).r = d__1, A(j,j).i = 0.;
} else {
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
d__1 = A(j,j).r;
A(j,j).r = d__1, A(j,j).i = 0.;
}
jx += *incx;
@@ -308,118 +277,88 @@
/* Form A when A is stored in the lower triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
- i__3 = j;
if (X(j).r != 0. || X(j).i != 0. || (Y(j).r != 0. ||
Y(j).i != 0.)) {
d_cnjg(&z__2, &Y(j));
z__1.r = alpha->r * z__2.r - alpha->i * z__2.i, z__1.i =
alpha->r * z__2.i + alpha->i * z__2.r;
temp1.r = z__1.r, temp1.i = z__1.i;
- i__2 = j;
z__2.r = alpha->r * X(j).r - alpha->i * X(j).i,
z__2.i = alpha->r * X(j).i + alpha->i * X(j)
.r;
d_cnjg(&z__1, &z__2);
temp2.r = z__1.r, temp2.i = z__1.i;
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
- i__4 = j;
z__2.r = X(j).r * temp1.r - X(j).i * temp1.i,
z__2.i = X(j).r * temp1.i + X(j).i *
temp1.r;
- i__5 = j;
z__3.r = Y(j).r * temp2.r - Y(j).i * temp2.i,
z__3.i = Y(j).r * temp2.i + Y(j).i *
temp2.r;
z__1.r = z__2.r + z__3.r, z__1.i = z__2.i + z__3.i;
d__1 = A(j,j).r + z__1.r;
A(j,j).r = d__1, A(j,j).i = 0.;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = i;
z__3.r = X(i).r * temp1.r - X(i).i * temp1.i,
z__3.i = X(i).r * temp1.i + X(i).i *
temp1.r;
z__2.r = A(i,j).r + z__3.r, z__2.i = A(i,j).i +
z__3.i;
- i__6 = i;
z__4.r = Y(i).r * temp2.r - Y(i).i * temp2.i,
z__4.i = Y(i).r * temp2.i + Y(i).i *
temp2.r;
z__1.r = z__2.r + z__4.r, z__1.i = z__2.i + z__4.i;
A(i,j).r = z__1.r, A(i,j).i = z__1.i;
/* L50: */
}
} else {
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
d__1 = A(j,j).r;
A(j,j).r = d__1, A(j,j).i = 0.;
}
/* L60: */
}
} else {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
- i__3 = jy;
if (X(jx).r != 0. || X(jx).i != 0. || (Y(jy).r != 0. ||
Y(jy).i != 0.)) {
d_cnjg(&z__2, &Y(jy));
z__1.r = alpha->r * z__2.r - alpha->i * z__2.i, z__1.i =
alpha->r * z__2.i + alpha->i * z__2.r;
temp1.r = z__1.r, temp1.i = z__1.i;
- i__2 = jx;
z__2.r = alpha->r * X(jx).r - alpha->i * X(jx).i,
z__2.i = alpha->r * X(jx).i + alpha->i * X(jx)
.r;
d_cnjg(&z__1, &z__2);
temp2.r = z__1.r, temp2.i = z__1.i;
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
- i__4 = jx;
z__2.r = X(jx).r * temp1.r - X(jx).i * temp1.i,
z__2.i = X(jx).r * temp1.i + X(jx).i *
temp1.r;
- i__5 = jy;
z__3.r = Y(jy).r * temp2.r - Y(jy).i * temp2.i,
z__3.i = Y(jy).r * temp2.i + Y(jy).i *
temp2.r;
z__1.r = z__2.r + z__3.r, z__1.i = z__2.i + z__3.i;
d__1 = A(j,j).r + z__1.r;
A(j,j).r = d__1, A(j,j).i = 0.;
ix = jx;
iy = jy;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
iy += *incy;
- i__3 = i + j * a_dim1;
- i__4 = i + j * a_dim1;
- i__5 = ix;
z__3.r = X(ix).r * temp1.r - X(ix).i * temp1.i,
z__3.i = X(ix).r * temp1.i + X(ix).i *
temp1.r;
z__2.r = A(i,j).r + z__3.r, z__2.i = A(i,j).i +
z__3.i;
- i__6 = iy;
z__4.r = Y(iy).r * temp2.r - Y(iy).i * temp2.i,
z__4.i = Y(iy).r * temp2.i + Y(iy).i *
temp2.r;
z__1.r = z__2.r + z__4.r, z__1.i = z__2.i + z__4.i;
A(i,j).r = z__1.r, A(i,j).i = z__1.i;
/* L70: */
}
} else {
- i__2 = j + j * a_dim1;
- i__3 = j + j * a_dim1;
d__1 = A(j,j).r;
A(j,j).r = d__1, A(j,j).i = 0.;
}
jx += *incx;
Index: SuperLU_3.0/CBLAS/ztrsv.c
===================================================================
--- SuperLU_3.0.orig/CBLAS/ztrsv.c
+++ SuperLU_3.0/CBLAS/ztrsv.c
@@ -11,9 +11,8 @@
{
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5;
doublecomplex z__1, z__2, z__3;
/* Builtin functions */
void z_div(doublecomplex *, doublecomplex *, doublecomplex *), d_cnjg(
@@ -187,21 +186,15 @@
if (lsame_(uplo, "U")) {
if (*incx == 1) {
for (j = *n; j >= 1; --j) {
- i__1 = j;
if (X(j).r != 0. || X(j).i != 0.) {
if (nounit) {
- i__1 = j;
z_div(&z__1, &X(j), &A(j,j));
X(j).r = z__1.r, X(j).i = z__1.i;
}
- i__1 = j;
temp.r = X(j).r, temp.i = X(j).i;
for (i = j - 1; i >= 1; --i) {
- i__1 = i;
- i__2 = i;
- i__3 = i + j * a_dim1;
z__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
z__2.i = temp.r * A(i,j).i + temp.i * A(i,j).r;
z__1.r = X(i).r - z__2.r, z__1.i = X(i).i -
z__2.i;
@@ -213,23 +206,17 @@
}
} else {
jx = kx + (*n - 1) * *incx;
for (j = *n; j >= 1; --j) {
- i__1 = jx;
if (X(jx).r != 0. || X(jx).i != 0.) {
if (nounit) {
- i__1 = jx;
z_div(&z__1, &X(jx), &A(j,j));
X(jx).r = z__1.r, X(jx).i = z__1.i;
}
- i__1 = jx;
temp.r = X(jx).r, temp.i = X(jx).i;
ix = jx;
for (i = j - 1; i >= 1; --i) {
ix -= *incx;
- i__1 = ix;
- i__2 = ix;
- i__3 = i + j * a_dim1;
z__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
z__2.i = temp.r * A(i,j).i + temp.i * A(i,j).r;
z__1.r = X(ix).r - z__2.r, z__1.i = X(ix).i -
z__2.i;
@@ -242,24 +229,16 @@
}
}
} else {
if (*incx == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
if (X(j).r != 0. || X(j).i != 0.) {
if (nounit) {
- i__2 = j;
z_div(&z__1, &X(j), &A(j,j));
X(j).r = z__1.r, X(j).i = z__1.i;
}
- i__2 = j;
temp.r = X(j).r, temp.i = X(j).i;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
z__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
z__2.i = temp.r * A(i,j).i + temp.i * A(i,j).r;
z__1.r = X(i).r - z__2.r, z__1.i = X(i).i -
z__2.i;
@@ -270,26 +249,18 @@
/* L60: */
}
} else {
jx = kx;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
if (X(jx).r != 0. || X(jx).i != 0.) {
if (nounit) {
- i__2 = jx;
z_div(&z__1, &X(jx), &A(j,j));
X(jx).r = z__1.r, X(jx).i = z__1.i;
}
- i__2 = jx;
temp.r = X(jx).r, temp.i = X(jx).i;
ix = jx;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
- i__3 = ix;
- i__4 = ix;
- i__5 = i + j * a_dim1;
z__2.r = temp.r * A(i,j).r - temp.i * A(i,j).i,
z__2.i = temp.r * A(i,j).i + temp.i * A(i,j).r;
z__1.r = X(ix).r - z__2.r, z__1.i = X(ix).i -
z__2.i;
@@ -307,17 +278,12 @@
/* Form x := inv( A' )*x or x := inv( conjg( A' ) )*x. */
if (lsame_(uplo, "U")) {
if (*incx == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
temp.r = X(j).r, temp.i = X(j).i;
if (noconj) {
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = i;
z__2.r = A(i,j).r * X(i).r - A(i,j).i * X(
i).i, z__2.i = A(i,j).r * X(i).i +
A(i,j).i * X(i).r;
z__1.r = temp.r - z__2.r, z__1.i = temp.i -
@@ -329,12 +295,10 @@
z_div(&z__1, &temp, &A(j,j));
temp.r = z__1.r, temp.i = z__1.i;
}
} else {
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
d_cnjg(&z__3, &A(i,j));
- i__3 = i;
z__2.r = z__3.r * X(i).r - z__3.i * X(i).i,
z__2.i = z__3.r * X(i).i + z__3.i * X(
i).r;
z__1.r = temp.r - z__2.r, z__1.i = temp.i -
@@ -347,24 +311,18 @@
z_div(&z__1, &temp, &z__2);
temp.r = z__1.r, temp.i = z__1.i;
}
}
- i__2 = j;
X(j).r = temp.r, X(j).i = temp.i;
/* L110: */
}
} else {
jx = kx;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
ix = kx;
- i__2 = jx;
temp.r = X(jx).r, temp.i = X(jx).i;
if (noconj) {
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i + j * a_dim1;
- i__4 = ix;
z__2.r = A(i,j).r * X(ix).r - A(i,j).i * X(
ix).i, z__2.i = A(i,j).r * X(ix).i +
A(i,j).i * X(ix).r;
z__1.r = temp.r - z__2.r, z__1.i = temp.i -
@@ -377,12 +335,10 @@
z_div(&z__1, &temp, &A(j,j));
temp.r = z__1.r, temp.i = z__1.i;
}
} else {
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
d_cnjg(&z__3, &A(i,j));
- i__3 = ix;
z__2.r = z__3.r * X(ix).r - z__3.i * X(ix).i,
z__2.i = z__3.r * X(ix).i + z__3.i * X(
ix).r;
z__1.r = temp.r - z__2.r, z__1.i = temp.i -
@@ -396,24 +352,19 @@
z_div(&z__1, &temp, &z__2);
temp.r = z__1.r, temp.i = z__1.i;
}
}
- i__2 = jx;
X(jx).r = temp.r, X(jx).i = temp.i;
jx += *incx;
/* L140: */
}
}
} else {
if (*incx == 1) {
for (j = *n; j >= 1; --j) {
- i__1 = j;
temp.r = X(j).r, temp.i = X(j).i;
if (noconj) {
- i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
- i__2 = i + j * a_dim1;
- i__3 = i;
z__2.r = A(i,j).r * X(i).r - A(i,j).i * X(
i).i, z__2.i = A(i,j).r * X(i).i +
A(i,j).i * X(i).r;
z__1.r = temp.r - z__2.r, z__1.i = temp.i -
@@ -425,12 +376,10 @@
z_div(&z__1, &temp, &A(j,j));
temp.r = z__1.r, temp.i = z__1.i;
}
} else {
- i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
d_cnjg(&z__3, &A(i,j));
- i__2 = i;
z__2.r = z__3.r * X(i).r - z__3.i * X(i).i,
z__2.i = z__3.r * X(i).i + z__3.i * X(
i).r;
z__1.r = temp.r - z__2.r, z__1.i = temp.i -
@@ -443,24 +392,19 @@
z_div(&z__1, &temp, &z__2);
temp.r = z__1.r, temp.i = z__1.i;
}
}
- i__1 = j;
X(j).r = temp.r, X(j).i = temp.i;
/* L170: */
}
} else {
kx += (*n - 1) * *incx;
jx = kx;
for (j = *n; j >= 1; --j) {
ix = kx;
- i__1 = jx;
temp.r = X(jx).r, temp.i = X(jx).i;
if (noconj) {
- i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
- i__2 = i + j * a_dim1;
- i__3 = ix;
z__2.r = A(i,j).r * X(ix).r - A(i,j).i * X(
ix).i, z__2.i = A(i,j).r * X(ix).i +
A(i,j).i * X(ix).r;
z__1.r = temp.r - z__2.r, z__1.i = temp.i -
@@ -473,12 +417,10 @@
z_div(&z__1, &temp, &A(j,j));
temp.r = z__1.r, temp.i = z__1.i;
}
} else {
- i__1 = j + 1;
for (i = *n; i >= j+1; --i) {
d_cnjg(&z__3, &A(i,j));
- i__2 = ix;
z__2.r = z__3.r * X(ix).r - z__3.i * X(ix).i,
z__2.i = z__3.r * X(ix).i + z__3.i * X(
ix).r;
z__1.r = temp.r - z__2.r, z__1.i = temp.i -
@@ -492,9 +434,8 @@
z_div(&z__1, &temp, &z__2);
temp.r = z__1.r, temp.i = z__1.i;
}
}
- i__1 = jx;
X(jx).r = temp.r, X(jx).i = temp.i;
jx -= *incx;
/* L200: */
}
Index: SuperLU_3.0/SRC/sutil.c
===================================================================
--- SuperLU_3.0.orig/SRC/sutil.c
+++ SuperLU_3.0/SRC/sutil.c
@@ -266,14 +266,15 @@ sPrint_SuperNode_Matrix(char *what, Supe
void
sPrint_Dense_Matrix(char *what, SuperMatrix *A)
{
DNformat *Astore;
- register int i, j, lda = Astore->lda;
+ register int i, j, lda;
float *dp;
printf("\nDense matrix %s:\n", what);
printf("Stype %d, Dtype %d, Mtype %d\n", A->Stype,A->Dtype,A->Mtype);
Astore = (DNformat *) A->Store;
+ lda = Astore->lda;
dp = (float *) Astore->nzval;
printf("nrow %d, ncol %d, lda %d\n", A->nrow,A->ncol,lda);
printf("\nnzval: ");
for (j = 0; j < A->ncol; ++j) {
Index: SuperLU_3.0/SRC/cutil.c
===================================================================
--- SuperLU_3.0.orig/SRC/cutil.c
+++ SuperLU_3.0/SRC/cutil.c
@@ -266,14 +266,15 @@ cPrint_SuperNode_Matrix(char *what, Supe
void
cPrint_Dense_Matrix(char *what, SuperMatrix *A)
{
DNformat *Astore;
- register int i, j, lda = Astore->lda;
+ register int i, j, lda;
float *dp;
printf("\nDense matrix %s:\n", what);
printf("Stype %d, Dtype %d, Mtype %d\n", A->Stype,A->Dtype,A->Mtype);
Astore = (DNformat *) A->Store;
+ lda = Astore->lda;
dp = (float *) Astore->nzval;
printf("nrow %d, ncol %d, lda %d\n", A->nrow,A->ncol,lda);
printf("\nnzval: ");
for (j = 0; j < A->ncol; ++j) {
Index: SuperLU_3.0/SRC/dutil.c
===================================================================
--- SuperLU_3.0.orig/SRC/dutil.c
+++ SuperLU_3.0/SRC/dutil.c
@@ -266,14 +266,15 @@ dPrint_SuperNode_Matrix(char *what, Supe
void
dPrint_Dense_Matrix(char *what, SuperMatrix *A)
{
DNformat *Astore;
- register int i, j, lda = Astore->lda;
+ register int i, j, lda;
double *dp;
printf("\nDense matrix %s:\n", what);
printf("Stype %d, Dtype %d, Mtype %d\n", A->Stype,A->Dtype,A->Mtype);
Astore = (DNformat *) A->Store;
+ lda = Astore->lda;
dp = (double *) Astore->nzval;
printf("nrow %d, ncol %d, lda %d\n", A->nrow,A->ncol,lda);
printf("\nnzval: ");
for (j = 0; j < A->ncol; ++j) {
Index: SuperLU_3.0/SRC/zutil.c
===================================================================
--- SuperLU_3.0.orig/SRC/zutil.c
+++ SuperLU_3.0/SRC/zutil.c
@@ -266,14 +266,15 @@ zPrint_SuperNode_Matrix(char *what, Supe
void
zPrint_Dense_Matrix(char *what, SuperMatrix *A)
{
DNformat *Astore;
- register int i, j, lda = Astore->lda;
+ register int i, j, lda;
double *dp;
printf("\nDense matrix %s:\n", what);
printf("Stype %d, Dtype %d, Mtype %d\n", A->Stype,A->Dtype,A->Mtype);
Astore = (DNformat *) A->Store;
+ lda = Astore->lda;
dp = (double *) Astore->nzval;
printf("nrow %d, ncol %d, lda %d\n", A->nrow,A->ncol,lda);
printf("\nnzval: ");
for (j = 0; j < A->ncol; ++j) {
Index: SuperLU_3.0/TESTING/MATGEN/claset.c
===================================================================
--- SuperLU_3.0.orig/TESTING/MATGEN/claset.c
+++ SuperLU_3.0/TESTING/MATGEN/claset.c
@@ -55,9 +55,8 @@
Parameter adjustments
Function Body */
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3;
/* Local variables */
static integer i, j;
extern logical lsame_(char *, char *);
@@ -70,23 +69,17 @@
/* Set the diagonal to BETA and the strictly upper triangular
part of the array to ALPHA. */
- i__1 = *n;
for (j = 2; j <= *n; ++j) {
/* Computing MIN */
- i__3 = j - 1;
- i__2 = min(i__3,*m);
for (i = 1; i <= min(j-1,*m); ++i) {
- i__3 = i + j * a_dim1;
A(i,j).r = alpha->r, A(i,j).i = alpha->i;
/* L10: */
}
/* L20: */
}
- i__1 = min(*n,*m);
for (i = 1; i <= min(*n,*m); ++i) {
- i__2 = i + i * a_dim1;
A(i,i).r = beta->r, A(i,i).i = beta->i;
/* L30: */
}
@@ -95,21 +88,16 @@
/* Set the diagonal to BETA and the strictly lower triangular
part of the array to ALPHA. */
- i__1 = min(*m,*n);
for (j = 1; j <= min(*m,*n); ++j) {
- i__2 = *m;
for (i = j + 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
A(i,j).r = alpha->r, A(i,j).i = alpha->i;
/* L40: */
}
/* L50: */
}
- i__1 = min(*n,*m);
for (i = 1; i <= min(*n,*m); ++i) {
- i__2 = i + i * a_dim1;
A(i,i).r = beta->r, A(i,i).i = beta->i;
/* L60: */
}
@@ -117,21 +105,16 @@
/* Set the array to BETA on the diagonal and ALPHA on the
offdiagonal. */
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
A(i,j).r = alpha->r, A(i,j).i = alpha->i;
/* L70: */
}
/* L80: */
}
- i__1 = min(*m,*n);
for (i = 1; i <= min(*m,*n); ++i) {
- i__2 = i + i * a_dim1;
A(i,i).r = beta->r, A(i,i).i = beta->i;
/* L90: */
}
}
Index: SuperLU_3.0/TESTING/MATGEN/csymv.c
===================================================================
--- SuperLU_3.0.orig/TESTING/MATGEN/csymv.c
+++ SuperLU_3.0/TESTING/MATGEN/csymv.c
@@ -101,9 +101,8 @@
Parameter adjustments
Function Body */
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5;
complex q__1, q__2, q__3, q__4;
/* Local variables */
static integer info;
static complex temp1, temp2;
@@ -163,19 +162,14 @@
if (beta->r != 1.f || beta->i != 0.f) {
if (*incy == 1) {
if (beta->r == 0.f && beta->i == 0.f) {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
Y(i).r = 0.f, Y(i).i = 0.f;
/* L10: */
}
} else {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
- i__3 = i;
q__1.r = beta->r * Y(i).r - beta->i * Y(i).i,
q__1.i = beta->r * Y(i).i + beta->i * Y(i)
.r;
Y(i).r = q__1.r, Y(i).i = q__1.i;
@@ -184,20 +178,15 @@
}
} else {
iy = ky;
if (beta->r == 0.f && beta->i == 0.f) {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
Y(iy).r = 0.f, Y(iy).i = 0.f;
iy += *incy;
/* L30: */
}
} else {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
- i__3 = iy;
q__1.r = beta->r * Y(iy).r - beta->i * Y(iy).i,
q__1.i = beta->r * Y(iy).i + beta->i * Y(iy)
.r;
Y(iy).r = q__1.r, Y(iy).i = q__1.i;
@@ -214,37 +203,26 @@
/* Form y when A is stored in upper triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
q__1.r = alpha->r * X(j).r - alpha->i * X(j).i, q__1.i =
alpha->r * X(j).i + alpha->i * X(j).r;
temp1.r = q__1.r, temp1.i = q__1.i;
temp2.r = 0.f, temp2.i = 0.f;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
q__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
q__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
q__1.r = Y(i).r + q__2.r, q__1.i = Y(i).i + q__2.i;
Y(i).r = q__1.r, Y(i).i = q__1.i;
- i__3 = i + j * a_dim1;
- i__4 = i;
q__2.r = A(i,j).r * X(i).r - A(i,j).i * X(i).i,
q__2.i = A(i,j).r * X(i).i + A(i,j).i * X(
i).r;
q__1.r = temp2.r + q__2.r, q__1.i = temp2.i + q__2.i;
temp2.r = q__1.r, temp2.i = q__1.i;
/* L50: */
}
- i__2 = j;
- i__3 = j;
- i__4 = j + j * a_dim1;
q__3.r = temp1.r * A(j,j).r - temp1.i * A(j,j).i, q__3.i =
temp1.r * A(j,j).i + temp1.i * A(j,j).r;
q__2.r = Y(j).r + q__3.r, q__2.i = Y(j).i + q__3.i;
q__4.r = alpha->r * temp2.r - alpha->i * temp2.i, q__4.i =
@@ -255,29 +233,21 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
q__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i, q__1.i =
alpha->r * X(jx).i + alpha->i * X(jx).r;
temp1.r = q__1.r, temp1.i = q__1.i;
temp2.r = 0.f, temp2.i = 0.f;
ix = kx;
iy = ky;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = iy;
- i__4 = iy;
- i__5 = i + j * a_dim1;
q__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
q__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
q__1.r = Y(iy).r + q__2.r, q__1.i = Y(iy).i + q__2.i;
Y(iy).r = q__1.r, Y(iy).i = q__1.i;
- i__3 = i + j * a_dim1;
- i__4 = ix;
q__2.r = A(i,j).r * X(ix).r - A(i,j).i * X(ix).i,
q__2.i = A(i,j).r * X(ix).i + A(i,j).i * X(
ix).r;
q__1.r = temp2.r + q__2.r, q__1.i = temp2.i + q__2.i;
@@ -285,11 +255,8 @@
ix += *incx;
iy += *incy;
/* L70: */
}
- i__2 = jy;
- i__3 = jy;
- i__4 = j + j * a_dim1;
q__3.r = temp1.r * A(j,j).r - temp1.i * A(j,j).i, q__3.i =
temp1.r * A(j,j).i + temp1.i * A(j,j).r;
q__2.r = Y(jy).r + q__3.r, q__2.i = Y(jy).i + q__3.i;
q__4.r = alpha->r * temp2.r - alpha->i * temp2.i, q__4.i =
@@ -305,43 +272,30 @@
/* Form y when A is stored in lower triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
q__1.r = alpha->r * X(j).r - alpha->i * X(j).i, q__1.i =
alpha->r * X(j).i + alpha->i * X(j).r;
temp1.r = q__1.r, temp1.i = q__1.i;
temp2.r = 0.f, temp2.i = 0.f;
- i__2 = j;
- i__3 = j;
- i__4 = j + j * a_dim1;
q__2.r = temp1.r * A(j,j).r - temp1.i * A(j,j).i, q__2.i =
temp1.r * A(j,j).i + temp1.i * A(j,j).r;
q__1.r = Y(j).r + q__2.r, q__1.i = Y(j).i + q__2.i;
Y(j).r = q__1.r, Y(j).i = q__1.i;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
q__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
q__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
q__1.r = Y(i).r + q__2.r, q__1.i = Y(i).i + q__2.i;
Y(i).r = q__1.r, Y(i).i = q__1.i;
- i__3 = i + j * a_dim1;
- i__4 = i;
q__2.r = A(i,j).r * X(i).r - A(i,j).i * X(i).i,
q__2.i = A(i,j).r * X(i).i + A(i,j).i * X(
i).r;
q__1.r = temp2.r + q__2.r, q__1.i = temp2.i + q__2.i;
temp2.r = q__1.r, temp2.i = q__1.i;
/* L90: */
}
- i__2 = j;
- i__3 = j;
q__2.r = alpha->r * temp2.r - alpha->i * temp2.i, q__2.i =
alpha->r * temp2.i + alpha->i * temp2.r;
q__1.r = Y(j).r + q__2.r, q__1.i = Y(j).i + q__2.i;
Y(j).r = q__1.r, Y(j).i = q__1.i;
@@ -349,47 +303,34 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
q__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i, q__1.i =
alpha->r * X(jx).i + alpha->i * X(jx).r;
temp1.r = q__1.r, temp1.i = q__1.i;
temp2.r = 0.f, temp2.i = 0.f;
- i__2 = jy;
- i__3 = jy;
- i__4 = j + j * a_dim1;
q__2.r = temp1.r * A(j,j).r - temp1.i * A(j,j).i, q__2.i =
temp1.r * A(j,j).i + temp1.i * A(j,j).r;
q__1.r = Y(jy).r + q__2.r, q__1.i = Y(jy).i + q__2.i;
Y(jy).r = q__1.r, Y(jy).i = q__1.i;
ix = jx;
iy = jy;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
iy += *incy;
- i__3 = iy;
- i__4 = iy;
- i__5 = i + j * a_dim1;
q__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
q__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
q__1.r = Y(iy).r + q__2.r, q__1.i = Y(iy).i + q__2.i;
Y(iy).r = q__1.r, Y(iy).i = q__1.i;
- i__3 = i + j * a_dim1;
- i__4 = ix;
q__2.r = A(i,j).r * X(ix).r - A(i,j).i * X(ix).i,
q__2.i = A(i,j).r * X(ix).i + A(i,j).i * X(
ix).r;
q__1.r = temp2.r + q__2.r, q__1.i = temp2.i + q__2.i;
temp2.r = q__1.r, temp2.i = q__1.i;
/* L110: */
}
- i__2 = jy;
- i__3 = jy;
q__2.r = alpha->r * temp2.r - alpha->i * temp2.i, q__2.i =
alpha->r * temp2.i + alpha->i * temp2.r;
q__1.r = Y(jy).r + q__2.r, q__1.i = Y(jy).i + q__2.i;
Y(jy).r = q__1.r, Y(jy).i = q__1.i;
Index: SuperLU_3.0/TESTING/MATGEN/zlaset.c
===================================================================
--- SuperLU_3.0.orig/TESTING/MATGEN/zlaset.c
+++ SuperLU_3.0/TESTING/MATGEN/zlaset.c
@@ -56,9 +56,8 @@
Parameter adjustments
Function Body */
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3;
/* Local variables */
static integer i, j;
extern logical lsame_(char *, char *);
@@ -71,23 +70,17 @@
/* Set the diagonal to BETA and the strictly upper triangular
part of the array to ALPHA. */
- i__1 = *n;
for (j = 2; j <= *n; ++j) {
/* Computing MIN */
- i__3 = j - 1;
- i__2 = min(i__3,*m);
for (i = 1; i <= min(j-1,*m); ++i) {
- i__3 = i + j * a_dim1;
A(i,j).r = alpha->r, A(i,j).i = alpha->i;
/* L10: */
}
/* L20: */
}
- i__1 = min(*n,*m);
for (i = 1; i <= min(*n,*m); ++i) {
- i__2 = i + i * a_dim1;
A(i,i).r = beta->r, A(i,i).i = beta->i;
/* L30: */
}
@@ -96,21 +89,16 @@
/* Set the diagonal to BETA and the strictly lower triangular
part of the array to ALPHA. */
- i__1 = min(*m,*n);
for (j = 1; j <= min(*m,*n); ++j) {
- i__2 = *m;
for (i = j + 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
A(i,j).r = alpha->r, A(i,j).i = alpha->i;
/* L40: */
}
/* L50: */
}
- i__1 = min(*n,*m);
for (i = 1; i <= min(*n,*m); ++i) {
- i__2 = i + i * a_dim1;
A(i,i).r = beta->r, A(i,i).i = beta->i;
/* L60: */
}
@@ -118,21 +106,16 @@
/* Set the array to BETA on the diagonal and ALPHA on the
offdiagonal. */
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = *m;
for (i = 1; i <= *m; ++i) {
- i__3 = i + j * a_dim1;
A(i,j).r = alpha->r, A(i,j).i = alpha->i;
/* L70: */
}
/* L80: */
}
- i__1 = min(*m,*n);
for (i = 1; i <= min(*m,*n); ++i) {
- i__2 = i + i * a_dim1;
A(i,i).r = beta->r, A(i,i).i = beta->i;
/* L90: */
}
}
Index: SuperLU_3.0/TESTING/MATGEN/zsymv.c
===================================================================
--- SuperLU_3.0.orig/TESTING/MATGEN/zsymv.c
+++ SuperLU_3.0/TESTING/MATGEN/zsymv.c
@@ -101,9 +101,8 @@
Parameter adjustments
Function Body */
/* System generated locals */
- integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5;
doublecomplex z__1, z__2, z__3, z__4;
/* Local variables */
static integer info;
static doublecomplex temp1, temp2;
@@ -163,19 +162,14 @@
if (beta->r != 1. || beta->i != 0.) {
if (*incy == 1) {
if (beta->r == 0. && beta->i == 0.) {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
Y(i).r = 0., Y(i).i = 0.;
/* L10: */
}
} else {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = i;
- i__3 = i;
z__1.r = beta->r * Y(i).r - beta->i * Y(i).i,
z__1.i = beta->r * Y(i).i + beta->i * Y(i)
.r;
Y(i).r = z__1.r, Y(i).i = z__1.i;
@@ -184,20 +178,15 @@
}
} else {
iy = ky;
if (beta->r == 0. && beta->i == 0.) {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
Y(iy).r = 0., Y(iy).i = 0.;
iy += *incy;
/* L30: */
}
} else {
- i__1 = *n;
for (i = 1; i <= *n; ++i) {
- i__2 = iy;
- i__3 = iy;
z__1.r = beta->r * Y(iy).r - beta->i * Y(iy).i,
z__1.i = beta->r * Y(iy).i + beta->i * Y(iy)
.r;
Y(iy).r = z__1.r, Y(iy).i = z__1.i;
@@ -214,37 +203,26 @@
/* Form y when A is stored in upper triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
z__1.r = alpha->r * X(j).r - alpha->i * X(j).i, z__1.i =
alpha->r * X(j).i + alpha->i * X(j).r;
temp1.r = z__1.r, temp1.i = z__1.i;
temp2.r = 0., temp2.i = 0.;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
z__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
z__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
z__1.r = Y(i).r + z__2.r, z__1.i = Y(i).i + z__2.i;
Y(i).r = z__1.r, Y(i).i = z__1.i;
- i__3 = i + j * a_dim1;
- i__4 = i;
z__2.r = A(i,j).r * X(i).r - A(i,j).i * X(i).i,
z__2.i = A(i,j).r * X(i).i + A(i,j).i * X(
i).r;
z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i;
temp2.r = z__1.r, temp2.i = z__1.i;
/* L50: */
}
- i__2 = j;
- i__3 = j;
- i__4 = j + j * a_dim1;
z__3.r = temp1.r * A(j,j).r - temp1.i * A(j,j).i, z__3.i =
temp1.r * A(j,j).i + temp1.i * A(j,j).r;
z__2.r = Y(j).r + z__3.r, z__2.i = Y(j).i + z__3.i;
z__4.r = alpha->r * temp2.r - alpha->i * temp2.i, z__4.i =
@@ -255,29 +233,21 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
z__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i, z__1.i =
alpha->r * X(jx).i + alpha->i * X(jx).r;
temp1.r = z__1.r, temp1.i = z__1.i;
temp2.r = 0., temp2.i = 0.;
ix = kx;
iy = ky;
- i__2 = j - 1;
for (i = 1; i <= j-1; ++i) {
- i__3 = iy;
- i__4 = iy;
- i__5 = i + j * a_dim1;
z__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
z__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
z__1.r = Y(iy).r + z__2.r, z__1.i = Y(iy).i + z__2.i;
Y(iy).r = z__1.r, Y(iy).i = z__1.i;
- i__3 = i + j * a_dim1;
- i__4 = ix;
z__2.r = A(i,j).r * X(ix).r - A(i,j).i * X(ix).i,
z__2.i = A(i,j).r * X(ix).i + A(i,j).i * X(
ix).r;
z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i;
@@ -285,11 +255,8 @@
ix += *incx;
iy += *incy;
/* L70: */
}
- i__2 = jy;
- i__3 = jy;
- i__4 = j + j * a_dim1;
z__3.r = temp1.r * A(j,j).r - temp1.i * A(j,j).i, z__3.i =
temp1.r * A(j,j).i + temp1.i * A(j,j).r;
z__2.r = Y(jy).r + z__3.r, z__2.i = Y(jy).i + z__3.i;
z__4.r = alpha->r * temp2.r - alpha->i * temp2.i, z__4.i =
@@ -305,43 +272,30 @@
/* Form y when A is stored in lower triangle. */
if (*incx == 1 && *incy == 1) {
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = j;
z__1.r = alpha->r * X(j).r - alpha->i * X(j).i, z__1.i =
alpha->r * X(j).i + alpha->i * X(j).r;
temp1.r = z__1.r, temp1.i = z__1.i;
temp2.r = 0., temp2.i = 0.;
- i__2 = j;
- i__3 = j;
- i__4 = j + j * a_dim1;
z__2.r = temp1.r * A(j,j).r - temp1.i * A(j,j).i, z__2.i =
temp1.r * A(j,j).i + temp1.i * A(j,j).r;
z__1.r = Y(j).r + z__2.r, z__1.i = Y(j).i + z__2.i;
Y(j).r = z__1.r, Y(j).i = z__1.i;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
- i__3 = i;
- i__4 = i;
- i__5 = i + j * a_dim1;
z__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
z__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
z__1.r = Y(i).r + z__2.r, z__1.i = Y(i).i + z__2.i;
Y(i).r = z__1.r, Y(i).i = z__1.i;
- i__3 = i + j * a_dim1;
- i__4 = i;
z__2.r = A(i,j).r * X(i).r - A(i,j).i * X(i).i,
z__2.i = A(i,j).r * X(i).i + A(i,j).i * X(
i).r;
z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i;
temp2.r = z__1.r, temp2.i = z__1.i;
/* L90: */
}
- i__2 = j;
- i__3 = j;
z__2.r = alpha->r * temp2.r - alpha->i * temp2.i, z__2.i =
alpha->r * temp2.i + alpha->i * temp2.r;
z__1.r = Y(j).r + z__2.r, z__1.i = Y(j).i + z__2.i;
Y(j).r = z__1.r, Y(j).i = z__1.i;
@@ -349,47 +303,34 @@
}
} else {
jx = kx;
jy = ky;
- i__1 = *n;
for (j = 1; j <= *n; ++j) {
- i__2 = jx;
z__1.r = alpha->r * X(jx).r - alpha->i * X(jx).i, z__1.i =
alpha->r * X(jx).i + alpha->i * X(jx).r;
temp1.r = z__1.r, temp1.i = z__1.i;
temp2.r = 0., temp2.i = 0.;
- i__2 = jy;
- i__3 = jy;
- i__4 = j + j * a_dim1;
z__2.r = temp1.r * A(j,j).r - temp1.i * A(j,j).i, z__2.i =
temp1.r * A(j,j).i + temp1.i * A(j,j).r;
z__1.r = Y(jy).r + z__2.r, z__1.i = Y(jy).i + z__2.i;
Y(jy).r = z__1.r, Y(jy).i = z__1.i;
ix = jx;
iy = jy;
- i__2 = *n;
for (i = j + 1; i <= *n; ++i) {
ix += *incx;
iy += *incy;
- i__3 = iy;
- i__4 = iy;
- i__5 = i + j * a_dim1;
z__2.r = temp1.r * A(i,j).r - temp1.i * A(i,j).i,
z__2.i = temp1.r * A(i,j).i + temp1.i * A(i,j)
.r;
z__1.r = Y(iy).r + z__2.r, z__1.i = Y(iy).i + z__2.i;
Y(iy).r = z__1.r, Y(iy).i = z__1.i;
- i__3 = i + j * a_dim1;
- i__4 = ix;
z__2.r = A(i,j).r * X(ix).r - A(i,j).i * X(ix).i,
z__2.i = A(i,j).r * X(ix).i + A(i,j).i * X(
ix).r;
z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i;
temp2.r = z__1.r, temp2.i = z__1.i;
/* L110: */
}
- i__2 = jy;
- i__3 = jy;
z__2.r = alpha->r * temp2.r - alpha->i * temp2.i, z__2.i =
alpha->r * temp2.i + alpha->i * temp2.r;
z__1.r = Y(jy).r + z__2.r, z__1.i = Y(jy).i + z__2.i;
Y(jy).r = z__1.r, Y(jy).i = z__1.i;
Reference in New Issue View Git Blame Copy Permalink