Esope source files

   1 : C HMC1OC    SOURCE    CHAT      05/01/13    00:24:55     5004           
   2 : 
   3 :       SUBROUTINE HESSMAC1O2C (X,NDIMX,AMAT,NDIMA,NMODEL,DELTAX,TYPX)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       integer ndimx,ndima,nmodel,i,j
   6 :       real*8    x(ndimx),Amat(ndima,ndima),deltax,dx,typx(ndimx)
   7 :       complex*16 auxx(8),auxvecm(8)
   8 :       do i=1,8
   9 :          auxx(i)=cmplx(0.D0,0.D0)
  10 :          auxvecm(i)=cmplx(0.D0,0.D0)
  11 :       enddo
  12 :       do i=1,ndimx
  13 :          auxx(i)=x(i)
  14 :       enddo
  15 :       do i=1,ndimx
  16 :          dx=deltax*max(abs(x(i)),typx(i))
  17 :          auxx(i)=cmplx(x(i),dx)
  17 :          call c_vecflMAC(auxx,ndimx,auxvecm,nmodel)
  18 :          auxx(i)=x(i)
  19 :          do j=1,ndimx
  20 :             Amat(j,i)=imag(auxvecm(j))/dx
  21 :          enddo
  22 :       enddo
  23 :       return
  24 :       end
  25 :  
  26 :  
  27 :