Skip to content
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

Checking IDs of adjacent entities in the mixed mesh test #18

Open
joshia5 opened this issue Jul 1, 2021 · 0 comments
Open

Checking IDs of adjacent entities in the mixed mesh test #18

joshia5 opened this issue Jul 1, 2021 · 0 comments

Comments

@joshia5
Copy link
Collaborator

joshia5 commented Jul 1, 2021

In the unit test for mixed mesh, the IDs of adjacent entities returned
from query APIs are checked against known values.

omega_h/src/mixed_test.cpp

Lines 132 to 153 in 466e470

OMEGA_H_CHECK(tet2vtx.ab2b == LOs({0, 1, 8, 2}));
OMEGA_H_CHECK(hex2vtx.ab2b == LOs({4, 5, 9, 11, 7, 6, 1, 0}));
OMEGA_H_CHECK(wedge2vtx.ab2b == LOs({11, 9, 10, 0, 1, 8}));
OMEGA_H_CHECK(pyram2vtx.ab2b == LOs({8, 1, 9, 10, 3}));
OMEGA_H_CHECK(tet2edge.ab2b == LOs({0, 22, 20, 1, 3, 6}));
OMEGA_H_CHECK(hex2edge.ab2b == LOs({13, 9, 18, 12, 14, 15, 23, 19, 16, 5, 0, 2}));
OMEGA_H_CHECK(wedge2edge.ab2b == LOs({18, 17, 11, 19, 23, 21, 0, 22, 20}));
OMEGA_H_CHECK(pyram2edge.ab2b == LOs({22, 23, 17, 21, 7, 4, 8, 10}));
}
else if (num_vertex == 9) {
//pyram on hex
OMEGA_H_CHECK(hex2edge.ab2b == LOs({1, 7, 4, 0, 2, 10, 8, 5, 12, 13, 14, 15}));
OMEGA_H_CHECK(hex2vtx.ab2b == LOs({0, 3, 2, 1, 5, 8, 7, 6}));
OMEGA_H_CHECK(pyram2edge.ab2b == LOs({12, 13, 14, 15, 3, 11, 9, 6}));
OMEGA_H_CHECK(pyram2vtx.ab2b == LOs({5, 8, 7, 6, 4}));
}
else if (num_vertex == 7) {
//tet on wedge
OMEGA_H_CHECK(tet2edge.ab2b == LOs({9, 10, 11, 3, 8, 6}));
OMEGA_H_CHECK(tet2vtx.ab2b == LOs({4, 6, 5, 3}));
OMEGA_H_CHECK(wedge2edge.ab2b == LOs({1, 4, 0, 2, 7, 5, 9, 10, 11}));
OMEGA_H_CHECK(wedge2vtx.ab2b == LOs({0, 2, 1, 4, 6, 5}));

An open issue is to define a new OMEGA_H_VERIFY in Omega_h_fail.hpp which
prints the values of adjacent entity IDs incase of a mis-match, before failing.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant