diff options
author | deaglanhalligan <deaglanhalligan@8a072113-8704-0410-8d35-dd094bca7971> | 2009-02-11 06:06:54 +0000 |
---|---|---|
committer | deaglanhalligan <deaglanhalligan@8a072113-8704-0410-8d35-dd094bca7971> | 2009-02-11 06:06:54 +0000 |
commit | 1f0792306441cdcd175199ad6720671df249dc47 (patch) | |
tree | 3a1c16a498b614128bc3977416ee7b3090d03dae /SRC/zsyrfsx.f | |
parent | e9454e224e5fa4d61076493be3916a64f274a5ae (diff) |
Added Jasons patches to bring work and rwork down to 2*N in the complex refinement routines.
Diffstat (limited to 'SRC/zsyrfsx.f')
-rw-r--r-- | SRC/zsyrfsx.f | 9 |
1 files changed, 5 insertions, 4 deletions
diff --git a/SRC/zsyrfsx.f b/SRC/zsyrfsx.f index 2d204ab2..e08f1784 100644 --- a/SRC/zsyrfsx.f +++ b/SRC/zsyrfsx.f @@ -268,9 +268,9 @@ * is true, 0.0 is false. * Default: 1.0 (attempt componentwise convergence) * -* WORK (workspace) DOUBLE PRECISION array, dimension (4*N) +* WORK (workspace) COMPLEX*16 array, dimension (2*N) * -* IWORK (workspace) INTEGER array, dimension (N) +* RWORK (workspace) DOUBLE PRECISION array, dimension (2*N) * * INFO (output) INTEGER * = 0: Successful exit. The solution to every right-hand side is @@ -334,7 +334,7 @@ EXTERNAL XERBLA, ZSYCON, ZLA_SYRFSX_EXTENDED * .. * .. Intrinsic Functions .. - INTRINSIC MAX, SQRT + INTRINSIC MAX, SQRT, TRANSFER * .. * .. External Functions .. EXTERNAL LSAME, BLAS_FPINFO_X, ILATRANS, ILAPREC @@ -472,7 +472,8 @@ CALL ZLA_SYRFSX_EXTENDED( PREC_TYPE, UPLO, N, $ NRHS, A, LDA, AF, LDAF, IPIV, RCEQU, S, B, $ LDB, X, LDX, BERR, N_NORMS, ERR_BNDS_NORM, ERR_BNDS_COMP, - $ WORK(N+1), RWORK, WORK(2*N+1), WORK(1), RCOND, + $ WORK, RWORK, WORK(N+1), + $ TRANSFER (RWORK(1:2*N), (/ (ZERO, ZERO) /), N), RCOND, $ ITHRESH, RTHRESH, UNSTABLE_THRESH, IGNORE_CWISE, $ INFO ) END IF |