-
Notifications
You must be signed in to change notification settings - Fork 54
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Issue with t8_forest_leaf_face_neighbors() #1220
Comments
Got it! Wait for your good news! |
Hey @CFDER-HF, currently, most of us are on holiday and we had to postpone the discussion to Monday of the next week. Greetings, |
Hey @CFDER-HF, do you have a specific tree boundary, at which these problems appear in the example you have given? Also, could you clarify what you mean by not correct? Ideally, I would want to know for which refinement level you did this, and for which element the neighbors are not what you expected. Even better if you can tell me what you expected. Thank you for your help! |
Just putting down my first observation here: The cmesh that you provided seems to have all element vertices in an order, that does not fit to the expected orientation for t8_code cmesh elements.
Thus, the t8code internal enumeration of vertices will not correspond to your ordering of vertices, and the t8code internal ordering of faces might also not correspond to your ordering of faces. That might already explain the discrepancies that you are experiencing, but I'm looking forward to getting more detail |
Are you able to reorient the original mesh in gmsh and try it again? You could also try commenting out this part:
in the gmsh4 file reader. Then, all of your refined elements would be oriented negatively, but the face-neighbor relation should work. |
If the problem at the quad tri boundaries still occurs after reorienting, it might be that your mesh from gmsh is not connected at these boundaries, so the boundary vertices are duplicated and t8code cannot identify the face connection. Could you check at the left right border between quads and tris? You might have the gmsh connectivity there. The edges need to be the same transfinite line, so probably you have to split the upper and lower rectangular element into 3 parts, so that you can connect the middle parts to the triangle part. |
I have changed the orientation in GMSH and everything is going as expected. Thanks for your enthusiastic help! |
Just read the conversation here. @lukasdreyer: Is there now still a bug in t8code, i.e. the reorientation? |
Bug report
Describe the bug
I am working on flow past a cylinder, and following is the coarse mesh I use which is created by Gmsh.
I've noticed that the face neighbors returned by
t8_forest_leaf_face_neighbors()
appear to be incorrect when the neighboring cells belong to differenttree_id
. I'm unsure whether this issue stems from a potential bug in t8code related to the imported hybrid mesh file, or if I should be using different interface functions for this specific scenario. Could you clarify this for me?cylinder_hybrid1.msh.zip
To Reproduce
Which command did you run/which function did you call?
t8_forest_leaf_face_neighbors (problem->forest,itree,elem,&neighs,iface,&elem_data->dual_faces[iface],
&elem_data->num_neighbors[iface],&elem_data->neighs[iface],&neigh_scheme,1);
Estimated priority
Which of these is most applicable (remove the others):
"Priority: high" Should be solved as soon as possible
The text was updated successfully, but these errors were encountered: