masksolve.py - Solve a masked linear problem quickly

everest.masksolve.MaskSolve(A, b, w=5, progress=True, niter=None)

Finds the solution x to the linear problem

A x = b

for all contiguous w-sized masks applied to the rows and columns of A and to the entries of b.

Returns an array X of shape (N - w + 1, N - w), where the nth row is the solution to the equation

A[![n,n+w)] x = b[![n,n+w)]

where ![n,n+w) indicates that indices in the range [n,n+w) have been masked.