Skip to content

Normalize vectors passed to MatNullSpaceCreate()

Jose E. Roman requested to merge jose/maint/normalize-nullspace into maint

In GHEP, basis vectors are normalized with respect to the B-norm, including the constraints from the deflation space. Hence, we must normalize them explicitly so that MatNullSpaceCreate() does not complain.

@prj- I tried to adapt one of the test examples to reproduce the issue without the data files, but I always get 'Zero pivot in LU factorization'. Regarding the other issue, it appears only in master, but the second commit here should fix it.

Edited by Jose E. Roman

Merge request reports