Actually, this issue is now resolved. The problem was that the default linear solver was not converging. This issue can be resolved by using a stronger linear solver.
Thank you, Mike > On Apr 25, 2018, at 2:56 PM, Lee, Jae Ho <jaeho...@live.unc.edu> wrote: > > I tried the subdivision example provided in libMesh, and I had some questions > about the example that I wasn’t able to find an answer to by myself. > > In the example (that follows Cirak’s paper from 2000), I see that there are > the following steps before the subdivision mesh gets prepared: > > // Quadrisect the mesh triangles a few times to obtain a > > // finer mesh. Subdivision surface elements require the > > // refinement data to be removed afterward. > > MeshRefinement mesh_refinement(mesh); > > mesh_refinement.uniformly_refine(3); > > MeshTools::Modification::flatten(mesh); > > > Here, what’s the exact effect of the mesh refinement? I just tried increasing > the refinement to 5 or 6, and the solution started deviating quite a lot (by > orders of magnitudes) from the analytical solution. I was just wondering if > it’s due to just the nature of subdivision or due to the way libMesh > works/handles it? Is there a reason 3 is chosen for refinement? > > Thank you, > Mike > ------------------------------------------------------------------------------ > Check out the vibrant tech community on one of the world's most > engaging tech sites, Slashdot.org! http://sdm.link/slashdot > _______________________________________________ > Libmesh-users mailing list > Libmesh-users@lists.sourceforge.net > https://lists.sourceforge.net/lists/listinfo/libmesh-users ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
