Normalize vectors passed to MatNullSpaceCreate()
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