glenco / SLsimLib

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

Problem with making a PixelMap from a fits file #123

Closed rbmetcalf closed 7 years ago

rbmetcalf commented 7 years ago

From Giulia:

I didn’t have time last week to go on with those test with my source, but I looked into it in the last few days and I think I got why we could not get an image of the source. It seems that in pixelize.cpp when you define a pixel map from a fits file as

PixelMap::PixelMap( std::string fitsfilename /// file name of fits file to be read ,double my_res /// resolution (rad) of fits image if not given in fits file, use default or -1 otherwise ) {

then the variable “resolution” is only set equal to my_res for the case when the resolution has to be red from the fits file and is not given by hand. So then it was =0 when it computed rangeX and rangeY. Adding a resolution=my_res at the beginning now I can get an image of the source using the dummy lens!