| 1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071 |
- *DECK SSLUTI
- SUBROUTINE SSLUTI (N, B, X, NELT, IA, JA, A, ISYM, RWORK, IWORK)
- C***BEGIN PROLOGUE SSLUTI
- C***PURPOSE SLAP MTSOLV for LDU Factorization.
- C This routine acts as an interface between the SLAP generic
- C MTSOLV calling convention and the routine that actually
- C -T
- C computes (LDU) B = X.
- C***LIBRARY SLATEC (SLAP)
- C***CATEGORY D2E
- C***TYPE SINGLE PRECISION (SSLUTI-S, DSLUTI-D)
- C***KEYWORDS ITERATIVE PRECONDITION, LINEAR SYSTEM SOLVE, SLAP, SPARSE
- C***AUTHOR Greenbaum, Anne, (Courant Institute)
- C Seager, Mark K., (LLNL)
- C Lawrence Livermore National Laboratory
- C PO BOX 808, L-60
- C Livermore, CA 94550 (510) 423-3141
- C seager@llnl.gov
- C***DESCRIPTION
- C It is assumed that RWORK and IWORK have initialized with
- C the information required for SSLUI4:
- C IWORK(1) = Starting location of IL in IWORK.
- C IWORK(2) = Starting location of JL in IWORK.
- C IWORK(3) = Starting location of IU in IWORK.
- C IWORK(4) = Starting location of JU in IWORK.
- C IWORK(5) = Starting location of L in RWORK.
- C IWORK(6) = Starting location of DINV in RWORK.
- C IWORK(7) = Starting location of U in RWORK.
- C See the DESCRIPTION of SSLUI4 for details.
- C***REFERENCES (NONE)
- C***ROUTINES CALLED SSLUI4
- C***REVISION HISTORY (YYMMDD)
- C 871119 DATE WRITTEN
- C 881213 Previous REVISION DATE
- C 890915 Made changes requested at July 1989 CML Meeting. (MKS)
- C 890922 Numerous changes to prologue to make closer to SLATEC
- C standard. (FNF)
- C 890929 Numerous changes to reduce SP/DP differences. (FNF)
- C 910411 Prologue converted to Version 4.0 format. (BAB)
- C 920511 Added complete declaration section. (WRB)
- C 921113 Corrected C***CATEGORY line. (FNF)
- C 930701 Updated CATEGORY section. (FNF, WRB)
- C***END PROLOGUE SSLUTI
- C .. Scalar Arguments ..
- INTEGER ISYM, N, NELT
- C .. Array Arguments ..
- REAL A(N), B(N), RWORK(*), X(N)
- INTEGER IA(NELT), IWORK(10), JA(NELT)
- C .. Local Scalars ..
- INTEGER LOCDIN, LOCIL, LOCIU, LOCJL, LOCJU, LOCL, LOCU
- C .. External Subroutines ..
- EXTERNAL SSLUI4
- C***FIRST EXECUTABLE STATEMENT SSLUTI
- C
- C Pull out the pointers to the L, D and U matrices and call
- C the workhorse routine.
- C
- LOCIL = IWORK(1)
- LOCJL = IWORK(2)
- LOCIU = IWORK(3)
- LOCJU = IWORK(4)
- LOCL = IWORK(5)
- LOCDIN = IWORK(6)
- LOCU = IWORK(7)
- C
- CALL SSLUI4(N, B, X, IWORK(LOCIL), IWORK(LOCJL), RWORK(LOCL),
- $ RWORK(LOCDIN), IWORK(LOCIU), IWORK(LOCJU), RWORK(LOCU))
- C
- RETURN
- C------------- LAST LINE OF SSLUTI FOLLOWS ----------------------------
- END
|
