Numérotation des lignes :

C MTC21     SOURCE    CHAT      05/01/13    01:53:37     5004CCCC **********************************************************************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

© Cast3M 2003 - Tous droits réservés.
Mentions légales