Hi, I pulled and built libMesh after a long time (~2 years) recently and noticed that the new implementation of PetscDM does not work with DirichletBoundary class anymore. Anytime I try to set the dirichlet BCs for the problem when running with petscdm (i.e. using SNES types vinewtonrsls or vinewtonssls) it behaves as if the boundary value is zero! The problem disappears, if I not use PetscDM (use newtonls). The old implementation used to work with very well with dirichlet BCs and it’s such a shame if the new one cannot use that part of the package.
My questions is if this is a bug or am I supposed to do something differently for example when I assemble the residual and jacobian? Best, Ata _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users