Sources Esope | Cast3M

hmc1oc
C HMC1OC    SOURCE    CB215821  17/07/21    21:15:16     9513           
 
      SUBROUTINE HESSMAC1O2C (X,NDIMX,AMAT,NDIMA,NMODEL,DELTAX,TYPX)
      IMPLICIT INTEGER(I-N)
      integer ndimx,ndima,nmodel,i,j
      real*8    x(ndimx),Amat(ndima,ndima),deltax,dx,typx(ndimx)
      complex*16 auxx(8),auxvecm(8)
      do i=1,8
         auxx(i)=cmplx(0.D0,0.D0)
         auxvecm(i)=cmplx(0.D0,0.D0)
      enddo
      do i=1,ndimx
         auxx(i)=x(i)
      enddo
      do i=1,ndimx
         dx=deltax*max(abs(x(i)),typx(i))
         auxx(i)=cmplx(x(i),dx)
         call c_vecflMAC(auxx,ndimx,auxvecm,nmodel)
         auxx(i)=x(i)
         do j=1,ndimx
            Amat(j,i)=aimag(auxvecm(j))/dx
         enddo
      enddo
      return
      end