Composite image test. Here we will produce a composite lens as a sum of three lenses with the same velocity dispersions and different centroids.
In [1]:
%matplotlib inline
import evillens as evil
In [38]:
lens1 = evil.AnalyticSIELens(0.4,1.5)
lens2 = evil.AnalyticSIELens(0.4,1.5)
lens3 = evil.AnalyticSIELens(0.4,1.5)
lens1.build_kappa_map(sigma=200.0, r_c=0.1)
lens2.build_kappa_map(sigma=170.0, q = 0.45, centroid=[-1.51,0.51])
lens3.build_kappa_map(sigma =100.0, centroid=[1.01,1.01])
In [39]:
lens4 = lens1+lens2
lens5 = lens3+lens4
lens5.plot('kappa')
In [40]:
lens5.deflect()
In [41]:
lens5.source = evil.Source(1.5)
lens5.source.Read_Source_from("analytic_SIE_kappa.fits")
In [42]:
lens5.raytrace()
In [43]:
lens5.plot('non-lensed image')
In [44]:
lens5.plot('lensed image')
In [8]:
In [ ]: