def METHOD_NAME(self):
    for arr, shp in zip(self.w0s, self.hessian_return_shapes):
        if self.E_unit.is_equivalent(u.one) and hasattr(arr, 'pos') and \
                not arr.xyz.unit.is_equivalent(u.one):
            continue
        g = self.obj.hessian(arr)
        assert g.shape == shp
        # TODO: check return units