reserv.for

来自「著名的berkeley大学开发的结构动力非线性分析程序」· FOR 代码 · 共 31 行

C*****
      SUBROUTINE RESERV (A,AA,M,L,NF,KODST,NEQ,JCOL,K)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
C
      DIMENSION A(L),AA(L),M(NEQ)
C
      IF (KODST.EQ.0) GO TO 10
C
      CALL STORE (A,L,NF,K)
      REWIND NF
      GO TO 60
C
   10 IF (JCOL.GT.1) GO TO 30
      DO 20 I=1,L
   20 A(I)=AA(I)
      GO TO 60
   30 JDP=M(JCOL-1)
      JJ=1
      DO 50 J=JCOL,NEQ
      JD=M(J)
      JH=JD-JDP
      JDP=JD
      NE=MIN0(JH,JJ)
      DO 40 I=1,NE
      A(JD)=AA(JD)
   40 JD=JD-1
   50 JJ=JJ+1
C
   60 RETURN
      END