Skip to content

hessian

esta.xtbBag.hessian

reading full hessian from xtb hessian output file

get_full_hessian(fname, natom)

read full hessian from the xtb hessian output file

Parameters:

  • fname

    input file: str name containing hessian information

  • natom
  • natom

Returns:

  • hess ( array, rank 2d array of dim (3*natom, 3*natom) ) –
  • hessian matrix
Note

read hessian matrix from the xtb hessain file

the equivalent fortran90 code is:

DO i=1,natom*3 maxcol = 0 DO WHILE(maxcol .lt. natom*3) mincol = maxcol + 1 maxcol = MIN(maxcol+5,natom*3) READ(10, ) (hessian(j,i), j = mincol, maxcol) !write(,*) (hessian(j,i), j = mincol, maxcol) ENDDO ENDDO