MINICERN/packlib/kernlib/kernnum/f011fort/feqn.inc

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/02/15 17:49:03  mclareni
   6 * Kernlib
   7 *
   8 *
   9 *
  10 * feqn.inc
  11 *
  12           IF(IDIM .GE. N  .AND.  N .GT. 0  .AND.  K .GT. 0)  GOTO 210
  13           CALL TMPRNT(HNAME,N,IDIM,K)
  14           RETURN
  15  210      NXCH  =  IR(N)
  16           IF(NXCH .EQ. 0)  GOTO 220
  17           DO 212    M  =  1, NXCH
  18              IJ  =  IR(M)
  19              I   =  IJ / 4096
  20              J   =  MOD(IJ,4096)
  21              DO 211   L  =  1, K
  22                 TE      =  B(I,L)
  23                 B(I,L)  =  B(J,L)
  24                 B(J,L)  =  TE
  25  211            CONTINUE
  26  212         CONTINUE
  27  220      DO 221    L  =  1, K
  28              B(1,L)  =  A(1,1)*B(1,L)
  29  221         CONTINUE
  30           IF(N .EQ. 1)  GOTO 299
  31           DO 243    L  =  1, K
  32              DO 232   I  =  2, N
  33                 IM1  =  I-1
  34                 S21  =  - B(I,L)
  35                 DO 231   J  =  1, IM1
  36                    S21  =  DOTF(A(I,J),B(J,L),S21)
  37  231               CONTINUE
  38                 B(I,L)  =  - A(I,I)*S21
  39  232            CONTINUE
  40              NM1  =  N-1
  41              DO 242   I  =  1, NM1
  42                 NMI  =  N-I
  43                 S22  =  - B(NMI,L)
  44                 DO 241   J  =  1, I
  45                    NMJP1  =  N - J+1
  46                    S22    =  DOTF(A(NMI,NMJP1),B(NMJP1,L),S22)
  47  241               CONTINUE
  48                 B(NMI,L)  =  - S22
  49  242            CONTINUE
  50  243         CONTINUE
  51  299      CONTINUE