share/hompack/fortran/mfacds.f

   1       SUBROUTINE MFACDS(NN,Q,LENAA,MAXA)
   2 C
   3 C     SETS UP Q AS THE FACTORIZATION OF THE ENTIRE
   4 C     (NN+1) X (NN+1) MATRIX M.
   5 C
   6 C on input:
   7 C
   8 C NN  is the dimension of the symmetric matrix AA, the upper left corner
   9 C     of the augmented Jacobian matrix B.
  10 C
  11 C Q  contains AA in its first LENAA positions.
  12 C
  13 C LENAA  is the length of the one dimensional packed array AA.
  14 C
  15 C MAXA  is the integer array used for packed skyline storage.  It describes
  16 C       AA and M, the symmetric piece of B.
  17 C
  18 C on output:
  19 C
  20 C NN, LENAA, and MAXA are unchanged.
  21 C
  22 C Q  contains an approximate factorization of M, in packed skyline storage
  23 C    form.
  24 C
  25 C
  26 C Calls  GMFADS .
  27 C
  28       INTEGER I,IMAX,LENAA,LENQ,NN,MAXA(NN+2),NQ
  29       DOUBLE PRECISION Q(LENAA+NN+1)
  30 C
  31       NQ=NN+1
  32       IMAX=MAXA(NN+2)-LENAA-2
  33       LENQ=MAXA(NN+2)-1
  34 C
  35       DO 100 I=1,IMAX,1
  36          Q(LENAA+I)=0.0
  37 100   CONTINUE
  38       Q(LENQ)=1.0D0
  39 C
  40       CALL GMFADS(NQ,Q,LENQ,MAXA)
  41 C
  42       RETURN
  43       END