Sources Esope | Cast3M

mtc21
C MTC21     SOURCE    CHAT      05/01/13    01:53:37     5004
CCC
C **********************************************************************
CCC
      SUBROUTINE MTC21 (K,NDIMK,X,NDIMX,LAM,NDIMS,
     .                nmodel,nescri,ues,nnumer,deltax,kerre)
      IMPLICIT INTEGER(I-N)
      integer ndims,nmodel,nescri,nnumer,ndimx,ues
      integer ifi,jc,i,j,kerre,ndimk
      real*8    x(ndimx),lam,deltax,K(ndimk,ndimk),g
      real*8    Amat(64),Gmat(64),vecn(8),vecm(8),Avecm(8),vecntA(8)
      kerre=0
      do i=1,8
          vecntA(i)=0.D0
          Avecm(i)=0.D0
          vecn(i)=0.D0
          vecm(i)=0.D0
        enddo
      do i=1,64
          Amat(i)=0.D0
          Gmat(i)=0.D0
        enddo
      call MatGenHookinv(Gmat,ndimx,ndims)
      call VecFlMAC(x,ndimx,vecm,nmodel)
      call VecYiMAC(x,ndimx,vecn,nmodel)
      call HessMAC(x,ndimx,Amat,ndimx,nnumer,deltax,nmodel)
      do i=1,ndimx*ndimx
         Amat(i)=Gmat(i)+lam*Amat(i)
      enddo
      call DescLU(Amat,ndimx)
      call LUinv(Amat,Gmat,ndimx)
      g=0.D0
      do i=1,ndimx
         do j=1,ndimx
             g=g+vecn(i)*Gmat(i+ndimx*(j-1))*vecm(j)
         enddo
      enddo
      do i=1,ndims
         do j=1,ndimx
           Avecm(i)=Avecm(i)+Gmat(i+ndimx*(j-1))*vecm(j)
           vecntA(i)=vecntA(i)+vecn(j)*Gmat(j+ndimx*(i-1))
         enddo
      enddo
      do ifi=1,ndims
         do jc=1,ndims
            K(ifi,jc)=Gmat((jc-1)*ndimx+ifi)-Avecm(ifi)*vecntA(jc)/g
         enddo
      enddo
      return
      end