Sources Esope | Cast3M

jac_j2
C JAC_J2    SOURCE    CHAT      05/01/13    00:48:08     5004
CCC
C **********************************************************************
CCC
      SUBROUTINE JAC_J2 (KMAT,NDIMK,SIG,NDIMS,VINT,LAM,XMAT,
     .                   nescri,ues,kerre)
      IMPLICIT INTEGER(I-N)
      integer ndims,kerre,nescri,ues,ndimk,ndimx,ndimv,i,j
      real*8  sig(ndims),vint,Kmat(*),xmat(*),x(7),void,lam
      real*8            sigy0,kiso,siginf,velo
      common /j2j2data/ sigy0,kiso,siginf,velo
      kerre=0
      do i=1,7
        x(i)=0.D0
      enddo
      call carac_mate_j2(xmat)
      if (lam.le.0.D0) then
         do i=1,ndims
           do j=1,ndims
             Kmat(ndimk*(i-1)+j)=0.D0
           enddo
           Kmat(ndimk*(i-1)+i)=1.D0
         enddo
         return
      endif
      if ((kiso.eq.0.D0).and.(velo.eq.0.D0)) then
         ndimv=0
         call jac0J2(Kmat,ndimk,sig,ndims,ndimv,lam,
     .               nescri,ues,kerre)
         return
      else
         ndimx=ndims+1
         do i=1,ndims
              x(i)=sig(i)
         enddo
         x(ndimx)=vint
         void=0.D0
         call jac1J2(Kmat,ndimk,x,ndimx,lam,ndims,
     .               nescri,ues,kerre)
      endif
      return
      end