GNU/Linux |
CentOS 4.8 |
i386 |
![]() |
dlasq5(l) |
![]() |
DLASQ5 - compute one dqds transform in ping-pong form, one version for IEEE machines another for non IEEE machines
SUBROUTINE DLASQ5( |
I0, N0, Z, PP, TAU, DMIN, DMIN1, DMIN2, DN, DNM1, DNM2, IEEE ) | ||
LOGICAL |
IEEE | ||
INTEGER |
I0, N0, PP | ||
DOUBLE |
PRECISION DMIN, DMIN1, DMIN2, DN, DNM1, DNM2, TAU | ||
DOUBLE |
PRECISION Z( * ) |
DLASQ5 computes one dqds transform in ping-pong form, one version for IEEE machines another for non IEEE machines.
I0 (input) INTEGER
First index.
N0 (input) INTEGER
Last index.
Z (input) DOUBLE PRECISION array, dimension ( 4*N )
Z holds the qd array. EMIN is stored in Z(4*N0) to avoid an extra argument.
PP (input) INTEGER
PP=0 for ping, PP=1 for pong.
TAU (input) DOUBLE PRECISION
This is the shift.
DMIN (output) DOUBLE PRECISION
Minimum value of d.
DMIN1 (output) DOUBLE PRECISION Minimum value of d, excluding D( N0 ).
DMIN2 (output) DOUBLE PRECISION Minimum value of d, excluding D( N0 ) and D( N0-1 ).
DN (output) DOUBLE PRECISION
d(N0), the last value of d.
DNM1 (output) DOUBLE PRECISION
d(N0-1).
DNM2 (output) DOUBLE PRECISION
d(N0-2).
IEEE (input) LOGICAL
Flag for IEEE or non IEEE arithmetic.
![]() |
dlasq5(l) | ![]() |