search

cmlibs / zinc

Source code repository for OpenCMISS-Zinc
Mozilla Public License 2.0
15 stars 18 forks source link

Define faces failing with general linear map #104

Open rchristie opened 6 years ago

rchristie commented 6 years ago

Limitation is in function global_to_element_map_nodes which only reports the first node in a general linear map. This results in some faces using general linear maps being treated as collapsed, especially some cases of hanging nodes, or in other cases not matching each other correctly. A workaround that fixes some cases is to mix the order of terms in each general linear map to ensure all nodes are reported in the above function, and hence will match faces.

Note that the find faces algorithm currently assumes only that faces with the same nodes are the same, but in reality it should match the expression for the coordinate field on the face.

rchristie commented 1 year ago

This is limiting some uses of general linear maps for more efficient elements with fewer degrees of freedom as surfaces and lines are not correctly set up nor able to be drawn.