version 1.14, 2014/01/27 09:28:19
|
version 1.15, 2015/11/26 11:44:17
|
Line 91
|
Line 91
|
* ===================================================================== |
* ===================================================================== |
SUBROUTINE DLADIV( A, B, C, D, P, Q ) |
SUBROUTINE DLADIV( A, B, C, D, P, Q ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.5.0) -- |
* -- LAPACK auxiliary routine (version 3.6.0) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* January 2013 |
* January 2013 |
Line 177
|
Line 177
|
|
|
SUBROUTINE DLADIV1( A, B, C, D, P, Q ) |
SUBROUTINE DLADIV1( A, B, C, D, P, Q ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.5.0) -- |
* -- LAPACK auxiliary routine (version 3.6.0) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* January 2013 |
* January 2013 |
Line 215
|
Line 215
|
|
|
DOUBLE PRECISION FUNCTION DLADIV2( A, B, C, D, R, T ) |
DOUBLE PRECISION FUNCTION DLADIV2( A, B, C, D, R, T ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.5.0) -- |
* -- LAPACK auxiliary routine (version 3.6.0) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* January 2013 |
* January 2013 |
Line 237
|
Line 237
|
* |
* |
IF( R.NE.ZERO ) THEN |
IF( R.NE.ZERO ) THEN |
BR = B * R |
BR = B * R |
if( BR.NE.ZERO ) THEN |
IF( BR.NE.ZERO ) THEN |
DLADIV2 = (A + BR) * T |
DLADIV2 = (A + BR) * T |
ELSE |
ELSE |
DLADIV2 = A * T + (B * T) * R |
DLADIV2 = A * T + (B * T) * R |