]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Wolfgang's code removed again
authorGuido Kanschat <dr.guido.kanschat@gmail.com>
Mon, 24 Oct 2005 14:51:36 +0000 (14:51 +0000)
committerGuido Kanschat <dr.guido.kanschat@gmail.com>
Mon, 24 Oct 2005 14:51:36 +0000 (14:51 +0000)
git-svn-id: https://svn.dealii.org/trunk@11654 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_raviart_thomas_nodal.cc

index 93628517a552f4393c1fa5c4f0e5da569137f771..3c978fa9d21f48f158e0651dc7cbcffdcee36367 100644 (file)
@@ -377,20 +377,18 @@ FE_RaviartThomasNodal<dim>::initialize_support_points (const unsigned int deg)
 
   for (unsigned int d=0;d<dim;++d)
     {
-                                       // guard the following statement
-      Assert (dim <= 3, ExcNotImplemented());
-      const QAnisotropic<dim>
-        quadrature (dim == 1 ?
-                    QAnisotropic<dim>(high) :
-                    (dim == 2 ?
-                     QAnisotropic<dim>(((d==0) ? low : high),
-                                       ((d==1) ? low : high)) :
-                     QAnisotropic<dim>(((d==0) ? low : high),
-                                       ((d==1) ? low : high),
-                                       ((d==2) ? low : high))));
+      QAnisotropic<dim>* quadrature;
+      if (dim == 1) quadrature = new QAnisotropic<dim>(high);
+      if (dim == 2) quadrature = new QAnisotropic<dim>(((d==0) ? low : high),
+                                                      ((d==1) ? low : high));
+      if (dim == 3) quadrature = new QAnisotropic<dim>(((d==0) ? low : high),
+                                                      ((d==1) ? low : high),
+                                                      ((d==2) ? low : high));
+      Assert(dim<=3, ExcNotImplemented());
       
-      for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
-       this->generalized_support_points[current++] = quadrature.point(k);
+      for (unsigned int k=0;k<quadrature->n_quadrature_points;++k)
+       this->generalized_support_points[current++] = quadrature->point(k);
+      delete quadrature;
     }
   Assert (current == this->dofs_per_cell, 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.