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 dae76e0614f71569456e3839e069bd96869a0732..356f519d7694dbffa7610a05a29d864ec5e1cc5f 100644 (file)
--- 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 4e82d95c1234445a00355ce1a915227e3805f9df..c9e1a804286f6862c3cb370ab506f76339b4a048 100644 (file)
--- 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<dim, spacedim> fe_values(map_fe,
                                     fe_sys,
                                     quadrature_formula,
-                                    update_values | update_JxW_values);
+                                    update_values | update_JxW_values |
+                                      update_volume_elements);
 
   typename DoFHandler<dim, spacedim>::active_cell_iterator cell =
                                                              dof_sys