Sources Esope | Cast3M

jhcrhm
C JHCRHM    SOURCE    CHAT      05/01/13    00:49:33     5004
CCC
C **********************************************************************
CCC
      SUBROUTINE JAC_RHMC (KMAT,NDIMK,SIG,NDIMS,LAM,XMAT,
     .                     nescri,ues,kerre)
      IMPLICIT INTEGER(I-N)
      integer   ndims,kerre,nescri,ues,ndimk,ndimv,i,j
      real*8    sig(ndims),Kmat(*),xmat(*),lam
      kerre=0
c      it is not necessary: call carac_mate_rhmc(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
      ndimv=0
      call jac0RHMC(Kmat,ndimk,sig,ndims,ndimv,lam,
     .              nescri,ues,kerre)
      return
      end