]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix the assertion of MappingFEField's Jacobian determinant. 8828/head
authorDoug Shi-Dong <doug.shidong@gmail.com>
Sat, 21 Sep 2019 16:35:04 +0000 (12:35 -0400)
committerDoug Shi-Dong <doug.shidong@gmail.com>
Tue, 24 Sep 2019 02:18:45 +0000 (22:18 -0400)
Also add update flag to tests/mapping/mapping_fe_field_01.cc to check this assertion.

source/fe/mapping_fe_field.cc
tests/mappings/mapping_fe_field_01.cc

index dae76e0614f71569456e3839e069bd96869a0732..356f519d7694dbffa7610a05a29d864ec5e1cc5f 100644 (file)
@@ -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)
index 4e82d95c1234445a00355ce1a915227e3805f9df..c9e1a804286f6862c3cb370ab506f76339b4a048 100644 (file)
@@ -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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.