def METHOD_NAME(self):
    for calc in [{(1, 1): LennardJonesCut(1, 1, 3)}]:
        atoms = FaceCenteredCubic(
            'H', size=[2, 2, 2], latticeconstant=1.550)
        a = calculator.PairPotential(calc)
        atoms.calc = a
        H_numerical = numerical_hessian(atoms, d=1e-5, indices=None)
        H_numerical = H_numerical.todense()
        self.assertArrayAlmostEqual(
            np.sum(np.abs(H_numerical - H_numerical.T)), 0, tol=1e-5)