Like the one they have in PyAMG, it should use some sort of stencil to define a banded matrix.
Like the one they have in PyAMG, it should use some sort of stencil to define a banded matrix.