diff options
| author | jofret | 2008-04-11 09:46:18 +0000 |
|---|---|---|
| committer | jofret | 2008-04-11 09:46:18 +0000 |
| commit | c679afbd8d08c322d8323db5f57e0ab31db0cfca (patch) | |
| tree | b308487f2ba1252003bbd964a12437e70aa6c6b1 /src/lib/lapack/zladiv.f | |
| parent | dd0287279b0299e0446d6b9aa9b028814dd30137 (diff) | |
| download | scilab2c-c679afbd8d08c322d8323db5f57e0ab31db0cfca.tar.gz scilab2c-c679afbd8d08c322d8323db5f57e0ab31db0cfca.tar.bz2 scilab2c-c679afbd8d08c322d8323db5f57e0ab31db0cfca.zip | |
Adding LAPACK and compilation process
Diffstat (limited to 'src/lib/lapack/zladiv.f')
| -rw-r--r-- | src/lib/lapack/zladiv.f | 46 |
1 files changed, 46 insertions, 0 deletions
diff --git a/src/lib/lapack/zladiv.f b/src/lib/lapack/zladiv.f new file mode 100644 index 00000000..4a12055e --- /dev/null +++ b/src/lib/lapack/zladiv.f @@ -0,0 +1,46 @@ + COMPLEX*16 FUNCTION ZLADIV( X, Y ) +* +* -- LAPACK auxiliary routine (version 3.1) -- +* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. +* November 2006 +* +* .. Scalar Arguments .. + COMPLEX*16 X, Y +* .. +* +* Purpose +* ======= +* +* ZLADIV := X / Y, where X and Y are complex. The computation of X / Y +* will not overflow on an intermediary step unless the results +* overflows. +* +* Arguments +* ========= +* +* X (input) COMPLEX*16 +* Y (input) COMPLEX*16 +* The complex scalars X and Y. +* +* ===================================================================== +* +* .. Local Scalars .. + DOUBLE PRECISION ZI, ZR +* .. +* .. External Subroutines .. + EXTERNAL DLADIV +* .. +* .. Intrinsic Functions .. + INTRINSIC DBLE, DCMPLX, DIMAG +* .. +* .. Executable Statements .. +* + CALL DLADIV( DBLE( X ), DIMAG( X ), DBLE( Y ), DIMAG( Y ), ZR, + $ ZI ) + ZLADIV = DCMPLX( ZR, ZI ) +* + RETURN +* +* End of ZLADIV +* + END |