From: Doug Shi-Dong Date: Sat, 21 Sep 2019 16:35:04 +0000 (-0400) Subject: Fix the assertion of MappingFEField's Jacobian determinant. X-Git-Tag: v9.2.0-rc1~1042^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F8828%2Fhead;p=dealii.git Fix the assertion of MappingFEField's Jacobian determinant. Also add update flag to tests/mapping/mapping_fe_field_01.cc to check this assertion. --- diff --git a/source/fe/mapping_fe_field.cc b/source/fe/mapping_fe_field.cc index dae76e0614..356f519d76 100644 --- a/source/fe/mapping_fe_field.cc +++ b/source/fe/mapping_fe_field.cc @@ -838,7 +838,8 @@ namespace internal if (update_flags & update_volume_elements) { - AssertDimension(data.covariant.size(), data.volume_elements.size()); + AssertDimension(data.contravariant.size(), + data.volume_elements.size()); if (cell_similarity != CellSimilarity::translation) for (unsigned int point = 0; point < data.contravariant.size(); ++point) diff --git a/tests/mappings/mapping_fe_field_01.cc b/tests/mappings/mapping_fe_field_01.cc index 4e82d95c12..c9e1a80428 100644 --- a/tests/mappings/mapping_fe_field_01.cc +++ b/tests/mappings/mapping_fe_field_01.cc @@ -77,7 +77,8 @@ test(const unsigned int degree) FEValues fe_values(map_fe, fe_sys, quadrature_formula, - update_values | update_JxW_values); + update_values | update_JxW_values | + update_volume_elements); typename DoFHandler::active_cell_iterator cell = dof_sys