C MONDE1 SOURCE PV 09/11/05 21:15:21 6535 SUBROUTINE MONDE1 (IPPVV,VECTBB,VAL,IVPO, > NA,IPREL,MULRE,INC,IFIB) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) DIMENSION VECTBB(1),VAL(1),IVPO(1),IPPVV(1) IMOI1 = IVPO(2*IPPVV(NA+1)-1) DO 50 ILM=NA,1,-1 II=IPREL-1+ILM IDEB2=IPPVV(ILM)*2 IMOI2= IVPO(IDEB2-1) LLOM=IMOI1-IMOI2-1 IF (LLOM.GT.0) THEN IPOSM=IMOI2-1 VKON=VECTBB(II) IPLAC=IVPO(IDEB2)-1 IDEBZ=1 IPLAC2=II-LLOM-1 DO 2 IDEB3=IDEB2,IFIB,2 IMOI = IVPO ( IDEB3 +2) ILONZ=IMOI -IPLAC-IDEBZ IPLAC=IPLAC-IPLAC2 IDEBZC=IDEBZ+IPLAC2 DO 1 L=IDEBZC,MIN(IDEBZC+ILONZ,II)-1-3,4 VECTBB(L)=VECTBB(L)-VKON*VAL(IPLAC+L) VECTBB(L+1)=VECTBB(L+1)-VKON*VAL(IPLAC+L+1) VECTBB(L+2)=VECTBB(L+2)-VKON*VAL(IPLAC+L+2) VECTBB(L+3)=VECTBB(L+3)-VKON*VAL(IPLAC+L+3) 1 CONTINUE LDEB=L DO L=LDEB,MIN(IDEBZC+ILONZ,II)-1 VECTBB(L)=VECTBB(L)-VKON*VAL(IPLAC+L) ENDDO IF (IDEBZ.GE.LLOM) GOTO 3 IDEBZ=IVPO(IDEB3+1)-IPOSM IPLAC=IMOI-IDEBZ 2 CONTINUE 3 CONTINUE IF (MULRE.EQ.1) GOTO 60 IPLAC=IVPO(IDEB2)-1 IDEBZ=1 DO 12 IDEB3=IDEB2,IFIB,2 IMOI = IVPO ( IDEB3 +2) ILONZ=IMOI -IPLAC-IDEBZ IPLAC=IPLAC-IPLAC2 IDEBZC=IDEBZ+IPLAC2 LHAU=MIN(IDEBZC+ILONZ,II)-1 DO 15 K=INC,(MULRE-1)*INC,INC VKON=VECTBB(II+K) DO 11 L=IDEBZC,LHAU-3,4 J=L+K VECTBB(J)=VECTBB(J)-VKON*VAL(IPLAC+L) VECTBB(J+1)=VECTBB(J+1)-VKON*VAL(IPLAC+L+1) VECTBB(J+2)=VECTBB(J+2)-VKON*VAL(IPLAC+L+2) VECTBB(J+3)=VECTBB(J+3)-VKON*VAL(IPLAC+L+3) 11 CONTINUE ldeb=l DO L=ldeb,LHAU J=L+K VECTBB(J)=VECTBB(J)-VKON*VAL(IPLAC+L) ENDDO 15 CONTINUE IF (IDEBZ.GE.LLOM) GOTO 13 IDEBZ=IVPO(IDEB3+1)-IPOSM IPLAC=IMOI-IDEBZ 12 CONTINUE 13 CONTINUE 60 CONTINUE ENDIF IMOI1=IMOI2 50 CONTINUE RETURN END