Hi Corrado and FEniCS list, I'm just looking onto option of using semi-smooth Newton. It seems being implemented as SNESVINEWTONSSLS [1] in PETSc. As I understand it, it should work without specifying any bounds (letting them +/-INF). Nevertheless it raises
[0]PETSC ERROR: Object is in wrong state! [0]PETSC ERROR: This DM is incapable of computing variable bounds.! when calling solve() without bounds. Is it possible bug in PETScSNESSolver.cpp? Tested with PETSc 3.4.0. [1]http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/SNESVINEWTONSSLS.html Jan _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
