dslvd.f

Go to the documentation of this file.

C Work performed under the auspices of the U.S. Department of Energy
C by Lawrence Livermore National Laboratory under contract number
C W-7405-Eng-48.
C
      SUBROUTINE dslvd(NEQ,DELTA,WM,IWM)
C
C***BEGIN PROLOGUE  DSLVD
C***REFER TO  DDASPK
C***DATE WRITTEN   890101   (YYMMDD)
C***REVISION DATE  900926   (YYMMDD)
C***REVISION DATE  940701   (YYMMDD) (new LIPVT)
C
C-----------------------------------------------------------------------
C***DESCRIPTION
C
C     This routine manages the solution of the linear
C     system arising in the Newton iteration.
C     Real matrix information and real temporary storage
C     is stored in the array WM.
C     Integer matrix information is stored in the array IWM.
C     For a dense matrix, the LAPACK routine DGETRS is called.
C     For a banded matrix, the LAPACK routine DGBTRS is called.
C-----------------------------------------------------------------------
C***ROUTINES CALLED
C   DGETRS, DGBTRS
C
C***END PROLOGUE  DSLVD
C
C
      IMPLICIT DOUBLE PRECISION(a-h,o-z)
      dimension delta(*),wm(*),iwm(*)
C
      parameter(lml=1, lmu=2, lmtype=4, llciwp=30)
C
      lipvt = iwm(llciwp)
      mtype=iwm(lmtype)
      GO TO(100,100,300,400,400),mtype
C
C     Dense matrix.
C
100   CALL dgetrs('N', neq, 1, wm, neq, iwm(lipvt), delta, neq, inlpck)
      RETURN
C
C     Dummy section for MTYPE=3.
C
300   CONTINUE
      RETURN
C
C     Banded matrix.
C
400   meband=2*iwm(lml)+iwm(lmu)+1
      CALL dgbtrs('N', neq, iwm(lml), iwm(lmu), 1, wm, meband,
     *     iwm(lipvt), delta, neq, inlpck)
      RETURN
C
C------END OF SUBROUTINE DSLVD------------------------------------------
      END