If we do an inexact nullspace projection, it's possible that we lose orthogonality in Krylov subspace vectors in the truncated conjugate direction algorithm. To compensate, we add an extra stopping criteria for when we lose orthgonality. We need a unit to test to verify this functionality.
If we do an inexact nullspace projection, it's possible that we lose orthogonality in Krylov subspace vectors in the truncated conjugate direction algorithm. To compensate, we add an extra stopping criteria for when we lose orthgonality. We need a unit to test to verify this functionality.