> I was trying to find if there is a method that can provide the
> reference to the interior element, given a side on the boundary.
> 
> I tried using interior_parent but it returns null for a given side.
> Looking at the code, it appears to me that there might be a bug in the
> implementation since for a level-0 element, the higher-dimensional
> interior element is only asserted and never found. So, parent will
> always be null for level-0 calls. This should also trigger a
> segmentation fault in debug mode, although I have not yet tested this.
> 
> I might be wrong on this but do let me know if there is any alternate
> way to find the interior element of a side.

I'll look into the interior_parent() issue, but in my application I

(1) build a lower dimensional mesh by extracting all the sides on the
boundaries of interest, and then simply call

(2) interior_elem = boundary_elem->parent()

Does that work for you?


------------------------------------------------------------------------------

_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to