diff options
| -rw-r--r-- | SRC/zhetrd_hb2st.F | 23 |
1 files changed, 20 insertions, 3 deletions
diff --git a/SRC/zhetrd_hb2st.F b/SRC/zhetrd_hb2st.F index 71419481..8fed56e8 100644 --- a/SRC/zhetrd_hb2st.F +++ b/SRC/zhetrd_hb2st.F @@ -382,7 +382,10 @@ DO 40 I = 1, N-1 E( I ) = RZERO 40 CONTINUE - RETURN +* + HOUS( 1 ) = LHMIN + WORK( 1 ) = LWMIN + RETURN END IF * * Case KD=1: @@ -437,6 +440,20 @@ C CALL ZSCAL( N, TMP, Q( 1, I+1 ), 1 ) C END IF 70 CONTINUE ENDIF +#else + IF( UPPER ) THEN + DO 60 I = 1, N-1 + E( I ) = DBLE( AB( ABOFDPOS, I+1 ) ) + 60 CONTINUE + ELSE + DO 70 I = 1, N-1 + E( I ) = DBLE( AB( ABOFDPOS, I ) ) + 70 CONTINUE + ENDIF +#endif +* + HOUS( 1 ) = LHMIN + WORK( 1 ) = LWMIN RETURN END IF * @@ -473,7 +490,7 @@ C END IF THED = MIN( (STT + THGRSIZ -1), (N-1)) DO 110 I = STT, N-1 ED = MIN( I, THED ) - IF( STT.GT.ED ) GOTO 100 + IF( STT.GT.ED ) EXIT DO 120 M = 1, STEPERCOL ST = STT DO 130 SWEEPID = ST, ED @@ -537,7 +554,7 @@ C END IF #endif IF ( BLKLASTIND.GE.(N-1) ) THEN STT = STT + 1 - GOTO 130 + EXIT ENDIF 140 CONTINUE 130 CONTINUE |