From: wolf Date: Mon, 29 Sep 2003 17:33:59 +0000 (+0000) Subject: Reinstate previous order. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d26b0f2b2fe54263aab05d8d98eb5a4977728166;p=dealii-svn.git Reinstate previous order. git-svn-id: https://svn.dealii.org/trunk@8058 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/mapping_q1.cc b/deal.II/deal.II/source/fe/mapping_q1.cc index da3643d8a9..78823cd27a 100644 --- a/deal.II/deal.II/source/fe/mapping_q1.cc +++ b/deal.II/deal.II/source/fe/mapping_q1.cc @@ -467,12 +467,18 @@ MappingQ1::compute_fill (const typename DoFHandler::cell_iterator &cel // if necessary, recompute the // support points of the // transformation of this cell - if ((cell != data.cell_of_current_support_points) + // (note that we need to first + // check the triangulation pointer, + // since otherwise the second test + // might trigger an exception if + // the triangulations are not the + // same) + if ((data.mapping_support_points.size() == 0) || (&cell->get_triangulation() != &data.cell_of_current_support_points->get_triangulation()) || - (data.mapping_support_points.size() == 0)) + (cell != data.cell_of_current_support_points)) { compute_mapping_support_points(cell, data.mapping_support_points); data.cell_of_current_support_points = cell;