glenco / SLsimLib

Library for Gravitational Lensing Simulations
MIT License
2 stars 1 forks source link

Time-delays #156

Open rbmetcalf opened 4 years ago

rbmetcalf commented 4 years ago

@EnziWolfgang wrote:

attached are two time delay maps created with GLAMER. One of them is created for the analytical LensHaloPowerLaw halo(1.e+13,1,zl,0.5,0.7,20,0,Pseudo); :

The other one I created with LensHaloParticles::makeSIE("testSIE",zl,1e8,1e13,500,0.7,ran); LensHaloParticles halo("testSIE",zl,64,cosmo,rot,false,false,0); :

I think there might be something odd going on with the time delay calculation of the LensHaloParticles case, because qualitatively they look quite different.

The code snippet is also in the attachments:

The results obtained for the analytic halo is closer to what I obtain when I calculate time delays for an SIE with my own code.