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)