From: David Wells Date: Mon, 19 Jul 2021 20:53:52 +0000 (-0400) Subject: Use the correct reference cell in MappingFEField. X-Git-Tag: v9.4.0-rc1~1133^2~3 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=cc762798476b18728b5b9e3ecd0318472c15e6b8;p=dealii.git Use the correct reference cell in MappingFEField. --- diff --git a/source/fe/mapping_fe_field.cc b/source/fe/mapping_fe_field.cc index 5d0b581f68..b23c21fb18 100644 --- a/source/fe/mapping_fe_field.cc +++ b/source/fe/mapping_fe_field.cc @@ -638,8 +638,7 @@ MappingFEField::get_face_data( std::make_unique(euler_dof_handler->get_fe(), fe_mask); auto & data = dynamic_cast(*data_ptr); const Quadrature q( - QProjector::project_to_all_faces(ReferenceCells::get_hypercube(), - quadrature[0])); + QProjector::project_to_all_faces(reference_cell, quadrature[0])); this->compute_face_data(update_flags, q, quadrature[0].size(), data); return data_ptr; @@ -656,7 +655,7 @@ MappingFEField::get_subface_data( std::make_unique(euler_dof_handler->get_fe(), fe_mask); auto & data = dynamic_cast(*data_ptr); const Quadrature q(QProjector::project_to_all_subfaces( - ReferenceCells::get_hypercube(), quadrature)); + reference_cell, quadrature)); this->compute_face_data(update_flags, q, quadrature.size(), data); return data_ptr; @@ -1720,7 +1719,7 @@ MappingFEField::fill_fe_face_values( face_no, numbers::invalid_unsigned_int, QProjector::DataSetDescriptor::face( - ReferenceCells::get_hypercube(), + reference_cell, face_no, cell->face_orientation(face_no), cell->face_flip(face_no), @@ -1761,7 +1760,7 @@ MappingFEField::fill_fe_subface_values( face_no, numbers::invalid_unsigned_int, QProjector::DataSetDescriptor::subface( - ReferenceCells::get_hypercube(), + reference_cell, face_no, subface_no, cell->face_orientation(face_no),