From: Martin Kronbichler Date: Tue, 26 May 2009 10:10:02 +0000 (+0000) Subject: One trivial case where CellSimilarity can even be used for MappingQ (when degree... X-Git-Tag: v8.0.0~7661 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=78fdbb4e57234b81ca631dc021104a6e104ab048;p=dealii.git 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 --- 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