User-defined near null space about algebraicmultigrid.jl HOT 3 OPEN

Thank you for your kind words and for using the package!

Yes, I can have definitely have that be an optional argument with default as ones(T,n). Feel free to send a PR. It should be a quick change. I'll merge as soon as tests pass. Thank you!

from algebraicmultigrid.jl.

Yes, I can have definitely have that be an optional argument with default as ones(T,n). Feel free to send a PR. It should be a quick change.

Perhaps this is not sufficient. In general, the near null space can have several vectors, so size(B,2) can be > 1. I don't know if the code will support the case when B is a matrix. I'll try once I have some time.

from algebraicmultigrid.jl.

Came across the same issue. Note that pyamg.aggregation.smoothed_aggregation_solver allows near-nullspace candidates B of shape (N, K). As long as the SA setup logic follows PyAMG, passing a 2-D matrix B should be fine.

from algebraicmultigrid.jl.