Esope source files

   1 : C HMC1O1    SOURCE    CHAT      05/01/13    00:24:46     5004           
   2 : 
   3 :       SUBROUTINE HESSMAC1O1 (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
   7 :       real*8    auxx(8),vecm(8),dx,typx(ndimx)
   8 :       do i=1,8
   9 :           auxx(i)=0.D0
  10 :           vecm(i)=0.D0
  11 :         enddo
  11 :       call vecflMAC(x,ndimx,vecm,nmodel)
  12 :       do i=1,ndimx
  13 :          auxx(i)=x(i)
  14 :       enddo
  15 :       do i=1,ndimx
  16 :          dx=sign(deltax*max(abs(x(i)),typx(i)),x(i))
  17 :          auxx(i)=x(i)+dx
  17 :          call vecflMAC(auxx,ndimx,Amat(1,i),nmodel)
  18 :          auxx(i)=x(i)
  19 :          do j=1,ndimx
  20 :             Amat(j,i)=(Amat(j,i)-vecm(j))/dx
  21 :          enddo
  22 :       enddo
  23 :       return
  24 :       end
  25 :  
  26 :  
  27 :