Lenstool

I agree, but do I understand correctly that there are many loops over I.n_mult (using it as the number of multiple images)?
It seems easier to read this integer as a different variable name in r_image.c

In any case currently the code accepts :
0 = nothing
any number >0 but different from 2 : usual multiple image format
2 : "shape" multiple image format.

We should only accept 0, 1, and 2 and otherwise return a warning (format not supported) just for consistency.

Sorry there are kind of 2 issues mixed together (and also 3 if we count the documentation).

Thanks!

Johan

The format definition in f_shape is * flag=1 : ID RA DEC A B THETA Z MAG (default) * flag=2 : ID RA DEC A B THETA Z MAG VARE1 VARE2 (for WL or source with Sersic index) * flag=3 : ID RA DEC E1 E2 Z VARE1 VARE2 (for WL)

For multiple images, if n_mult=0, nothing happens, if n_mult=2, then we use flag=2 format, otherwise we use flag=1 format.

In any case, only n_mult systems are used.

It seems natural to me to specify the format of the input catalog in front of the catalog name, as an argument of the multifile keyword. Then, we can add another keyword in the image section with the number of multiple systems to take into account (the actual n_mult variable)

What do you think?

But flag=3 is (I believe) not coded? Only flag=2 is used (in o_global) for the format.

I think we agree that we need two variable names. I don't think specifying the n_mult in the parameter
file is very useful, so just storing the "file format integer" in a new variable is necessary, and also
if we want to have flag=3 used it needs to be coded?

flag=3 is coded, but not used for SL (it is used for WL).

So, shall we create a new variable in the image section, let's say I.multmode that contains the format of the input catalog?

Already, for WL we have I.statmode that contains the arclet catalog format.