offset_i,
ad_entity_types[offset_i]
.template get_combined_orientation<unsigned int>(
- make_array_view(ref_indices.begin(),
- ref_indices.begin() +
- ad_entity_types[offset_i].n_vertices()),
make_array_view(ad_entity_vertices[offset_i].begin(),
ad_entity_vertices[offset_i].begin() +
+ ad_entity_types[offset_i].n_vertices()),
+ make_array_view(ref_indices.begin(),
+ ref_indices.begin() +
ad_entity_types[offset_i].n_vertices())));
}
col_d[offset_i] = counter;
// face_orientation=true, face_rotation=true, face_flip=false
if (v0_equals({vertices_1[1], vertices_1[2], vertices_1[0]}))
- return 3;
+ return 5;
// face_orientation=true, face_rotation=false, face_flip=true
if (v0_equals({vertices_1[2], vertices_1[0], vertices_1[1]}))
- return 5;
+ return 3;
// face_orientation=false, face_rotation=false, face_flip=false
if (v0_equals({vertices_1[0], vertices_1[2], vertices_1[1]}))
CELLS 3 21
8 0 1 3 2 4 5 7 6
5 4 5 7 6 8
-5 3 7 5 1 9
+5 5 1 3 7 9
CELL_TYPES 3
12 14 14
// check orientation
deallog << "face orientation: " << orientation << ' '
<< int(cell->combined_face_orientation(0)) << ' ' << std::endl;
- AssertDimension(orientation,
- (cell->face_orientation(0) * 1 + cell->face_rotation(0) * 2 +
- cell->face_flip(0) * 4));
+ // AssertDimension(orientation,
+ // (cell->face_orientation(0) * 1 + cell->face_rotation(0) * 2
+ // +
+ // cell->face_flip(0) * 4));
// check vertices
{
const auto permuted = type.permute_according_orientation(origin, o);
const unsigned int origin_back =
- type.compute_orientation(permuted, origin);
+ type.compute_orientation(origin, permuted);
// AssertDimension(o, origin_back);
deallog << o << " " << origin_back << std::endl;
permuted[permutedi[i]] = origin[i];
const unsigned int origin_back =
- type.compute_orientation(origin, permuted);
+ type.compute_orientation(permuted, origin);
// AssertDimension(o, origin_back);
deallog << o << " " << origin_back << std::endl;