diff options
| author | jofret | 2009-04-28 07:17:00 +0000 |
|---|---|---|
| committer | jofret | 2009-04-28 07:17:00 +0000 |
| commit | 8c8d2f518968ce7057eec6aa5cd5aec8faab861a (patch) | |
| tree | 3dd1788b71d6a3ce2b73d2d475a3133580e17530 /src/lib/lapack/dladiv.f | |
| parent | 9f652ffc16a310ac6641a9766c5b9e2671e0e9cb (diff) | |
| download | scilab2c-8c8d2f518968ce7057eec6aa5cd5aec8faab861a.tar.gz scilab2c-8c8d2f518968ce7057eec6aa5cd5aec8faab861a.tar.bz2 scilab2c-8c8d2f518968ce7057eec6aa5cd5aec8faab861a.zip | |
Moving lapack to right place
Diffstat (limited to 'src/lib/lapack/dladiv.f')
| -rw-r--r-- | src/lib/lapack/dladiv.f | 62 |
1 files changed, 0 insertions, 62 deletions
diff --git a/src/lib/lapack/dladiv.f b/src/lib/lapack/dladiv.f deleted file mode 100644 index b6a74d1b..00000000 --- a/src/lib/lapack/dladiv.f +++ /dev/null @@ -1,62 +0,0 @@ - SUBROUTINE DLADIV( A, B, C, D, P, Q ) -* -* -- LAPACK auxiliary routine (version 3.1) -- -* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. -* November 2006 -* -* .. Scalar Arguments .. - DOUBLE PRECISION A, B, C, D, P, Q -* .. -* -* Purpose -* ======= -* -* DLADIV performs complex division in real arithmetic -* -* a + i*b -* p + i*q = --------- -* c + i*d -* -* The algorithm is due to Robert L. Smith and can be found -* in D. Knuth, The art of Computer Programming, Vol.2, p.195 -* -* Arguments -* ========= -* -* A (input) DOUBLE PRECISION -* B (input) DOUBLE PRECISION -* C (input) DOUBLE PRECISION -* D (input) DOUBLE PRECISION -* The scalars a, b, c, and d in the above expression. -* -* P (output) DOUBLE PRECISION -* Q (output) DOUBLE PRECISION -* The scalars p and q in the above expression. -* -* ===================================================================== -* -* .. Local Scalars .. - DOUBLE PRECISION E, F -* .. -* .. Intrinsic Functions .. - INTRINSIC ABS -* .. -* .. Executable Statements .. -* - IF( ABS( D ).LT.ABS( C ) ) THEN - E = D / C - F = C + D*E - P = ( A+B*E ) / F - Q = ( B-A*E ) / F - ELSE - E = C / D - F = D + C*E - P = ( B+A*E ) / F - Q = ( -A+B*E ) / F - END IF -* - RETURN -* -* End of DLADIV -* - END |