LEI: Run your code with option '-snes_converged_reason -ksp_converged_reason -ksp_monitor -snes_monitor' to check which solver leads to divergence. Then use option '-snes_view' to check what solver is used. Hong
> Dear all > > > > I am using SNES for solving a large nonlinear system for Material Point > Method. But it always diverged after server steps, I think there are some > bug in the computation of Jacobian matrix. Do anyone have practical > approach for identified these bugs in a large matrix? Currently, I saved > the matrix to a txt file and load it to Matlab for checking. Is there > function for checking Petsc Mat in Petsc itself? > > > > Cheers > > Lei >
