fsmp.py: 减去 A_vec 和 c_vec 的平均值，等价变换。 (!11) · Merge requests · XMASS / waveform

Benda Xu requested to merge beta-condition into master Aug 12, 2022

设 A_vec (简记 a) 和 c_vec (简记 c) 是 2 x l_wave 的二重向量。

当 A_vec[0] == A_vec[1] ，L = diag((-1, 1)) 时，a L a^T == 0。故可以给 A_vec 和 c_vec 的两个分量同时减去它们的平均值，使得总的二范数极小， L^{-1} + a^T c 的行列式尽可能由满秩的 L 主导。

Reference: metropolis_waveform (22) (23) 式，L 对应 \Lambda。 Closes: #32 (closed)

fsmp.py: 减去 A_vec 和 c_vec 的平均值，等价变换。