Mercurial > octave-nkf

comparison libcruft/odepack/prepj.f @ 4329:d53c33d93440

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2003-02-18 20:00:48 by jwe]
author jwe
date Tue, 18 Feb 2003 20:08:20 +0000
parents 44ed237bdc1e
children
comparison
equal deleted inserted replaced
4328:f7b63f362168 4329:d53c33d93440
27 C MITER = 1 OR 4, OR BY FINITE DIFFERENCING IF MITER = 2, 3, OR 5. 27 C MITER = 1 OR 4, OR BY FINITE DIFFERENCING IF MITER = 2, 3, OR 5.
28 C IF MITER = 3, A DIAGONAL APPROXIMATION TO J IS USED. 28 C IF MITER = 3, A DIAGONAL APPROXIMATION TO J IS USED.
29 C J IS STORED IN WM AND REPLACED BY P. IF MITER .NE. 3, P IS THEN 29 C J IS STORED IN WM AND REPLACED BY P. IF MITER .NE. 3, P IS THEN
30 C SUBJECTED TO LU DECOMPOSITION IN PREPARATION FOR LATER SOLUTION 30 C SUBJECTED TO LU DECOMPOSITION IN PREPARATION FOR LATER SOLUTION
31 C OF LINEAR SYSTEMS WITH P AS COEFFICIENT MATRIX. THIS IS DONE 31 C OF LINEAR SYSTEMS WITH P AS COEFFICIENT MATRIX. THIS IS DONE
32 C BY DGEFA IF MITER = 1 OR 2, AND BY DGBFA IF MITER = 4 OR 5. 32 C BY DGETRF IF MITER = 1 OR 2, AND BY DGBTRF IF MITER = 4 OR 5.
33 C 33 C
34 C IN ADDITION TO VARIABLES DESCRIBED PREVIOUSLY, COMMUNICATION 34 C IN ADDITION TO VARIABLES DESCRIBED PREVIOUSLY, COMMUNICATION
35 C WITH PREPJ USES THE FOLLOWING.. 35 C WITH PREPJ USES THE FOLLOWING..
36 C Y = ARRAY CONTAINING PREDICTED VALUES ON ENTRY. 36 C Y = ARRAY CONTAINING PREDICTED VALUES ON ENTRY.
37 C FTEM = WORK ARRAY OF LENGTH N (ACOR IN STODE). 37 C FTEM = WORK ARRAY OF LENGTH N (ACOR IN STODE).
93 NP1 = N + 1 93 NP1 = N + 1
94 DO 250 I = 1,N 94 DO 250 I = 1,N
95 WM(J) = WM(J) + 1.0D0 95 WM(J) = WM(J) + 1.0D0
96 250 J = J + NP1 96 250 J = J + NP1
97 C DO LU DECOMPOSITION ON P. -------------------------------------------- 97 C DO LU DECOMPOSITION ON P. --------------------------------------------
98 CALL DGEFA (WM(3), N, N, IWM(21), IER) 98 CALL DGETRF ( N, N, WM(3), N, IWM(21), IER)
99 IF (IER .NE. 0) IERPJ = 1 99 IF (IER .NE. 0) IERPJ = 1
100 RETURN 100 RETURN
101 C IF MITER = 3, CONSTRUCT A DIAGONAL APPROXIMATION TO J AND P. --------- 101 C IF MITER = 3, CONSTRUCT A DIAGONAL APPROXIMATION TO J AND P. ---------
102 300 WM(2) = HL0 102 300 WM(2) = HL0
103 R = EL0*0.1D0 103 R = EL0*0.1D0
168 570 II = MBAND + 2 168 570 II = MBAND + 2
169 DO 580 I = 1,N 169 DO 580 I = 1,N
170 WM(II) = WM(II) + 1.0D0 170 WM(II) = WM(II) + 1.0D0
171 580 II = II + MEBAND 171 580 II = II + MEBAND
172 C DO LU DECOMPOSITION OF P. -------------------------------------------- 172 C DO LU DECOMPOSITION OF P. --------------------------------------------
173 CALL DGBFA (WM(3), MEBAND, N, ML, MU, IWM(21), IER) 173 CALL DGBTRF ( N, N, ML, MU, WM(3), MEBAND, IWM(21), IER)
174 IF (IER .NE. 0) IERPJ = 1 174 IF (IER .NE. 0) IERPJ = 1
175 RETURN 175 RETURN
176 C----------------------- END OF SUBROUTINE PREPJ ----------------------- 176 C----------------------- END OF SUBROUTINE PREPJ -----------------------
177 END 177 END