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

Reply via email to