Mercurial > octave-nkf

comparison libcruft/daspk/dmatd.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 8389e78e67d4
children
comparison
equal deleted inserted replaced
4328:f7b63f362168 4329:d53c33d93440
54 C are used for communication between the 54 C are used for communication between the
55 C calling program and external user routines. 55 C calling program and external user routines.
56 C They are not altered by DMATD. 56 C They are not altered by DMATD.
57 C----------------------------------------------------------------------- 57 C-----------------------------------------------------------------------
58 C***ROUTINES CALLED 58 C***ROUTINES CALLED
59 C JACD, RES, DGEFA, DGBFA 59 C JACD, RES, DGETRF, DGBTRF
60 C 60 C
61 C***END PROLOGUE DMATD 61 C***END PROLOGUE DMATD
62 C 62 C
63 C 63 C
64 IMPLICIT DOUBLE PRECISION(A-H,O-Z) 64 IMPLICIT DOUBLE PRECISION(A-H,O-Z)
109 210 CONTINUE 109 210 CONTINUE
110 C 110 C
111 C 111 C
112 C Do dense-matrix LU decomposition on J. 112 C Do dense-matrix LU decomposition on J.
113 C 113 C
114 230 CALL DGEFA(WM,NEQ,NEQ,IWM(LIPVT),IER) 114 230 CALL DGETRF( NEQ, NEQ, WM, NEQ, IWM(LIPVT), IER)
115 RETURN 115 RETURN
116 C 116 C
117 C 117 C
118 C Dummy section for IWM(MTYPE)=3. 118 C Dummy section for IWM(MTYPE)=3.
119 C 119 C
173 540 CONTINUE 173 540 CONTINUE
174 C 174 C
175 C 175 C
176 C Do LU decomposition of banded J. 176 C Do LU decomposition of banded J.
177 C 177 C
178 550 CALL DGBFA (WM,MEBAND,NEQ,IWM(LML),IWM(LMU),IWM(LIPVT),IER) 178 550 CALL DGBTRF(NEQ, NEQ, IWM(LML), IWM(LMU), WM, MEBAND,
179 * IWM(LIPVT), IER)
179 RETURN 180 RETURN
180 C 181 C
181 C------END OF SUBROUTINE DMATD------------------------------------------ 182 C------END OF SUBROUTINE DMATD------------------------------------------
182 END 183 END