On Fri, Jun 21, 2019 at 9:36 AM 서승진 <tjtmdwls1...@gmail.com> wrote:
> I tried to use the smoother after reading your email, and it seems that > the LaplaceMeshSmoother is quitely similar with smooth function in > MeshTools::Modification namespace. > > After trying this, I found a new problem. > > As you said, boundary nodes are not moved. > > However, internal nodes escape from its domain, and as you said, it seems > that they ingnore the boundary nodes... > Are you still running the smoother for 1000 iterations as in the original code snippet you posted? That may be causing an issue, as I don't think I've ever tested Laplace smoothing for more than about 2 iterations. Note that even though boundary nodes are not moved, they should still be accounted for in the "averaging" process of the smoother, and therefore nodes should not be able to go outside the domain. Theoretically I think it's possible for Laplace mesh smoothing to send points outside the domain for non-convex geometries, though off the top of my head I can't really construct a case where this would happen for Quad4s... and I don't think you have a non-convex domain (?) so I'm not really clear on how your internal nodes are escaping. -- John _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users