Re: [Libmesh-users] reinitializing on refined mesh

2016-08-02 Thread Roy Stogner
On Tue, 2 Aug 2016, Harriet Li wrote: > MeshRefinement::refine_elements() to refine them, like in Misc Ex 4. > However, when I call EquationSystems::reinit(), I get > " > Assertion `!libmesh_assert_pass' failed. > " > with the backtrace > " > #0 libMesh::MeshRefinement::test_level_one (this=0x7f

[Libmesh-users] reinitializing on refined mesh

2016-08-02 Thread Harriet Li
Hello, I'm trying to solve a problem on a uniform mesh, refine part of it, and then solve the problem again on the new mesh, using the previous solution as an initial guess. To that end, I set FEMSystem::project_solution_on_reinit() = true and call EquationSystems::reinit() after each mesh refinem