Numérotation des lignes :

C JAC1J2    SOURCE    CHAT      05/01/13    00:47:59     5004CCCC **********************************************************************CCC      SUBROUTINE JAC1J2 (KMAT,NDIMK,X,NDIMX,LAM,NDIMS,     .                   nescri,ues,kerre)      IMPLICIT INTEGER(I-N)      integer ndims,nescri,ues,ndimv,ndimx,i,j,k,kerre,ndimk      real*8  x(ndimx),lam,void,Kmat(ndimk,ndimk),g      real*8  Amat(49),Gmat(49),vecn(7),vecm(7),Avecm(7),vecntA(7),     .        vecEm(7),Hmat(49)      kerre=0        ndimv=1      void =0.D0        do i=1,7         Avecm(i)=0.D0         vecntA(i)=0.D0         vecn(i)=0.D0         vecm(i)=0.D0         vecEm(i)=0.D0        enddo        do i=1,49         amat(i)=0.D0         gmat(i)=0.D0         hmat(i)=0.D0        enddoc m      call vflsigJ2(x,ndims,x(ndimx),ndimv,vecm)        call vflvarJ2(x,ndims,x(ndimx),ndimv,vecm(ndimx))c n      call vyisigJ2(x,ndims,x(ndimx),ndimv,vecn)        call vyivarJ2(x,ndims,x(ndimx),ndimv,vecn(ndimx))c H      call HessFlsigJ2(x,ndims,x(ndimx),ndimv,Hmat,ndimx) c G=E      call MatGenHook(Gmat,ndimx,ndims)c Em = E*m      do i=1,ndimx       do j=1,ndimx         vecEm(i)=vecEm(i)+Gmat(ndimx*(j-1)+i)*vecm(j)       enddo      enddoc A = I+l*E*H      do i=1,ndimx       do j=1,ndimx        do k=1,ndimx          Amat(ndimx*(j-1)+i)=Amat(ndimx*(j-1)+i)+     .                        Gmat(ndimx*(k-1)+i)*Hmat(ndimx*(j-1)+k)        enddo        Amat(ndimx*(j-1)+i)=lam*Amat(ndimx*(j-1)+i)       enddo       Amat(ndimx*(i-1)+i)=1.D0+Amat(ndimx*(i-1)+i)      enddoc G = A-1      call DescLU(Amat,ndimx)      call LUinv(Amat,Gmat,ndimx)c g = nt*A-1*Em      g=0.D0      do i=1,ndimx       do j=1,ndimx          g=g+vecn(i)*Gmat(i+ndimx*(j-1))*vecEm(j)       enddo      enddoc G*Em; nt*G      do i=1,ndimx       do j=1,ndimx         Avecm(i)=Avecm(i)+Gmat(i+ndimx*(j-1))*vecEm(j)         vecntA(i)=vecntA(i)+vecn(j)*Gmat(j+ndimx*(i-1))       enddo      enddoc K = G-(G*m * nt*G)/g      do i=1,ndimx       do j=1,ndimx          Kmat(i,j)=Gmat((j-1)*ndimx+i)-Avecm(i)*vecntA(j)/g       enddo      enddo      return      end

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