def METHOD_NAME(self):
    ee, ff, vv = self.dp.eval(self.coords, self.box, self.atype, eval_fv=True)
    # check shape of the returns
    nframes = 1
    natoms = len(self.atype)
    self.assertEqual(ee.shape, (nframes,))
    self.assertEqual(ff.shape, (nframes, natoms, 3))
    self.assertEqual(vv.shape, (nframes, 9))
    self.assertEqual(self.expected_e.shape, (nframes,))
    self.assertEqual(self.expected_f.shape, (nframes * natoms * 3,))
    self.assertEqual(self.expected_v.shape, (nframes * 9,))
    # np.savetxt('ee.out', ee.reshape([1, -1]), delimiter=',')
    # np.savetxt('ff.out', ff.reshape([1, -1]), delimiter=',')
    # np.savetxt('vv.out', vv.reshape([1, -1]), delimiter=',')
    ee = ee.reshape([-1])
    ff = ff.reshape([-1])
    vv = vv.reshape([-1])
    np.testing.assert_almost_equal(ee, self.expected_e)
    np.testing.assert_almost_equal(ff, self.expected_f)
    np.testing.assert_almost_equal(vv, self.expected_v)