]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Undo previous accidental commit.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 25 Oct 2006 01:59:15 +0000 (01:59 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 25 Oct 2006 01:59:15 +0000 (01:59 +0000)
git-svn-id: https://svn.dealii.org/trunk@14083 0785d39b-7218-0410-832d-ea1e28bc413d

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

index fc284daadeccc05ffaa1ae11b50f5f03f2120b44..ceb3474d86e6c6029947ec35f0f742c933d82ca8 100644 (file)
@@ -1243,7 +1243,7 @@ FE_Q<dim>::initialize_embedding ()
                                                // value is somewhat
                                                // fuzzy, but it works
                                                // for
-                                               // 1e-14*(degree^dim)*dim,
+                                               // 1e-14*degree*dim,
                                                // which is kind of
                                                // reasonable given
                                                // that we compute the
@@ -1257,15 +1257,15 @@ FE_Q<dim>::initialize_embedding ()
                                                // growth in
                                                // degree*dim, times a
                                                // small constant.
-             if (std::fabs(cell_value) < 2e-14*std::pow(1.*this->degree,dim)*dim)
+             if (std::fabs(cell_value) < 2e-14*this->degree*dim)
                cell_interpolation(j, i) = 0.;
              else
                cell_interpolation(j, i) = cell_value;
 
-             if (std::fabs(subcell_value) < 2e-14*std::pow(1.*this->degree,dim)*dim)
+             if (std::fabs(subcell_value) < 2e-14*this->degree*dim)
                subcell_interpolation(j, i) = 0.;
              else
-               if (std::fabs(subcell_value-1) < 2e-14*std::pow(1.*this->degree,dim)*dim)
+               if (std::fabs(subcell_value-1) < 2e-14*this->degree*dim)
                  subcell_interpolation(j, i) = 1.;
                else                    
                                                   // we have put our
@@ -1317,7 +1317,7 @@ FE_Q<dim>::initialize_embedding ()
          double sum = 0;
          for (unsigned int col=0; col<this->dofs_per_cell; ++col)
            sum += this->prolongation[child](row,col);
-         Assert (std::fabs(sum-1.) < 2e-14*std::pow(1.*this->degree,dim)*dim,
+         Assert (std::fabs(sum-1.) < 2e-14*this->degree*dim,
                  ExcInternalError());
        }
     }

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.