I have a linear system of equations that I using fenics to solve in parallel with MPI. The solver works but when I increase the number of cores I begin to get an error:
*** Error: Unable to solve linear system using PETSc Krylov solver.
*** Reason: Solution failed to converge in 0 iterations (PETSc reason DIVERGED_PCSETUP_FAILED, residual norm ||r|| = 0.000000e+00).
*** Where: This error was encountered inside PETScKrylovSolver.cpp.
*** Process: 7
*** DOLFIN version: 2017.1.0
The solver works fine with one core so I am assuming it is an issue with how MPI splits up the problem. Any help would be greatly appreciated!