Sources Esope | Cast3M

graco13
C GRACO13   SOURCE    PV        19/01/30    21:15:06     10089          
      SUBROUTINE GRACO13(ilideb,posm,valm,vect1,nbinc)
*  descente de la resolution
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      integer ilideb(*),posm(*)
      real*8 valm(*),vect1(*)
      DO 2 idli=1,nbinc
**     p=0.d0
**     do 5 i=ilideb(idli),ilideb(idli+1)-2
**      inc=posm(i)
**      p=p+valm(i)*vect1(inc)
** 5   continue
**     vect1(idli)=vect1(idli)-p
      id=ilideb(idli)
      if=ilideb(idli+1)-2
      lon=if-id+1
      if (lon.gt.0)
     > vect1(idli)=vect1(idli)-ddotip(lon,valm(id),posm(id),vect1(1))
   2  continue
      RETURN
      END