Remove a check that was not only long satisfied, but also was in a function

(2d case) where it would never trigger (because it checked for dim==3).

git-svn-id: https://svn.dealii.org/trunk@9781 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/deal.II/source/fe/fe_q.cc b/deal.II/deal.II/source/fe/fe_q.cc
index f46af800504c11b8164b9d760a705df80c58a084..85a61d26e54021f712236a43ccecd85791c41168 100644 (file)
--- a/deal.II/deal.II/source/fe/fe_q.cc
+++ b/deal.II/deal.II/source/fe/fe_q.cc
@@ -1215,25 +1215,6 @@ FE_Q<2>::initialize_constraints ()
     .TableBase<2,double>::reinit (this->interface_constraints_size());
   subface_interpolation.Tmmult (this->interface_constraints,
                                face_interpolation);
-
-                                  // in 3d we still have the
-                                  // constraint matrices, so make the
-                                  // check
-//TODO[WB]: get rid of this, the check has worked a long time, so it's time to get rid of the file in which constraint_matrices are actually declared.  
-  if (dim == 3)
-    if (this->degree < Matrices::n_constraint_matrices+1)
-      {
-       FullMatrix<double> x;
-       x.TableBase<2,double>::reinit (this->interface_constraints_size());
-       x.fill (Matrices::constraint_matrices[this->degree-1]);
-
-       for (unsigned int i=0; i<x.m(); ++i)
-         for (unsigned int j=0; j<x.n(); ++j)
-           Assert (std::fabs (x(i,j) - this->interface_constraints(i,j))
-                    <
-                    1e-14,
-                   ExcInternalError());
-      }    
 }
 
 #endif