emzed.align.mz_align_routines module

emzed.align.mz_align_routines.irls_fit(A, b, p, N=10)

IRLS iterations for minimizing ||Ax-b||_p

iterative weighting maxtrix is

w_ii = (eps + resid_i)^[(p-2)/2]

see also https://en.wikipedia.org/wiki/Iteratively_reweighted_least_squares