Sources Esope | Cast3M

res2lv
C RES2LV    SOURCE    PV        22/04/19    16:18:10     11344          
 
      SUBROUTINE RESIDUO_2LEVELS(X,R,N,NDIMX,NMODEL,XTRI,SIG,VECM,LAM)
      IMPLICIT INTEGER(I-N)
      IMPLICIT real*8(a-h,o-z)
      integer n,i,j,ndims,ndimv,nmodel,ndimx
      dimension void(1)
      real*8  x(n),r(n),void,sig(4),vecm(4),xtri(ndimx),lam
      integer             augla
      real*8                    c
      common /auglagrang1/ augla
      common /auglagrang2/ c
      void=0.D0
      ndimv=1
      ndims=3
      call der_enerelas_dpral(x,sig,nmodel)
      if (ndimx.eq.3) then
       call yieldd(sig,ndims,void,ndimv,r(n),nmodel)
       call vflsig(sig,ndims,void,ndimv,vecm,nmodel)
      else if (ndimx.eq.4) then
       call yieldd(sig,ndims,x(ndimx),ndimv,r(n),nmodel)
       call vflsig(sig,ndims,x(ndimx),ndimv,vecm,nmodel)
       call vflvar(sig,ndims,x(ndimx),ndimv,vecm(ndimx),nmodel)
      endif
      do i=1,ndimx
         r(i)=x(i)-xtri(i)+ABS(lam+c*r(n))*vecm(i)
      enddo
      return
      end