diff options
| author | Jean-Marc Zucconi <jmz@FreeBSD.org> | 1996-06-07 15:39:54 +0000 |
|---|---|---|
| committer | Jean-Marc Zucconi <jmz@FreeBSD.org> | 1996-06-07 15:39:54 +0000 |
| commit | d2e95f05c163f196a666702adf838324fd6ea6dd (patch) | |
| tree | e4b08b49cd31550bb706b566347298b0112a0593 /math/lapack/files/patch-ad | |
| parent | Upgrade to 19.31 (diff) | |
LAPACK is a library of Fortran 77 subroutines for linear algebra.
Notes
Notes:
svn path=/head/; revision=3178
Diffstat (limited to 'math/lapack/files/patch-ad')
| -rw-r--r-- | math/lapack/files/patch-ad | 22 |
1 files changed, 22 insertions, 0 deletions
diff --git a/math/lapack/files/patch-ad b/math/lapack/files/patch-ad new file mode 100644 index 000000000000..d88b6f107fd7 --- /dev/null +++ b/math/lapack/files/patch-ad @@ -0,0 +1,22 @@ +*** BLAS/SRC/zrotg.f~ Fri Oct 30 21:09:39 1992 +--- BLAS/SRC/zrotg.f Mon Jul 3 01:28:54 1995 +*************** +*** 10,17 **** + go to 20 + 10 continue + scale = cdabs(ca) + cdabs(cb) +! norm = scale*dsqrt((cdabs(ca/dcmplx(scale,0.0d0)))**2 + +! * (cdabs(cb/dcmplx(scale,0.0d0)))**2) + alpha = ca /cdabs(ca) + c = cdabs(ca) / norm + s = alpha * dconjg(cb) / norm +--- 10,18 ---- + go to 20 + 10 continue + scale = cdabs(ca) + cdabs(cb) +! c = (cdabs(ca/dcmplx(scale,0.0d0)))**2 + +! * (cdabs(cb/dcmplx(scale,0.0d0)))**2 +! norm = scale*dsqrt(c) + alpha = ca /cdabs(ca) + c = cdabs(ca) / norm + s = alpha * dconjg(cb) / norm |