C RES2LV SOURCE CHAT 05/01/13 02:51:15 5004 SUBROUTINE RESIDUO_2LEVELS(X,R,N,NDIMX,NMODEL,XTRI,SIG,VECM,LAM) IMPLICIT INTEGER(I-N) integer n,i,j,ndims,ndimv,nmodel,ndimx 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)+dabs(lam+c*r(n))*vecm(i) enddo return end