Sources Esope | Cast3M

vysgmc
C VYSGMC    SOURCE    CHAT      05/01/13    04:11:11     5004
 
      SUBROUTINE VYISIGMAC (SIG,NDIMS,VAR,NDIMV,VECM,NMODEL)
      IMPLICIT INTEGER(I-N)
      integer  ndims,ndimv,nmodel,i
      real*8     sig(ndims),siginv(3),var(ndimv),vecm(ndims),
     .           x1(6),x2(6),x3(6),xp,xq,xt
      call InvariantesPQT(sig,ndims,siginv)
      call vyiinvMAC(siginv,var,ndimv,xp,xq,xt,nmodel)
      call DerInvariPQT(sig,x1,x2,x3,ndims)
      do i=1,ndims
         vecm(i)=xp*x1(i)+xq*x2(i)+xt*x3(i)
      enddo
      return
      end