1 : C JHCRHM SOURCE CHAT 05/01/13 00:49:33 5004 2 : CCC 3 : C ********************************************************************** 4 : CCC 5 : SUBROUTINE JAC_RHMC (KMAT,NDIMK,SIG,NDIMS,LAM,XMAT, 6 : . nescri,ues,kerre) 7 : IMPLICIT INTEGER(I-N) 8 : integer ndims,kerre,nescri,ues,ndimk,ndimv,i,j 9 : real*8 sig(ndims),Kmat(*),xmat(*),lam 10 : kerre=0 10 : c it is not necessary: call carac_mate_rhmc(xmat) 11 : if (lam.le.0.D0) then 12 : do i=1,ndims 13 : do j=1,ndims 14 : Kmat(ndimk*(i-1)+j)=0.D0 15 : enddo 16 : Kmat(ndimk*(i-1)+i)=1.D0 17 : enddo 18 : return 19 : endif 20 : ndimv=0 20 : call jac0RHMC(Kmat,ndimk,sig,ndims,ndimv,lam, 21 : . nescri,ues,kerre) 22 : return 23 : end 24 : 25 : 26 :
© Cast3M 2003 - All rights reserved.
Disclaimer