diff --git a/lapack-netlib/SRC/cla_heamv.c b/lapack-netlib/SRC/cla_heamv.c
index 59cee37371..37907dc94c 100644
--- a/lapack-netlib/SRC/cla_heamv.c
+++ b/lapack-netlib/SRC/cla_heamv.c
@@ -627,7 +627,7 @@ mv.f"> */
i__3 = j + i__ * a_dim1;
temp = (r__1 = a[i__3].r, abs(r__1)) + (r__2 = r_imag(
&a[j + i__ * a_dim1]), abs(r__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0.f && x[
i__3].i == 0.f) || temp == 0.f);
i__3 = jx;
@@ -640,7 +640,7 @@ mv.f"> */
i__3 = i__ + j * a_dim1;
temp = (r__1 = a[i__3].r, abs(r__1)) + (r__2 = r_imag(
&a[i__ + j * a_dim1]), abs(r__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0.f && x[
i__3].i == 0.f) || temp == 0.f);
i__3 = jx;
@@ -673,7 +673,7 @@ mv.f"> */
i__3 = i__ + j * a_dim1;
temp = (r__1 = a[i__3].r, abs(r__1)) + (r__2 = r_imag(
&a[i__ + j * a_dim1]), abs(r__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0.f && x[
i__3].i == 0.f) || temp == 0.f);
i__3 = jx;
@@ -686,7 +686,7 @@ mv.f"> */
i__3 = j + i__ * a_dim1;
temp = (r__1 = a[i__3].r, abs(r__1)) + (r__2 = r_imag(
&a[j + i__ * a_dim1]), abs(r__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0.f && x[
i__3].i == 0.f) || temp == 0.f);
i__3 = jx;
diff --git a/lapack-netlib/SRC/cla_heamv.f b/lapack-netlib/SRC/cla_heamv.f
index 6230de4af7..0f9c23b177 100644
--- a/lapack-netlib/SRC/cla_heamv.f
+++ b/lapack-netlib/SRC/cla_heamv.f
@@ -5,7 +5,6 @@
* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/
*
-*> \htmlonly
*> Download CLA_HEAMV + dependencies
*>
*> [TGZ]
@@ -13,7 +12,6 @@
*> [ZIP]
*>
*> [TXT]
-*> \endhtmlonly
*
* Definition:
* ===========
@@ -154,7 +152,7 @@
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \ingroup complexHEcomputational
+*> \ingroup la_heamv
*
*> \par Further Details:
* =====================
@@ -175,6 +173,7 @@
* =====================================================================
SUBROUTINE CLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
$ INCY )
+ IMPLICIT NONE
*
* -- LAPACK computational routine --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -262,7 +261,7 @@ SUBROUTINE CLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
* number of additions in each row.
*
SAFE1 = SLAMCH( 'Safe minimum' )
- SAFE1 = (N+1)*SAFE1
+ SAFE1 = REAL( N+1 )*SAFE1
*
* Form y := alpha*abs(A)*abs(x) + beta*abs(y).
*
@@ -356,7 +355,7 @@ SUBROUTINE CLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = CABS1( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -364,7 +363,7 @@ SUBROUTINE CLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = CABS1( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -392,7 +391,7 @@ SUBROUTINE CLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = CABS1( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -400,7 +399,7 @@ SUBROUTINE CLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = CABS1( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
diff --git a/lapack-netlib/SRC/cla_syamv.c b/lapack-netlib/SRC/cla_syamv.c
index bafa49f2d9..a9d763aee4 100644
--- a/lapack-netlib/SRC/cla_syamv.c
+++ b/lapack-netlib/SRC/cla_syamv.c
@@ -628,7 +628,7 @@ mv.f"> */
i__3 = j + i__ * a_dim1;
temp = (r__1 = a[i__3].r, abs(r__1)) + (r__2 = r_imag(
&a[j + i__ * a_dim1]), abs(r__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0.f && x[
i__3].i == 0.f) || temp == 0.f);
i__3 = jx;
@@ -641,7 +641,7 @@ mv.f"> */
i__3 = i__ + j * a_dim1;
temp = (r__1 = a[i__3].r, abs(r__1)) + (r__2 = r_imag(
&a[i__ + j * a_dim1]), abs(r__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0.f && x[
i__3].i == 0.f) || temp == 0.f);
i__3 = jx;
@@ -674,7 +674,7 @@ mv.f"> */
i__3 = i__ + j * a_dim1;
temp = (r__1 = a[i__3].r, abs(r__1)) + (r__2 = r_imag(
&a[i__ + j * a_dim1]), abs(r__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0.f && x[
i__3].i == 0.f) || temp == 0.f);
i__3 = jx;
@@ -687,7 +687,7 @@ mv.f"> */
i__3 = j + i__ * a_dim1;
temp = (r__1 = a[i__3].r, abs(r__1)) + (r__2 = r_imag(
&a[j + i__ * a_dim1]), abs(r__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0.f && x[
i__3].i == 0.f) || temp == 0.f);
i__3 = jx;
diff --git a/lapack-netlib/SRC/cla_syamv.f b/lapack-netlib/SRC/cla_syamv.f
index 6507495838..4115576ef8 100644
--- a/lapack-netlib/SRC/cla_syamv.f
+++ b/lapack-netlib/SRC/cla_syamv.f
@@ -5,7 +5,6 @@
* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/
*
-*> \htmlonly
*> Download CLA_SYAMV + dependencies
*>
*> [TGZ]
@@ -13,7 +12,6 @@
*> [ZIP]
*>
*> [TXT]
-*> \endhtmlonly
*
* Definition:
* ===========
@@ -155,7 +153,7 @@
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \ingroup complexSYcomputational
+*> \ingroup la_heamv
*
*> \par Further Details:
* =====================
@@ -176,6 +174,7 @@
* =====================================================================
SUBROUTINE CLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
$ INCY )
+ IMPLICIT NONE
*
* -- LAPACK computational routine --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -264,7 +263,7 @@ SUBROUTINE CLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
* number of additions in each row.
*
SAFE1 = SLAMCH( 'Safe minimum' )
- SAFE1 = (N+1)*SAFE1
+ SAFE1 = REAL( N+1 )*SAFE1
*
* Form y := alpha*abs(A)*abs(x) + beta*abs(y).
*
@@ -358,7 +357,7 @@ SUBROUTINE CLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = CABS1( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -366,7 +365,7 @@ SUBROUTINE CLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = CABS1( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -394,7 +393,7 @@ SUBROUTINE CLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = CABS1( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -402,7 +401,7 @@ SUBROUTINE CLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = CABS1( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
diff --git a/lapack-netlib/SRC/dla_syamv.c b/lapack-netlib/SRC/dla_syamv.c
index 0c994832ed..f8ee4f30ce 100644
--- a/lapack-netlib/SRC/dla_syamv.c
+++ b/lapack-netlib/SRC/dla_syamv.c
@@ -578,7 +578,7 @@ mv.f"> */
i__2 = i__;
for (j = 1; j <= i__2; ++j) {
temp = (d__1 = a[j + i__ * a_dim1], abs(d__1));
- symb_zero__ = symb_zero__ && (x[j] == 0. || temp ==
+ symb_zero__ = symb_zero__ && (x[jx] == 0. || temp ==
0.);
y[iy] += *alpha * (d__1 = x[jx], abs(d__1)) * temp;
jx += *incx;
@@ -586,7 +586,7 @@ mv.f"> */
i__2 = *n;
for (j = i__ + 1; j <= i__2; ++j) {
temp = (d__1 = a[i__ + j * a_dim1], abs(d__1));
- symb_zero__ = symb_zero__ && (x[j] == 0. || temp ==
+ symb_zero__ = symb_zero__ && (x[jx] == 0. || temp ==
0.);
y[iy] += *alpha * (d__1 = x[jx], abs(d__1)) * temp;
jx += *incx;
@@ -614,7 +614,7 @@ mv.f"> */
i__2 = i__;
for (j = 1; j <= i__2; ++j) {
temp = (d__1 = a[i__ + j * a_dim1], abs(d__1));
- symb_zero__ = symb_zero__ && (x[j] == 0. || temp ==
+ symb_zero__ = symb_zero__ && (x[jx] == 0. || temp ==
0.);
y[iy] += *alpha * (d__1 = x[jx], abs(d__1)) * temp;
jx += *incx;
@@ -622,7 +622,7 @@ mv.f"> */
i__2 = *n;
for (j = i__ + 1; j <= i__2; ++j) {
temp = (d__1 = a[j + i__ * a_dim1], abs(d__1));
- symb_zero__ = symb_zero__ && (x[j] == 0. || temp ==
+ symb_zero__ = symb_zero__ && (x[jx] == 0. || temp ==
0.);
y[iy] += *alpha * (d__1 = x[jx], abs(d__1)) * temp;
jx += *incx;
diff --git a/lapack-netlib/SRC/dla_syamv.f b/lapack-netlib/SRC/dla_syamv.f
index 3911e2ab2f..b0c273002d 100644
--- a/lapack-netlib/SRC/dla_syamv.f
+++ b/lapack-netlib/SRC/dla_syamv.f
@@ -5,7 +5,6 @@
* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/
*
-*> \htmlonly
*> Download DLA_SYAMV + dependencies
*>
*> [TGZ]
@@ -13,7 +12,6 @@
*> [ZIP]
*>
*> [TXT]
-*> \endhtmlonly
*
* Definition:
* ===========
@@ -153,7 +151,7 @@
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \ingroup doubleSYcomputational
+*> \ingroup la_heamv
*
*> \par Further Details:
* =====================
@@ -174,6 +172,7 @@
* =====================================================================
SUBROUTINE DLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
$ INCY )
+ IMPLICIT NONE
*
* -- LAPACK computational routine --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -347,7 +346,7 @@ SUBROUTINE DLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = ABS( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*ABS( X( JX ) )*TEMP
JX = JX + INCX
@@ -355,7 +354,7 @@ SUBROUTINE DLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = ABS( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*ABS( X( JX ) )*TEMP
JX = JX + INCX
@@ -383,7 +382,7 @@ SUBROUTINE DLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = ABS( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*ABS( X( JX ) )*TEMP
JX = JX + INCX
@@ -391,7 +390,7 @@ SUBROUTINE DLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = ABS( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*ABS( X( JX ) )*TEMP
JX = JX + INCX
diff --git a/lapack-netlib/SRC/sla_syamv.c b/lapack-netlib/SRC/sla_syamv.c
index c8cad82aed..a91d7558b0 100644
--- a/lapack-netlib/SRC/sla_syamv.c
+++ b/lapack-netlib/SRC/sla_syamv.c
@@ -823,7 +823,7 @@ mv.f"> */
i__2 = i__;
for (j = 1; j <= i__2; ++j) {
temp = (r__1 = a[j + i__ * a_dim1], abs(r__1));
- symb_zero__ = symb_zero__ && (x[j] == 0.f || temp ==
+ symb_zero__ = symb_zero__ && (x[jx] == 0.f || temp ==
0.f);
y[iy] += *alpha * (r__1 = x[jx], abs(r__1)) * temp;
jx += *incx;
@@ -831,7 +831,7 @@ mv.f"> */
i__2 = *n;
for (j = i__ + 1; j <= i__2; ++j) {
temp = (r__1 = a[i__ + j * a_dim1], abs(r__1));
- symb_zero__ = symb_zero__ && (x[j] == 0.f || temp ==
+ symb_zero__ = symb_zero__ && (x[jx] == 0.f || temp ==
0.f);
y[iy] += *alpha * (r__1 = x[jx], abs(r__1)) * temp;
jx += *incx;
@@ -859,7 +859,7 @@ mv.f"> */
i__2 = i__;
for (j = 1; j <= i__2; ++j) {
temp = (r__1 = a[i__ + j * a_dim1], abs(r__1));
- symb_zero__ = symb_zero__ && (x[j] == 0.f || temp ==
+ symb_zero__ = symb_zero__ && (x[jx] == 0.f || temp ==
0.f);
y[iy] += *alpha * (r__1 = x[jx], abs(r__1)) * temp;
jx += *incx;
@@ -867,7 +867,7 @@ mv.f"> */
i__2 = *n;
for (j = i__ + 1; j <= i__2; ++j) {
temp = (r__1 = a[j + i__ * a_dim1], abs(r__1));
- symb_zero__ = symb_zero__ && (x[j] == 0.f || temp ==
+ symb_zero__ = symb_zero__ && (x[jx] == 0.f || temp ==
0.f);
y[iy] += *alpha * (r__1 = x[jx], abs(r__1)) * temp;
jx += *incx;
diff --git a/lapack-netlib/SRC/sla_syamv.f b/lapack-netlib/SRC/sla_syamv.f
index a44cc6e0d2..58783a2424 100644
--- a/lapack-netlib/SRC/sla_syamv.f
+++ b/lapack-netlib/SRC/sla_syamv.f
@@ -5,7 +5,6 @@
* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/
*
-*> \htmlonly
*> Download SLA_SYAMV + dependencies
*>
*> [TGZ]
@@ -13,7 +12,6 @@
*> [ZIP]
*>
*> [TXT]
-*> \endhtmlonly
*
* Definition:
* ===========
@@ -153,7 +151,7 @@
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \ingroup realSYcomputational
+*> \ingroup la_heamv
*
*> \par Further Details:
* =====================
@@ -174,6 +172,7 @@
* =====================================================================
SUBROUTINE SLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
$ INCY )
+ IMPLICIT NONE
*
* -- LAPACK computational routine --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -253,7 +252,7 @@ SUBROUTINE SLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
* number of additions in each row.
*
SAFE1 = SLAMCH( 'Safe minimum' )
- SAFE1 = (N+1)*SAFE1
+ SAFE1 = REAL( N+1 )*SAFE1
*
* Form y := alpha*abs(A)*abs(x) + beta*abs(y).
*
@@ -347,7 +346,7 @@ SUBROUTINE SLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = ABS( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*ABS( X( JX ) )*TEMP
JX = JX + INCX
@@ -355,7 +354,7 @@ SUBROUTINE SLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = ABS( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*ABS( X( JX ) )*TEMP
JX = JX + INCX
@@ -383,7 +382,7 @@ SUBROUTINE SLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = ABS( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*ABS( X( JX ) )*TEMP
JX = JX + INCX
@@ -391,7 +390,7 @@ SUBROUTINE SLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = ABS( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*ABS( X( JX ) )*TEMP
JX = JX + INCX
diff --git a/lapack-netlib/SRC/zla_heamv.c b/lapack-netlib/SRC/zla_heamv.c
index b6da294d6b..fece9464a9 100644
--- a/lapack-netlib/SRC/zla_heamv.c
+++ b/lapack-netlib/SRC/zla_heamv.c
@@ -620,7 +620,7 @@ mv.f"> */
i__3 = j + i__ * a_dim1;
temp = (d__1 = a[i__3].r, abs(d__1)) + (d__2 = d_imag(
&a[j + i__ * a_dim1]), abs(d__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0. && x[
i__3].i == 0.) || temp == 0.);
i__3 = jx;
@@ -633,7 +633,7 @@ mv.f"> */
i__3 = i__ + j * a_dim1;
temp = (d__1 = a[i__3].r, abs(d__1)) + (d__2 = d_imag(
&a[i__ + j * a_dim1]), abs(d__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0. && x[
i__3].i == 0.) || temp == 0.);
i__3 = jx;
@@ -666,7 +666,7 @@ mv.f"> */
i__3 = i__ + j * a_dim1;
temp = (d__1 = a[i__3].r, abs(d__1)) + (d__2 = d_imag(
&a[i__ + j * a_dim1]), abs(d__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0. && x[
i__3].i == 0.) || temp == 0.);
i__3 = jx;
@@ -679,7 +679,7 @@ mv.f"> */
i__3 = j + i__ * a_dim1;
temp = (d__1 = a[i__3].r, abs(d__1)) + (d__2 = d_imag(
&a[j + i__ * a_dim1]), abs(d__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0. && x[
i__3].i == 0.) || temp == 0.);
i__3 = jx;
diff --git a/lapack-netlib/SRC/zla_heamv.f b/lapack-netlib/SRC/zla_heamv.f
index 66cece241c..64dd5fb90e 100644
--- a/lapack-netlib/SRC/zla_heamv.f
+++ b/lapack-netlib/SRC/zla_heamv.f
@@ -5,7 +5,6 @@
* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/
*
-*> \htmlonly
*> Download ZLA_HEAMV + dependencies
*>
*> [TGZ]
@@ -13,7 +12,6 @@
*> [ZIP]
*>
*> [TXT]
-*> \endhtmlonly
*
* Definition:
* ===========
@@ -154,7 +152,7 @@
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \ingroup complex16HEcomputational
+*> \ingroup la_heamv
*
*> \par Further Details:
* =====================
@@ -175,6 +173,7 @@
* =====================================================================
SUBROUTINE ZLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
$ INCY )
+ IMPLICIT NONE
*
* -- LAPACK computational routine --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -356,7 +355,7 @@ SUBROUTINE ZLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = CABS1( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -364,7 +363,7 @@ SUBROUTINE ZLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = CABS1( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -392,7 +391,7 @@ SUBROUTINE ZLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = CABS1( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -400,7 +399,7 @@ SUBROUTINE ZLA_HEAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = CABS1( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
diff --git a/lapack-netlib/SRC/zla_syamv.c b/lapack-netlib/SRC/zla_syamv.c
index 21333fffbc..80c9086a33 100644
--- a/lapack-netlib/SRC/zla_syamv.c
+++ b/lapack-netlib/SRC/zla_syamv.c
@@ -631,7 +631,7 @@ mv.f"> */
i__3 = j + i__ * a_dim1;
temp = (d__1 = a[i__3].r, abs(d__1)) + (d__2 = d_imag(
&a[j + i__ * a_dim1]), abs(d__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0. && x[
i__3].i == 0.) || temp == 0.);
i__3 = jx;
@@ -644,7 +644,7 @@ mv.f"> */
i__3 = i__ + j * a_dim1;
temp = (d__1 = a[i__3].r, abs(d__1)) + (d__2 = d_imag(
&a[i__ + j * a_dim1]), abs(d__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0. && x[
i__3].i == 0.) || temp == 0.);
i__3 = jx;
@@ -677,7 +677,7 @@ mv.f"> */
i__3 = i__ + j * a_dim1;
temp = (d__1 = a[i__3].r, abs(d__1)) + (d__2 = d_imag(
&a[i__ + j * a_dim1]), abs(d__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0. && x[
i__3].i == 0.) || temp == 0.);
i__3 = jx;
@@ -690,7 +690,7 @@ mv.f"> */
i__3 = j + i__ * a_dim1;
temp = (d__1 = a[i__3].r, abs(d__1)) + (d__2 = d_imag(
&a[j + i__ * a_dim1]), abs(d__2));
- i__3 = j;
+ i__3 = jx;
symb_zero__ = symb_zero__ && ((x[i__3].r == 0. && x[
i__3].i == 0.) || temp == 0.);
i__3 = jx;
diff --git a/lapack-netlib/SRC/zla_syamv.f b/lapack-netlib/SRC/zla_syamv.f
index d00406fb57..9b8f726ba1 100644
--- a/lapack-netlib/SRC/zla_syamv.f
+++ b/lapack-netlib/SRC/zla_syamv.f
@@ -5,7 +5,6 @@
* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/
*
-*> \htmlonly
*> Download ZLA_SYAMV + dependencies
*>
*> [TGZ]
@@ -13,7 +12,6 @@
*> [ZIP]
*>
*> [TXT]
-*> \endhtmlonly
*
* Definition:
* ===========
@@ -155,7 +153,7 @@
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \ingroup complex16SYcomputational
+*> \ingroup la_heamv
*
*> \par Further Details:
* =====================
@@ -176,6 +174,7 @@
* =====================================================================
SUBROUTINE ZLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
$ INCY )
+ IMPLICIT NONE
*
* -- LAPACK computational routine --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -358,7 +357,7 @@ SUBROUTINE ZLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = CABS1( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -366,7 +365,7 @@ SUBROUTINE ZLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = CABS1( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -394,7 +393,7 @@ SUBROUTINE ZLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = 1, I
TEMP = CABS1( A( I, J ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX
@@ -402,7 +401,7 @@ SUBROUTINE ZLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y,
DO J = I+1, N
TEMP = CABS1( A( J, I ) )
SYMB_ZERO = SYMB_ZERO .AND.
- $ ( X( J ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
+ $ ( X( JX ) .EQ. ZERO .OR. TEMP .EQ. ZERO )
Y( IY ) = Y( IY ) + ALPHA*CABS1( X( JX ) )*TEMP
JX = JX + INCX