def test_raise(self): with self.assertRaises(ValueError): analysis = LensAnalysis(kwargs_model={'lens_model_list': ['SIS']}) analysis.multi_gaussian_lens(kwargs_lens=[{'theta_E'}]) with self.assertRaises(ValueError): analysis = LensAnalysis( kwargs_model={'lens_light_model_list': ['GAUSSIAN']}) analysis.flux_components(kwargs_light=[{}], n_grid=400, delta_grid=0.01, deltaPix=1., type="wrong")
def test_flux_components(self): phi, q = -0.37221683730659516, 0.70799587973181288 e1, e2 = param_util.phi_q2_ellipticity(phi, q) phi2, q2 = 0.14944144075912402, 0.4105628122365978 e12, e22 = param_util.phi_q2_ellipticity(phi2, q2) kwargs_profile = [{ 'Rs': 0.16350224766074103, 'e1': e12, 'e2': e22, 'center_x': -0.019983826426838536, 'center_y': 0.90000011282957304, 'amp': 1.3168943578511678 }, { 'Rs': 0.29187068596715743, 'e1': e1, 'e2': e2, 'center_x': 0.020568531548241405, 'center_y': 0.036038490364800925, 'Ra': 0.020000382843298824, 'amp': 85.948773973262391 }] kwargs_options = { 'lens_model_list': ['SPEP'], 'lens_model_internal_bool': [True], 'lens_light_model_internal_bool': [True, True], 'lens_light_model_list': ['HERNQUIST_ELLIPSE', 'PJAFFE_ELLIPSE'] } lensAnalysis = LensAnalysis(kwargs_options) flux_list, R_h_list = lensAnalysis.flux_components(kwargs_profile, n_grid=400, delta_grid=0.01, deltaPix=1., type="lens") assert len(flux_list) == 2 npt.assert_almost_equal(flux_list[0], 0.23898248741810812, decimal=8) npt.assert_almost_equal(flux_list[1], 3.0565768930826662, decimal=8) kwargs_profile = [{'amp': 1.}] kwargs_options = { 'source_light_model_list': ['UNIFORM'], 'lens_model_list': [] } lensAnalysis = LensAnalysis(kwargs_options) flux_list, R_h_list = lensAnalysis.flux_components(kwargs_profile, n_grid=400, delta_grid=0.01, deltaPix=1., type="source") assert len(flux_list) == 1 npt.assert_almost_equal(flux_list[0], 16, decimal=8)
def test_flux_components(self): kwargs_profile = [{ 'Rs': 0.16350224766074103, 'q': 0.4105628122365978, 'center_x': -0.019983826426838536, 'center_y': 0.90000011282957304, 'phi_G': 0.14944144075912402, 'sigma0': 1.3168943578511678 }, { 'Rs': 0.29187068596715743, 'q': 0.70799587973181288, 'center_x': -0.01, 'center_y': 0.9, 'Ra': 0.020000382843298824, 'phi_G': -0.37221683730659516, 'sigma0': 85.948773973262391 }] kwargs_options = { 'lens_model_list': ['SPEP'], 'lens_model_internal_bool': [True], 'lens_light_model_internal_bool': [True, True], 'lens_light_model_list': ['HERNQUIST_ELLIPSE', 'PJAFFE_ELLIPSE'] } lensAnalysis = LensAnalysis(kwargs_options) flux_list, R_h_list = lensAnalysis.flux_components(kwargs_profile, n_grid=400, delta_grid=0.01, deltaPix=1., type="lens") assert len(flux_list) == 2 npt.assert_almost_equal(flux_list[0], 0.23898248741810812, decimal=8) npt.assert_almost_equal(flux_list[1], 3.0565768930826662, decimal=8) kwargs_profile = [{'mean': 1.}] kwargs_options = { 'lens_light_model_list': ['UNIFORM'], 'lens_model_list': ['NONE'] } lensAnalysis = LensAnalysis(kwargs_options) flux_list, R_h_list = lensAnalysis.flux_components(kwargs_profile, n_grid=400, delta_grid=0.01, deltaPix=1., type="lens") assert len(flux_list) == 1 npt.assert_almost_equal(flux_list[0], 16, decimal=8)