I want so solve a system with a singular matrix (just the laplace discretized with the fem using pure Neumann boundary conditions) with BoomerAMG. Is there any way to do it directly in BoomerAMG/PETSc, or must I change the matrix to make it nonsingular? Thomas