Hi all, I use PETSC's nonelinear solver to do emiconductor simulation. I found a strange problem with trust region solvers. It usually works well if LU method is used for inner linear solver. But when krylov method i.e. GMRES, is employed, it often stops due to pnorm is small, but still remains large error of function norm. How to solve this problem? LU method is much slower than i.e. GMRES.
G.D. -- USTC Alumni Email System
