These RDTs are computed there:
|
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.
These RDTs are computed there:
MAD-X/src/twiss.f90
Line 2524 in 88211ac
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.