92
92
KrylovSolver solver(A, "gmres");
94
94
// Create vector that spans null space (normalised)
96
96
V.dofmap()->set(*null_space_ptr, sqrt(1.0/null_space_ptr->size()));
98
98
null_space_basis.push_back(null_space_ptr);
100
100
// Create null space basis object and attach to Krylov solver