From 8624a33ccc12dff6f6c4f92992ca5636af1576a6 Mon Sep 17 00:00:00 2001 From: Devin Matthews Date: Thu, 25 Feb 2016 13:51:26 -0600 Subject: [PATCH] Fix remaining f2c conflicts. --- frame/compat/f2c/bla_gbmv.c | 78 +++++++-------- frame/compat/f2c/bla_hbmv.c | 134 +++++++++++++------------- frame/compat/f2c/bla_hpmv.c | 134 +++++++++++++------------- frame/compat/f2c/bla_hpr.c | 72 +++++++------- frame/compat/f2c/bla_hpr2.c | 122 ++++++++++++------------ frame/compat/f2c/bla_rot.c | 36 +++---- frame/compat/f2c/bla_rotg.c | 16 ++-- frame/compat/f2c/bla_tbmv.c | 152 +++++++++++++++--------------- frame/compat/f2c/bla_tbsv.c | 128 ++++++++++++------------- frame/compat/f2c/bla_tpmv.c | 152 +++++++++++++++--------------- frame/compat/f2c/bla_tpsv.c | 128 ++++++++++++------------- frame/compat/f2c/util/bla_z_abs.c | 2 +- 12 files changed, 577 insertions(+), 577 deletions(-) diff --git a/frame/compat/f2c/bla_gbmv.c b/frame/compat/f2c/bla_gbmv.c index 9c19ec734..0c66786d7 100644 --- a/frame/compat/f2c/bla_gbmv.c +++ b/frame/compat/f2c/bla_gbmv.c @@ -1382,7 +1382,7 @@ /* 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.))) { return 0; } @@ -1415,9 +1415,9 @@ /* 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 (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { + if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) { i__1 = leny; for (i__ = 1; i__ <= i__1; ++i__) { i__2 = i__; @@ -1429,14 +1429,14 @@ for (i__ = 1; i__ <= i__1; ++i__) { i__2 = 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_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); + 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_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); /* L20: */ } } } else { iy = ky; - if (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { + if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) { i__1 = leny; for (i__ = 1; i__ <= i__1; ++i__) { i__2 = iy; @@ -1449,15 +1449,15 @@ for (i__ = 1; i__ <= i__1; ++i__) { i__2 = 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_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); + 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_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); iy += *incy; /* L40: */ } } } } - if (bli_zbla_real(*alpha) == 0. && bli_zimag(*alpha) == 0.) { + if (bli_zreal(*alpha) == 0. && bli_zimag(*alpha) == 0.) { return 0; } kup1 = *ku + 1; @@ -1470,10 +1470,10 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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; - 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); k = kup1 - j; /* Computing MAX */ i__2 = 1, i__3 = j - *ku; @@ -1484,9 +1484,9 @@ i__2 = i__; i__3 = i__; 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_zbla_real(y[i__3]) + bli_zbla_real(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(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_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); /* L50: */ } } @@ -1497,10 +1497,10 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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; - 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); iy = ky; k = kup1 - j; /* Computing MAX */ @@ -1512,9 +1512,9 @@ i__4 = iy; i__2 = iy; 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_zbla_real(y[i__2]) + bli_zbla_real(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(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_zreal(y[i__2]) + bli_zreal(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] ); iy += *incy; /* L70: */ } @@ -1545,9 +1545,9 @@ for (i__ = f2c_max(i__3,i__4); i__ <= i__2; ++i__) { i__3 = k + i__ + j * a_dim1; 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_zbla_real(temp) + bli_zbla_real(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(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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L90: */ } } else { @@ -1559,17 +1559,17 @@ for (i__ = f2c_max(i__2,i__3); i__ <= i__4; ++i__) { bla_d_cnjg(&z__3, &a[k + i__ + j * a_dim1]); 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_zbla_real(temp) + bli_zbla_real(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__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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L100: */ } } i__4 = 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_zbla_real(y[i__2]) + bli_zbla_real(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(*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_zreal(y[i__2]) + bli_zreal(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] ); jy += *incy; /* L110: */ } @@ -1588,9 +1588,9 @@ for (i__ = f2c_max(i__4,i__2); i__ <= i__3; ++i__) { i__4 = k + i__ + j * a_dim1; 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_zbla_real(temp) + bli_zbla_real(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(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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix += *incx; /* L120: */ } @@ -1603,18 +1603,18 @@ for (i__ = f2c_max(i__3,i__4); i__ <= i__2; ++i__) { bla_d_cnjg(&z__3, &a[k + i__ + j * a_dim1]); 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_zbla_real(temp) + bli_zbla_real(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__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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix += *incx; /* L130: */ } } i__2 = 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_zbla_real(y[i__3]) + bli_zbla_real(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(*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_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); jy += *incy; if (j > *ku) { kx += *incx; diff --git a/frame/compat/f2c/bla_hbmv.c b/frame/compat/f2c/bla_hbmv.c index 125f058ad..33ee628fb 100644 --- a/frame/compat/f2c/bla_hbmv.c +++ b/frame/compat/f2c/bla_hbmv.c @@ -671,7 +671,7 @@ /* 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.))) { return 0; } @@ -694,9 +694,9 @@ /* 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 (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { + if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) { i__1 = *n; for (i__ = 1; i__ <= i__1; ++i__) { i__2 = i__; @@ -708,14 +708,14 @@ for (i__ = 1; i__ <= i__1; ++i__) { i__2 = 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_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); + 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_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); /* L20: */ } } } else { iy = ky; - if (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { + if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) { i__1 = *n; for (i__ = 1; i__ <= i__1; ++i__) { i__2 = iy; @@ -728,15 +728,15 @@ for (i__ = 1; i__ <= i__1; ++i__) { i__2 = 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_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); + 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_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); iy += *incy; /* L40: */ } } } } - if (bli_zbla_real(*alpha) == 0. && bli_zimag(*alpha) == 0.) { + if (bli_zreal(*alpha) == 0. && bli_zimag(*alpha) == 0.) { return 0; } if (PASTEF770(lsame)(uplo, "U", (ftnlen)1, (ftnlen)1)) { @@ -748,8 +748,8 @@ i__1 = *n; for (j = 1; j <= i__1; ++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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (0.), (0.), temp2 ); l = kplus1 - j; /* Computing MAX */ @@ -759,25 +759,25 @@ i__2 = i__; i__3 = i__; 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_zbla_real(y[i__3]) + bli_zbla_real(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(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_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp2) + bli_zbla_real(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__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_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 ); /* L50: */ } i__4 = j; i__2 = j; i__3 = kplus1 + j * a_dim1; - d__1 = bli_zbla_real(a[i__3]); - bli_zsets( (d__1 * bli_zbla_real(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_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_zbla_real(z__2) + bli_zbla_real(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] ); + d__1 = bli_zreal(a[i__3]); + bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__3 ); + 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_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_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] ); /* L60: */ } } else { @@ -786,8 +786,8 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (0.), (0.), temp2 ); ix = kx; iy = ky; @@ -799,14 +799,14 @@ i__4 = iy; i__2 = iy; 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_zbla_real(y[i__2]) + bli_zbla_real(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(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_zreal(y[i__2]) + bli_zreal(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] ); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp2) + bli_zbla_real(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__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_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 ); ix += *incx; iy += *incy; /* L70: */ @@ -814,12 +814,12 @@ i__3 = jy; i__4 = jy; i__2 = kplus1 + j * a_dim1; - d__1 = bli_zbla_real(a[i__2]); - bli_zsets( (d__1 * bli_zbla_real(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_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_zbla_real(z__2) + bli_zbla_real(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] ); + d__1 = bli_zreal(a[i__2]); + bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__3 ); + 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_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_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] ); jx += *incx; jy += *incy; if (j > *k) { @@ -837,16 +837,16 @@ i__1 = *n; for (j = 1; j <= i__1; ++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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (0.), (0.), temp2 ); i__3 = j; i__4 = j; i__2 = j * a_dim1 + 1; - d__1 = bli_zbla_real(a[i__2]); - bli_zsets( (d__1 * bli_zbla_real(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_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); + d__1 = bli_zreal(a[i__2]); + bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__2 ); + 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_zreal(z__1)), (bli_zimag(z__1)), y[i__3] ); l = 1 - j; /* Computing MIN */ i__4 = *n, i__2 = j + *k; @@ -855,21 +855,21 @@ i__4 = i__; i__2 = i__; 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_zbla_real(y[i__2]) + bli_zbla_real(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(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_zreal(y[i__2]) + bli_zreal(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] ); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp2) + bli_zbla_real(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__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_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 ); /* L90: */ } i__3 = 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_zbla_real(y[i__4]) + bli_zbla_real(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(*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_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] ); /* L100: */ } } else { @@ -878,16 +878,16 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (0.), (0.), temp2 ); i__3 = jy; i__4 = jy; i__2 = j * a_dim1 + 1; - d__1 = bli_zbla_real(a[i__2]); - bli_zsets( (d__1 * bli_zbla_real(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_zbla_real(z__1)), (bli_zimag(z__1)), y[i__3] ); + d__1 = bli_zreal(a[i__2]); + bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__2 ); + 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_zreal(z__1)), (bli_zimag(z__1)), y[i__3] ); l = 1 - j; ix = jx; iy = jy; @@ -900,21 +900,21 @@ i__4 = iy; i__2 = iy; 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_zbla_real(y[i__2]) + bli_zbla_real(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(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_zreal(y[i__2]) + bli_zreal(z__2)), (bli_zimag(y[i__2]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__4] ); bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp2) + bli_zbla_real(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__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_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 ); /* L110: */ } i__3 = 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_zbla_real(y[i__4]) + bli_zbla_real(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(*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_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] ); jx += *incx; jy += *incy; /* L120: */ diff --git a/frame/compat/f2c/bla_hpmv.c b/frame/compat/f2c/bla_hpmv.c index 3fc294f9c..e43508fd6 100644 --- a/frame/compat/f2c/bla_hpmv.c +++ b/frame/compat/f2c/bla_hpmv.c @@ -583,7 +583,7 @@ /* 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.))) { return 0; } @@ -606,9 +606,9 @@ /* 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 (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { + if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) { i__1 = *n; for (i__ = 1; i__ <= i__1; ++i__) { i__2 = i__; @@ -620,14 +620,14 @@ for (i__ = 1; i__ <= i__1; ++i__) { i__2 = 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_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); + 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_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); /* L20: */ } } } else { iy = ky; - if (bli_zbla_real(*beta) == 0. && bli_zimag(*beta) == 0.) { + if (bli_zreal(*beta) == 0. && bli_zimag(*beta) == 0.) { i__1 = *n; for (i__ = 1; i__ <= i__1; ++i__) { i__2 = iy; @@ -640,15 +640,15 @@ for (i__ = 1; i__ <= i__1; ++i__) { i__2 = 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_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); + 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_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); iy += *incy; /* L40: */ } } } } - if (bli_zbla_real(*alpha) == 0. && bli_zimag(*alpha) == 0.) { + if (bli_zreal(*alpha) == 0. && bli_zimag(*alpha) == 0.) { return 0; } kk = 1; @@ -660,8 +660,8 @@ i__1 = *n; for (j = 1; j <= i__1; ++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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (0.), (0.), temp2 ); k = kk; i__2 = j - 1; @@ -669,26 +669,26 @@ i__3 = i__; i__4 = i__; 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_zbla_real(y[i__4]) + bli_zbla_real(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(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_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] ); bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp2) + bli_zbla_real(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__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_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 ); ++k; /* L50: */ } i__2 = j; i__3 = j; i__4 = kk + j - 1; - d__1 = bli_zbla_real(ap[i__4]); - bli_zsets( (d__1 * bli_zbla_real(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_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_zbla_real(z__2) + bli_zbla_real(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] ); + d__1 = bli_zreal(ap[i__4]); + bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__3 ); + 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_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_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); kk += j; /* L60: */ } @@ -698,8 +698,8 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (0.), (0.), temp2 ); ix = kx; iy = ky; @@ -708,14 +708,14 @@ i__3 = iy; i__4 = iy; 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_zbla_real(y[i__4]) + bli_zbla_real(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(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_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] ); bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp2) + bli_zbla_real(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__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_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 ); ix += *incx; iy += *incy; /* L70: */ @@ -723,12 +723,12 @@ i__2 = jy; i__3 = jy; i__4 = kk + j - 1; - d__1 = bli_zbla_real(ap[i__4]); - bli_zsets( (d__1 * bli_zbla_real(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_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_zbla_real(z__2) + bli_zbla_real(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] ); + d__1 = bli_zreal(ap[i__4]); + bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__3 ); + 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_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_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); jx += *incx; jy += *incy; kk += j; @@ -743,38 +743,38 @@ i__1 = *n; for (j = 1; j <= i__1; ++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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (0.), (0.), temp2 ); i__2 = j; i__3 = j; i__4 = kk; - d__1 = bli_zbla_real(ap[i__4]); - bli_zsets( (d__1 * bli_zbla_real(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_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); + d__1 = bli_zreal(ap[i__4]); + bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__2 ); + 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_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); k = kk + 1; i__2 = *n; for (i__ = j + 1; i__ <= i__2; ++i__) { i__3 = i__; i__4 = i__; 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_zbla_real(y[i__4]) + bli_zbla_real(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(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_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] ); bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp2) + bli_zbla_real(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__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_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 ); ++k; /* L90: */ } i__2 = 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_zbla_real(y[i__3]) + bli_zbla_real(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(*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_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); kk += *n - j + 1; /* L100: */ } @@ -784,16 +784,16 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); bli_zsets( (0.), (0.), temp2 ); i__2 = jy; i__3 = jy; i__4 = kk; - d__1 = bli_zbla_real(ap[i__4]); - bli_zsets( (d__1 * bli_zbla_real(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_zbla_real(z__1)), (bli_zimag(z__1)), y[i__2] ); + d__1 = bli_zreal(ap[i__4]); + bli_zsets( (d__1 * bli_zreal(temp1)), (d__1 * bli_zimag(temp1)), z__2 ); + 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_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); ix = jx; iy = jy; i__2 = kk + *n - j; @@ -803,21 +803,21 @@ i__3 = iy; i__4 = iy; 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_zbla_real(y[i__4]) + bli_zbla_real(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(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_zreal(y[i__4]) + bli_zreal(z__2)), (bli_zimag(y[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__3] ); bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp2) + bli_zbla_real(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__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_zreal(temp2) + bli_zreal(z__2)), (bli_zimag(temp2) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp2 ); /* L110: */ } i__2 = 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_zbla_real(y[i__3]) + bli_zbla_real(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(*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_zreal(y[i__3]) + bli_zreal(z__2)), (bli_zimag(y[i__3]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), y[i__2] ); jx += *incx; jy += *incy; kk += *n - j + 1; diff --git a/frame/compat/f2c/bla_hpr.c b/frame/compat/f2c/bla_hpr.c index 2e4148b54..7d110170c 100644 --- a/frame/compat/f2c/bla_hpr.c +++ b/frame/compat/f2c/bla_hpr.c @@ -506,32 +506,32 @@ i__1 = *n; for (j = 1; j <= i__1; ++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]); - bli_zsets( (*alpha * bli_zbla_real(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); - bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + bli_zsets( (*alpha * bli_zreal(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); k = kk; i__2 = j - 1; for (i__ = 1; i__ <= i__2; ++i__) { i__3 = k; i__4 = k; 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_zbla_real(ap[i__4]) + bli_zbla_real(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(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_zreal(ap[i__4]) + bli_zreal(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] ); ++k; /* L10: */ } i__2 = kk + j - 1; i__3 = kk + j - 1; 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 ); - d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(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_zreal(ap[i__3]) + bli_zreal(z__1); bli_zsets( (d__1), (0.), ap[i__2] ); } else { i__2 = 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] ); } kk += j; @@ -542,32 +542,32 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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]); - bli_zsets( (*alpha * bli_zbla_real(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); - bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + bli_zsets( (*alpha * bli_zreal(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix = kx; i__2 = kk + j - 2; for (k = kk; k <= i__2; ++k) { i__3 = k; i__4 = k; 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_zbla_real(ap[i__4]) + bli_zbla_real(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(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_zreal(ap[i__4]) + bli_zreal(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] ); ix += *incx; /* L30: */ } i__2 = kk + j - 1; i__3 = kk + j - 1; 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 ); - d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(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_zreal(ap[i__3]) + bli_zreal(z__1); bli_zsets( (d__1), (0.), ap[i__2] ); } else { i__2 = 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] ); } jx += *incx; @@ -583,15 +583,15 @@ i__1 = *n; for (j = 1; j <= i__1; ++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]); - bli_zsets( (*alpha * bli_zbla_real(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); - bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + bli_zsets( (*alpha * bli_zreal(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); i__2 = kk; i__3 = kk; 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 ); - d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(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_zreal(ap[i__3]) + bli_zreal(z__1); bli_zsets( (d__1), (0.), ap[i__2] ); k = kk + 1; i__2 = *n; @@ -599,16 +599,16 @@ i__3 = k; i__4 = k; 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_zbla_real(ap[i__4]) + bli_zbla_real(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(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_zreal(ap[i__4]) + bli_zreal(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] ); ++k; /* L50: */ } } else { i__2 = 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] ); } kk = kk + *n - j + 1; @@ -619,15 +619,15 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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]); - bli_zsets( (*alpha * bli_zbla_real(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); - bli_zsets( (bli_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + bli_zsets( (*alpha * bli_zreal(z__2)), (*alpha * bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); i__2 = kk; i__3 = kk; 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 ); - d__1 = bli_zbla_real(ap[i__3]) + bli_zbla_real(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_zreal(ap[i__3]) + bli_zreal(z__1); bli_zsets( (d__1), (0.), ap[i__2] ); ix = jx; i__2 = kk + *n - j; @@ -636,15 +636,15 @@ i__3 = k; i__4 = k; 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_zbla_real(ap[i__4]) + bli_zbla_real(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(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_zreal(ap[i__4]) + bli_zreal(z__2)), (bli_zimag(ap[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] ); /* L70: */ } } else { i__2 = 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] ); } jx += *incx; diff --git a/frame/compat/f2c/bla_hpr2.c b/frame/compat/f2c/bla_hpr2.c index 104b67703..c1563fc56 100644 --- a/frame/compat/f2c/bla_hpr2.c +++ b/frame/compat/f2c/bla_hpr2.c @@ -572,7 +572,7 @@ /* 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; } @@ -607,43 +607,43 @@ for (j = 1; j <= i__1; ++j) { i__2 = 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.)) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); 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); - 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; i__2 = j - 1; for (i__ = 1; i__ <= i__2; ++i__) { i__3 = k; i__4 = k; 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_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(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_zreal(ap[i__4]) + bli_zreal(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 ); 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_zbla_real(z__2) + bli_zbla_real(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(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_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] ); ++k; /* L10: */ } i__2 = kk + j - 1; i__3 = kk + j - 1; 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; - 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_zbla_real(z__2) + bli_zbla_real(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); + 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_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); + d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1); bli_zsets( (d__1), (0.), ap[i__2] ); } else { i__2 = 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] ); } kk += j; @@ -654,15 +654,15 @@ for (j = 1; j <= i__1; ++j) { i__2 = jx; 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.)) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); 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); - 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; iy = ky; i__2 = kk + j - 2; @@ -670,12 +670,12 @@ i__3 = k; i__4 = k; 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_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(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_zreal(ap[i__4]) + bli_zreal(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 ); 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_zbla_real(z__2) + bli_zbla_real(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(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_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] ); ix += *incx; iy += *incy; /* L30: */ @@ -683,16 +683,16 @@ i__2 = kk + j - 1; i__3 = kk + j - 1; 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; - 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_zbla_real(z__2) + bli_zbla_real(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); + 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_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); + d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1); bli_zsets( (d__1), (0.), ap[i__2] ); } else { i__2 = 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] ); } jx += *incx; @@ -710,23 +710,23 @@ for (j = 1; j <= i__1; ++j) { i__2 = 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.)) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); 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); - 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__3 = kk; 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; - 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_zbla_real(z__2) + bli_zbla_real(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); + 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_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); + d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1); bli_zsets( (d__1), (0.), ap[i__2] ); k = kk + 1; i__2 = *n; @@ -734,19 +734,19 @@ i__3 = k; i__4 = k; 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_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(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_zreal(ap[i__4]) + bli_zreal(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 ); 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_zbla_real(z__2) + bli_zbla_real(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(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_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] ); ++k; /* L50: */ } } else { i__2 = 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] ); } kk = kk + *n - j + 1; @@ -757,23 +757,23 @@ for (j = 1; j <= i__1; ++j) { i__2 = jx; 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.)) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp1 ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp1 ); 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); - 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__3 = kk; 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; - 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_zbla_real(z__2) + bli_zbla_real(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); + 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_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); + d__1 = bli_zreal(ap[i__3]) + bli_zreal(z__1); bli_zsets( (d__1), (0.), ap[i__2] ); ix = jx; iy = jy; @@ -784,18 +784,18 @@ i__3 = k; i__4 = k; 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_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(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_zreal(ap[i__4]) + bli_zreal(z__3)), (bli_zimag(ap[i__4]) + bli_zimag(z__3)), z__2 ); 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_zbla_real(z__2) + bli_zbla_real(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(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_zreal(z__2) + bli_zreal(z__4)), (bli_zimag(z__2) + bli_zimag(z__4)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ap[i__3] ); /* L70: */ } } else { i__2 = 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] ); } jx += *incx; diff --git a/frame/compat/f2c/bla_rot.c b/frame/compat/f2c/bla_rot.c index c81619a96..18618ee5e 100644 --- a/frame/compat/f2c/bla_rot.c +++ b/frame/compat/f2c/bla_rot.c @@ -313,20 +313,20 @@ L20: i__1 = *n; for (i__ = 1; i__ <= i__1; ++i__) { 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; - bli_zsets( (*s * bli_zbla_real(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_zbla_real(z__1)), (bli_zimag(z__1)), ztemp ); + bli_zsets( (*s * bli_zreal(zy[i__3])), (*s * bli_zimag(zy[i__3])), z__3 ); + bli_zsets( (bli_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ztemp ); i__2 = 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; - bli_zsets( (*s * bli_zbla_real(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_zbla_real(z__1)), (bli_zimag(z__1)), zy[i__2] ); + bli_zsets( (*s * bli_zreal(zx[i__4])), (*s * bli_zimag(zx[i__4])), z__3 ); + bli_zsets( (bli_zreal(z__2) - bli_zreal(z__3)), (bli_zimag(z__2) - bli_zimag(z__3)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), zy[i__2] ); 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; iy += *incy; /* L10: */ @@ -339,20 +339,20 @@ L20: i__1 = *n; for (i__ = 1; i__ <= i__1; ++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__; - bli_zsets( (*s * bli_zbla_real(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_zbla_real(z__1)), (bli_zimag(z__1)), ztemp ); + bli_zsets( (*s * bli_zreal(zy[i__3])), (*s * bli_zimag(zy[i__3])), z__3 ); + bli_zsets( (bli_zreal(z__2) + bli_zreal(z__3)), (bli_zimag(z__2) + bli_zimag(z__3)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), ztemp ); i__2 = 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__; - bli_zsets( (*s * bli_zbla_real(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_zbla_real(z__1)), (bli_zimag(z__1)), zy[i__2] ); + bli_zsets( (*s * bli_zreal(zx[i__4])), (*s * bli_zimag(zx[i__4])), z__3 ); + bli_zsets( (bli_zreal(z__2) - bli_zreal(z__3)), (bli_zimag(z__2) - bli_zimag(z__3)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), zy[i__2] ); 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: */ } return 0; diff --git a/frame/compat/f2c/bla_rotg.c b/frame/compat/f2c/bla_rotg.c index de0db5c48..490b01a59 100644 --- a/frame/compat/f2c/bla_rotg.c +++ b/frame/compat/f2c/bla_rotg.c @@ -233,7 +233,7 @@ L20: } *c__ = 0.; 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; L10: scale = bla_z_abs(ca) + bla_z_abs(cb); @@ -247,15 +247,15 @@ L10: d__2 = bla_z_abs(&z__3); norm = scale * sqrt(d__1 * d__1 + d__2 * d__2); d__1 = bla_z_abs(ca); - bli_zsets( (bli_zbla_real(*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(*ca) / d__1), (bli_zimag(*ca) / d__1), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), alpha ); *c__ = bla_z_abs(ca) / norm; 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_zbla_real(z__2) / norm), (bli_zimag(z__2) / norm), z__1 ); - bli_zsets( bli_zbla_real(z__1), bli_zimag(z__1), *s ); - bli_zsets( (norm * bli_zbla_real(alpha)), (norm * bli_zimag(alpha)), z__1 ); - bli_zsets( bli_zbla_real(z__1), bli_zimag(z__1), *ca ); + 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_zreal(z__2) / norm), (bli_zimag(z__2) / norm), z__1 ); + bli_zsets( bli_zreal(z__1), bli_zimag(z__1), *s ); + bli_zsets( (norm * bli_zreal(alpha)), (norm * bli_zimag(alpha)), z__1 ); + bli_zsets( bli_zreal(z__1), bli_zimag(z__1), *ca ); L20: return 0; } /* zrotg_ */ diff --git a/frame/compat/f2c/bla_tbmv.c b/frame/compat/f2c/bla_tbmv.c index d7b09b697..82d963348 100644 --- a/frame/compat/f2c/bla_tbmv.c +++ b/frame/compat/f2c/bla_tbmv.c @@ -1656,9 +1656,9 @@ i__1 = *n; for (j = 1; j <= i__1; ++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; - 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; /* Computing MAX */ i__2 = 1, i__3 = j - *k; @@ -1667,17 +1667,17 @@ i__2 = i__; i__3 = i__; 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_zbla_real(x[i__3]) + bli_zbla_real(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(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_zreal(x[i__3]) + bli_zreal(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] ); /* L10: */ } if (nounit) { i__4 = j; i__2 = j; 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_zbla_real(z__1)), (bli_zimag(z__1)), x[i__4] ); + 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_zreal(z__1)), (bli_zimag(z__1)), x[i__4] ); } } /* L20: */ @@ -1687,9 +1687,9 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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; - 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; l = kplus1 - j; /* Computing MAX */ @@ -1699,9 +1699,9 @@ i__4 = ix; i__2 = ix; 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_zbla_real(x[i__2]) + bli_zbla_real(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(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_zreal(x[i__2]) + bli_zreal(z__2)), (bli_zimag(x[i__2]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__4] ); ix += *incx; /* L30: */ } @@ -1709,8 +1709,8 @@ i__3 = jx; i__4 = jx; 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_zbla_real(z__1)), (bli_zimag(z__1)), x[i__3] ); + 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_zreal(z__1)), (bli_zimag(z__1)), x[i__3] ); } } jx += *incx; @@ -1724,9 +1724,9 @@ if (*incx == 1) { for (j = *n; j >= 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; - 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; /* Computing MIN */ i__1 = *n, i__3 = j + *k; @@ -1735,17 +1735,17 @@ i__1 = i__; i__3 = i__; 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_zbla_real(x[i__3]) + bli_zbla_real(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(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_zreal(x[i__3]) + bli_zreal(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__1] ); /* L50: */ } if (nounit) { i__4 = j; i__1 = j; 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_zbla_real(z__1)), (bli_zimag(z__1)), x[i__4] ); + 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_zreal(z__1)), (bli_zimag(z__1)), x[i__4] ); } } /* L60: */ @@ -1755,9 +1755,9 @@ jx = kx; for (j = *n; j >= 1; --j) { 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; - 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; l = 1 - j; /* Computing MIN */ @@ -1767,9 +1767,9 @@ i__4 = ix; i__1 = ix; 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_zbla_real(x[i__1]) + bli_zbla_real(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(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_zreal(x[i__1]) + bli_zreal(z__2)), (bli_zimag(x[i__1]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__4] ); ix -= *incx; /* L70: */ } @@ -1777,8 +1777,8 @@ i__3 = jx; i__4 = jx; 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_zbla_real(z__1)), (bli_zimag(z__1)), x[i__3] ); + 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_zreal(z__1)), (bli_zimag(z__1)), x[i__3] ); } } jx -= *incx; @@ -1798,13 +1798,13 @@ if (*incx == 1) { for (j = *n; j >= 1; --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; if (noconj) { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } /* Computing MAX */ i__4 = 1, i__1 = j - *k; @@ -1812,16 +1812,16 @@ for (i__ = j - 1; i__ >= i__3; --i__) { i__4 = l + i__ + j * a_dim1; 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_zbla_real(temp) + bli_zbla_real(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(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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L90: */ } } else { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } /* Computing MAX */ i__4 = 1, i__1 = j - *k; @@ -1829,14 +1829,14 @@ for (i__ = j - 1; i__ >= i__3; --i__) { bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp) + bli_zbla_real(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__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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L100: */ } } 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: */ } } else { @@ -1844,15 +1844,15 @@ jx = kx; for (j = *n; j >= 1; --j) { 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; ix = kx; l = kplus1 - j; if (noconj) { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } /* Computing MAX */ i__4 = 1, i__1 = j - *k; @@ -1860,17 +1860,17 @@ for (i__ = j - 1; i__ >= i__3; --i__) { i__4 = l + i__ + j * a_dim1; 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_zbla_real(temp) + bli_zbla_real(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(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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix -= *incx; /* L120: */ } } else { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } /* Computing MAX */ i__4 = 1, i__1 = j - *k; @@ -1878,15 +1878,15 @@ for (i__ = j - 1; i__ >= i__3; --i__) { bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp) + bli_zbla_real(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__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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix -= *incx; /* L130: */ } } 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; /* L140: */ } @@ -1896,13 +1896,13 @@ i__3 = *n; for (j = 1; j <= i__3; ++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; if (noconj) { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } /* Computing MIN */ i__1 = *n, i__2 = j + *k; @@ -1910,16 +1910,16 @@ for (i__ = j + 1; i__ <= i__4; ++i__) { i__1 = l + i__ + j * a_dim1; 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_zbla_real(temp) + bli_zbla_real(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(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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L150: */ } } else { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } /* Computing MIN */ i__1 = *n, i__2 = j + *k; @@ -1927,14 +1927,14 @@ for (i__ = j + 1; i__ <= i__4; ++i__) { bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp) + bli_zbla_real(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__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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L160: */ } } 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: */ } } else { @@ -1942,15 +1942,15 @@ i__3 = *n; for (j = 1; j <= i__3; ++j) { 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; ix = kx; l = 1 - j; if (noconj) { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } /* Computing MIN */ i__1 = *n, i__2 = j + *k; @@ -1958,17 +1958,17 @@ for (i__ = j + 1; i__ <= i__4; ++i__) { i__1 = l + i__ + j * a_dim1; 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_zbla_real(temp) + bli_zbla_real(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(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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix += *incx; /* L180: */ } } else { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } /* Computing MIN */ i__1 = *n, i__2 = j + *k; @@ -1976,15 +1976,15 @@ for (i__ = j + 1; i__ <= i__4; ++i__) { bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp) + bli_zbla_real(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__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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix += *incx; /* L190: */ } } 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; /* L200: */ } diff --git a/frame/compat/f2c/bla_tbsv.c b/frame/compat/f2c/bla_tbsv.c index e3b2541a2..848c19aaf 100644 --- a/frame/compat/f2c/bla_tbsv.c +++ b/frame/compat/f2c/bla_tbsv.c @@ -1660,15 +1660,15 @@ if (*incx == 1) { for (j = *n; j >= 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; if (nounit) { i__1 = j; 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; - 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 */ i__2 = 1, i__3 = j - *k; i__1 = f2c_max(i__2,i__3); @@ -1676,9 +1676,9 @@ i__2 = i__; i__3 = i__; 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_zbla_real(x[i__3]) - bli_zbla_real(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(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_zreal(x[i__3]) - bli_zreal(z__2)), (bli_zimag(x[i__3]) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] ); /* L10: */ } } @@ -1690,16 +1690,16 @@ for (j = *n; j >= 1; --j) { kx -= *incx; 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; l = kplus1 - j; if (nounit) { i__1 = jx; 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; - 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 */ i__2 = 1, i__3 = j - *k; i__1 = f2c_max(i__2,i__3); @@ -1707,9 +1707,9 @@ i__2 = ix; i__3 = ix; 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_zbla_real(x[i__3]) - bli_zbla_real(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(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_zreal(x[i__3]) - bli_zreal(z__2)), (bli_zimag(x[i__3]) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] ); ix -= *incx; /* L30: */ } @@ -1723,15 +1723,15 @@ i__1 = *n; for (j = 1; j <= i__1; ++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; if (nounit) { i__2 = j; 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; - 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 */ i__3 = *n, i__4 = j + *k; i__2 = f2c_min(i__3,i__4); @@ -1739,9 +1739,9 @@ i__3 = i__; i__4 = i__; 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_zbla_real(x[i__4]) - bli_zbla_real(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(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_zreal(x[i__4]) - bli_zreal(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] ); /* L50: */ } } @@ -1753,16 +1753,16 @@ for (j = 1; j <= i__1; ++j) { kx += *incx; 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; l = 1 - j; if (nounit) { i__2 = jx; 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; - 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 */ i__3 = *n, i__4 = j + *k; i__2 = f2c_min(i__3,i__4); @@ -1770,9 +1770,9 @@ i__3 = ix; i__4 = ix; 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_zbla_real(x[i__4]) - bli_zbla_real(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(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_zreal(x[i__4]) - bli_zreal(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] ); ix += *incx; /* L70: */ } @@ -1792,7 +1792,7 @@ i__1 = *n; for (j = 1; j <= i__1; ++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; if (noconj) { /* Computing MAX */ @@ -1801,14 +1801,14 @@ for (i__ = f2c_max(i__2,i__3); i__ <= i__4; ++i__) { i__2 = l + i__ + j * a_dim1; 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_zbla_real(temp) - bli_zbla_real(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(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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L90: */ } if (nounit) { 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 { /* Computing MAX */ @@ -1817,19 +1817,19 @@ for (i__ = f2c_max(i__4,i__2); i__ <= i__3; ++i__) { bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp) - bli_zbla_real(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__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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L100: */ } if (nounit) { bla_d_cnjg(&z__2, &a[kplus1 + j * a_dim1]); 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; - bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__3] ); + bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__3] ); /* L110: */ } } else { @@ -1837,7 +1837,7 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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; l = kplus1 - j; if (noconj) { @@ -1847,15 +1847,15 @@ for (i__ = f2c_max(i__3,i__4); i__ <= i__2; ++i__) { i__3 = l + i__ + j * a_dim1; 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_zbla_real(temp) - bli_zbla_real(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(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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix += *incx; /* L120: */ } if (nounit) { 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 { /* Computing MAX */ @@ -1864,20 +1864,20 @@ for (i__ = f2c_max(i__2,i__3); i__ <= i__4; ++i__) { bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp) - bli_zbla_real(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__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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix += *incx; /* L130: */ } if (nounit) { bla_d_cnjg(&z__2, &a[kplus1 + j * a_dim1]); 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; - 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; if (j > *k) { kx += *incx; @@ -1889,7 +1889,7 @@ if (*incx == 1) { for (j = *n; j >= 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; if (noconj) { /* Computing MIN */ @@ -1898,14 +1898,14 @@ for (i__ = f2c_min(i__1,i__4); i__ >= i__2; --i__) { i__1 = l + i__ + j * a_dim1; 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_zbla_real(temp) - bli_zbla_real(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(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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L150: */ } if (nounit) { 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 { /* Computing MIN */ @@ -1914,19 +1914,19 @@ for (i__ = f2c_min(i__2,i__1); i__ >= i__4; --i__) { bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp) - bli_zbla_real(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__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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L160: */ } if (nounit) { bla_d_cnjg(&z__2, &a[j * a_dim1 + 1]); 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; - bli_zsets( (bli_zbla_real(temp)), (bli_zimag(temp)), x[i__4] ); + bli_zsets( (bli_zreal(temp)), (bli_zimag(temp)), x[i__4] ); /* L170: */ } } else { @@ -1934,7 +1934,7 @@ jx = kx; for (j = *n; j >= 1; --j) { 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; l = 1 - j; if (noconj) { @@ -1944,15 +1944,15 @@ for (i__ = f2c_min(i__4,i__2); i__ >= i__1; --i__) { i__4 = l + i__ + j * a_dim1; 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_zbla_real(temp) - bli_zbla_real(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(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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix -= *incx; /* L180: */ } if (nounit) { 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 { /* Computing MIN */ @@ -1961,20 +1961,20 @@ for (i__ = f2c_min(i__1,i__4); i__ >= i__2; --i__) { bla_d_cnjg(&z__3, &a[l + i__ + j * a_dim1]); 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_zbla_real(temp) - bli_zbla_real(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__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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix -= *incx; /* L190: */ } if (nounit) { bla_d_cnjg(&z__2, &a[j * a_dim1 + 1]); 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; - 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; if (*n - j >= *k) { kx -= *incx; diff --git a/frame/compat/f2c/bla_tpmv.c b/frame/compat/f2c/bla_tpmv.c index cc897d1af..5cafe7941 100644 --- a/frame/compat/f2c/bla_tpmv.c +++ b/frame/compat/f2c/bla_tpmv.c @@ -1415,18 +1415,18 @@ i__1 = *n; for (j = 1; j <= i__1; ++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; - 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; i__2 = j - 1; for (i__ = 1; i__ <= i__2; ++i__) { i__3 = i__; i__4 = i__; 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_zbla_real(x[i__4]) + bli_zbla_real(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(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_zreal(x[i__4]) + bli_zreal(z__2)), (bli_zimag(x[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] ); ++k; /* L10: */ } @@ -1434,8 +1434,8 @@ i__2 = j; i__3 = j; 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_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); + 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_zreal(z__1)), (bli_zimag(z__1)), x[i__2] ); } } kk += j; @@ -1446,18 +1446,18 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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; - 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; i__2 = kk + j - 2; for (k = kk; k <= i__2; ++k) { i__3 = ix; i__4 = ix; 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_zbla_real(x[i__4]) + bli_zbla_real(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(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_zreal(x[i__4]) + bli_zreal(z__2)), (bli_zimag(x[i__4]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] ); ix += *incx; /* L30: */ } @@ -1465,8 +1465,8 @@ i__2 = jx; i__3 = jx; 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_zbla_real(z__1)), (bli_zimag(z__1)), x[i__2] ); + 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_zreal(z__1)), (bli_zimag(z__1)), x[i__2] ); } } jx += *incx; @@ -1479,18 +1479,18 @@ if (*incx == 1) { for (j = *n; j >= 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; - 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; i__1 = j + 1; for (i__ = *n; i__ >= i__1; --i__) { i__2 = i__; i__3 = i__; 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_zbla_real(x[i__3]) + bli_zbla_real(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(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_zreal(x[i__3]) + bli_zreal(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] ); --k; /* L50: */ } @@ -1498,8 +1498,8 @@ i__1 = j; i__2 = 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_zbla_real(z__1)), (bli_zimag(z__1)), x[i__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_zreal(z__1)), (bli_zimag(z__1)), x[i__1] ); } } kk -= *n - j + 1; @@ -1510,18 +1510,18 @@ jx = kx; for (j = *n; j >= 1; --j) { 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; - 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; i__1 = kk - (*n - (j + 1)); for (k = kk; k >= i__1; --k) { i__2 = ix; i__3 = ix; 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_zbla_real(x[i__3]) + bli_zbla_real(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(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_zreal(x[i__3]) + bli_zreal(z__2)), (bli_zimag(x[i__3]) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] ); ix -= *incx; /* L70: */ } @@ -1529,8 +1529,8 @@ i__1 = jx; i__2 = jx; 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_zbla_real(z__1)), (bli_zimag(z__1)), x[i__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_zreal(z__1)), (bli_zimag(z__1)), x[i__1] ); } } jx -= *incx; @@ -1548,41 +1548,41 @@ if (*incx == 1) { for (j = *n; j >= 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; if (noconj) { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } for (i__ = j - 1; i__ >= 1; --i__) { i__1 = k; 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_zbla_real(temp) + bli_zbla_real(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(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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); --k; /* L90: */ } } else { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } for (i__ = j - 1; i__ >= 1; --i__) { bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp) + bli_zbla_real(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__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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); --k; /* L100: */ } } 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; /* L110: */ } @@ -1590,43 +1590,43 @@ jx = kx + (*n - 1) * *incx; for (j = *n; j >= 1; --j) { 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; if (noconj) { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } i__1 = kk - j + 1; for (k = kk - 1; k >= i__1; --k) { ix -= *incx; i__2 = k; 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_zbla_real(temp) + bli_zbla_real(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(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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L120: */ } } else { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } i__1 = kk - j + 1; for (k = kk - 1; k >= i__1; --k) { ix -= *incx; bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp) + bli_zbla_real(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__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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L130: */ } } 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; kk -= j; /* L140: */ @@ -1638,43 +1638,43 @@ i__1 = *n; for (j = 1; j <= i__1; ++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; if (noconj) { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } i__2 = *n; for (i__ = j + 1; i__ <= i__2; ++i__) { i__3 = k; 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_zbla_real(temp) + bli_zbla_real(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(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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ++k; /* L150: */ } } else { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } i__2 = *n; for (i__ = j + 1; i__ <= i__2; ++i__) { bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp) + bli_zbla_real(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__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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ++k; /* L160: */ } } 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; /* L170: */ } @@ -1683,43 +1683,43 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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; if (noconj) { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } i__2 = kk + *n - j; for (k = kk + 1; k <= i__2; ++k) { ix += *incx; i__3 = k; 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_zbla_real(temp) + bli_zbla_real(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(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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L180: */ } } else { if (nounit) { 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_zbla_real(z__1)), (bli_zimag(z__1)), temp ); + 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_zreal(z__1)), (bli_zimag(z__1)), temp ); } i__2 = kk + *n - j; for (k = kk + 1; k <= i__2; ++k) { ix += *incx; bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp) + bli_zbla_real(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__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_zreal(temp) + bli_zreal(z__2)), (bli_zimag(temp) + bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); /* L190: */ } } 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; kk += *n - j + 1; /* L200: */ diff --git a/frame/compat/f2c/bla_tpsv.c b/frame/compat/f2c/bla_tpsv.c index 7421e6631..2ff23575d 100644 --- a/frame/compat/f2c/bla_tpsv.c +++ b/frame/compat/f2c/bla_tpsv.c @@ -1416,22 +1416,22 @@ if (*incx == 1) { for (j = *n; j >= 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) { i__1 = j; 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; - 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; for (i__ = j - 1; i__ >= 1; --i__) { i__1 = i__; i__2 = i__; 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_zbla_real(x[i__2]) - bli_zbla_real(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(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_zreal(x[i__2]) - bli_zreal(z__2)), (bli_zimag(x[i__2]) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__1] ); --k; /* L10: */ } @@ -1443,14 +1443,14 @@ jx = kx + (*n - 1) * *incx; for (j = *n; j >= 1; --j) { 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) { i__1 = jx; 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; - 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; i__1 = kk - j + 1; for (k = kk - 1; k >= i__1; --k) { @@ -1458,9 +1458,9 @@ i__2 = ix; i__3 = ix; 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_zbla_real(x[i__3]) - bli_zbla_real(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(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_zreal(x[i__3]) - bli_zreal(z__2)), (bli_zimag(x[i__3]) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__2] ); /* L30: */ } } @@ -1475,23 +1475,23 @@ i__1 = *n; for (j = 1; j <= i__1; ++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) { i__2 = j; 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; - 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; i__2 = *n; for (i__ = j + 1; i__ <= i__2; ++i__) { i__3 = i__; i__4 = i__; 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_zbla_real(x[i__4]) - bli_zbla_real(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(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_zreal(x[i__4]) - bli_zreal(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] ); ++k; /* L50: */ } @@ -1504,14 +1504,14 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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) { i__2 = jx; 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; - 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; i__2 = kk + *n - j; for (k = kk + 1; k <= i__2; ++k) { @@ -1519,9 +1519,9 @@ i__3 = ix; i__4 = ix; 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_zbla_real(x[i__4]) - bli_zbla_real(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(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_zreal(x[i__4]) - bli_zreal(z__2)), (bli_zimag(x[i__4]) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), x[i__3] ); /* L70: */ } } @@ -1541,42 +1541,42 @@ i__1 = *n; for (j = 1; j <= i__1; ++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; if (noconj) { i__2 = j - 1; for (i__ = 1; i__ <= i__2; ++i__) { i__3 = k; 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_zbla_real(temp) - bli_zbla_real(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(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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ++k; /* L90: */ } if (nounit) { 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 { i__2 = j - 1; for (i__ = 1; i__ <= i__2; ++i__) { bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp) - bli_zbla_real(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__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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ++k; /* L100: */ } if (nounit) { bla_d_cnjg(&z__2, &ap[kk + j - 1]); 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; - 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; /* L110: */ } @@ -1585,42 +1585,42 @@ i__1 = *n; for (j = 1; j <= i__1; ++j) { 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; if (noconj) { i__2 = kk + j - 2; for (k = kk; k <= i__2; ++k) { i__3 = k; 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_zbla_real(temp) - bli_zbla_real(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(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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix += *incx; /* L120: */ } if (nounit) { 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 { i__2 = kk + j - 2; for (k = kk; k <= i__2; ++k) { bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp) - bli_zbla_real(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__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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix += *incx; /* L130: */ } if (nounit) { bla_d_cnjg(&z__2, &ap[kk + j - 1]); 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; - 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; kk += j; /* L140: */ @@ -1631,42 +1631,42 @@ if (*incx == 1) { for (j = *n; j >= 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; if (noconj) { i__1 = j + 1; for (i__ = *n; i__ >= i__1; --i__) { i__2 = k; 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_zbla_real(temp) - bli_zbla_real(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(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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); --k; /* L150: */ } if (nounit) { 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 { i__1 = j + 1; for (i__ = *n; i__ >= i__1; --i__) { bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp) - bli_zbla_real(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__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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); --k; /* L160: */ } if (nounit) { bla_d_cnjg(&z__2, &ap[kk - *n + j]); 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; - 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; /* L170: */ } @@ -1675,42 +1675,42 @@ jx = kx; for (j = *n; j >= 1; --j) { 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; if (noconj) { i__1 = kk - (*n - (j + 1)); for (k = kk; k >= i__1; --k) { i__2 = k; 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_zbla_real(temp) - bli_zbla_real(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(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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix -= *incx; /* L180: */ } if (nounit) { 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 { i__1 = kk - (*n - (j + 1)); for (k = kk; k >= i__1; --k) { bla_d_cnjg(&z__3, &ap[k]); 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_zbla_real(temp) - bli_zbla_real(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__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_zreal(temp) - bli_zreal(z__2)), (bli_zimag(temp) - bli_zimag(z__2)), z__1 ); + bli_zsets( (bli_zreal(z__1)), (bli_zimag(z__1)), temp ); ix -= *incx; /* L190: */ } if (nounit) { bla_d_cnjg(&z__2, &ap[kk - *n + j]); 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; - 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; kk -= *n - j + 1; /* L200: */ diff --git a/frame/compat/f2c/util/bla_z_abs.c b/frame/compat/f2c/util/bla_z_abs.c index 5a9cec965..0ad781aee 100644 --- a/frame/compat/f2c/util/bla_z_abs.c +++ b/frame/compat/f2c/util/bla_z_abs.c @@ -38,7 +38,7 @@ 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 ) ) ); }