*DECK SSLUI SUBROUTINE SSLUI (N, B, X, NELT, IA, JA, A, ISYM, RWORK, IWORK) C***BEGIN PROLOGUE SSLUI C***PURPOSE SLAP MSOLVE for LDU Factorization. C This routine acts as an interface between the SLAP generic C MSOLVE calling convention and the routine that actually C -1 C computes (LDU) B = X. C***LIBRARY SLATEC (SLAP) C***CATEGORY D2E C***TYPE SINGLE PRECISION (SSLUI-S, DSLUI-D) C***KEYWORDS ITERATIVE PRECONDITION, NON-SYMMETRIC LINEAR SYSTEM SOLVE, C 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 SSLUI2: 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 SSLUI2 for details. C***REFERENCES (NONE) C***ROUTINES CALLED SSLUI2 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 SSLUI C .. Scalar Arguments .. INTEGER ISYM, N, NELT C .. Array Arguments .. REAL A(NELT), 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 SSLUI2 C***FIRST EXECUTABLE STATEMENT SSLUI C C Pull out the locations of the arrays holding the ILU C factorization. C LOCIL = IWORK(1) LOCJL = IWORK(2) LOCIU = IWORK(3) LOCJU = IWORK(4) LOCL = IWORK(5) LOCDIN = IWORK(6) LOCU = IWORK(7) C C Solve the system LUx = b CALL SSLUI2(N, B, X, IWORK(LOCIL), IWORK(LOCJL), RWORK(LOCL), $ RWORK(LOCDIN), IWORK(LOCIU), IWORK(LOCJU), RWORK(LOCU) ) C RETURN C------------- LAST LINE OF SSLUI FOLLOWS ---------------------------- END