Esope source files

   1 : C HMC1O2    SOURCE    CHAT      05/01/13    00:24:50     5004           
   2 : 
   3 :       SUBROUTINE HESSMAC1O2 (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),auxvecm(8),dx,typx(ndimx)
   8 :       do i=1,8
   9 :           auxx(i)=0.D0
  10 :           auxvecm(i)=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)=x(i)+dx
  17 :          call vecflMAC(auxx,ndimx,Amat(1,i),nmodel)
  18 :          auxx(i)=x(i)-dx
  18 :          call vecflMAC(auxx,ndimx,auxvecm,nmodel)
  19 :          auxx(i)=x(i)
  20 :          do j=1,ndimx
  21 :             Amat(j,i)=(Amat(j,i)-auxvecm(j))/(2.D0*dx)
  22 :          enddo
  23 :       enddo
  24 :       return
  25 :       end
  26 :  
  27 :  
  28 :