How to implement problem with only Neumann BCs, no Dirichlet

Simply send in an empty list bcs=[]
or not send in any argument, i.e. fem.petsc.NonlinearProblem(F, N) as it has a default value:

1 Like