On Mon, May 16, 2011 at 16:50, Tian(ICT) <rongtian at ncic.ac.cn> wrote:
> ./pgfem -pc_type none does not work either, leading to negative jacobian of > element. > During finite differencing or in the line search? Perhaps a trust region would be better for this problem. What is the output for -mat_view, are all the output the jacobian, or do > they inlcude preconditioning matrix? > It includes all matrices. You can call MatSetOptionsPrefix(B,"jacpre_") and then run with -jacpre_mat_view to see only that matrix. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110516/205d7a39/attachment.htm>
