const EPS_L1: f64 = 1e-5;
ε in the L1 weights 1/sqrt(|wG·∇χ|² + ε) (chi_sep_MEDI_func’s eps__).
ε
1/sqrt(|wG·∇χ|² + ε)
chi_sep_MEDI_func
eps__