f1001 and f1010 in twiss table

[ldeniau]

These RDTs are computed there:

MAD-X/src/twiss.f90

Line 2524 in 88211ac

f1001 = cmplx((Ctmp(1,2) - Ctmp(2,1)), Ctmp(1,1) + Ctmp(2,2))/ 4d0 / gamma

And it would be useful to have them in the twiss table for coupling calculation...

Alternatively, it would be more MAD-X oriented to save the coupling matrix Cbar (aka Ctmp in the code), since f1001 and f1010 can be easily recalculated from it.