public interface latdf
LATDF computes the contribution to the reciprocal Dif-estimate
by solving for x in Z * x = b, where b is chosen such that the norm
of x is as large as possible. It is assumed that LU decomposition
of Z has been computed by CGETC2. On entry RHS = f holds the
contribution from earlier solved sub-systems, and on return RHS = x.
The factorization of Z returned by CGETC2 has the form
Z = P * L * U * Q, where P and Q are permutation matrices. L is lower
triangular with unit diagonal elements and U is upper triangular.
Subroutines
Arguments
| Type |
Intent | Optional | Attributes |
|
Name |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ijob |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
n |
|
|
complex(kind=sp),
|
intent(inout) |
|
|
:: |
z(ldz,*) |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ldz |
|
|
complex(kind=sp),
|
intent(inout) |
|
|
:: |
rhs(*) |
|
|
real(kind=sp),
|
intent(inout) |
|
|
:: |
rdsum |
|
|
real(kind=sp),
|
intent(inout) |
|
|
:: |
rdscal |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ipiv(*) |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
jpiv(*) |
|
Arguments
| Type |
Intent | Optional | Attributes |
|
Name |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ijob |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
n |
|
|
real(kind=dp),
|
intent(inout) |
|
|
:: |
z(ldz,*) |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ldz |
|
|
real(kind=dp),
|
intent(inout) |
|
|
:: |
rhs(*) |
|
|
real(kind=dp),
|
intent(inout) |
|
|
:: |
rdsum |
|
|
real(kind=dp),
|
intent(inout) |
|
|
:: |
rdscal |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ipiv(*) |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
jpiv(*) |
|
Arguments
| Type |
Intent | Optional | Attributes |
|
Name |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ijob |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
n |
|
|
real(kind=sp),
|
intent(inout) |
|
|
:: |
z(ldz,*) |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ldz |
|
|
real(kind=sp),
|
intent(inout) |
|
|
:: |
rhs(*) |
|
|
real(kind=sp),
|
intent(inout) |
|
|
:: |
rdsum |
|
|
real(kind=sp),
|
intent(inout) |
|
|
:: |
rdscal |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ipiv(*) |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
jpiv(*) |
|
Arguments
| Type |
Intent | Optional | Attributes |
|
Name |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ijob |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
n |
|
|
complex(kind=dp),
|
intent(inout) |
|
|
:: |
z(ldz,*) |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ldz |
|
|
complex(kind=dp),
|
intent(inout) |
|
|
:: |
rhs(*) |
|
|
real(kind=dp),
|
intent(inout) |
|
|
:: |
rdsum |
|
|
real(kind=dp),
|
intent(inout) |
|
|
:: |
rdscal |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
ipiv(*) |
|
|
integer(kind=ilp),
|
intent(in) |
|
|
:: |
jpiv(*) |
|
Module Procedures