C----------------------------------------------------------------------- C IMSL Name: LFTXG/DLFTXG (Single/Double precision version) C C Purpose: Compute the LU factorization of a real general sparse C matrix. C C Usage: CALL LFTXG (N, NZ, A, IROW, JCOL, IPARAM, RPARAM, C NFAC, NL, FAC, IRFAC, JCFAC, IPVT, JPVT) C C Example 1: INTEGER N, NZ PARAMETER (N=6, NZ=15) INTEGER IPARAM(6), IROW(NZ), JCOL(NZ), NFAC, NL, & IRFAC(3*NZ), JCFAC(3*NZ), IPVT(N), JPVT(N) REAL RPARAM(5), A(NZ), FAC(3*NZ) C DATA A/6., 10., 15., -3., 10., -1., -1., -3., -5., 1., 10., -1., & -2., -1., -2./ DATA IROW/6, 2, 3, 2, 4, 4, 5, 5, 5, 5, 1, 6, 6, 2, 4/ DATA JCOL/6, 2, 3, 3, 4, 5, 1, 6, 4, 5, 1, 1, 2, 4, 1/ C NFAC = 3*NZ C Use default options IPARAM(1) = 0 CALL LFTXG (N, NZ, A, IROW, JCOL, IPARAM, RPARAM, NFAC, NL, & FAC, IRFAC, JCFAC, IPVT, JPVT) C CALL WRRRN (' fac ', 1, NFAC, FAC, 1, 0) CALL WRIRN (' irfac ', 1, NFAC, IRFAC, 1, 0) CALL WRIRN (' jcfac ', 1, NFAC, JCFAC, 1, 0) CALL WRIRN (' p ', 1, N, IPVT, 1, 0) CALL WRIRN (' q ', 1, N, JPVT, 1, 0) C END C C fac C 1 2 3 4 5 6 7 8 9 10 C -0.10 -5.00 -0.20 -0.10 -0.10 -1.00 -0.20 4.90 -5.10 1.00 C C 11 12 13 14 15 16 C -1.00 30.00 6.00 -2.00 10.00 15.00 C C irfac C 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 C 3 4 4 5 5 6 6 6 5 5 4 4 3 3 2 1 C C jcfac C 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 C 2 3 1 4 2 5 2 6 6 5 6 4 4 3 2 1 C C p C 1 2 3 4 5 6 C 3 1 6 2 5 4 C C q C 1 2 3 4 5 6 C 3 1 2 6 5 4