pmelchior / skylens Goto Github PK

need the shapelet library for complex galaxy morphologies. requires port from numla to tmv for vector/matrix operations.

Switching on OpenMP works fine with the new boost::geometry tree implementation, but the LensingLayer creates a race condition that leads to horizontal streaks of positive flux from each GalaxyModel.

Moreover, when several GalaxyLayers are present, the iterative nature of the switch on/off for each lens layer (needed to compute the correct angle for the entire background layer) is not parallelizable. Solution unclear!

filter, sed, and cosmology should be extracted from libastrocpp and made available. Ideally, the CCfits dependency could be removed entirely by replacing it with cfitsio functions.

We need to write a cluster member class ported from the f90 version of sky lens

SQLite certainly needed, WCS not really. If the libraries are needed, remove the compilation flags, e.g. -DHAS_SQLiteDB from Makefile and code segments.

If the deflection angle field in LensLayer is smaller than the requested FoV, the fluxes at the boundary pixels are discontinuous. The only way to make sure this doesn't happen is to use large enough maps for the deflection angles, but then they are getting huge, while they contain very little information on such large scales.

Some physically motivated way to extrapolate would be advised.

All units are currently calculated from left-lower corner, which requires POS_X/Y in lens config and puts the reference also in that corner in the resulting FITS file. Also, in the global config file, POINTING_X/Y refers to the ll corner.

It's more natural if the origin was in the center of the FoV, deviations could still be done simply with the POS_X/Y and POINTING_X/Y parameters.

Requires:

• optional reading of POS_X/Y in LensLayer
• careful check for consistency in programs, Observation and potentially some other classes.