dev-aym branch

[aymgal]

That's the current status of my fork, mostly with fixes for Pyhon 3.

Below are some questions about the code.

[aymgal]

• This line in Solve.SLIT(...) (also in SLIT_MCA(...)), why is the argument PSFconj overwritten here ?
  

  SLIT/SLIT/Solve.py

  Line 53 in ce591d5

  PSFconj = PSF.T

• Can you explain in more details the behaviour of this soft thresholding operator Soft(...) ? it seems to do special things at specific starlet scales ?
  

  SLIT/SLIT/tools.py

  Line 128 in ce591d5

  def Soft(X, level, k, supp =1, Kill = 0):

[aymgal]

Thanks 😉
Extra question : why do you normalise by onelens here ?

SLIT/SLIT/Lens.py

Line 176 in ce591d5

return Image/onelens

[herjy]

To answer your questions in order:

• PSFconj should be removed from the arguments list.
• In soft thresholding, I raise the threshold by 1 in the first scale. That's more of a familly recipe, but you will notice that in astronomical images, you will want to cut extra power in the high frequencies.
• If I recal, this is used to compute the noise levels only. onelens will be set to 1 in all other cases.