home *** CD-ROM | disk | FTP | other *** search
/ Geek Gadgets 1 / ADE-1.bin / ade-dist / octave-1.1.1p1-src.tgz / tar.out / fsf / octave / libcruft / dassl / ddaslv.f < prev    next >
Text File  |  1996-09-28  |  2KB  |  61 lines
  1.       SUBROUTINE DDASLV (NEQ, DELTA, WM, IWM)
  2. C***BEGIN PROLOGUE  DDASLV
  3. C***SUBSIDIARY
  4. C***PURPOSE  Linear system solver for DDASSL.
  5. C***LIBRARY   SLATEC (DASSL)
  6. C***TYPE      DOUBLE PRECISION (SDASLV-S, DDASLV-D)
  7. C***AUTHOR  PETZOLD, LINDA R., (LLNL)
  8. C***DESCRIPTION
  9. C-----------------------------------------------------------------------
  10. C     THIS ROUTINE MANAGES THE SOLUTION OF THE LINEAR
  11. C     SYSTEM ARISING IN THE NEWTON ITERATION.
  12. C     MATRICES AND REAL TEMPORARY STORAGE AND
  13. C     REAL INFORMATION ARE STORED IN THE ARRAY WM.
  14. C     INTEGER MATRIX INFORMATION IS STORED IN
  15. C     THE ARRAY IWM.
  16. C     FOR A DENSE MATRIX, THE LINPACK ROUTINE
  17. C     DGESL IS CALLED.
  18. C     FOR A BANDED MATRIX,THE LINPACK ROUTINE
  19. C     DGBSL IS CALLED.
  20. C-----------------------------------------------------------------------
  21. C***ROUTINES CALLED  DGBSL, DGESL
  22. C***REVISION HISTORY  (YYMMDD)
  23. C   830315  DATE WRITTEN
  24. C   901009  Finished conversion to SLATEC 4.0 format (F.N.Fritsch)
  25. C   901019  Merged changes made by C. Ulrich with SLATEC 4.0 format.
  26. C   901026  Added explicit declarations for all variables and minor
  27. C           cosmetic changes to prologue.  (FNF)
  28. C***END PROLOGUE  DDASLV
  29. C
  30.       INTEGER  NEQ, IWM(*)
  31.       DOUBLE PRECISION  DELTA(*), WM(*)
  32. C
  33.       EXTERNAL  DGBSL, DGESL
  34. C
  35.       INTEGER  LIPVT, LML, LMU, LMTYPE, MEBAND, MTYPE, NPD
  36.       PARAMETER (NPD=1)
  37.       PARAMETER (LML=1)
  38.       PARAMETER (LMU=2)
  39.       PARAMETER (LMTYPE=4)
  40.       PARAMETER (LIPVT=21)
  41. C
  42. C***FIRST EXECUTABLE STATEMENT  DDASLV
  43.       MTYPE=IWM(LMTYPE)
  44.       GO TO(100,100,300,400,400),MTYPE
  45. C
  46. C     DENSE MATRIX
  47. 100   CALL DGESL(WM(NPD),NEQ,NEQ,IWM(LIPVT),DELTA,0)
  48.       RETURN
  49. C
  50. C     DUMMY SECTION FOR MTYPE=3
  51. 300   CONTINUE
  52.       RETURN
  53. C
  54. C     BANDED MATRIX
  55. 400   MEBAND=2*IWM(LML)+IWM(LMU)+1
  56.       CALL DGBSL(WM(NPD),MEBAND,NEQ,IWM(LML),
  57.      *  IWM(LMU),IWM(LIPVT),DELTA,0)
  58.       RETURN
  59. C------END OF SUBROUTINE DDASLV------
  60.       END
  61.