share/odepack/fortran/dslsbt.f

   1 *DECK DSLSBT
   2       SUBROUTINE DSLSBT (WM, IWM, X, TEM)
   3       INTEGER IWM
   4       INTEGER LBLOX, LPB, LPC, MB, NB
   5       DOUBLE PRECISION WM, X, TEM
   6       DIMENSION WM(*), IWM(*), X(*), TEM(*)
   7 C-----------------------------------------------------------------------
   8 C This routine acts as an interface between the core integrator
   9 C routine and the DSOLBT routine for the solution of the linear system
  10 C arising from chord iteration.
  11 C Communication with DSLSBT uses the following variables:
  12 C WM    = real work space containing the LU decomposition,
  13 C         starting at WM(3).
  14 C IWM   = integer work space containing pivot information, starting at
  15 C         IWM(21).  IWM also contains block structure parameters
  16 C         MB = IWM(1) and NB = IWM(2).
  17 C X     = the right-hand side vector on input, and the solution vector
  18 C         on output, of length N.
  19 C TEM   = vector of work space of length N, not used in this version.
  20 C-----------------------------------------------------------------------
  21       MB = IWM(1)
  22       NB = IWM(2)
  23       LBLOX = MB*MB*NB
  24       LPB = 3 + LBLOX
  25       LPC = LPB + LBLOX
  26       CALL DSOLBT (MB, NB, WM(3), WM(LPB), WM(LPC), X, IWM(21))
  27       RETURN
  28 C----------------------- End of Subroutine DSLSBT ----------------------
  29       END