Adapted code to use the '.solve' function of the N-matrix objects. this allows using it with PTA models containing ECORR noise
Former code calculated and subtracted the object XNTSigmaTNy twice, leading to too small values of the Fe-statistic. This mistake is solved. Now the code matches the calculations outlined in Ellis, Siemens, Creighton 2012.
Note: The updated version is an EPTA contribution to IPTA DR3. It is used in EPTA and MPTA data analyses.
Adapted code to use the '.solve' function of the
N-matrix objects. this allows using it with PTA models containingECORRnoiseFormer code calculated and subtracted the object
XNTSigmaTNytwice, leading to too small values of the Fe-statistic. This mistake is solved. Now the code matches the calculations outlined in Ellis, Siemens, Creighton 2012.Note: The updated version is an EPTA contribution to IPTA DR3. It is used in EPTA and MPTA data analyses.