def test_flexion(self):
lensModel = LensModel(lens_model_list=['FLEXION'])
g1, g2, g3, g4 = 0.01, 0.02, 0.03, 0.04
kwargs = [{'g1': g1, 'g2': g2, 'g3': g3, 'g4': g4}]
f_xxx, f_xxy, f_xyy, f_yyy = lensModel.flexion(x=1., y=1., kwargs=kwargs, diff=0.0001)
npt.assert_almost_equal(f_xxx, g1, decimal=2)
npt.assert_almost_equal(f_xxy, g2, decimal=2)
npt.assert_almost_equal(f_xyy, g3, decimal=2)
npt.assert_almost_equal(f_yyy, g4, decimal=2)
f_xxx, f_xxy, f_xyy, f_yyy = lensModel.flexion(x=1., y=1., kwargs=kwargs, diff=0.0001, hessian_diff=0.001)
npt.assert_almost_equal(f_xxx, g1, decimal=2)
npt.assert_almost_equal(f_xxy, g2, decimal=2)
npt.assert_almost_equal(f_xyy, g3, decimal=2)
npt.assert_almost_equal(f_yyy, g4, decimal=2)
def test_flexion(self):
x = np.array(0)
y = np.array(2)
flex = LensModel(['FLEXION'])
f_xxx, f_xxy, f_xyy, f_yyy = flex.flexion(x, y, [self.kwargs_lens])
npt.assert_almost_equal(f_xxx, self.kwargs_lens['g1'], decimal=9)
npt.assert_almost_equal(f_xxy, self.kwargs_lens['g2'], decimal=9)
npt.assert_almost_equal(f_xyy, self.kwargs_lens['g3'], decimal=9)
npt.assert_almost_equal(f_yyy, self.kwargs_lens['g4'], decimal=9)
def test_flexion(self):
x = np.array(0)
y = np.array(2)
flex = LensModel(['FLEXIONFG'])
f_xxx, f_xxy, f_xyy, f_yyy = flex.flexion(x, y, [self.kwargs_lens])
_g1, _g2, _g3, _g4 = self.flex.transform_fg(**self.kwargs_lens)
npt.assert_almost_equal(f_xxx, _g1, decimal=9)
npt.assert_almost_equal(f_xxy, _g2, decimal=9)
npt.assert_almost_equal(f_xyy, _g3, decimal=9)
npt.assert_almost_equal(f_yyy, _g4, decimal=9)
