amd/blis
1
0
Fork 0
mirror of https://github.com/amd/blis.git synced 2026-05-11 01:30:00 +00:00
Code Issues Packages Projects Releases Wiki Activity

Fix remaining f2c conflicts.

Browse Source
This commit is contained in:
Devin Matthews
2016-02-25 13:51:26 -06:00
parent 372eef0b6c
commit 8624a33ccc
12 changed files with 577 additions and 577 deletions
Download Patch File Download Diff File Expand all files Collapse all files

78
frame/compat/f2c/bla_gbmv.c
View File

@@ -1382,7 +1382,7 @@
/* Quick return if possible. */ /* Quick return if possible. */
if (*m == 0 || *n == 0 || (bli_zbla_real(*alpha) == 0. && bli_zimag(*alpha) == 0. && (bli_zbla_real(*beta) == if (*m == 0 || *n == 0 || (bli_zreal(*alpha) == 0. && bli_zimag(*alpha) == 0. && (bli_zreal(*beta) ==
1. && bli_zimag(*beta) == 0.))) { 1. && bli_zimag(*beta) == 0.))) {
return 0; return 0;
} }
@@ -1415,9 +1415,9 @@
/* First form y := beta*y. */ /* First form y := beta*y. */
if (bli_zbla_real(*beta) != 1. || bli_zimag(*beta) != 0.) { if (bli_zreal(*beta) != 1. || bli_zimag(*beta) != 0.) {
if (*incy == 1) { if (*incy == 1) {
if (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) {
i__1 = leny; i__1 = leny;
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = i__; i__2 = i__;
@@ -1429,14 +1429,14 @@
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = i__; i__2 = i__;
i__3 = i__; i__3 = i__;
bli_zsets( (bli_zbla_real(*beta) * bli_zbla_real(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zbla_real(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zbla_real(y[i__3])), z__1 ); bli_zsets( (bli_zreal(*beta) * bli_zreal(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zreal(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zreal(y[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
/* L20: */ /* L20: */
} }
} }
} else { } else {
iy = ky; iy = ky;
if (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) {
i__1 = leny; i__1 = leny;
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = iy; i__2 = iy;
@@ -1449,15 +1449,15 @@
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = iy; i__2 = iy;
i__3 = iy; i__3 = iy;
bli_zsets( (bli_zbla_real(*beta) * bli_zbla_real(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zbla_real(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zbla_real(y[i__3])), z__1 ); bli_zsets( (bli_zreal(*beta) * bli_zreal(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zreal(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zreal(y[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
iy += *incy; iy += *incy;
/* L40: */ /* L40: */
} }
} }
} }
} }
if (bli_zbla_real(*alpha) == 0. && bli_zimag(*alpha) == 0.) { if (bli_zreal(*alpha) == 0. && bli_zimag(*alpha) == 0.) {
return 0; return 0;
} }
kup1 = *ku + 1; kup1 = *ku + 1;
@@ -1470,10 +1470,10 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zbla_real(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zbla_real(x[i__2])), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zreal(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zreal(x[i__2])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
k = kup1 - j; k = kup1 - j;
/* Computing MAX */ /* Computing MAX */
i__2 = 1, i__3 = j - *ku; i__2 = 1, i__3 = j - *ku;
@@ -1484,9 +1484,9 @@
i__2 = i__; i__2 = i__;
i__3 = i__; i__3 = i__;
i__5 = k + i__ + j * a_dim1; i__5 = k + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zbla_real(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zbla_real(a[i__5])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zreal(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zreal(a[i__5])), z__2 );
bli_zsets( (bli_zbla_real(y[i__3]) + bli_zbla_real(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
/* L50: */ /* L50: */
} }
} }
@@ -1497,10 +1497,10 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__4 = jx; i__4 = jx;
if (bli_zbla_real(x[i__4]) != 0. || bli_zimag(x[i__4]) != 0.) { if (bli_zreal(x[i__4]) != 0. || bli_zimag(x[i__4]) != 0.) {
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__4]) - bli_zimag(*alpha) * bli_zimag(x[i__4])), (bli_zbla_real(*alpha) * bli_zimag(x[i__4]) + bli_zimag(*alpha) * bli_zbla_real(x[i__4])), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__4]) - bli_zimag(*alpha) * bli_zimag(x[i__4])), (bli_zreal(*alpha) * bli_zimag(x[i__4]) + bli_zimag(*alpha) * bli_zreal(x[i__4])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
iy = ky; iy = ky;
k = kup1 - j; k = kup1 - j;
/* Computing MAX */ /* Computing MAX */
@@ -1512,9 +1512,9 @@
i__4 = iy; i__4 = iy;
i__2 = iy; i__2 = iy;
i__5 = k + i__ + j * a_dim1; i__5 = k + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zbla_real(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zbla_real(a[i__5])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zreal(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zreal(a[i__5])), z__2 );
bli_zsets( (bli_zbla_real(y[i__2]) + bli_zbla_real(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__2]) + bli_zreal(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__4] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] );
iy += *incy; iy += *incy;
/* L70: */ /* L70: */
} }
@@ -1545,9 +1545,9 @@
for (i__ = f2c_max(i__3,i__4); i__ <= i__2; ++i__) { for (i__ = f2c_max(i__3,i__4); i__ <= i__2; ++i__) {
i__3 = k + i__ + j * a_dim1; i__3 = k + i__ + j * a_dim1;
i__4 = i__; i__4 = i__;
bli_zsets( (bli_zbla_real(a[i__3]) * bli_zbla_real(x[i__4]) - bli_zimag(a[i__3]) * bli_zimag(x[i__4])), (bli_zbla_real(a[i__3]) * bli_zimag(x[i__4]) + bli_zimag(a[i__3]) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(a[i__3]) * bli_zreal(x[i__4]) - bli_zimag(a[i__3]) * bli_zimag(x[i__4])), (bli_zreal(a[i__3]) * bli_zimag(x[i__4]) + bli_zimag(a[i__3]) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L90: */ /* L90: */
} }
} else { } else {
@@ -1559,17 +1559,17 @@
for (i__ = f2c_max(i__2,i__3); i__ <= i__4; ++i__) { for (i__ = f2c_max(i__2,i__3); i__ <= i__4; ++i__) {
bla_d_cnjg(&z__3, &a[k + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[k + i__ + j * a_dim1]);
i__2 = i__; i__2 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zbla_real(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zreal(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L100: */ /* L100: */
} }
} }
i__4 = jy; i__4 = jy;
i__2 = jy; i__2 = jy;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(temp) - bli_zimag(*alpha) * bli_zimag(temp)), (bli_zbla_real(*alpha) * bli_zimag(temp) + bli_zimag(*alpha) * bli_zbla_real(temp)), z__2 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(temp) - bli_zimag(*alpha) * bli_zimag(temp)), (bli_zreal(*alpha) * bli_zimag(temp) + bli_zimag(*alpha) * bli_zreal(temp)), z__2 );
bli_zsets( (bli_zbla_real(y[i__2]) + bli_zbla_real(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__2]) + bli_zreal(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__4] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] );
jy += *incy; jy += *incy;
/* L110: */ /* L110: */
} }
@@ -1588,9 +1588,9 @@
for (i__ = f2c_max(i__4,i__2); i__ <= i__3; ++i__) { for (i__ = f2c_max(i__4,i__2); i__ <= i__3; ++i__) {
i__4 = k + i__ + j * a_dim1; i__4 = k + i__ + j * a_dim1;
i__2 = ix; i__2 = ix;
bli_zsets( (bli_zbla_real(a[i__4]) * bli_zbla_real(x[i__2]) - bli_zimag(a[i__4]) * bli_zimag(x[i__2])), (bli_zbla_real(a[i__4]) * bli_zimag(x[i__2]) + bli_zimag(a[i__4]) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(a[i__4]) * bli_zreal(x[i__2]) - bli_zimag(a[i__4]) * bli_zimag(x[i__2])), (bli_zreal(a[i__4]) * bli_zimag(x[i__2]) + bli_zimag(a[i__4]) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix += *incx; ix += *incx;
/* L120: */ /* L120: */
} }
@@ -1603,18 +1603,18 @@
for (i__ = f2c_max(i__3,i__4); i__ <= i__2; ++i__) { for (i__ = f2c_max(i__3,i__4); i__ <= i__2; ++i__) {
bla_d_cnjg(&z__3, &a[k + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[k + i__ + j * a_dim1]);
i__3 = ix; i__3 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zbla_real(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zreal(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix += *incx; ix += *incx;
/* L130: */ /* L130: */
} }
} }
i__2 = jy; i__2 = jy;
i__3 = jy; i__3 = jy;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(temp) - bli_zimag(*alpha) * bli_zimag(temp)), (bli_zbla_real(*alpha) * bli_zimag(temp) + bli_zimag(*alpha) * bli_zbla_real(temp)), z__2 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(temp) - bli_zimag(*alpha) * bli_zimag(temp)), (bli_zreal(*alpha) * bli_zimag(temp) + bli_zimag(*alpha) * bli_zreal(temp)), z__2 );
bli_zsets( (bli_zbla_real(y[i__3]) + bli_zbla_real(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
jy += *incy; jy += *incy;
if (j > *ku) { if (j > *ku) {
kx += *incx; kx += *incx;

134
frame/compat/f2c/bla_hbmv.c
View File

@@ -671,7 +671,7 @@
/* Quick return if possible. */ /* Quick return if possible. */
if (*n == 0 || (bli_zbla_real(*alpha) == 0. && bli_zimag(*alpha) == 0. && (bli_zbla_real(*beta) == 1. && if (*n == 0 || (bli_zreal(*alpha) == 0. && bli_zimag(*alpha) == 0. && (bli_zreal(*beta) == 1. &&
bli_zimag(*beta) == 0.))) { bli_zimag(*beta) == 0.))) {
return 0; return 0;
} }
@@ -694,9 +694,9 @@
/* First form y := beta*y. */ /* First form y := beta*y. */
if (bli_zbla_real(*beta) != 1. || bli_zimag(*beta) != 0.) { if (bli_zreal(*beta) != 1. || bli_zimag(*beta) != 0.) {
if (*incy == 1) { if (*incy == 1) {
if (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) {
i__1 = *n; i__1 = *n;
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = i__; i__2 = i__;
@@ -708,14 +708,14 @@
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = i__; i__2 = i__;
i__3 = i__; i__3 = i__;
bli_zsets( (bli_zbla_real(*beta) * bli_zbla_real(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zbla_real(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zbla_real(y[i__3])), z__1 ); bli_zsets( (bli_zreal(*beta) * bli_zreal(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zreal(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zreal(y[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
/* L20: */ /* L20: */
} }
} }
} else { } else {
iy = ky; iy = ky;
if (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) {
i__1 = *n; i__1 = *n;
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = iy; i__2 = iy;
@@ -728,15 +728,15 @@
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = iy; i__2 = iy;
i__3 = iy; i__3 = iy;
bli_zsets( (bli_zbla_real(*beta) * bli_zbla_real(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zbla_real(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zbla_real(y[i__3])), z__1 ); bli_zsets( (bli_zreal(*beta) * bli_zreal(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zreal(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zreal(y[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
iy += *incy; iy += *incy;
/* L40: */ /* L40: */
} }
} }
} }
} }
if (bli_zbla_real(*alpha) == 0. && bli_zimag(*alpha) == 0.) { if (bli_zreal(*alpha) == 0. && bli_zimag(*alpha) == 0.) {
return 0; return 0;
} }
if (PASTEF770(lsame)(uplo, "U", (ftnlen)1, (ftnlen)1)) { if (PASTEF770(lsame)(uplo, "U", (ftnlen)1, (ftnlen)1)) {
@@ -748,8 +748,8 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zbla_real(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zbla_real(x[i__2])), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zreal(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zreal(x[i__2])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
bli_zsets( (0.), (0.), temp2 ); bli_zsets( (0.), (0.), temp2 );
l = kplus1 - j; l = kplus1 - j;
/* Computing MAX */ /* Computing MAX */
@@ -759,25 +759,25 @@
i__2 = i__; i__2 = i__;
i__3 = i__; i__3 = i__;
i__5 = l + i__ + j * a_dim1; i__5 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp1) * bli_zbla_real(a[i__5]) - bli_zimag(temp1) * bli_zimag(a[i__5])), (bli_zbla_real(temp1) * bli_zimag(a[i__5]) + bli_zimag(temp1) * bli_zbla_real(a[i__5])), z__2 ); bli_zsets( (bli_zreal(temp1) * bli_zreal(a[i__5]) - bli_zimag(temp1) * bli_zimag(a[i__5])), (bli_zreal(temp1) * bli_zimag(a[i__5]) + bli_zimag(temp1) * bli_zreal(a[i__5])), z__2 );
bli_zsets( (bli_zbla_real(y[i__3]) + bli_zbla_real(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__2 = i__; i__2 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zbla_real(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zreal(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp2) + bli_zbla_real(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
/* L50: */ /* L50: */
} }
i__4 = j; i__4 = j;
i__2 = j; i__2 = j;
i__3 = kplus1 + j * a_dim1; i__3 = kplus1 + j * a_dim1;
d__1 = bli_zbla_real(a[i__3]); d__1 = bli_zreal(a[i__3]);
bli_zsets( (d__1 * bli_zbla_real(temp1)), (d__1 * bli_zimag(temp1)), z__3 ); bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__3 );
bli_zsets( (bli_zbla_real(y[i__2]) + bli_zbla_real(z__3)), (bli_zimag(y[i__2]) + bli_zimag(z__3)), z__2 ); bli_zsets( (bli_zreal(y[i__2]) + bli_zreal(z__3)), (bli_zimag(y[i__2]) + bli_zimag(z__3)), z__2 );
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zbla_real(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zbla_real(temp2)), z__4 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zreal(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zreal(temp2)), z__4 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__4] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] );
/* L60: */ /* L60: */
} }
} else { } else {
@@ -786,8 +786,8 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__4]) - bli_zimag(*alpha) * bli_zimag(x[i__4])), (bli_zbla_real(*alpha) * bli_zimag(x[i__4]) + bli_zimag(*alpha) * bli_zbla_real(x[i__4])), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__4]) - bli_zimag(*alpha) * bli_zimag(x[i__4])), (bli_zreal(*alpha) * bli_zimag(x[i__4]) + bli_zimag(*alpha) * bli_zreal(x[i__4])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
bli_zsets( (0.), (0.), temp2 ); bli_zsets( (0.), (0.), temp2 );
ix = kx; ix = kx;
iy = ky; iy = ky;
@@ -799,14 +799,14 @@
i__4 = iy; i__4 = iy;
i__2 = iy; i__2 = iy;
i__5 = l + i__ + j * a_dim1; i__5 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp1) * bli_zbla_real(a[i__5]) - bli_zimag(temp1) * bli_zimag(a[i__5])), (bli_zbla_real(temp1) * bli_zimag(a[i__5]) + bli_zimag(temp1) * bli_zbla_real(a[i__5])), z__2 ); bli_zsets( (bli_zreal(temp1) * bli_zreal(a[i__5]) - bli_zimag(temp1) * bli_zimag(a[i__5])), (bli_zreal(temp1) * bli_zimag(a[i__5]) + bli_zimag(temp1) * bli_zreal(a[i__5])), z__2 );
bli_zsets( (bli_zbla_real(y[i__2]) + bli_zbla_real(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__2]) + bli_zreal(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__4] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] );
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__4 = ix; i__4 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zbla_real(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zreal(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp2) + bli_zbla_real(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
ix += *incx; ix += *incx;
iy += *incy; iy += *incy;
/* L70: */ /* L70: */
@@ -814,12 +814,12 @@
i__3 = jy; i__3 = jy;
i__4 = jy; i__4 = jy;
i__2 = kplus1 + j * a_dim1; i__2 = kplus1 + j * a_dim1;
d__1 = bli_zbla_real(a[i__2]); d__1 = bli_zreal(a[i__2]);
bli_zsets( (d__1 * bli_zbla_real(temp1)), (d__1 * bli_zimag(temp1)), z__3 ); bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__3 );
bli_zsets( (bli_zbla_real(y[i__4]) + bli_zbla_real(z__3)), (bli_zimag(y[i__4]) + bli_zimag(z__3)), z__2 ); bli_zsets( (bli_zreal(y[i__4]) + bli_zreal(z__3)), (bli_zimag(y[i__4]) + bli_zimag(z__3)), z__2 );
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zbla_real(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zbla_real(temp2)), z__4 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zreal(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zreal(temp2)), z__4 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] );
jx += *incx; jx += *incx;
jy += *incy; jy += *incy;
if (j > *k) { if (j > *k) {
@@ -837,16 +837,16 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__3 = j; i__3 = j;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__3]) - bli_zimag(*alpha) * bli_zimag(x[i__3])), (bli_zbla_real(*alpha) * bli_zimag(x[i__3]) + bli_zimag(*alpha) * bli_zbla_real(x[i__3])), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__3]) - bli_zimag(*alpha) * bli_zimag(x[i__3])), (bli_zreal(*alpha) * bli_zimag(x[i__3]) + bli_zimag(*alpha) * bli_zreal(x[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
bli_zsets( (0.), (0.), temp2 ); bli_zsets( (0.), (0.), temp2 );
i__3 = j; i__3 = j;
i__4 = j; i__4 = j;
i__2 = j * a_dim1 + 1; i__2 = j * a_dim1 + 1;
d__1 = bli_zbla_real(a[i__2]); d__1 = bli_zreal(a[i__2]);
bli_zsets( (d__1 * bli_zbla_real(temp1)), (d__1 * bli_zimag(temp1)), z__2 ); bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__2 );
bli_zsets( (bli_zbla_real(y[i__4]) + bli_zbla_real(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] );
l = 1 - j; l = 1 - j;
/* Computing MIN */ /* Computing MIN */
i__4 = *n, i__2 = j + *k; i__4 = *n, i__2 = j + *k;
@@ -855,21 +855,21 @@
i__4 = i__; i__4 = i__;
i__2 = i__; i__2 = i__;
i__5 = l + i__ + j * a_dim1; i__5 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp1) * bli_zbla_real(a[i__5]) - bli_zimag(temp1) * bli_zimag(a[i__5])), (bli_zbla_real(temp1) * bli_zimag(a[i__5]) + bli_zimag(temp1) * bli_zbla_real(a[i__5])), z__2 ); bli_zsets( (bli_zreal(temp1) * bli_zreal(a[i__5]) - bli_zimag(temp1) * bli_zimag(a[i__5])), (bli_zreal(temp1) * bli_zimag(a[i__5]) + bli_zimag(temp1) * bli_zreal(a[i__5])), z__2 );
bli_zsets( (bli_zbla_real(y[i__2]) + bli_zbla_real(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__2]) + bli_zreal(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__4] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] );
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__4 = i__; i__4 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zbla_real(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zreal(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp2) + bli_zbla_real(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
/* L90: */ /* L90: */
} }
i__3 = j; i__3 = j;
i__4 = j; i__4 = j;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zbla_real(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zbla_real(temp2)), z__2 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zreal(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zreal(temp2)), z__2 );
bli_zsets( (bli_zbla_real(y[i__4]) + bli_zbla_real(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] );
/* L100: */ /* L100: */
} }
} else { } else {
@@ -878,16 +878,16 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__3 = jx; i__3 = jx;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__3]) - bli_zimag(*alpha) * bli_zimag(x[i__3])), (bli_zbla_real(*alpha) * bli_zimag(x[i__3]) + bli_zimag(*alpha) * bli_zbla_real(x[i__3])), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__3]) - bli_zimag(*alpha) * bli_zimag(x[i__3])), (bli_zreal(*alpha) * bli_zimag(x[i__3]) + bli_zimag(*alpha) * bli_zreal(x[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
bli_zsets( (0.), (0.), temp2 ); bli_zsets( (0.), (0.), temp2 );
i__3 = jy; i__3 = jy;
i__4 = jy; i__4 = jy;
i__2 = j * a_dim1 + 1; i__2 = j * a_dim1 + 1;
d__1 = bli_zbla_real(a[i__2]); d__1 = bli_zreal(a[i__2]);
bli_zsets( (d__1 * bli_zbla_real(temp1)), (d__1 * bli_zimag(temp1)), z__2 ); bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__2 );
bli_zsets( (bli_zbla_real(y[i__4]) + bli_zbla_real(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] );
l = 1 - j; l = 1 - j;
ix = jx; ix = jx;
iy = jy; iy = jy;
@@ -900,21 +900,21 @@
i__4 = iy; i__4 = iy;
i__2 = iy; i__2 = iy;
i__5 = l + i__ + j * a_dim1; i__5 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp1) * bli_zbla_real(a[i__5]) - bli_zimag(temp1) * bli_zimag(a[i__5])), (bli_zbla_real(temp1) * bli_zimag(a[i__5]) + bli_zimag(temp1) * bli_zbla_real(a[i__5])), z__2 ); bli_zsets( (bli_zreal(temp1) * bli_zreal(a[i__5]) - bli_zimag(temp1) * bli_zimag(a[i__5])), (bli_zreal(temp1) * bli_zimag(a[i__5]) + bli_zimag(temp1) * bli_zreal(a[i__5])), z__2 );
bli_zsets( (bli_zbla_real(y[i__2]) + bli_zbla_real(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__2]) + bli_zreal(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__4] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] );
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__4 = ix; i__4 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zbla_real(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zreal(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp2) + bli_zbla_real(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
/* L110: */ /* L110: */
} }
i__3 = jy; i__3 = jy;
i__4 = jy; i__4 = jy;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zbla_real(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zbla_real(temp2)), z__2 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zreal(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zreal(temp2)), z__2 );
bli_zsets( (bli_zbla_real(y[i__4]) + bli_zbla_real(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] );
jx += *incx; jx += *incx;
jy += *incy; jy += *incy;
/* L120: */ /* L120: */

134
frame/compat/f2c/bla_hpmv.c
View File

@@ -583,7 +583,7 @@
/* Quick return if possible. */ /* Quick return if possible. */
if (*n == 0 || (bli_zbla_real(*alpha) == 0. && bli_zimag(*alpha) == 0. && (bli_zbla_real(*beta) == 1. && if (*n == 0 || (bli_zreal(*alpha) == 0. && bli_zimag(*alpha) == 0. && (bli_zreal(*beta) == 1. &&
bli_zimag(*beta) == 0.))) { bli_zimag(*beta) == 0.))) {
return 0; return 0;
} }
@@ -606,9 +606,9 @@
/* First form y := beta*y. */ /* First form y := beta*y. */
if (bli_zbla_real(*beta) != 1. || bli_zimag(*beta) != 0.) { if (bli_zreal(*beta) != 1. || bli_zimag(*beta) != 0.) {
if (*incy == 1) { if (*incy == 1) {
if (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) {
i__1 = *n; i__1 = *n;
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = i__; i__2 = i__;
@@ -620,14 +620,14 @@
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = i__; i__2 = i__;
i__3 = i__; i__3 = i__;
bli_zsets( (bli_zbla_real(*beta) * bli_zbla_real(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zbla_real(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zbla_real(y[i__3])), z__1 ); bli_zsets( (bli_zreal(*beta) * bli_zreal(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zreal(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zreal(y[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
/* L20: */ /* L20: */
} }
} }
} else { } else {
iy = ky; iy = ky;
if (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) {
i__1 = *n; i__1 = *n;
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = iy; i__2 = iy;
@@ -640,15 +640,15 @@
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = iy; i__2 = iy;
i__3 = iy; i__3 = iy;
bli_zsets( (bli_zbla_real(*beta) * bli_zbla_real(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zbla_real(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zbla_real(y[i__3])), z__1 ); bli_zsets( (bli_zreal(*beta) * bli_zreal(y[i__3]) - bli_zimag(*beta) * bli_zimag(y[i__3])), (bli_zreal(*beta) * bli_zimag(y[i__3]) + bli_zimag(*beta) * bli_zreal(y[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
iy += *incy; iy += *incy;
/* L40: */ /* L40: */
} }
} }
} }
} }
if (bli_zbla_real(*alpha) == 0. && bli_zimag(*alpha) == 0.) { if (bli_zreal(*alpha) == 0. && bli_zimag(*alpha) == 0.) {
return 0; return 0;
} }
kk = 1; kk = 1;
@@ -660,8 +660,8 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zbla_real(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zbla_real(x[i__2])), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zreal(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zreal(x[i__2])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
bli_zsets( (0.), (0.), temp2 ); bli_zsets( (0.), (0.), temp2 );
k = kk; k = kk;
i__2 = j - 1; i__2 = j - 1;
@@ -669,26 +669,26 @@
i__3 = i__; i__3 = i__;
i__4 = i__; i__4 = i__;
i__5 = k; i__5 = k;
bli_zsets( (bli_zbla_real(temp1) * bli_zbla_real(ap[i__5]) - bli_zimag(temp1) * bli_zimag(ap[i__5])), (bli_zbla_real(temp1) * bli_zimag(ap[i__5]) + bli_zimag(temp1) * bli_zbla_real(ap[i__5])), z__2 ); bli_zsets( (bli_zreal(temp1) * bli_zreal(ap[i__5]) - bli_zimag(temp1) * bli_zimag(ap[i__5])), (bli_zreal(temp1) * bli_zimag(ap[i__5]) + bli_zimag(temp1) * bli_zreal(ap[i__5])), z__2 );
bli_zsets( (bli_zbla_real(y[i__4]) + bli_zbla_real(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] );
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__3 = i__; i__3 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zbla_real(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zreal(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp2) + bli_zbla_real(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
++k; ++k;
/* L50: */ /* L50: */
} }
i__2 = j; i__2 = j;
i__3 = j; i__3 = j;
i__4 = kk + j - 1; i__4 = kk + j - 1;
d__1 = bli_zbla_real(ap[i__4]); d__1 = bli_zreal(ap[i__4]);
bli_zsets( (d__1 * bli_zbla_real(temp1)), (d__1 * bli_zimag(temp1)), z__3 ); bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__3 );
bli_zsets( (bli_zbla_real(y[i__3]) + bli_zbla_real(z__3)), (bli_zimag(y[i__3]) + bli_zimag(z__3)), z__2 ); bli_zsets( (bli_zreal(y[i__3]) + bli_zreal(z__3)), (bli_zimag(y[i__3]) + bli_zimag(z__3)), z__2 );
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zbla_real(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zbla_real(temp2)), z__4 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zreal(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zreal(temp2)), z__4 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
kk += j; kk += j;
/* L60: */ /* L60: */
} }
@@ -698,8 +698,8 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zbla_real(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zbla_real(x[i__2])), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zreal(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zreal(x[i__2])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
bli_zsets( (0.), (0.), temp2 ); bli_zsets( (0.), (0.), temp2 );
ix = kx; ix = kx;
iy = ky; iy = ky;
@@ -708,14 +708,14 @@
i__3 = iy; i__3 = iy;
i__4 = iy; i__4 = iy;
i__5 = k; i__5 = k;
bli_zsets( (bli_zbla_real(temp1) * bli_zbla_real(ap[i__5]) - bli_zimag(temp1) * bli_zimag(ap[i__5])), (bli_zbla_real(temp1) * bli_zimag(ap[i__5]) + bli_zimag(temp1) * bli_zbla_real(ap[i__5])), z__2 ); bli_zsets( (bli_zreal(temp1) * bli_zreal(ap[i__5]) - bli_zimag(temp1) * bli_zimag(ap[i__5])), (bli_zreal(temp1) * bli_zimag(ap[i__5]) + bli_zimag(temp1) * bli_zreal(ap[i__5])), z__2 );
bli_zsets( (bli_zbla_real(y[i__4]) + bli_zbla_real(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] );
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__3 = ix; i__3 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zbla_real(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zreal(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp2) + bli_zbla_real(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
ix += *incx; ix += *incx;
iy += *incy; iy += *incy;
/* L70: */ /* L70: */
@@ -723,12 +723,12 @@
i__2 = jy; i__2 = jy;
i__3 = jy; i__3 = jy;
i__4 = kk + j - 1; i__4 = kk + j - 1;
d__1 = bli_zbla_real(ap[i__4]); d__1 = bli_zreal(ap[i__4]);
bli_zsets( (d__1 * bli_zbla_real(temp1)), (d__1 * bli_zimag(temp1)), z__3 ); bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__3 );
bli_zsets( (bli_zbla_real(y[i__3]) + bli_zbla_real(z__3)), (bli_zimag(y[i__3]) + bli_zimag(z__3)), z__2 ); bli_zsets( (bli_zreal(y[i__3]) + bli_zreal(z__3)), (bli_zimag(y[i__3]) + bli_zimag(z__3)), z__2 );
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zbla_real(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zbla_real(temp2)), z__4 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zreal(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zreal(temp2)), z__4 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
jx += *incx; jx += *incx;
jy += *incy; jy += *incy;
kk += j; kk += j;
@@ -743,38 +743,38 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zbla_real(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zbla_real(x[i__2])), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zreal(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zreal(x[i__2])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
bli_zsets( (0.), (0.), temp2 ); bli_zsets( (0.), (0.), temp2 );
i__2 = j; i__2 = j;
i__3 = j; i__3 = j;
i__4 = kk; i__4 = kk;
d__1 = bli_zbla_real(ap[i__4]); d__1 = bli_zreal(ap[i__4]);
bli_zsets( (d__1 * bli_zbla_real(temp1)), (d__1 * bli_zimag(temp1)), z__2 ); bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__2 );
bli_zsets( (bli_zbla_real(y[i__3]) + bli_zbla_real(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
k = kk + 1; k = kk + 1;
i__2 = *n; i__2 = *n;
for (i__ = j + 1; i__ <= i__2; ++i__) { for (i__ = j + 1; i__ <= i__2; ++i__) {
i__3 = i__; i__3 = i__;
i__4 = i__; i__4 = i__;
i__5 = k; i__5 = k;
bli_zsets( (bli_zbla_real(temp1) * bli_zbla_real(ap[i__5]) - bli_zimag(temp1) * bli_zimag(ap[i__5])), (bli_zbla_real(temp1) * bli_zimag(ap[i__5]) + bli_zimag(temp1) * bli_zbla_real(ap[i__5])), z__2 ); bli_zsets( (bli_zreal(temp1) * bli_zreal(ap[i__5]) - bli_zimag(temp1) * bli_zimag(ap[i__5])), (bli_zreal(temp1) * bli_zimag(ap[i__5]) + bli_zimag(temp1) * bli_zreal(ap[i__5])), z__2 );
bli_zsets( (bli_zbla_real(y[i__4]) + bli_zbla_real(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] );
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__3 = i__; i__3 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zbla_real(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zreal(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp2) + bli_zbla_real(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
++k; ++k;
/* L90: */ /* L90: */
} }
i__2 = j; i__2 = j;
i__3 = j; i__3 = j;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zbla_real(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zbla_real(temp2)), z__2 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zreal(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zreal(temp2)), z__2 );
bli_zsets( (bli_zbla_real(y[i__3]) + bli_zbla_real(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
kk += *n - j + 1; kk += *n - j + 1;
/* L100: */ /* L100: */
} }
@@ -784,16 +784,16 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zbla_real(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zbla_real(x[i__2])), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zreal(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zreal(x[i__2])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
bli_zsets( (0.), (0.), temp2 ); bli_zsets( (0.), (0.), temp2 );
i__2 = jy; i__2 = jy;
i__3 = jy; i__3 = jy;
i__4 = kk; i__4 = kk;
d__1 = bli_zbla_real(ap[i__4]); d__1 = bli_zreal(ap[i__4]);
bli_zsets( (d__1 * bli_zbla_real(temp1)), (d__1 * bli_zimag(temp1)), z__2 ); bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__2 );
bli_zsets( (bli_zbla_real(y[i__3]) + bli_zbla_real(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
ix = jx; ix = jx;
iy = jy; iy = jy;
i__2 = kk + *n - j; i__2 = kk + *n - j;
@@ -803,21 +803,21 @@
i__3 = iy; i__3 = iy;
i__4 = iy; i__4 = iy;
i__5 = k; i__5 = k;
bli_zsets( (bli_zbla_real(temp1) * bli_zbla_real(ap[i__5]) - bli_zimag(temp1) * bli_zimag(ap[i__5])), (bli_zbla_real(temp1) * bli_zimag(ap[i__5]) + bli_zimag(temp1) * bli_zbla_real(ap[i__5])), z__2 ); bli_zsets( (bli_zreal(temp1) * bli_zreal(ap[i__5]) - bli_zimag(temp1) * bli_zimag(ap[i__5])), (bli_zreal(temp1) * bli_zimag(ap[i__5]) + bli_zimag(temp1) * bli_zreal(ap[i__5])), z__2 );
bli_zsets( (bli_zbla_real(y[i__4]) + bli_zbla_real(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] );
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__3 = ix; i__3 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zbla_real(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zreal(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp2) + bli_zbla_real(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
/* L110: */ /* L110: */
} }
i__2 = jy; i__2 = jy;
i__3 = jy; i__3 = jy;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zbla_real(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zbla_real(temp2)), z__2 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(temp2) - bli_zimag(*alpha) * bli_zimag(temp2)), (bli_zreal(*alpha) * bli_zimag(temp2) + bli_zimag(*alpha) * bli_zreal(temp2)), z__2 );
bli_zsets( (bli_zbla_real(y[i__3]) + bli_zbla_real(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] );
jx += *incx; jx += *incx;
jy += *incy; jy += *incy;
kk += *n - j + 1; kk += *n - j + 1;

72
frame/compat/f2c/bla_hpr.c
View File

@@ -506,32 +506,32 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
bla_d_cnjg(&z__2, &x[j]); bla_d_cnjg(&z__2, &x[j]);
bli_zsets( (*alpha * bli_zbla_real(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); bli_zsets( (*alpha * bli_zreal(z__2)), (*alpha * bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
k = kk; k = kk;
i__2 = j - 1; i__2 = j - 1;
for (i__ = 1; i__ <= i__2; ++i__) { for (i__ = 1; i__ <= i__2; ++i__) {
i__3 = k; i__3 = k;
i__4 = k; i__4 = k;
i__5 = i__; i__5 = i__;
bli_zsets( (bli_zbla_real(x[i__5]) * bli_zbla_real(temp) - bli_zimag(x[i__5]) * bli_zimag(temp)), (bli_zbla_real(x[i__5]) * bli_zimag(temp) + bli_zimag(x[i__5]) * bli_zbla_real(temp)), z__2 ); bli_zsets( (bli_zreal(x[i__5]) * bli_zreal(temp) - bli_zimag(x[i__5]) * bli_zimag(temp)), (bli_zreal(x[i__5]) * bli_zimag(temp) + bli_zimag(x[i__5]) * bli_zreal(temp)), z__2 );
bli_zsets( (bli_zbla_real(ap[i__4]) + bli_zbla_real(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(ap[i__4]) + bli_zreal(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), ap[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] );
++k; ++k;
/* L10: */ /* L10: */
} }
i__2 = kk + j - 1; i__2 = kk + j - 1;
i__3 = kk + j - 1; i__3 = kk + j - 1;
i__4 = j; i__4 = j;
bli_zsets( (bli_zbla_real(x[i__4]) * bli_zbla_real(temp) - bli_zimag(x[i__4]) * bli_zimag(temp)), (bli_zbla_real(x[i__4]) * bli_zimag(temp) + bli_zimag(x[i__4]) * bli_zbla_real(temp)), z__1 ); bli_zsets( (bli_zreal(x[i__4]) * bli_zreal(temp) - bli_zimag(x[i__4]) * bli_zimag(temp)), (bli_zreal(x[i__4]) * bli_zimag(temp) + bli_zimag(x[i__4]) * bli_zreal(temp)), z__1 );
d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(z__1); d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} else { } else {
i__2 = kk + j - 1; i__2 = kk + j - 1;
i__3 = kk + j - 1; i__3 = kk + j - 1;
d__1 = bli_zbla_real(ap[i__3]); d__1 = bli_zreal(ap[i__3]);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} }
kk += j; kk += j;
@@ -542,32 +542,32 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
bla_d_cnjg(&z__2, &x[jx]); bla_d_cnjg(&z__2, &x[jx]);
bli_zsets( (*alpha * bli_zbla_real(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); bli_zsets( (*alpha * bli_zreal(z__2)), (*alpha * bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix = kx; ix = kx;
i__2 = kk + j - 2; i__2 = kk + j - 2;
for (k = kk; k <= i__2; ++k) { for (k = kk; k <= i__2; ++k) {
i__3 = k; i__3 = k;
i__4 = k; i__4 = k;
i__5 = ix; i__5 = ix;
bli_zsets( (bli_zbla_real(x[i__5]) * bli_zbla_real(temp) - bli_zimag(x[i__5]) * bli_zimag(temp)), (bli_zbla_real(x[i__5]) * bli_zimag(temp) + bli_zimag(x[i__5]) * bli_zbla_real(temp)), z__2 ); bli_zsets( (bli_zreal(x[i__5]) * bli_zreal(temp) - bli_zimag(x[i__5]) * bli_zimag(temp)), (bli_zreal(x[i__5]) * bli_zimag(temp) + bli_zimag(x[i__5]) * bli_zreal(temp)), z__2 );
bli_zsets( (bli_zbla_real(ap[i__4]) + bli_zbla_real(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(ap[i__4]) + bli_zreal(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), ap[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] );
ix += *incx; ix += *incx;
/* L30: */ /* L30: */
} }
i__2 = kk + j - 1; i__2 = kk + j - 1;
i__3 = kk + j - 1; i__3 = kk + j - 1;
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(x[i__4]) * bli_zbla_real(temp) - bli_zimag(x[i__4]) * bli_zimag(temp)), (bli_zbla_real(x[i__4]) * bli_zimag(temp) + bli_zimag(x[i__4]) * bli_zbla_real(temp)), z__1 ); bli_zsets( (bli_zreal(x[i__4]) * bli_zreal(temp) - bli_zimag(x[i__4]) * bli_zimag(temp)), (bli_zreal(x[i__4]) * bli_zimag(temp) + bli_zimag(x[i__4]) * bli_zreal(temp)), z__1 );
d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(z__1); d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} else { } else {
i__2 = kk + j - 1; i__2 = kk + j - 1;
i__3 = kk + j - 1; i__3 = kk + j - 1;
d__1 = bli_zbla_real(ap[i__3]); d__1 = bli_zreal(ap[i__3]);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} }
jx += *incx; jx += *incx;
@@ -583,15 +583,15 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
bla_d_cnjg(&z__2, &x[j]); bla_d_cnjg(&z__2, &x[j]);
bli_zsets( (*alpha * bli_zbla_real(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); bli_zsets( (*alpha * bli_zreal(z__2)), (*alpha * bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
i__2 = kk; i__2 = kk;
i__3 = kk; i__3 = kk;
i__4 = j; i__4 = j;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(x[i__4]) - bli_zimag(temp) * bli_zimag(x[i__4])), (bli_zbla_real(temp) * bli_zimag(x[i__4]) + bli_zimag(temp) * bli_zbla_real(x[i__4])), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(x[i__4]) - bli_zimag(temp) * bli_zimag(x[i__4])), (bli_zreal(temp) * bli_zimag(x[i__4]) + bli_zimag(temp) * bli_zreal(x[i__4])), z__1 );
d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(z__1); d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
k = kk + 1; k = kk + 1;
i__2 = *n; i__2 = *n;
@@ -599,16 +599,16 @@
i__3 = k; i__3 = k;
i__4 = k; i__4 = k;
i__5 = i__; i__5 = i__;
bli_zsets( (bli_zbla_real(x[i__5]) * bli_zbla_real(temp) - bli_zimag(x[i__5]) * bli_zimag(temp)), (bli_zbla_real(x[i__5]) * bli_zimag(temp) + bli_zimag(x[i__5]) * bli_zbla_real(temp)), z__2 ); bli_zsets( (bli_zreal(x[i__5]) * bli_zreal(temp) - bli_zimag(x[i__5]) * bli_zimag(temp)), (bli_zreal(x[i__5]) * bli_zimag(temp) + bli_zimag(x[i__5]) * bli_zreal(temp)), z__2 );
bli_zsets( (bli_zbla_real(ap[i__4]) + bli_zbla_real(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(ap[i__4]) + bli_zreal(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), ap[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] );
++k; ++k;
/* L50: */ /* L50: */
} }
} else { } else {
i__2 = kk; i__2 = kk;
i__3 = kk; i__3 = kk;
d__1 = bli_zbla_real(ap[i__3]); d__1 = bli_zreal(ap[i__3]);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} }
kk = kk + *n - j + 1; kk = kk + *n - j + 1;
@@ -619,15 +619,15 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
bla_d_cnjg(&z__2, &x[jx]); bla_d_cnjg(&z__2, &x[jx]);
bli_zsets( (*alpha * bli_zbla_real(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); bli_zsets( (*alpha * bli_zreal(z__2)), (*alpha * bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
i__2 = kk; i__2 = kk;
i__3 = kk; i__3 = kk;
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(x[i__4]) - bli_zimag(temp) * bli_zimag(x[i__4])), (bli_zbla_real(temp) * bli_zimag(x[i__4]) + bli_zimag(temp) * bli_zbla_real(x[i__4])), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(x[i__4]) - bli_zimag(temp) * bli_zimag(x[i__4])), (bli_zreal(temp) * bli_zimag(x[i__4]) + bli_zimag(temp) * bli_zreal(x[i__4])), z__1 );
d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(z__1); d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
ix = jx; ix = jx;
i__2 = kk + *n - j; i__2 = kk + *n - j;
@@ -636,15 +636,15 @@
i__3 = k; i__3 = k;
i__4 = k; i__4 = k;
i__5 = ix; i__5 = ix;
bli_zsets( (bli_zbla_real(x[i__5]) * bli_zbla_real(temp) - bli_zimag(x[i__5]) * bli_zimag(temp)), (bli_zbla_real(x[i__5]) * bli_zimag(temp) + bli_zimag(x[i__5]) * bli_zbla_real(temp)), z__2 ); bli_zsets( (bli_zreal(x[i__5]) * bli_zreal(temp) - bli_zimag(x[i__5]) * bli_zimag(temp)), (bli_zreal(x[i__5]) * bli_zimag(temp) + bli_zimag(x[i__5]) * bli_zreal(temp)), z__2 );
bli_zsets( (bli_zbla_real(ap[i__4]) + bli_zbla_real(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(ap[i__4]) + bli_zreal(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), ap[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] );
/* L70: */ /* L70: */
} }
} else { } else {
i__2 = kk; i__2 = kk;
i__3 = kk; i__3 = kk;
d__1 = bli_zbla_real(ap[i__3]); d__1 = bli_zreal(ap[i__3]);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} }
jx += *incx; jx += *incx;

122
frame/compat/f2c/bla_hpr2.c
View File

@@ -572,7 +572,7 @@
/* Quick return if possible. */ /* Quick return if possible. */
if (*n == 0 || (bli_zbla_real(*alpha) == 0. && bli_zimag(*alpha) == 0.)) { if (*n == 0 || (bli_zreal(*alpha) == 0. && bli_zimag(*alpha) == 0.)) {
return 0; return 0;
} }
@@ -607,43 +607,43 @@
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
i__3 = j; i__3 = j;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0. || (bli_zbla_real(y[i__3]) != 0. || if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0. || (bli_zreal(y[i__3]) != 0. ||
bli_zimag(y[i__3]) != 0.)) { bli_zimag(y[i__3]) != 0.)) {
bla_d_cnjg(&z__2, &y[j]); bla_d_cnjg(&z__2, &y[j]);
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(z__2) - bli_zimag(*alpha) * bli_zimag(z__2)), (bli_zbla_real(*alpha) * bli_zimag(z__2) + bli_zimag(*alpha) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(z__2) - bli_zimag(*alpha) * bli_zimag(z__2)), (bli_zreal(*alpha) * bli_zimag(z__2) + bli_zimag(*alpha) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zbla_real(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zreal(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zreal(x[i__2])), z__2 );
bla_d_cnjg(&z__1, &z__2); bla_d_cnjg(&z__1, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
k = kk; k = kk;
i__2 = j - 1; i__2 = j - 1;
for (i__ = 1; i__ <= i__2; ++i__) { for (i__ = 1; i__ <= i__2; ++i__) {
i__3 = k; i__3 = k;
i__4 = k; i__4 = k;
i__5 = i__; i__5 = i__;
bli_zsets( (bli_zbla_real(x[i__5]) * bli_zbla_real(temp1) - bli_zimag(x[i__5]) * bli_zimag(temp1)), (bli_zbla_real(x[i__5]) * bli_zimag(temp1) + bli_zimag(x[i__5]) * bli_zbla_real(temp1)), z__3 ); bli_zsets( (bli_zreal(x[i__5]) * bli_zreal(temp1) - bli_zimag(x[i__5]) * bli_zimag(temp1)), (bli_zreal(x[i__5]) * bli_zimag(temp1) + bli_zimag(x[i__5]) * bli_zreal(temp1)), z__3 );
bli_zsets( (bli_zbla_real(ap[i__4]) + bli_zbla_real(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 ); bli_zsets( (bli_zreal(ap[i__4]) + bli_zreal(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 );
i__6 = i__; i__6 = i__;
bli_zsets( (bli_zbla_real(y[i__6]) * bli_zbla_real(temp2) - bli_zimag(y[i__6]) * bli_zimag(temp2)), (bli_zbla_real(y[i__6]) * bli_zimag(temp2) + bli_zimag(y[i__6]) * bli_zbla_real(temp2)), z__4 ); bli_zsets( (bli_zreal(y[i__6]) * bli_zreal(temp2) - bli_zimag(y[i__6]) * bli_zimag(temp2)), (bli_zreal(y[i__6]) * bli_zimag(temp2) + bli_zimag(y[i__6]) * bli_zreal(temp2)), z__4 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), ap[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] );
++k; ++k;
/* L10: */ /* L10: */
} }
i__2 = kk + j - 1; i__2 = kk + j - 1;
i__3 = kk + j - 1; i__3 = kk + j - 1;
i__4 = j; i__4 = j;
bli_zsets( (bli_zbla_real(x[i__4]) * bli_zbla_real(temp1) - bli_zimag(x[i__4]) * bli_zimag(temp1)), (bli_zbla_real(x[i__4]) * bli_zimag(temp1) + bli_zimag(x[i__4]) * bli_zbla_real(temp1)), z__2 ); bli_zsets( (bli_zreal(x[i__4]) * bli_zreal(temp1) - bli_zimag(x[i__4]) * bli_zimag(temp1)), (bli_zreal(x[i__4]) * bli_zimag(temp1) + bli_zimag(x[i__4]) * bli_zreal(temp1)), z__2 );
i__5 = j; i__5 = j;
bli_zsets( (bli_zbla_real(y[i__5]) * bli_zbla_real(temp2) - bli_zimag(y[i__5]) * bli_zimag(temp2)), (bli_zbla_real(y[i__5]) * bli_zimag(temp2) + bli_zimag(y[i__5]) * bli_zbla_real(temp2)), z__3 ); bli_zsets( (bli_zreal(y[i__5]) * bli_zreal(temp2) - bli_zimag(y[i__5]) * bli_zimag(temp2)), (bli_zreal(y[i__5]) * bli_zimag(temp2) + bli_zimag(y[i__5]) * bli_zreal(temp2)), z__3 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 );
d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(z__1); d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} else { } else {
i__2 = kk + j - 1; i__2 = kk + j - 1;
i__3 = kk + j - 1; i__3 = kk + j - 1;
d__1 = bli_zbla_real(ap[i__3]); d__1 = bli_zreal(ap[i__3]);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} }
kk += j; kk += j;
@@ -654,15 +654,15 @@
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
i__3 = jy; i__3 = jy;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0. || (bli_zbla_real(y[i__3]) != 0. || if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0. || (bli_zreal(y[i__3]) != 0. ||
bli_zimag(y[i__3]) != 0.)) { bli_zimag(y[i__3]) != 0.)) {
bla_d_cnjg(&z__2, &y[jy]); bla_d_cnjg(&z__2, &y[jy]);
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(z__2) - bli_zimag(*alpha) * bli_zimag(z__2)), (bli_zbla_real(*alpha) * bli_zimag(z__2) + bli_zimag(*alpha) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(z__2) - bli_zimag(*alpha) * bli_zimag(z__2)), (bli_zreal(*alpha) * bli_zimag(z__2) + bli_zimag(*alpha) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zbla_real(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zreal(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zreal(x[i__2])), z__2 );
bla_d_cnjg(&z__1, &z__2); bla_d_cnjg(&z__1, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
ix = kx; ix = kx;
iy = ky; iy = ky;
i__2 = kk + j - 2; i__2 = kk + j - 2;
@@ -670,12 +670,12 @@
i__3 = k; i__3 = k;
i__4 = k; i__4 = k;
i__5 = ix; i__5 = ix;
bli_zsets( (bli_zbla_real(x[i__5]) * bli_zbla_real(temp1) - bli_zimag(x[i__5]) * bli_zimag(temp1)), (bli_zbla_real(x[i__5]) * bli_zimag(temp1) + bli_zimag(x[i__5]) * bli_zbla_real(temp1)), z__3 ); bli_zsets( (bli_zreal(x[i__5]) * bli_zreal(temp1) - bli_zimag(x[i__5]) * bli_zimag(temp1)), (bli_zreal(x[i__5]) * bli_zimag(temp1) + bli_zimag(x[i__5]) * bli_zreal(temp1)), z__3 );
bli_zsets( (bli_zbla_real(ap[i__4]) + bli_zbla_real(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 ); bli_zsets( (bli_zreal(ap[i__4]) + bli_zreal(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 );
i__6 = iy; i__6 = iy;
bli_zsets( (bli_zbla_real(y[i__6]) * bli_zbla_real(temp2) - bli_zimag(y[i__6]) * bli_zimag(temp2)), (bli_zbla_real(y[i__6]) * bli_zimag(temp2) + bli_zimag(y[i__6]) * bli_zbla_real(temp2)), z__4 ); bli_zsets( (bli_zreal(y[i__6]) * bli_zreal(temp2) - bli_zimag(y[i__6]) * bli_zimag(temp2)), (bli_zreal(y[i__6]) * bli_zimag(temp2) + bli_zimag(y[i__6]) * bli_zreal(temp2)), z__4 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), ap[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] );
ix += *incx; ix += *incx;
iy += *incy; iy += *incy;
/* L30: */ /* L30: */
@@ -683,16 +683,16 @@
i__2 = kk + j - 1; i__2 = kk + j - 1;
i__3 = kk + j - 1; i__3 = kk + j - 1;
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(x[i__4]) * bli_zbla_real(temp1) - bli_zimag(x[i__4]) * bli_zimag(temp1)), (bli_zbla_real(x[i__4]) * bli_zimag(temp1) + bli_zimag(x[i__4]) * bli_zbla_real(temp1)), z__2 ); bli_zsets( (bli_zreal(x[i__4]) * bli_zreal(temp1) - bli_zimag(x[i__4]) * bli_zimag(temp1)), (bli_zreal(x[i__4]) * bli_zimag(temp1) + bli_zimag(x[i__4]) * bli_zreal(temp1)), z__2 );
i__5 = jy; i__5 = jy;
bli_zsets( (bli_zbla_real(y[i__5]) * bli_zbla_real(temp2) - bli_zimag(y[i__5]) * bli_zimag(temp2)), (bli_zbla_real(y[i__5]) * bli_zimag(temp2) + bli_zimag(y[i__5]) * bli_zbla_real(temp2)), z__3 ); bli_zsets( (bli_zreal(y[i__5]) * bli_zreal(temp2) - bli_zimag(y[i__5]) * bli_zimag(temp2)), (bli_zreal(y[i__5]) * bli_zimag(temp2) + bli_zimag(y[i__5]) * bli_zreal(temp2)), z__3 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 );
d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(z__1); d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} else { } else {
i__2 = kk + j - 1; i__2 = kk + j - 1;
i__3 = kk + j - 1; i__3 = kk + j - 1;
d__1 = bli_zbla_real(ap[i__3]); d__1 = bli_zreal(ap[i__3]);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} }
jx += *incx; jx += *incx;
@@ -710,23 +710,23 @@
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
i__3 = j; i__3 = j;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0. || (bli_zbla_real(y[i__3]) != 0. || if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0. || (bli_zreal(y[i__3]) != 0. ||
bli_zimag(y[i__3]) != 0.)) { bli_zimag(y[i__3]) != 0.)) {
bla_d_cnjg(&z__2, &y[j]); bla_d_cnjg(&z__2, &y[j]);
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(z__2) - bli_zimag(*alpha) * bli_zimag(z__2)), (bli_zbla_real(*alpha) * bli_zimag(z__2) + bli_zimag(*alpha) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(z__2) - bli_zimag(*alpha) * bli_zimag(z__2)), (bli_zreal(*alpha) * bli_zimag(z__2) + bli_zimag(*alpha) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zbla_real(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zreal(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zreal(x[i__2])), z__2 );
bla_d_cnjg(&z__1, &z__2); bla_d_cnjg(&z__1, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
i__2 = kk; i__2 = kk;
i__3 = kk; i__3 = kk;
i__4 = j; i__4 = j;
bli_zsets( (bli_zbla_real(x[i__4]) * bli_zbla_real(temp1) - bli_zimag(x[i__4]) * bli_zimag(temp1)), (bli_zbla_real(x[i__4]) * bli_zimag(temp1) + bli_zimag(x[i__4]) * bli_zbla_real(temp1)), z__2 ); bli_zsets( (bli_zreal(x[i__4]) * bli_zreal(temp1) - bli_zimag(x[i__4]) * bli_zimag(temp1)), (bli_zreal(x[i__4]) * bli_zimag(temp1) + bli_zimag(x[i__4]) * bli_zreal(temp1)), z__2 );
i__5 = j; i__5 = j;
bli_zsets( (bli_zbla_real(y[i__5]) * bli_zbla_real(temp2) - bli_zimag(y[i__5]) * bli_zimag(temp2)), (bli_zbla_real(y[i__5]) * bli_zimag(temp2) + bli_zimag(y[i__5]) * bli_zbla_real(temp2)), z__3 ); bli_zsets( (bli_zreal(y[i__5]) * bli_zreal(temp2) - bli_zimag(y[i__5]) * bli_zimag(temp2)), (bli_zreal(y[i__5]) * bli_zimag(temp2) + bli_zimag(y[i__5]) * bli_zreal(temp2)), z__3 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 );
d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(z__1); d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
k = kk + 1; k = kk + 1;
i__2 = *n; i__2 = *n;
@@ -734,19 +734,19 @@
i__3 = k; i__3 = k;
i__4 = k; i__4 = k;
i__5 = i__; i__5 = i__;
bli_zsets( (bli_zbla_real(x[i__5]) * bli_zbla_real(temp1) - bli_zimag(x[i__5]) * bli_zimag(temp1)), (bli_zbla_real(x[i__5]) * bli_zimag(temp1) + bli_zimag(x[i__5]) * bli_zbla_real(temp1)), z__3 ); bli_zsets( (bli_zreal(x[i__5]) * bli_zreal(temp1) - bli_zimag(x[i__5]) * bli_zimag(temp1)), (bli_zreal(x[i__5]) * bli_zimag(temp1) + bli_zimag(x[i__5]) * bli_zreal(temp1)), z__3 );
bli_zsets( (bli_zbla_real(ap[i__4]) + bli_zbla_real(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 ); bli_zsets( (bli_zreal(ap[i__4]) + bli_zreal(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 );
i__6 = i__; i__6 = i__;
bli_zsets( (bli_zbla_real(y[i__6]) * bli_zbla_real(temp2) - bli_zimag(y[i__6]) * bli_zimag(temp2)), (bli_zbla_real(y[i__6]) * bli_zimag(temp2) + bli_zimag(y[i__6]) * bli_zbla_real(temp2)), z__4 ); bli_zsets( (bli_zreal(y[i__6]) * bli_zreal(temp2) - bli_zimag(y[i__6]) * bli_zimag(temp2)), (bli_zreal(y[i__6]) * bli_zimag(temp2) + bli_zimag(y[i__6]) * bli_zreal(temp2)), z__4 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), ap[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] );
++k; ++k;
/* L50: */ /* L50: */
} }
} else { } else {
i__2 = kk; i__2 = kk;
i__3 = kk; i__3 = kk;
d__1 = bli_zbla_real(ap[i__3]); d__1 = bli_zreal(ap[i__3]);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} }
kk = kk + *n - j + 1; kk = kk + *n - j + 1;
@@ -757,23 +757,23 @@
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
i__3 = jy; i__3 = jy;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0. || (bli_zbla_real(y[i__3]) != 0. || if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0. || (bli_zreal(y[i__3]) != 0. ||
bli_zimag(y[i__3]) != 0.)) { bli_zimag(y[i__3]) != 0.)) {
bla_d_cnjg(&z__2, &y[jy]); bla_d_cnjg(&z__2, &y[jy]);
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(z__2) - bli_zimag(*alpha) * bli_zimag(z__2)), (bli_zbla_real(*alpha) * bli_zimag(z__2) + bli_zimag(*alpha) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(z__2) - bli_zimag(*alpha) * bli_zimag(z__2)), (bli_zreal(*alpha) * bli_zimag(z__2) + bli_zimag(*alpha) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp1 );
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(*alpha) * bli_zbla_real(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zbla_real(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(*alpha) * bli_zreal(x[i__2]) - bli_zimag(*alpha) * bli_zimag(x[i__2])), (bli_zreal(*alpha) * bli_zimag(x[i__2]) + bli_zimag(*alpha) * bli_zreal(x[i__2])), z__2 );
bla_d_cnjg(&z__1, &z__2); bla_d_cnjg(&z__1, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp2 ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 );
i__2 = kk; i__2 = kk;
i__3 = kk; i__3 = kk;
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(x[i__4]) * bli_zbla_real(temp1) - bli_zimag(x[i__4]) * bli_zimag(temp1)), (bli_zbla_real(x[i__4]) * bli_zimag(temp1) + bli_zimag(x[i__4]) * bli_zbla_real(temp1)), z__2 ); bli_zsets( (bli_zreal(x[i__4]) * bli_zreal(temp1) - bli_zimag(x[i__4]) * bli_zimag(temp1)), (bli_zreal(x[i__4]) * bli_zimag(temp1) + bli_zimag(x[i__4]) * bli_zreal(temp1)), z__2 );
i__5 = jy; i__5 = jy;
bli_zsets( (bli_zbla_real(y[i__5]) * bli_zbla_real(temp2) - bli_zimag(y[i__5]) * bli_zimag(temp2)), (bli_zbla_real(y[i__5]) * bli_zimag(temp2) + bli_zimag(y[i__5]) * bli_zbla_real(temp2)), z__3 ); bli_zsets( (bli_zreal(y[i__5]) * bli_zreal(temp2) - bli_zimag(y[i__5]) * bli_zimag(temp2)), (bli_zreal(y[i__5]) * bli_zimag(temp2) + bli_zimag(y[i__5]) * bli_zreal(temp2)), z__3 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 );
d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(z__1); d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
ix = jx; ix = jx;
iy = jy; iy = jy;
@@ -784,18 +784,18 @@
i__3 = k; i__3 = k;
i__4 = k; i__4 = k;
i__5 = ix; i__5 = ix;
bli_zsets( (bli_zbla_real(x[i__5]) * bli_zbla_real(temp1) - bli_zimag(x[i__5]) * bli_zimag(temp1)), (bli_zbla_real(x[i__5]) * bli_zimag(temp1) + bli_zimag(x[i__5]) * bli_zbla_real(temp1)), z__3 ); bli_zsets( (bli_zreal(x[i__5]) * bli_zreal(temp1) - bli_zimag(x[i__5]) * bli_zimag(temp1)), (bli_zreal(x[i__5]) * bli_zimag(temp1) + bli_zimag(x[i__5]) * bli_zreal(temp1)), z__3 );
bli_zsets( (bli_zbla_real(ap[i__4]) + bli_zbla_real(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 ); bli_zsets( (bli_zreal(ap[i__4]) + bli_zreal(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 );
i__6 = iy; i__6 = iy;
bli_zsets( (bli_zbla_real(y[i__6]) * bli_zbla_real(temp2) - bli_zimag(y[i__6]) * bli_zimag(temp2)), (bli_zbla_real(y[i__6]) * bli_zimag(temp2) + bli_zimag(y[i__6]) * bli_zbla_real(temp2)), z__4 ); bli_zsets( (bli_zreal(y[i__6]) * bli_zreal(temp2) - bli_zimag(y[i__6]) * bli_zimag(temp2)), (bli_zreal(y[i__6]) * bli_zimag(temp2) + bli_zimag(y[i__6]) * bli_zreal(temp2)), z__4 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), ap[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] );
/* L70: */ /* L70: */
} }
} else { } else {
i__2 = kk; i__2 = kk;
i__3 = kk; i__3 = kk;
d__1 = bli_zbla_real(ap[i__3]); d__1 = bli_zreal(ap[i__3]);
bli_zsets( (d__1), (0.), ap[i__2] ); bli_zsets( (d__1), (0.), ap[i__2] );
} }
jx += *incx; jx += *incx;

36
frame/compat/f2c/bla_rot.c
View File

@@ -313,20 +313,20 @@ L20:
i__1 = *n; i__1 = *n;
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = ix; i__2 = ix;
bli_zsets( (*c__ * bli_zbla_real(zx[i__2])), (*c__ * bli_zimag(zx[i__2])), z__2 ); bli_zsets( (*c__ * bli_zreal(zx[i__2])), (*c__ * bli_zimag(zx[i__2])), z__2 );
i__3 = iy; i__3 = iy;
bli_zsets( (*s * bli_zbla_real(zy[i__3])), (*s * bli_zimag(zy[i__3])), z__3 ); bli_zsets( (*s * bli_zreal(zy[i__3])), (*s * bli_zimag(zy[i__3])), z__3 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), ztemp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ztemp );
i__2 = iy; i__2 = iy;
i__3 = iy; i__3 = iy;
bli_zsets( (*c__ * bli_zbla_real(zy[i__3])), (*c__ * bli_zimag(zy[i__3])), z__2 ); bli_zsets( (*c__ * bli_zreal(zy[i__3])), (*c__ * bli_zimag(zy[i__3])), z__2 );
i__4 = ix; i__4 = ix;
bli_zsets( (*s * bli_zbla_real(zx[i__4])), (*s * bli_zimag(zx[i__4])), z__3 ); bli_zsets( (*s * bli_zreal(zx[i__4])), (*s * bli_zimag(zx[i__4])), z__3 );
bli_zsets( (bli_zbla_real(z__2) - bli_zbla_real(z__3)), (bli_zimag(z__2) - bli_zimag(z__3)), z__1 ); bli_zsets( (bli_zreal(z__2) - bli_zreal(z__3)), (bli_zimag(z__2) - bli_zimag(z__3)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), zy[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), zy[i__2] );
i__2 = ix; i__2 = ix;
bli_zsets( (bli_zbla_real(ztemp)), (bli_zimag(ztemp)), zx[i__2] ); bli_zsets( (bli_zreal(ztemp)), (bli_zimag(ztemp)), zx[i__2] );
ix += *incx; ix += *incx;
iy += *incy; iy += *incy;
/* L10: */ /* L10: */
@@ -339,20 +339,20 @@ L20:
i__1 = *n; i__1 = *n;
for (i__ = 1; i__ <= i__1; ++i__) { for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = i__; i__2 = i__;
bli_zsets( (*c__ * bli_zbla_real(zx[i__2])), (*c__ * bli_zimag(zx[i__2])), z__2 ); bli_zsets( (*c__ * bli_zreal(zx[i__2])), (*c__ * bli_zimag(zx[i__2])), z__2 );
i__3 = i__; i__3 = i__;
bli_zsets( (*s * bli_zbla_real(zy[i__3])), (*s * bli_zimag(zy[i__3])), z__3 ); bli_zsets( (*s * bli_zreal(zy[i__3])), (*s * bli_zimag(zy[i__3])), z__3 );
bli_zsets( (bli_zbla_real(z__2) + bli_zbla_real(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); bli_zsets( (bli_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), ztemp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ztemp );
i__2 = i__; i__2 = i__;
i__3 = i__; i__3 = i__;
bli_zsets( (*c__ * bli_zbla_real(zy[i__3])), (*c__ * bli_zimag(zy[i__3])), z__2 ); bli_zsets( (*c__ * bli_zreal(zy[i__3])), (*c__ * bli_zimag(zy[i__3])), z__2 );
i__4 = i__; i__4 = i__;
bli_zsets( (*s * bli_zbla_real(zx[i__4])), (*s * bli_zimag(zx[i__4])), z__3 ); bli_zsets( (*s * bli_zreal(zx[i__4])), (*s * bli_zimag(zx[i__4])), z__3 );
bli_zsets( (bli_zbla_real(z__2) - bli_zbla_real(z__3)), (bli_zimag(z__2) - bli_zimag(z__3)), z__1 ); bli_zsets( (bli_zreal(z__2) - bli_zreal(z__3)), (bli_zimag(z__2) - bli_zimag(z__3)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), zy[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), zy[i__2] );
i__2 = i__; i__2 = i__;
bli_zsets( (bli_zbla_real(ztemp)), (bli_zimag(ztemp)), zx[i__2] ); bli_zsets( (bli_zreal(ztemp)), (bli_zimag(ztemp)), zx[i__2] );
/* L30: */ /* L30: */
} }
return 0; return 0;

16
frame/compat/f2c/bla_rotg.c
View File

@@ -233,7 +233,7 @@ L20:
} }
*c__ = 0.; *c__ = 0.;
bli_zsets( 1., 0., *s ); bli_zsets( 1., 0., *s );
bli_zsets( bli_zbla_real(*cb), bli_zimag(*cb), *ca ); bli_zsets( bli_zreal(*cb), bli_zimag(*cb), *ca );
goto L20; goto L20;
L10: L10:
scale = bla_z_abs(ca) + bla_z_abs(cb); scale = bla_z_abs(ca) + bla_z_abs(cb);
@@ -247,15 +247,15 @@ L10:
d__2 = bla_z_abs(&z__3); d__2 = bla_z_abs(&z__3);
norm = scale * sqrt(d__1 * d__1 + d__2 * d__2); norm = scale * sqrt(d__1 * d__1 + d__2 * d__2);
d__1 = bla_z_abs(ca); d__1 = bla_z_abs(ca);
bli_zsets( (bli_zbla_real(*ca) / d__1), (bli_zimag(*ca) / d__1), z__1 ); bli_zsets( (bli_zreal(*ca) / d__1), (bli_zimag(*ca) / d__1), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), alpha ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), alpha );
*c__ = bla_z_abs(ca) / norm; *c__ = bla_z_abs(ca) / norm;
bla_d_cnjg(&z__3, cb); bla_d_cnjg(&z__3, cb);
bli_zsets( (bli_zbla_real(alpha) * bli_zbla_real(z__3) - bli_zimag(alpha) * bli_zimag(z__3)), (bli_zbla_real(alpha) * bli_zimag(z__3) + bli_zimag(alpha) * bli_zbla_real(z__3)), z__2 ); bli_zsets( (bli_zreal(alpha) * bli_zreal(z__3) - bli_zimag(alpha) * bli_zimag(z__3)), (bli_zreal(alpha) * bli_zimag(z__3) + bli_zimag(alpha) * bli_zreal(z__3)), z__2 );
bli_zsets( (bli_zbla_real(z__2) / norm), (bli_zimag(z__2) / norm), z__1 ); bli_zsets( (bli_zreal(z__2) / norm), (bli_zimag(z__2) / norm), z__1 );
bli_zsets( bli_zbla_real(z__1), bli_zimag(z__1), *s ); bli_zsets( bli_zreal(z__1), bli_zimag(z__1), *s );
bli_zsets( (norm * bli_zbla_real(alpha)), (norm * bli_zimag(alpha)), z__1 ); bli_zsets( (norm * bli_zreal(alpha)), (norm * bli_zimag(alpha)), z__1 );
bli_zsets( bli_zbla_real(z__1), bli_zimag(z__1), *ca ); bli_zsets( bli_zreal(z__1), bli_zimag(z__1), *ca );
L20: L20:
return 0; return 0;
} /* zrotg_ */ } /* zrotg_ */

152
frame/compat/f2c/bla_tbmv.c
View File

@@ -1656,9 +1656,9 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
l = kplus1 - j; l = kplus1 - j;
/* Computing MAX */ /* Computing MAX */
i__2 = 1, i__3 = j - *k; i__2 = 1, i__3 = j - *k;
@@ -1667,17 +1667,17 @@
i__2 = i__; i__2 = i__;
i__3 = i__; i__3 = i__;
i__5 = l + i__ + j * a_dim1; i__5 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zbla_real(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zbla_real(a[i__5])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zreal(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zreal(a[i__5])), z__2 );
bli_zsets( (bli_zbla_real(x[i__3]) + bli_zbla_real(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__3]) + bli_zreal(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
/* L10: */ /* L10: */
} }
if (nounit) { if (nounit) {
i__4 = j; i__4 = j;
i__2 = j; i__2 = j;
i__3 = kplus1 + j * a_dim1; i__3 = kplus1 + j * a_dim1;
bli_zsets( (bli_zbla_real(x[i__2]) * bli_zbla_real(a[i__3]) - bli_zimag(x[i__2]) * bli_zimag(a[i__3])), (bli_zbla_real(x[i__2]) * bli_zimag(a[i__3]) + bli_zimag(x[i__2]) * bli_zbla_real(a[i__3])), z__1 ); bli_zsets( (bli_zreal(x[i__2]) * bli_zreal(a[i__3]) - bli_zimag(x[i__2]) * bli_zimag(a[i__3])), (bli_zreal(x[i__2]) * bli_zimag(a[i__3]) + bli_zimag(x[i__2]) * bli_zreal(a[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__4] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__4] );
} }
} }
/* L20: */ /* L20: */
@@ -1687,9 +1687,9 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__4 = jx; i__4 = jx;
if (bli_zbla_real(x[i__4]) != 0. || bli_zimag(x[i__4]) != 0.) { if (bli_zreal(x[i__4]) != 0. || bli_zimag(x[i__4]) != 0.) {
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(x[i__4])), (bli_zimag(x[i__4])), temp ); bli_zsets( (bli_zreal(x[i__4])), (bli_zimag(x[i__4])), temp );
ix = kx; ix = kx;
l = kplus1 - j; l = kplus1 - j;
/* Computing MAX */ /* Computing MAX */
@@ -1699,9 +1699,9 @@
i__4 = ix; i__4 = ix;
i__2 = ix; i__2 = ix;
i__5 = l + i__ + j * a_dim1; i__5 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zbla_real(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zbla_real(a[i__5])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zreal(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zreal(a[i__5])), z__2 );
bli_zsets( (bli_zbla_real(x[i__2]) + bli_zbla_real(z__2)), (bli_zimag(x[i__2]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__2]) + bli_zreal(z__2)), (bli_zimag(x[i__2]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__4] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__4] );
ix += *incx; ix += *incx;
/* L30: */ /* L30: */
} }
@@ -1709,8 +1709,8 @@
i__3 = jx; i__3 = jx;
i__4 = jx; i__4 = jx;
i__2 = kplus1 + j * a_dim1; i__2 = kplus1 + j * a_dim1;
bli_zsets( (bli_zbla_real(x[i__4]) * bli_zbla_real(a[i__2]) - bli_zimag(x[i__4]) * bli_zimag(a[i__2])), (bli_zbla_real(x[i__4]) * bli_zimag(a[i__2]) + bli_zimag(x[i__4]) * bli_zbla_real(a[i__2])), z__1 ); bli_zsets( (bli_zreal(x[i__4]) * bli_zreal(a[i__2]) - bli_zimag(x[i__4]) * bli_zimag(a[i__2])), (bli_zreal(x[i__4]) * bli_zimag(a[i__2]) + bli_zimag(x[i__4]) * bli_zreal(a[i__2])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] );
} }
} }
jx += *incx; jx += *incx;
@@ -1724,9 +1724,9 @@
if (*incx == 1) { if (*incx == 1) {
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = j; i__1 = j;
if (bli_zbla_real(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) { if (bli_zreal(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) {
i__1 = j; i__1 = j;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
l = 1 - j; l = 1 - j;
/* Computing MIN */ /* Computing MIN */
i__1 = *n, i__3 = j + *k; i__1 = *n, i__3 = j + *k;
@@ -1735,17 +1735,17 @@
i__1 = i__; i__1 = i__;
i__3 = i__; i__3 = i__;
i__2 = l + i__ + j * a_dim1; i__2 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__2]) - bli_zimag(temp) * bli_zimag(a[i__2])), (bli_zbla_real(temp) * bli_zimag(a[i__2]) + bli_zimag(temp) * bli_zbla_real(a[i__2])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__2]) - bli_zimag(temp) * bli_zimag(a[i__2])), (bli_zreal(temp) * bli_zimag(a[i__2]) + bli_zimag(temp) * bli_zreal(a[i__2])), z__2 );
bli_zsets( (bli_zbla_real(x[i__3]) + bli_zbla_real(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__3]) + bli_zreal(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__1] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__1] );
/* L50: */ /* L50: */
} }
if (nounit) { if (nounit) {
i__4 = j; i__4 = j;
i__1 = j; i__1 = j;
i__3 = j * a_dim1 + 1; i__3 = j * a_dim1 + 1;
bli_zsets( (bli_zbla_real(x[i__1]) * bli_zbla_real(a[i__3]) - bli_zimag(x[i__1]) * bli_zimag(a[i__3])), (bli_zbla_real(x[i__1]) * bli_zimag(a[i__3]) + bli_zimag(x[i__1]) * bli_zbla_real(a[i__3])), z__1 ); bli_zsets( (bli_zreal(x[i__1]) * bli_zreal(a[i__3]) - bli_zimag(x[i__1]) * bli_zimag(a[i__3])), (bli_zreal(x[i__1]) * bli_zimag(a[i__3]) + bli_zimag(x[i__1]) * bli_zreal(a[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__4] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__4] );
} }
} }
/* L60: */ /* L60: */
@@ -1755,9 +1755,9 @@
jx = kx; jx = kx;
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__4 = jx; i__4 = jx;
if (bli_zbla_real(x[i__4]) != 0. || bli_zimag(x[i__4]) != 0.) { if (bli_zreal(x[i__4]) != 0. || bli_zimag(x[i__4]) != 0.) {
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(x[i__4])), (bli_zimag(x[i__4])), temp ); bli_zsets( (bli_zreal(x[i__4])), (bli_zimag(x[i__4])), temp );
ix = kx; ix = kx;
l = 1 - j; l = 1 - j;
/* Computing MIN */ /* Computing MIN */
@@ -1767,9 +1767,9 @@
i__4 = ix; i__4 = ix;
i__1 = ix; i__1 = ix;
i__2 = l + i__ + j * a_dim1; i__2 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__2]) - bli_zimag(temp) * bli_zimag(a[i__2])), (bli_zbla_real(temp) * bli_zimag(a[i__2]) + bli_zimag(temp) * bli_zbla_real(a[i__2])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__2]) - bli_zimag(temp) * bli_zimag(a[i__2])), (bli_zreal(temp) * bli_zimag(a[i__2]) + bli_zimag(temp) * bli_zreal(a[i__2])), z__2 );
bli_zsets( (bli_zbla_real(x[i__1]) + bli_zbla_real(z__2)), (bli_zimag(x[i__1]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__1]) + bli_zreal(z__2)), (bli_zimag(x[i__1]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__4] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__4] );
ix -= *incx; ix -= *incx;
/* L70: */ /* L70: */
} }
@@ -1777,8 +1777,8 @@
i__3 = jx; i__3 = jx;
i__4 = jx; i__4 = jx;
i__1 = j * a_dim1 + 1; i__1 = j * a_dim1 + 1;
bli_zsets( (bli_zbla_real(x[i__4]) * bli_zbla_real(a[i__1]) - bli_zimag(x[i__4]) * bli_zimag(a[i__1])), (bli_zbla_real(x[i__4]) * bli_zimag(a[i__1]) + bli_zimag(x[i__4]) * bli_zbla_real(a[i__1])), z__1 ); bli_zsets( (bli_zreal(x[i__4]) * bli_zreal(a[i__1]) - bli_zimag(x[i__4]) * bli_zimag(a[i__1])), (bli_zreal(x[i__4]) * bli_zimag(a[i__1]) + bli_zimag(x[i__4]) * bli_zreal(a[i__1])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] );
} }
} }
jx -= *incx; jx -= *incx;
@@ -1798,13 +1798,13 @@
if (*incx == 1) { if (*incx == 1) {
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__3 = j; i__3 = j;
bli_zsets( (bli_zbla_real(x[i__3])), (bli_zimag(x[i__3])), temp ); bli_zsets( (bli_zreal(x[i__3])), (bli_zimag(x[i__3])), temp );
l = kplus1 - j; l = kplus1 - j;
if (noconj) { if (noconj) {
if (nounit) { if (nounit) {
i__3 = kplus1 + j * a_dim1; i__3 = kplus1 + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__3]) - bli_zimag(temp) * bli_zimag(a[i__3])), (bli_zbla_real(temp) * bli_zimag(a[i__3]) + bli_zimag(temp) * bli_zbla_real(a[i__3])), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__3]) - bli_zimag(temp) * bli_zimag(a[i__3])), (bli_zreal(temp) * bli_zimag(a[i__3]) + bli_zimag(temp) * bli_zreal(a[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
/* Computing MAX */ /* Computing MAX */
i__4 = 1, i__1 = j - *k; i__4 = 1, i__1 = j - *k;
@@ -1812,16 +1812,16 @@
for (i__ = j - 1; i__ >= i__3; --i__) { for (i__ = j - 1; i__ >= i__3; --i__) {
i__4 = l + i__ + j * a_dim1; i__4 = l + i__ + j * a_dim1;
i__1 = i__; i__1 = i__;
bli_zsets( (bli_zbla_real(a[i__4]) * bli_zbla_real(x[i__1]) - bli_zimag(a[i__4]) * bli_zimag(x[i__1])), (bli_zbla_real(a[i__4]) * bli_zimag(x[i__1]) + bli_zimag(a[i__4]) * bli_zbla_real(x[i__1])), z__2 ); bli_zsets( (bli_zreal(a[i__4]) * bli_zreal(x[i__1]) - bli_zimag(a[i__4]) * bli_zimag(x[i__1])), (bli_zreal(a[i__4]) * bli_zimag(x[i__1]) + bli_zimag(a[i__4]) * bli_zreal(x[i__1])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L90: */ /* L90: */
} }
} else { } else {
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &a[kplus1 + j * a_dim1]); bla_d_cnjg(&z__2, &a[kplus1 + j * a_dim1]);
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zbla_real(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zreal(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
/* Computing MAX */ /* Computing MAX */
i__4 = 1, i__1 = j - *k; i__4 = 1, i__1 = j - *k;
@@ -1829,14 +1829,14 @@
for (i__ = j - 1; i__ >= i__3; --i__) { for (i__ = j - 1; i__ >= i__3; --i__) {
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__4 = i__; i__4 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zbla_real(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zreal(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L100: */ /* L100: */
} }
} }
i__3 = j; i__3 = j;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__3] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__3] );
/* L110: */ /* L110: */
} }
} else { } else {
@@ -1844,15 +1844,15 @@
jx = kx; jx = kx;
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__3 = jx; i__3 = jx;
bli_zsets( (bli_zbla_real(x[i__3])), (bli_zimag(x[i__3])), temp ); bli_zsets( (bli_zreal(x[i__3])), (bli_zimag(x[i__3])), temp );
kx -= *incx; kx -= *incx;
ix = kx; ix = kx;
l = kplus1 - j; l = kplus1 - j;
if (noconj) { if (noconj) {
if (nounit) { if (nounit) {
i__3 = kplus1 + j * a_dim1; i__3 = kplus1 + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__3]) - bli_zimag(temp) * bli_zimag(a[i__3])), (bli_zbla_real(temp) * bli_zimag(a[i__3]) + bli_zimag(temp) * bli_zbla_real(a[i__3])), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__3]) - bli_zimag(temp) * bli_zimag(a[i__3])), (bli_zreal(temp) * bli_zimag(a[i__3]) + bli_zimag(temp) * bli_zreal(a[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
/* Computing MAX */ /* Computing MAX */
i__4 = 1, i__1 = j - *k; i__4 = 1, i__1 = j - *k;
@@ -1860,17 +1860,17 @@
for (i__ = j - 1; i__ >= i__3; --i__) { for (i__ = j - 1; i__ >= i__3; --i__) {
i__4 = l + i__ + j * a_dim1; i__4 = l + i__ + j * a_dim1;
i__1 = ix; i__1 = ix;
bli_zsets( (bli_zbla_real(a[i__4]) * bli_zbla_real(x[i__1]) - bli_zimag(a[i__4]) * bli_zimag(x[i__1])), (bli_zbla_real(a[i__4]) * bli_zimag(x[i__1]) + bli_zimag(a[i__4]) * bli_zbla_real(x[i__1])), z__2 ); bli_zsets( (bli_zreal(a[i__4]) * bli_zreal(x[i__1]) - bli_zimag(a[i__4]) * bli_zimag(x[i__1])), (bli_zreal(a[i__4]) * bli_zimag(x[i__1]) + bli_zimag(a[i__4]) * bli_zreal(x[i__1])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix -= *incx; ix -= *incx;
/* L120: */ /* L120: */
} }
} else { } else {
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &a[kplus1 + j * a_dim1]); bla_d_cnjg(&z__2, &a[kplus1 + j * a_dim1]);
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zbla_real(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zreal(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
/* Computing MAX */ /* Computing MAX */
i__4 = 1, i__1 = j - *k; i__4 = 1, i__1 = j - *k;
@@ -1878,15 +1878,15 @@
for (i__ = j - 1; i__ >= i__3; --i__) { for (i__ = j - 1; i__ >= i__3; --i__) {
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__4 = ix; i__4 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zbla_real(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zreal(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix -= *incx; ix -= *incx;
/* L130: */ /* L130: */
} }
} }
i__3 = jx; i__3 = jx;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__3] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__3] );
jx -= *incx; jx -= *incx;
/* L140: */ /* L140: */
} }
@@ -1896,13 +1896,13 @@
i__3 = *n; i__3 = *n;
for (j = 1; j <= i__3; ++j) { for (j = 1; j <= i__3; ++j) {
i__4 = j; i__4 = j;
bli_zsets( (bli_zbla_real(x[i__4])), (bli_zimag(x[i__4])), temp ); bli_zsets( (bli_zreal(x[i__4])), (bli_zimag(x[i__4])), temp );
l = 1 - j; l = 1 - j;
if (noconj) { if (noconj) {
if (nounit) { if (nounit) {
i__4 = j * a_dim1 + 1; i__4 = j * a_dim1 + 1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__4]) - bli_zimag(temp) * bli_zimag(a[i__4])), (bli_zbla_real(temp) * bli_zimag(a[i__4]) + bli_zimag(temp) * bli_zbla_real(a[i__4])), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__4]) - bli_zimag(temp) * bli_zimag(a[i__4])), (bli_zreal(temp) * bli_zimag(a[i__4]) + bli_zimag(temp) * bli_zreal(a[i__4])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
/* Computing MIN */ /* Computing MIN */
i__1 = *n, i__2 = j + *k; i__1 = *n, i__2 = j + *k;
@@ -1910,16 +1910,16 @@
for (i__ = j + 1; i__ <= i__4; ++i__) { for (i__ = j + 1; i__ <= i__4; ++i__) {
i__1 = l + i__ + j * a_dim1; i__1 = l + i__ + j * a_dim1;
i__2 = i__; i__2 = i__;
bli_zsets( (bli_zbla_real(a[i__1]) * bli_zbla_real(x[i__2]) - bli_zimag(a[i__1]) * bli_zimag(x[i__2])), (bli_zbla_real(a[i__1]) * bli_zimag(x[i__2]) + bli_zimag(a[i__1]) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(a[i__1]) * bli_zreal(x[i__2]) - bli_zimag(a[i__1]) * bli_zimag(x[i__2])), (bli_zreal(a[i__1]) * bli_zimag(x[i__2]) + bli_zimag(a[i__1]) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L150: */ /* L150: */
} }
} else { } else {
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &a[j * a_dim1 + 1]); bla_d_cnjg(&z__2, &a[j * a_dim1 + 1]);
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zbla_real(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zreal(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
/* Computing MIN */ /* Computing MIN */
i__1 = *n, i__2 = j + *k; i__1 = *n, i__2 = j + *k;
@@ -1927,14 +1927,14 @@
for (i__ = j + 1; i__ <= i__4; ++i__) { for (i__ = j + 1; i__ <= i__4; ++i__) {
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__1 = i__; i__1 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__1]) - bli_zimag(z__3) * bli_zimag(x[i__1])), (bli_zbla_real(z__3) * bli_zimag(x[i__1]) + bli_zimag(z__3) * bli_zbla_real(x[i__1])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__1]) - bli_zimag(z__3) * bli_zimag(x[i__1])), (bli_zreal(z__3) * bli_zimag(x[i__1]) + bli_zimag(z__3) * bli_zreal(x[i__1])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L160: */ /* L160: */
} }
} }
i__4 = j; i__4 = j;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__4] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__4] );
/* L170: */ /* L170: */
} }
} else { } else {
@@ -1942,15 +1942,15 @@
i__3 = *n; i__3 = *n;
for (j = 1; j <= i__3; ++j) { for (j = 1; j <= i__3; ++j) {
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(x[i__4])), (bli_zimag(x[i__4])), temp ); bli_zsets( (bli_zreal(x[i__4])), (bli_zimag(x[i__4])), temp );
kx += *incx; kx += *incx;
ix = kx; ix = kx;
l = 1 - j; l = 1 - j;
if (noconj) { if (noconj) {
if (nounit) { if (nounit) {
i__4 = j * a_dim1 + 1; i__4 = j * a_dim1 + 1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__4]) - bli_zimag(temp) * bli_zimag(a[i__4])), (bli_zbla_real(temp) * bli_zimag(a[i__4]) + bli_zimag(temp) * bli_zbla_real(a[i__4])), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__4]) - bli_zimag(temp) * bli_zimag(a[i__4])), (bli_zreal(temp) * bli_zimag(a[i__4]) + bli_zimag(temp) * bli_zreal(a[i__4])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
/* Computing MIN */ /* Computing MIN */
i__1 = *n, i__2 = j + *k; i__1 = *n, i__2 = j + *k;
@@ -1958,17 +1958,17 @@
for (i__ = j + 1; i__ <= i__4; ++i__) { for (i__ = j + 1; i__ <= i__4; ++i__) {
i__1 = l + i__ + j * a_dim1; i__1 = l + i__ + j * a_dim1;
i__2 = ix; i__2 = ix;
bli_zsets( (bli_zbla_real(a[i__1]) * bli_zbla_real(x[i__2]) - bli_zimag(a[i__1]) * bli_zimag(x[i__2])), (bli_zbla_real(a[i__1]) * bli_zimag(x[i__2]) + bli_zimag(a[i__1]) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(a[i__1]) * bli_zreal(x[i__2]) - bli_zimag(a[i__1]) * bli_zimag(x[i__2])), (bli_zreal(a[i__1]) * bli_zimag(x[i__2]) + bli_zimag(a[i__1]) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix += *incx; ix += *incx;
/* L180: */ /* L180: */
} }
} else { } else {
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &a[j * a_dim1 + 1]); bla_d_cnjg(&z__2, &a[j * a_dim1 + 1]);
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zbla_real(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zreal(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
/* Computing MIN */ /* Computing MIN */
i__1 = *n, i__2 = j + *k; i__1 = *n, i__2 = j + *k;
@@ -1976,15 +1976,15 @@
for (i__ = j + 1; i__ <= i__4; ++i__) { for (i__ = j + 1; i__ <= i__4; ++i__) {
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__1 = ix; i__1 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__1]) - bli_zimag(z__3) * bli_zimag(x[i__1])), (bli_zbla_real(z__3) * bli_zimag(x[i__1]) + bli_zimag(z__3) * bli_zbla_real(x[i__1])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__1]) - bli_zimag(z__3) * bli_zimag(x[i__1])), (bli_zreal(z__3) * bli_zimag(x[i__1]) + bli_zimag(z__3) * bli_zreal(x[i__1])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix += *incx; ix += *incx;
/* L190: */ /* L190: */
} }
} }
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__4] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__4] );
jx += *incx; jx += *incx;
/* L200: */ /* L200: */
} }

128
frame/compat/f2c/bla_tbsv.c
View File

@@ -1660,15 +1660,15 @@
if (*incx == 1) { if (*incx == 1) {
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = j; i__1 = j;
if (bli_zbla_real(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) { if (bli_zreal(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) {
l = kplus1 - j; l = kplus1 - j;
if (nounit) { if (nounit) {
i__1 = j; i__1 = j;
bla_z_div(&z__1, &x[j], &a[kplus1 + j * a_dim1]); bla_z_div(&z__1, &x[j], &a[kplus1 + j * a_dim1]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__1] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__1] );
} }
i__1 = j; i__1 = j;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
/* Computing MAX */ /* Computing MAX */
i__2 = 1, i__3 = j - *k; i__2 = 1, i__3 = j - *k;
i__1 = f2c_max(i__2,i__3); i__1 = f2c_max(i__2,i__3);
@@ -1676,9 +1676,9 @@
i__2 = i__; i__2 = i__;
i__3 = i__; i__3 = i__;
i__4 = l + i__ + j * a_dim1; i__4 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__4]) - bli_zimag(temp) * bli_zimag(a[i__4])), (bli_zbla_real(temp) * bli_zimag(a[i__4]) + bli_zimag(temp) * bli_zbla_real(a[i__4])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__4]) - bli_zimag(temp) * bli_zimag(a[i__4])), (bli_zreal(temp) * bli_zimag(a[i__4]) + bli_zimag(temp) * bli_zreal(a[i__4])), z__2 );
bli_zsets( (bli_zbla_real(x[i__3]) - bli_zbla_real(z__2)), (bli_zimag(x[i__3]) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__3]) - bli_zreal(z__2)), (bli_zimag(x[i__3]) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
/* L10: */ /* L10: */
} }
} }
@@ -1690,16 +1690,16 @@
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
kx -= *incx; kx -= *incx;
i__1 = jx; i__1 = jx;
if (bli_zbla_real(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) { if (bli_zreal(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) {
ix = kx; ix = kx;
l = kplus1 - j; l = kplus1 - j;
if (nounit) { if (nounit) {
i__1 = jx; i__1 = jx;
bla_z_div(&z__1, &x[jx], &a[kplus1 + j * a_dim1]); bla_z_div(&z__1, &x[jx], &a[kplus1 + j * a_dim1]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__1] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__1] );
} }
i__1 = jx; i__1 = jx;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
/* Computing MAX */ /* Computing MAX */
i__2 = 1, i__3 = j - *k; i__2 = 1, i__3 = j - *k;
i__1 = f2c_max(i__2,i__3); i__1 = f2c_max(i__2,i__3);
@@ -1707,9 +1707,9 @@
i__2 = ix; i__2 = ix;
i__3 = ix; i__3 = ix;
i__4 = l + i__ + j * a_dim1; i__4 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__4]) - bli_zimag(temp) * bli_zimag(a[i__4])), (bli_zbla_real(temp) * bli_zimag(a[i__4]) + bli_zimag(temp) * bli_zbla_real(a[i__4])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__4]) - bli_zimag(temp) * bli_zimag(a[i__4])), (bli_zreal(temp) * bli_zimag(a[i__4]) + bli_zimag(temp) * bli_zreal(a[i__4])), z__2 );
bli_zsets( (bli_zbla_real(x[i__3]) - bli_zbla_real(z__2)), (bli_zimag(x[i__3]) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__3]) - bli_zreal(z__2)), (bli_zimag(x[i__3]) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
ix -= *incx; ix -= *incx;
/* L30: */ /* L30: */
} }
@@ -1723,15 +1723,15 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
l = 1 - j; l = 1 - j;
if (nounit) { if (nounit) {
i__2 = j; i__2 = j;
bla_z_div(&z__1, &x[j], &a[j * a_dim1 + 1]); bla_z_div(&z__1, &x[j], &a[j * a_dim1 + 1]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
} }
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
/* Computing MIN */ /* Computing MIN */
i__3 = *n, i__4 = j + *k; i__3 = *n, i__4 = j + *k;
i__2 = f2c_min(i__3,i__4); i__2 = f2c_min(i__3,i__4);
@@ -1739,9 +1739,9 @@
i__3 = i__; i__3 = i__;
i__4 = i__; i__4 = i__;
i__5 = l + i__ + j * a_dim1; i__5 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zbla_real(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zbla_real(a[i__5])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zreal(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zreal(a[i__5])), z__2 );
bli_zsets( (bli_zbla_real(x[i__4]) - bli_zbla_real(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__4]) - bli_zreal(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] );
/* L50: */ /* L50: */
} }
} }
@@ -1753,16 +1753,16 @@
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
kx += *incx; kx += *incx;
i__2 = jx; i__2 = jx;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
ix = kx; ix = kx;
l = 1 - j; l = 1 - j;
if (nounit) { if (nounit) {
i__2 = jx; i__2 = jx;
bla_z_div(&z__1, &x[jx], &a[j * a_dim1 + 1]); bla_z_div(&z__1, &x[jx], &a[j * a_dim1 + 1]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
} }
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
/* Computing MIN */ /* Computing MIN */
i__3 = *n, i__4 = j + *k; i__3 = *n, i__4 = j + *k;
i__2 = f2c_min(i__3,i__4); i__2 = f2c_min(i__3,i__4);
@@ -1770,9 +1770,9 @@
i__3 = ix; i__3 = ix;
i__4 = ix; i__4 = ix;
i__5 = l + i__ + j * a_dim1; i__5 = l + i__ + j * a_dim1;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zbla_real(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zbla_real(a[i__5])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(a[i__5]) - bli_zimag(temp) * bli_zimag(a[i__5])), (bli_zreal(temp) * bli_zimag(a[i__5]) + bli_zimag(temp) * bli_zreal(a[i__5])), z__2 );
bli_zsets( (bli_zbla_real(x[i__4]) - bli_zbla_real(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__4]) - bli_zreal(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] );
ix += *incx; ix += *incx;
/* L70: */ /* L70: */
} }
@@ -1792,7 +1792,7 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
l = kplus1 - j; l = kplus1 - j;
if (noconj) { if (noconj) {
/* Computing MAX */ /* Computing MAX */
@@ -1801,14 +1801,14 @@
for (i__ = f2c_max(i__2,i__3); i__ <= i__4; ++i__) { for (i__ = f2c_max(i__2,i__3); i__ <= i__4; ++i__) {
i__2 = l + i__ + j * a_dim1; i__2 = l + i__ + j * a_dim1;
i__3 = i__; i__3 = i__;
bli_zsets( (bli_zbla_real(a[i__2]) * bli_zbla_real(x[i__3]) - bli_zimag(a[i__2]) * bli_zimag(x[i__3])), (bli_zbla_real(a[i__2]) * bli_zimag(x[i__3]) + bli_zimag(a[i__2]) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(a[i__2]) * bli_zreal(x[i__3]) - bli_zimag(a[i__2]) * bli_zimag(x[i__3])), (bli_zreal(a[i__2]) * bli_zimag(x[i__3]) + bli_zimag(a[i__2]) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L90: */ /* L90: */
} }
if (nounit) { if (nounit) {
bla_z_div(&z__1, &temp, &a[kplus1 + j * a_dim1]); bla_z_div(&z__1, &temp, &a[kplus1 + j * a_dim1]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} else { } else {
/* Computing MAX */ /* Computing MAX */
@@ -1817,19 +1817,19 @@
for (i__ = f2c_max(i__4,i__2); i__ <= i__3; ++i__) { for (i__ = f2c_max(i__4,i__2); i__ <= i__3; ++i__) {
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__4 = i__; i__4 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zbla_real(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__4]) - bli_zimag(z__3) * bli_zimag(x[i__4])), (bli_zreal(z__3) * bli_zimag(x[i__4]) + bli_zimag(z__3) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L100: */ /* L100: */
} }
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &a[kplus1 + j * a_dim1]); bla_d_cnjg(&z__2, &a[kplus1 + j * a_dim1]);
bla_z_div(&z__1, &temp, &z__2); bla_z_div(&z__1, &temp, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} }
i__3 = j; i__3 = j;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__3] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__3] );
/* L110: */ /* L110: */
} }
} else { } else {
@@ -1837,7 +1837,7 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__3 = jx; i__3 = jx;
bli_zsets( (bli_zbla_real(x[i__3])), (bli_zimag(x[i__3])), temp ); bli_zsets( (bli_zreal(x[i__3])), (bli_zimag(x[i__3])), temp );
ix = kx; ix = kx;
l = kplus1 - j; l = kplus1 - j;
if (noconj) { if (noconj) {
@@ -1847,15 +1847,15 @@
for (i__ = f2c_max(i__3,i__4); i__ <= i__2; ++i__) { for (i__ = f2c_max(i__3,i__4); i__ <= i__2; ++i__) {
i__3 = l + i__ + j * a_dim1; i__3 = l + i__ + j * a_dim1;
i__4 = ix; i__4 = ix;
bli_zsets( (bli_zbla_real(a[i__3]) * bli_zbla_real(x[i__4]) - bli_zimag(a[i__3]) * bli_zimag(x[i__4])), (bli_zbla_real(a[i__3]) * bli_zimag(x[i__4]) + bli_zimag(a[i__3]) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(a[i__3]) * bli_zreal(x[i__4]) - bli_zimag(a[i__3]) * bli_zimag(x[i__4])), (bli_zreal(a[i__3]) * bli_zimag(x[i__4]) + bli_zimag(a[i__3]) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix += *incx; ix += *incx;
/* L120: */ /* L120: */
} }
if (nounit) { if (nounit) {
bla_z_div(&z__1, &temp, &a[kplus1 + j * a_dim1]); bla_z_div(&z__1, &temp, &a[kplus1 + j * a_dim1]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} else { } else {
/* Computing MAX */ /* Computing MAX */
@@ -1864,20 +1864,20 @@
for (i__ = f2c_max(i__2,i__3); i__ <= i__4; ++i__) { for (i__ = f2c_max(i__2,i__3); i__ <= i__4; ++i__) {
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__2 = ix; i__2 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zbla_real(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zreal(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix += *incx; ix += *incx;
/* L130: */ /* L130: */
} }
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &a[kplus1 + j * a_dim1]); bla_d_cnjg(&z__2, &a[kplus1 + j * a_dim1]);
bla_z_div(&z__1, &temp, &z__2); bla_z_div(&z__1, &temp, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} }
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__4] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__4] );
jx += *incx; jx += *incx;
if (j > *k) { if (j > *k) {
kx += *incx; kx += *incx;
@@ -1889,7 +1889,7 @@
if (*incx == 1) { if (*incx == 1) {
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = j; i__1 = j;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
l = 1 - j; l = 1 - j;
if (noconj) { if (noconj) {
/* Computing MIN */ /* Computing MIN */
@@ -1898,14 +1898,14 @@
for (i__ = f2c_min(i__1,i__4); i__ >= i__2; --i__) { for (i__ = f2c_min(i__1,i__4); i__ >= i__2; --i__) {
i__1 = l + i__ + j * a_dim1; i__1 = l + i__ + j * a_dim1;
i__4 = i__; i__4 = i__;
bli_zsets( (bli_zbla_real(a[i__1]) * bli_zbla_real(x[i__4]) - bli_zimag(a[i__1]) * bli_zimag(x[i__4])), (bli_zbla_real(a[i__1]) * bli_zimag(x[i__4]) + bli_zimag(a[i__1]) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(a[i__1]) * bli_zreal(x[i__4]) - bli_zimag(a[i__1]) * bli_zimag(x[i__4])), (bli_zreal(a[i__1]) * bli_zimag(x[i__4]) + bli_zimag(a[i__1]) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L150: */ /* L150: */
} }
if (nounit) { if (nounit) {
bla_z_div(&z__1, &temp, &a[j * a_dim1 + 1]); bla_z_div(&z__1, &temp, &a[j * a_dim1 + 1]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} else { } else {
/* Computing MIN */ /* Computing MIN */
@@ -1914,19 +1914,19 @@
for (i__ = f2c_min(i__2,i__1); i__ >= i__4; --i__) { for (i__ = f2c_min(i__2,i__1); i__ >= i__4; --i__) {
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__2 = i__; i__2 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zbla_real(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zreal(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L160: */ /* L160: */
} }
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &a[j * a_dim1 + 1]); bla_d_cnjg(&z__2, &a[j * a_dim1 + 1]);
bla_z_div(&z__1, &temp, &z__2); bla_z_div(&z__1, &temp, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} }
i__4 = j; i__4 = j;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__4] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__4] );
/* L170: */ /* L170: */
} }
} else { } else {
@@ -1934,7 +1934,7 @@
jx = kx; jx = kx;
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__4 = jx; i__4 = jx;
bli_zsets( (bli_zbla_real(x[i__4])), (bli_zimag(x[i__4])), temp ); bli_zsets( (bli_zreal(x[i__4])), (bli_zimag(x[i__4])), temp );
ix = kx; ix = kx;
l = 1 - j; l = 1 - j;
if (noconj) { if (noconj) {
@@ -1944,15 +1944,15 @@
for (i__ = f2c_min(i__4,i__2); i__ >= i__1; --i__) { for (i__ = f2c_min(i__4,i__2); i__ >= i__1; --i__) {
i__4 = l + i__ + j * a_dim1; i__4 = l + i__ + j * a_dim1;
i__2 = ix; i__2 = ix;
bli_zsets( (bli_zbla_real(a[i__4]) * bli_zbla_real(x[i__2]) - bli_zimag(a[i__4]) * bli_zimag(x[i__2])), (bli_zbla_real(a[i__4]) * bli_zimag(x[i__2]) + bli_zimag(a[i__4]) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(a[i__4]) * bli_zreal(x[i__2]) - bli_zimag(a[i__4]) * bli_zimag(x[i__2])), (bli_zreal(a[i__4]) * bli_zimag(x[i__2]) + bli_zimag(a[i__4]) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix -= *incx; ix -= *incx;
/* L180: */ /* L180: */
} }
if (nounit) { if (nounit) {
bla_z_div(&z__1, &temp, &a[j * a_dim1 + 1]); bla_z_div(&z__1, &temp, &a[j * a_dim1 + 1]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} else { } else {
/* Computing MIN */ /* Computing MIN */
@@ -1961,20 +1961,20 @@
for (i__ = f2c_min(i__1,i__4); i__ >= i__2; --i__) { for (i__ = f2c_min(i__1,i__4); i__ >= i__2; --i__) {
bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]);
i__1 = ix; i__1 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__1]) - bli_zimag(z__3) * bli_zimag(x[i__1])), (bli_zbla_real(z__3) * bli_zimag(x[i__1]) + bli_zimag(z__3) * bli_zbla_real(x[i__1])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__1]) - bli_zimag(z__3) * bli_zimag(x[i__1])), (bli_zreal(z__3) * bli_zimag(x[i__1]) + bli_zimag(z__3) * bli_zreal(x[i__1])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix -= *incx; ix -= *incx;
/* L190: */ /* L190: */
} }
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &a[j * a_dim1 + 1]); bla_d_cnjg(&z__2, &a[j * a_dim1 + 1]);
bla_z_div(&z__1, &temp, &z__2); bla_z_div(&z__1, &temp, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} }
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__2] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__2] );
jx -= *incx; jx -= *incx;
if (*n - j >= *k) { if (*n - j >= *k) {
kx -= *incx; kx -= *incx;

152
frame/compat/f2c/bla_tpmv.c
View File

@@ -1415,18 +1415,18 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
k = kk; k = kk;
i__2 = j - 1; i__2 = j - 1;
for (i__ = 1; i__ <= i__2; ++i__) { for (i__ = 1; i__ <= i__2; ++i__) {
i__3 = i__; i__3 = i__;
i__4 = i__; i__4 = i__;
i__5 = k; i__5 = k;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__5]) - bli_zimag(temp) * bli_zimag(ap[i__5])), (bli_zbla_real(temp) * bli_zimag(ap[i__5]) + bli_zimag(temp) * bli_zbla_real(ap[i__5])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__5]) - bli_zimag(temp) * bli_zimag(ap[i__5])), (bli_zreal(temp) * bli_zimag(ap[i__5]) + bli_zimag(temp) * bli_zreal(ap[i__5])), z__2 );
bli_zsets( (bli_zbla_real(x[i__4]) + bli_zbla_real(z__2)), (bli_zimag(x[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__4]) + bli_zreal(z__2)), (bli_zimag(x[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] );
++k; ++k;
/* L10: */ /* L10: */
} }
@@ -1434,8 +1434,8 @@
i__2 = j; i__2 = j;
i__3 = j; i__3 = j;
i__4 = kk + j - 1; i__4 = kk + j - 1;
bli_zsets( (bli_zbla_real(x[i__3]) * bli_zbla_real(ap[i__4]) - bli_zimag(x[i__3]) * bli_zimag(ap[i__4])), (bli_zbla_real(x[i__3]) * bli_zimag(ap[i__4]) + bli_zimag(x[i__3]) * bli_zbla_real(ap[i__4])), z__1 ); bli_zsets( (bli_zreal(x[i__3]) * bli_zreal(ap[i__4]) - bli_zimag(x[i__3]) * bli_zimag(ap[i__4])), (bli_zreal(x[i__3]) * bli_zimag(ap[i__4]) + bli_zimag(x[i__3]) * bli_zreal(ap[i__4])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
} }
} }
kk += j; kk += j;
@@ -1446,18 +1446,18 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
ix = kx; ix = kx;
i__2 = kk + j - 2; i__2 = kk + j - 2;
for (k = kk; k <= i__2; ++k) { for (k = kk; k <= i__2; ++k) {
i__3 = ix; i__3 = ix;
i__4 = ix; i__4 = ix;
i__5 = k; i__5 = k;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__5]) - bli_zimag(temp) * bli_zimag(ap[i__5])), (bli_zbla_real(temp) * bli_zimag(ap[i__5]) + bli_zimag(temp) * bli_zbla_real(ap[i__5])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__5]) - bli_zimag(temp) * bli_zimag(ap[i__5])), (bli_zreal(temp) * bli_zimag(ap[i__5]) + bli_zimag(temp) * bli_zreal(ap[i__5])), z__2 );
bli_zsets( (bli_zbla_real(x[i__4]) + bli_zbla_real(z__2)), (bli_zimag(x[i__4]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__4]) + bli_zreal(z__2)), (bli_zimag(x[i__4]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] );
ix += *incx; ix += *incx;
/* L30: */ /* L30: */
} }
@@ -1465,8 +1465,8 @@
i__2 = jx; i__2 = jx;
i__3 = jx; i__3 = jx;
i__4 = kk + j - 1; i__4 = kk + j - 1;
bli_zsets( (bli_zbla_real(x[i__3]) * bli_zbla_real(ap[i__4]) - bli_zimag(x[i__3]) * bli_zimag(ap[i__4])), (bli_zbla_real(x[i__3]) * bli_zimag(ap[i__4]) + bli_zimag(x[i__3]) * bli_zbla_real(ap[i__4])), z__1 ); bli_zsets( (bli_zreal(x[i__3]) * bli_zreal(ap[i__4]) - bli_zimag(x[i__3]) * bli_zimag(ap[i__4])), (bli_zreal(x[i__3]) * bli_zimag(ap[i__4]) + bli_zimag(x[i__3]) * bli_zreal(ap[i__4])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
} }
} }
jx += *incx; jx += *incx;
@@ -1479,18 +1479,18 @@
if (*incx == 1) { if (*incx == 1) {
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = j; i__1 = j;
if (bli_zbla_real(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) { if (bli_zreal(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) {
i__1 = j; i__1 = j;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
k = kk; k = kk;
i__1 = j + 1; i__1 = j + 1;
for (i__ = *n; i__ >= i__1; --i__) { for (i__ = *n; i__ >= i__1; --i__) {
i__2 = i__; i__2 = i__;
i__3 = i__; i__3 = i__;
i__4 = k; i__4 = k;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__4]) - bli_zimag(temp) * bli_zimag(ap[i__4])), (bli_zbla_real(temp) * bli_zimag(ap[i__4]) + bli_zimag(temp) * bli_zbla_real(ap[i__4])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__4]) - bli_zimag(temp) * bli_zimag(ap[i__4])), (bli_zreal(temp) * bli_zimag(ap[i__4]) + bli_zimag(temp) * bli_zreal(ap[i__4])), z__2 );
bli_zsets( (bli_zbla_real(x[i__3]) + bli_zbla_real(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__3]) + bli_zreal(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
--k; --k;
/* L50: */ /* L50: */
} }
@@ -1498,8 +1498,8 @@
i__1 = j; i__1 = j;
i__2 = j; i__2 = j;
i__3 = kk - *n + j; i__3 = kk - *n + j;
bli_zsets( (bli_zbla_real(x[i__2]) * bli_zbla_real(ap[i__3]) - bli_zimag(x[i__2]) * bli_zimag(ap[i__3])), (bli_zbla_real(x[i__2]) * bli_zimag(ap[i__3]) + bli_zimag(x[i__2]) * bli_zbla_real(ap[i__3])), z__1 ); bli_zsets( (bli_zreal(x[i__2]) * bli_zreal(ap[i__3]) - bli_zimag(x[i__2]) * bli_zimag(ap[i__3])), (bli_zreal(x[i__2]) * bli_zimag(ap[i__3]) + bli_zimag(x[i__2]) * bli_zreal(ap[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__1] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__1] );
} }
} }
kk -= *n - j + 1; kk -= *n - j + 1;
@@ -1510,18 +1510,18 @@
jx = kx; jx = kx;
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = jx; i__1 = jx;
if (bli_zbla_real(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) { if (bli_zreal(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) {
i__1 = jx; i__1 = jx;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
ix = kx; ix = kx;
i__1 = kk - (*n - (j + 1)); i__1 = kk - (*n - (j + 1));
for (k = kk; k >= i__1; --k) { for (k = kk; k >= i__1; --k) {
i__2 = ix; i__2 = ix;
i__3 = ix; i__3 = ix;
i__4 = k; i__4 = k;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__4]) - bli_zimag(temp) * bli_zimag(ap[i__4])), (bli_zbla_real(temp) * bli_zimag(ap[i__4]) + bli_zimag(temp) * bli_zbla_real(ap[i__4])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__4]) - bli_zimag(temp) * bli_zimag(ap[i__4])), (bli_zreal(temp) * bli_zimag(ap[i__4]) + bli_zimag(temp) * bli_zreal(ap[i__4])), z__2 );
bli_zsets( (bli_zbla_real(x[i__3]) + bli_zbla_real(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__3]) + bli_zreal(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
ix -= *incx; ix -= *incx;
/* L70: */ /* L70: */
} }
@@ -1529,8 +1529,8 @@
i__1 = jx; i__1 = jx;
i__2 = jx; i__2 = jx;
i__3 = kk - *n + j; i__3 = kk - *n + j;
bli_zsets( (bli_zbla_real(x[i__2]) * bli_zbla_real(ap[i__3]) - bli_zimag(x[i__2]) * bli_zimag(ap[i__3])), (bli_zbla_real(x[i__2]) * bli_zimag(ap[i__3]) + bli_zimag(x[i__2]) * bli_zbla_real(ap[i__3])), z__1 ); bli_zsets( (bli_zreal(x[i__2]) * bli_zreal(ap[i__3]) - bli_zimag(x[i__2]) * bli_zimag(ap[i__3])), (bli_zreal(x[i__2]) * bli_zimag(ap[i__3]) + bli_zimag(x[i__2]) * bli_zreal(ap[i__3])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__1] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__1] );
} }
} }
jx -= *incx; jx -= *incx;
@@ -1548,41 +1548,41 @@
if (*incx == 1) { if (*incx == 1) {
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = j; i__1 = j;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
k = kk - 1; k = kk - 1;
if (noconj) { if (noconj) {
if (nounit) { if (nounit) {
i__1 = kk; i__1 = kk;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__1]) - bli_zimag(temp) * bli_zimag(ap[i__1])), (bli_zbla_real(temp) * bli_zimag(ap[i__1]) + bli_zimag(temp) * bli_zbla_real(ap[i__1])), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__1]) - bli_zimag(temp) * bli_zimag(ap[i__1])), (bli_zreal(temp) * bli_zimag(ap[i__1]) + bli_zimag(temp) * bli_zreal(ap[i__1])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
for (i__ = j - 1; i__ >= 1; --i__) { for (i__ = j - 1; i__ >= 1; --i__) {
i__1 = k; i__1 = k;
i__2 = i__; i__2 = i__;
bli_zsets( (bli_zbla_real(ap[i__1]) * bli_zbla_real(x[i__2]) - bli_zimag(ap[i__1]) * bli_zimag(x[i__2])), (bli_zbla_real(ap[i__1]) * bli_zimag(x[i__2]) + bli_zimag(ap[i__1]) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(ap[i__1]) * bli_zreal(x[i__2]) - bli_zimag(ap[i__1]) * bli_zimag(x[i__2])), (bli_zreal(ap[i__1]) * bli_zimag(x[i__2]) + bli_zimag(ap[i__1]) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
--k; --k;
/* L90: */ /* L90: */
} }
} else { } else {
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &ap[kk]); bla_d_cnjg(&z__2, &ap[kk]);
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zbla_real(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zreal(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
for (i__ = j - 1; i__ >= 1; --i__) { for (i__ = j - 1; i__ >= 1; --i__) {
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__1 = i__; i__1 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__1]) - bli_zimag(z__3) * bli_zimag(x[i__1])), (bli_zbla_real(z__3) * bli_zimag(x[i__1]) + bli_zimag(z__3) * bli_zbla_real(x[i__1])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__1]) - bli_zimag(z__3) * bli_zimag(x[i__1])), (bli_zreal(z__3) * bli_zimag(x[i__1]) + bli_zimag(z__3) * bli_zreal(x[i__1])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
--k; --k;
/* L100: */ /* L100: */
} }
} }
i__1 = j; i__1 = j;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__1] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__1] );
kk -= j; kk -= j;
/* L110: */ /* L110: */
} }
@@ -1590,43 +1590,43 @@
jx = kx + (*n - 1) * *incx; jx = kx + (*n - 1) * *incx;
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = jx; i__1 = jx;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
ix = jx; ix = jx;
if (noconj) { if (noconj) {
if (nounit) { if (nounit) {
i__1 = kk; i__1 = kk;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__1]) - bli_zimag(temp) * bli_zimag(ap[i__1])), (bli_zbla_real(temp) * bli_zimag(ap[i__1]) + bli_zimag(temp) * bli_zbla_real(ap[i__1])), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__1]) - bli_zimag(temp) * bli_zimag(ap[i__1])), (bli_zreal(temp) * bli_zimag(ap[i__1]) + bli_zimag(temp) * bli_zreal(ap[i__1])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
i__1 = kk - j + 1; i__1 = kk - j + 1;
for (k = kk - 1; k >= i__1; --k) { for (k = kk - 1; k >= i__1; --k) {
ix -= *incx; ix -= *incx;
i__2 = k; i__2 = k;
i__3 = ix; i__3 = ix;
bli_zsets( (bli_zbla_real(ap[i__2]) * bli_zbla_real(x[i__3]) - bli_zimag(ap[i__2]) * bli_zimag(x[i__3])), (bli_zbla_real(ap[i__2]) * bli_zimag(x[i__3]) + bli_zimag(ap[i__2]) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(ap[i__2]) * bli_zreal(x[i__3]) - bli_zimag(ap[i__2]) * bli_zimag(x[i__3])), (bli_zreal(ap[i__2]) * bli_zimag(x[i__3]) + bli_zimag(ap[i__2]) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L120: */ /* L120: */
} }
} else { } else {
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &ap[kk]); bla_d_cnjg(&z__2, &ap[kk]);
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zbla_real(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zreal(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
i__1 = kk - j + 1; i__1 = kk - j + 1;
for (k = kk - 1; k >= i__1; --k) { for (k = kk - 1; k >= i__1; --k) {
ix -= *incx; ix -= *incx;
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__2 = ix; i__2 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zbla_real(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zreal(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L130: */ /* L130: */
} }
} }
i__1 = jx; i__1 = jx;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__1] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__1] );
jx -= *incx; jx -= *incx;
kk -= j; kk -= j;
/* L140: */ /* L140: */
@@ -1638,43 +1638,43 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
k = kk + 1; k = kk + 1;
if (noconj) { if (noconj) {
if (nounit) { if (nounit) {
i__2 = kk; i__2 = kk;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__2]) - bli_zimag(temp) * bli_zimag(ap[i__2])), (bli_zbla_real(temp) * bli_zimag(ap[i__2]) + bli_zimag(temp) * bli_zbla_real(ap[i__2])), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__2]) - bli_zimag(temp) * bli_zimag(ap[i__2])), (bli_zreal(temp) * bli_zimag(ap[i__2]) + bli_zimag(temp) * bli_zreal(ap[i__2])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
i__2 = *n; i__2 = *n;
for (i__ = j + 1; i__ <= i__2; ++i__) { for (i__ = j + 1; i__ <= i__2; ++i__) {
i__3 = k; i__3 = k;
i__4 = i__; i__4 = i__;
bli_zsets( (bli_zbla_real(ap[i__3]) * bli_zbla_real(x[i__4]) - bli_zimag(ap[i__3]) * bli_zimag(x[i__4])), (bli_zbla_real(ap[i__3]) * bli_zimag(x[i__4]) + bli_zimag(ap[i__3]) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(ap[i__3]) * bli_zreal(x[i__4]) - bli_zimag(ap[i__3]) * bli_zimag(x[i__4])), (bli_zreal(ap[i__3]) * bli_zimag(x[i__4]) + bli_zimag(ap[i__3]) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
++k; ++k;
/* L150: */ /* L150: */
} }
} else { } else {
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &ap[kk]); bla_d_cnjg(&z__2, &ap[kk]);
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zbla_real(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zreal(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
i__2 = *n; i__2 = *n;
for (i__ = j + 1; i__ <= i__2; ++i__) { for (i__ = j + 1; i__ <= i__2; ++i__) {
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__3 = i__; i__3 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zbla_real(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zreal(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
++k; ++k;
/* L160: */ /* L160: */
} }
} }
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__2] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__2] );
kk += *n - j + 1; kk += *n - j + 1;
/* L170: */ /* L170: */
} }
@@ -1683,43 +1683,43 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
ix = jx; ix = jx;
if (noconj) { if (noconj) {
if (nounit) { if (nounit) {
i__2 = kk; i__2 = kk;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__2]) - bli_zimag(temp) * bli_zimag(ap[i__2])), (bli_zbla_real(temp) * bli_zimag(ap[i__2]) + bli_zimag(temp) * bli_zbla_real(ap[i__2])), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__2]) - bli_zimag(temp) * bli_zimag(ap[i__2])), (bli_zreal(temp) * bli_zimag(ap[i__2]) + bli_zimag(temp) * bli_zreal(ap[i__2])), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
i__2 = kk + *n - j; i__2 = kk + *n - j;
for (k = kk + 1; k <= i__2; ++k) { for (k = kk + 1; k <= i__2; ++k) {
ix += *incx; ix += *incx;
i__3 = k; i__3 = k;
i__4 = ix; i__4 = ix;
bli_zsets( (bli_zbla_real(ap[i__3]) * bli_zbla_real(x[i__4]) - bli_zimag(ap[i__3]) * bli_zimag(x[i__4])), (bli_zbla_real(ap[i__3]) * bli_zimag(x[i__4]) + bli_zimag(ap[i__3]) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(ap[i__3]) * bli_zreal(x[i__4]) - bli_zimag(ap[i__3]) * bli_zimag(x[i__4])), (bli_zreal(ap[i__3]) * bli_zimag(x[i__4]) + bli_zimag(ap[i__3]) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L180: */ /* L180: */
} }
} else { } else {
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &ap[kk]); bla_d_cnjg(&z__2, &ap[kk]);
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zbla_real(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zbla_real(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) * bli_zreal(z__2) - bli_zimag(temp) * bli_zimag(z__2)), (bli_zreal(temp) * bli_zimag(z__2) + bli_zimag(temp) * bli_zreal(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
i__2 = kk + *n - j; i__2 = kk + *n - j;
for (k = kk + 1; k <= i__2; ++k) { for (k = kk + 1; k <= i__2; ++k) {
ix += *incx; ix += *incx;
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__3 = ix; i__3 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zbla_real(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zreal(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp) + bli_zbla_real(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
/* L190: */ /* L190: */
} }
} }
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__2] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__2] );
jx += *incx; jx += *incx;
kk += *n - j + 1; kk += *n - j + 1;
/* L200: */ /* L200: */

128
frame/compat/f2c/bla_tpsv.c
View File

@@ -1416,22 +1416,22 @@
if (*incx == 1) { if (*incx == 1) {
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = j; i__1 = j;
if (bli_zbla_real(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) { if (bli_zreal(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) {
if (nounit) { if (nounit) {
i__1 = j; i__1 = j;
bla_z_div(&z__1, &x[j], &ap[kk]); bla_z_div(&z__1, &x[j], &ap[kk]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__1] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__1] );
} }
i__1 = j; i__1 = j;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
k = kk - 1; k = kk - 1;
for (i__ = j - 1; i__ >= 1; --i__) { for (i__ = j - 1; i__ >= 1; --i__) {
i__1 = i__; i__1 = i__;
i__2 = i__; i__2 = i__;
i__3 = k; i__3 = k;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__3]) - bli_zimag(temp) * bli_zimag(ap[i__3])), (bli_zbla_real(temp) * bli_zimag(ap[i__3]) + bli_zimag(temp) * bli_zbla_real(ap[i__3])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__3]) - bli_zimag(temp) * bli_zimag(ap[i__3])), (bli_zreal(temp) * bli_zimag(ap[i__3]) + bli_zimag(temp) * bli_zreal(ap[i__3])), z__2 );
bli_zsets( (bli_zbla_real(x[i__2]) - bli_zbla_real(z__2)), (bli_zimag(x[i__2]) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__2]) - bli_zreal(z__2)), (bli_zimag(x[i__2]) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__1] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__1] );
--k; --k;
/* L10: */ /* L10: */
} }
@@ -1443,14 +1443,14 @@
jx = kx + (*n - 1) * *incx; jx = kx + (*n - 1) * *incx;
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = jx; i__1 = jx;
if (bli_zbla_real(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) { if (bli_zreal(x[i__1]) != 0. || bli_zimag(x[i__1]) != 0.) {
if (nounit) { if (nounit) {
i__1 = jx; i__1 = jx;
bla_z_div(&z__1, &x[jx], &ap[kk]); bla_z_div(&z__1, &x[jx], &ap[kk]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__1] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__1] );
} }
i__1 = jx; i__1 = jx;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
ix = jx; ix = jx;
i__1 = kk - j + 1; i__1 = kk - j + 1;
for (k = kk - 1; k >= i__1; --k) { for (k = kk - 1; k >= i__1; --k) {
@@ -1458,9 +1458,9 @@
i__2 = ix; i__2 = ix;
i__3 = ix; i__3 = ix;
i__4 = k; i__4 = k;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__4]) - bli_zimag(temp) * bli_zimag(ap[i__4])), (bli_zbla_real(temp) * bli_zimag(ap[i__4]) + bli_zimag(temp) * bli_zbla_real(ap[i__4])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__4]) - bli_zimag(temp) * bli_zimag(ap[i__4])), (bli_zreal(temp) * bli_zimag(ap[i__4]) + bli_zimag(temp) * bli_zreal(ap[i__4])), z__2 );
bli_zsets( (bli_zbla_real(x[i__3]) - bli_zbla_real(z__2)), (bli_zimag(x[i__3]) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__3]) - bli_zreal(z__2)), (bli_zimag(x[i__3]) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
/* L30: */ /* L30: */
} }
} }
@@ -1475,23 +1475,23 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
if (nounit) { if (nounit) {
i__2 = j; i__2 = j;
bla_z_div(&z__1, &x[j], &ap[kk]); bla_z_div(&z__1, &x[j], &ap[kk]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
} }
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
k = kk + 1; k = kk + 1;
i__2 = *n; i__2 = *n;
for (i__ = j + 1; i__ <= i__2; ++i__) { for (i__ = j + 1; i__ <= i__2; ++i__) {
i__3 = i__; i__3 = i__;
i__4 = i__; i__4 = i__;
i__5 = k; i__5 = k;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__5]) - bli_zimag(temp) * bli_zimag(ap[i__5])), (bli_zbla_real(temp) * bli_zimag(ap[i__5]) + bli_zimag(temp) * bli_zbla_real(ap[i__5])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__5]) - bli_zimag(temp) * bli_zimag(ap[i__5])), (bli_zreal(temp) * bli_zimag(ap[i__5]) + bli_zimag(temp) * bli_zreal(ap[i__5])), z__2 );
bli_zsets( (bli_zbla_real(x[i__4]) - bli_zbla_real(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__4]) - bli_zreal(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] );
++k; ++k;
/* L50: */ /* L50: */
} }
@@ -1504,14 +1504,14 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
if (bli_zbla_real(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) { if (bli_zreal(x[i__2]) != 0. || bli_zimag(x[i__2]) != 0.) {
if (nounit) { if (nounit) {
i__2 = jx; i__2 = jx;
bla_z_div(&z__1, &x[jx], &ap[kk]); bla_z_div(&z__1, &x[jx], &ap[kk]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] );
} }
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
ix = jx; ix = jx;
i__2 = kk + *n - j; i__2 = kk + *n - j;
for (k = kk + 1; k <= i__2; ++k) { for (k = kk + 1; k <= i__2; ++k) {
@@ -1519,9 +1519,9 @@
i__3 = ix; i__3 = ix;
i__4 = ix; i__4 = ix;
i__5 = k; i__5 = k;
bli_zsets( (bli_zbla_real(temp) * bli_zbla_real(ap[i__5]) - bli_zimag(temp) * bli_zimag(ap[i__5])), (bli_zbla_real(temp) * bli_zimag(ap[i__5]) + bli_zimag(temp) * bli_zbla_real(ap[i__5])), z__2 ); bli_zsets( (bli_zreal(temp) * bli_zreal(ap[i__5]) - bli_zimag(temp) * bli_zimag(ap[i__5])), (bli_zreal(temp) * bli_zimag(ap[i__5]) + bli_zimag(temp) * bli_zreal(ap[i__5])), z__2 );
bli_zsets( (bli_zbla_real(x[i__4]) - bli_zbla_real(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(x[i__4]) - bli_zreal(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), x[i__3] ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] );
/* L70: */ /* L70: */
} }
} }
@@ -1541,42 +1541,42 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
k = kk; k = kk;
if (noconj) { if (noconj) {
i__2 = j - 1; i__2 = j - 1;
for (i__ = 1; i__ <= i__2; ++i__) { for (i__ = 1; i__ <= i__2; ++i__) {
i__3 = k; i__3 = k;
i__4 = i__; i__4 = i__;
bli_zsets( (bli_zbla_real(ap[i__3]) * bli_zbla_real(x[i__4]) - bli_zimag(ap[i__3]) * bli_zimag(x[i__4])), (bli_zbla_real(ap[i__3]) * bli_zimag(x[i__4]) + bli_zimag(ap[i__3]) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(ap[i__3]) * bli_zreal(x[i__4]) - bli_zimag(ap[i__3]) * bli_zimag(x[i__4])), (bli_zreal(ap[i__3]) * bli_zimag(x[i__4]) + bli_zimag(ap[i__3]) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
++k; ++k;
/* L90: */ /* L90: */
} }
if (nounit) { if (nounit) {
bla_z_div(&z__1, &temp, &ap[kk + j - 1]); bla_z_div(&z__1, &temp, &ap[kk + j - 1]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} else { } else {
i__2 = j - 1; i__2 = j - 1;
for (i__ = 1; i__ <= i__2; ++i__) { for (i__ = 1; i__ <= i__2; ++i__) {
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__3 = i__; i__3 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zbla_real(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zreal(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
++k; ++k;
/* L100: */ /* L100: */
} }
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &ap[kk + j - 1]); bla_d_cnjg(&z__2, &ap[kk + j - 1]);
bla_z_div(&z__1, &temp, &z__2); bla_z_div(&z__1, &temp, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} }
i__2 = j; i__2 = j;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__2] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__2] );
kk += j; kk += j;
/* L110: */ /* L110: */
} }
@@ -1585,42 +1585,42 @@
i__1 = *n; i__1 = *n;
for (j = 1; j <= i__1; ++j) { for (j = 1; j <= i__1; ++j) {
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(x[i__2])), (bli_zimag(x[i__2])), temp ); bli_zsets( (bli_zreal(x[i__2])), (bli_zimag(x[i__2])), temp );
ix = kx; ix = kx;
if (noconj) { if (noconj) {
i__2 = kk + j - 2; i__2 = kk + j - 2;
for (k = kk; k <= i__2; ++k) { for (k = kk; k <= i__2; ++k) {
i__3 = k; i__3 = k;
i__4 = ix; i__4 = ix;
bli_zsets( (bli_zbla_real(ap[i__3]) * bli_zbla_real(x[i__4]) - bli_zimag(ap[i__3]) * bli_zimag(x[i__4])), (bli_zbla_real(ap[i__3]) * bli_zimag(x[i__4]) + bli_zimag(ap[i__3]) * bli_zbla_real(x[i__4])), z__2 ); bli_zsets( (bli_zreal(ap[i__3]) * bli_zreal(x[i__4]) - bli_zimag(ap[i__3]) * bli_zimag(x[i__4])), (bli_zreal(ap[i__3]) * bli_zimag(x[i__4]) + bli_zimag(ap[i__3]) * bli_zreal(x[i__4])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix += *incx; ix += *incx;
/* L120: */ /* L120: */
} }
if (nounit) { if (nounit) {
bla_z_div(&z__1, &temp, &ap[kk + j - 1]); bla_z_div(&z__1, &temp, &ap[kk + j - 1]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} else { } else {
i__2 = kk + j - 2; i__2 = kk + j - 2;
for (k = kk; k <= i__2; ++k) { for (k = kk; k <= i__2; ++k) {
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__3 = ix; i__3 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zbla_real(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__3]) - bli_zimag(z__3) * bli_zimag(x[i__3])), (bli_zreal(z__3) * bli_zimag(x[i__3]) + bli_zimag(z__3) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix += *incx; ix += *incx;
/* L130: */ /* L130: */
} }
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &ap[kk + j - 1]); bla_d_cnjg(&z__2, &ap[kk + j - 1]);
bla_z_div(&z__1, &temp, &z__2); bla_z_div(&z__1, &temp, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} }
i__2 = jx; i__2 = jx;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__2] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__2] );
jx += *incx; jx += *incx;
kk += j; kk += j;
/* L140: */ /* L140: */
@@ -1631,42 +1631,42 @@
if (*incx == 1) { if (*incx == 1) {
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = j; i__1 = j;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
k = kk; k = kk;
if (noconj) { if (noconj) {
i__1 = j + 1; i__1 = j + 1;
for (i__ = *n; i__ >= i__1; --i__) { for (i__ = *n; i__ >= i__1; --i__) {
i__2 = k; i__2 = k;
i__3 = i__; i__3 = i__;
bli_zsets( (bli_zbla_real(ap[i__2]) * bli_zbla_real(x[i__3]) - bli_zimag(ap[i__2]) * bli_zimag(x[i__3])), (bli_zbla_real(ap[i__2]) * bli_zimag(x[i__3]) + bli_zimag(ap[i__2]) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(ap[i__2]) * bli_zreal(x[i__3]) - bli_zimag(ap[i__2]) * bli_zimag(x[i__3])), (bli_zreal(ap[i__2]) * bli_zimag(x[i__3]) + bli_zimag(ap[i__2]) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
--k; --k;
/* L150: */ /* L150: */
} }
if (nounit) { if (nounit) {
bla_z_div(&z__1, &temp, &ap[kk - *n + j]); bla_z_div(&z__1, &temp, &ap[kk - *n + j]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} else { } else {
i__1 = j + 1; i__1 = j + 1;
for (i__ = *n; i__ >= i__1; --i__) { for (i__ = *n; i__ >= i__1; --i__) {
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__2 = i__; i__2 = i__;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zbla_real(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zreal(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
--k; --k;
/* L160: */ /* L160: */
} }
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &ap[kk - *n + j]); bla_d_cnjg(&z__2, &ap[kk - *n + j]);
bla_z_div(&z__1, &temp, &z__2); bla_z_div(&z__1, &temp, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} }
i__1 = j; i__1 = j;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__1] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__1] );
kk -= *n - j + 1; kk -= *n - j + 1;
/* L170: */ /* L170: */
} }
@@ -1675,42 +1675,42 @@
jx = kx; jx = kx;
for (j = *n; j >= 1; --j) { for (j = *n; j >= 1; --j) {
i__1 = jx; i__1 = jx;
bli_zsets( (bli_zbla_real(x[i__1])), (bli_zimag(x[i__1])), temp ); bli_zsets( (bli_zreal(x[i__1])), (bli_zimag(x[i__1])), temp );
ix = kx; ix = kx;
if (noconj) { if (noconj) {
i__1 = kk - (*n - (j + 1)); i__1 = kk - (*n - (j + 1));
for (k = kk; k >= i__1; --k) { for (k = kk; k >= i__1; --k) {
i__2 = k; i__2 = k;
i__3 = ix; i__3 = ix;
bli_zsets( (bli_zbla_real(ap[i__2]) * bli_zbla_real(x[i__3]) - bli_zimag(ap[i__2]) * bli_zimag(x[i__3])), (bli_zbla_real(ap[i__2]) * bli_zimag(x[i__3]) + bli_zimag(ap[i__2]) * bli_zbla_real(x[i__3])), z__2 ); bli_zsets( (bli_zreal(ap[i__2]) * bli_zreal(x[i__3]) - bli_zimag(ap[i__2]) * bli_zimag(x[i__3])), (bli_zreal(ap[i__2]) * bli_zimag(x[i__3]) + bli_zimag(ap[i__2]) * bli_zreal(x[i__3])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix -= *incx; ix -= *incx;
/* L180: */ /* L180: */
} }
if (nounit) { if (nounit) {
bla_z_div(&z__1, &temp, &ap[kk - *n + j]); bla_z_div(&z__1, &temp, &ap[kk - *n + j]);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} else { } else {
i__1 = kk - (*n - (j + 1)); i__1 = kk - (*n - (j + 1));
for (k = kk; k >= i__1; --k) { for (k = kk; k >= i__1; --k) {
bla_d_cnjg(&z__3, &ap[k]); bla_d_cnjg(&z__3, &ap[k]);
i__2 = ix; i__2 = ix;
bli_zsets( (bli_zbla_real(z__3) * bli_zbla_real(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zbla_real(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zbla_real(x[i__2])), z__2 ); bli_zsets( (bli_zreal(z__3) * bli_zreal(x[i__2]) - bli_zimag(z__3) * bli_zimag(x[i__2])), (bli_zreal(z__3) * bli_zimag(x[i__2]) + bli_zimag(z__3) * bli_zreal(x[i__2])), z__2 );
bli_zsets( (bli_zbla_real(temp) - bli_zbla_real(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); bli_zsets( (bli_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 );
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
ix -= *incx; ix -= *incx;
/* L190: */ /* L190: */
} }
if (nounit) { if (nounit) {
bla_d_cnjg(&z__2, &ap[kk - *n + j]); bla_d_cnjg(&z__2, &ap[kk - *n + j]);
bla_z_div(&z__1, &temp, &z__2); bla_z_div(&z__1, &temp, &z__2);
bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp );
} }
} }
i__1 = jx; i__1 = jx;
bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__1] ); bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__1] );
jx -= *incx; jx -= *incx;
kk -= *n - j + 1; kk -= *n - j + 1;
/* L200: */ /* L200: */

2
frame/compat/f2c/util/bla_z_abs.c
View File

@@ -38,7 +38,7 @@
double bla_z_abs(bla_dcomplex *z) double bla_z_abs(bla_dcomplex *z)
{ {
return( bla_f__cabs( bli_zbla_real( *z ), return( bla_f__cabs( bli_zreal( *z ),
bli_zimag( *z ) ) ); bli_zimag( *z ) ) );
} }