pele.utils.neighbor_list.NeighborListSubsetBuild.NumericalHessian

NeighborListSubsetBuild.NumericalHessian(coords, eps=1e-06)

return the Hessian matrix of second derivatives computed numerically

this takes 2*len(coords) calls to getGradient