search

sblauth / cashocs

computational adjoint-based shape optimization and optimal control software for python
https://cashocs.readthedocs.io/
GNU General Public License v3.0
53 stars 10 forks source link

[Development] SNES Solver may not work properly for newer PETSc versions #462

Closed sblauth closed 4 months ago

sblauth commented 4 months ago

See the bug described here https://fenicsproject.discourse.group/t/issue-with-petsc-snes-and-lu-solver/15198/3

A possible remedy would be to follow the approach used in https://bitbucket.org/fenics-project/dolfin/src/1c52e837eb54cc34627f90bde254be4aa8a2ae17/dolfin/nls/PETScSNESSolver.cpp#lines-277 and to call the formJacobian method one additional time, which fixes the bug.