From 78fdbb4e57234b81ca631dc021104a6e104ab048 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Tue, 26 May 2009 10:10:02 +0000 Subject: [PATCH] One trivial case where CellSimilarity can even be used for MappingQ (when degree is one). git-svn-id: https://svn.dealii.org/trunk@18876 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/mapping_q.cc | 3 ++- deal.II/deal.II/source/fe/mapping_q1.cc | 10 ++++------ 2 files changed, 6 insertions(+), 7 deletions(-) diff --git a/deal.II/deal.II/source/fe/mapping_q.cc b/deal.II/deal.II/source/fe/mapping_q.cc index 0cb43f683a..6aea6118d9 100644 --- a/deal.II/deal.II/source/fe/mapping_q.cc +++ b/deal.II/deal.II/source/fe/mapping_q.cc @@ -351,7 +351,8 @@ MappingQ::fill_fe_values ( else { p_data=&data; - cell_similarity = CellSimilarity::invalid_next_cell; + if (get_degree() > 1) + cell_similarity = CellSimilarity::invalid_next_cell; } MappingQ1::fill_fe_values(cell, q, *p_data, diff --git a/deal.II/deal.II/source/fe/mapping_q1.cc b/deal.II/deal.II/source/fe/mapping_q1.cc index 9309dca69f..b2f07faa87 100644 --- a/deal.II/deal.II/source/fe/mapping_q1.cc +++ b/deal.II/deal.II/source/fe/mapping_q1.cc @@ -718,7 +718,6 @@ MappingQ1::compute_fill (const typename Triangulation::fill_fe_values ( Assert( !(update_flags & update_cell_normal_vectors ) || (cell_normal_vectors.size() == n_q_points), ExcDimensionMismatch(cell_normal_vectors.size(), n_q_points)); + if (cell_similarity != CellSimilarity::translation) for (unsigned int point=0; point::fill_fe_values ( Assert (jacobians.size() == n_q_points, ExcDimensionMismatch(jacobians.size(), n_q_points)); if (cell_similarity != CellSimilarity::translation) - for (unsigned int point=0; point::fill_fe_values ( ExcDimensionMismatch(inverse_jacobians.size(), n_q_points)); if (cell_similarity != CellSimilarity::translation) for (unsigned int point=0; point