On Fri, Apr 24, 2015 at 1:13 AM, David Knezevic <[email protected]> wrote:
> On Fri, Apr 24, 2015 at 12:34 AM, Roy Stogner <[email protected]> > wrote: > >> >> On Fri, 24 Apr 2015, David Knezevic wrote: >> >> If anyone has any thoughts on what might be wrong here, I'd be interested >>> to hear. >>> >> >> Some interaction with the new multi-dimension mesh stuff? First the >> mesh does have 1D elements, so that's one of its active dimensions, >> then it doesn't, so that isn't. FEMContext tries to build FE objects >> for every active dimension and is some how "out of sync". >> > > > Yeah, I guess the issue is that the mesh still thinks it has 1D elements > even after they're all deleted. You can see this from the second > mesh.print_info() in the example I attached, which is called after the > deletion. It still says: > > Mesh Information: > elem_dimensions()={1, 3} > If I add _elem_dims.clear(); at the start of MeshBase::cache_elem_dims(), and I call cache_elem_dims() after the element deletion, then it works. Any objection to me adding _elem_dims.clear(); at the start of MeshBase::cache_elem_dims()? If not, then I'll make a PR for that. David ------------------------------------------------------------------------------ One dashboard for servers and applications across Physical-Virtual-Cloud Widest out-of-the-box monitoring support with 50+ applications Performance metrics, stats and reports that give you Actionable Insights Deep dive visibility with transaction tracing using APM Insight. http://ad.doubleclick.net/ddm/clk/290420510;117567292;y _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
