I have a couple of questions. Following the creation of a custom NewtonSolver as described here , is it possible to do a similar setup for a custom PETScSNESSolver.? Peeking at the underlying PETScSNESSolver.cpp I don’t see a solver_setup method (Am I missing anything)? In that case, would simply …

Hi, Did you manage to solve this problem? I am trying to setup a custom PETScSNESSolver where I supply the null space for the linear system. I am able to do this with a NewtonSolver using this code, however I have not managed to convert this to PETScSNESSolver. Cheers, K

Indeedn, it is much easier to directly deal with petsc4py.PETSc.SNES directly instead of the wrapper, and that is what I ended up doing.

Hi, Thanks for your reply. Currently within the NewtonSolver code the nullspace is attached in the solver block, but I’m not sure how to insert the nullspace within the snes solver. NewtonSolver syntax: class CustomSolver(NewtonSolver): def __init__(self, mesh, nullspace): NewtonSolv…

Creating a custom PETScSNESSolver

Linear Algebra

dokken April 24, 2023, 5:03pm 3

You could just use petsc4py directly, see Using petsc4py.PETSc.SNES directly - #12 by nate

1 Like

Topic		Replies	Views
Using petsc4py.PETSc.SNES directly Linear Algebra	23	4817	June 13, 2024
Implementing the Newton Method per hand Linear Algebra	5	3293	February 20, 2024
Issue with PETSc SNES and LU Solver General	5	226	July 11, 2024
Custom NewtonSolver using the mixed-dimensional branch and PETScNestMatrix	5	1246	April 3, 2020
PETSC solver Fenics 2018	2	1974	February 14, 2019

Creating a custom PETScSNESSolver

Related topics