]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Remove a check that was not only long satisfied, but also was in a function
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 17 Nov 2004 03:39:54 +0000 (03:39 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 17 Nov 2004 03:39:54 +0000 (03:39 +0000)
(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

deal.II/deal.II/source/fe/fe_q.cc

index f46af800504c11b8164b9d760a705df80c58a084..85a61d26e54021f712236a43ccecd85791c41168 100644 (file)
@@ -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

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.