def hessian(self):
result = np.zeros((self.numc, <NUM_LIT:3>, self.numc, <NUM_LIT:3>), float)<EOL>for index1 in range(self.numc):<EOL><INDENT>for index2 in range(self.numc):<EOL><INDENT>result[index1, :, index2, :] = self.hessian_component(index1, index2)<EOL><DEDENT><DEDENT>return result<EOL>
Compute the hessian of the energy