]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
defer "not implemented" to underlying function
authorkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 7 Dec 2011 16:33:42 +0000 (16:33 +0000)
committerkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 7 Dec 2011 16:33:42 +0000 (16:33 +0000)
git-svn-id: https://svn.dealii.org/trunk@24795 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/include/deal.II/fe/fe_poly.templates.h

index f24c672e14a75d3c3d1a6a83cb56a324272403ac..26288815024e032c6cb757e9bc63af36d5264ce0 100644 (file)
@@ -242,7 +242,7 @@ inline
 void
 FE_Poly<TensorProductPolynomials<1>,1,2>::fill_fe_values
   (const Mapping<1,2>                        &mapping,
-   const Triangulation<1,2>::cell_iterator   &/*cell*/,
+   const Triangulation<1,2>::cell_iterator   &cell,
    const Quadrature<1>                       &quadrature,
    Mapping<1,2>::InternalDataBase            &mapping_data,
    Mapping<1,2>::InternalDataBase            &fedata,
@@ -270,11 +270,8 @@ FE_Poly<TensorProductPolynomials<1>,1,2>::fill_fe_values
     }
 
   if (flags & update_hessians && cell_similarity != CellSimilarity::translation)
-    {
-      AssertThrow(false, ExcNotImplemented());
-/*       this->compute_2nd (mapping, cell, QProjector<dim>::DataSetDescriptor::cell(), */
-/*                      mapping_data, fe_data, data); */
-    }
+    this->compute_2nd (mapping, cell, QProjector<1>::DataSetDescriptor::cell(),
+                      mapping_data, fe_data, data);
 }
 
 
@@ -284,7 +281,7 @@ inline
 void
 FE_Poly<TensorProductPolynomials<2>,2,3>::fill_fe_values
   (const Mapping<2,3>                      &mapping,
-   const Triangulation<2,3>::cell_iterator &/* cell */,
+   const Triangulation<2,3>::cell_iterator &cell,
    const Quadrature<2>                     &quadrature,
    Mapping<2,3>::InternalDataBase          &mapping_data,
    Mapping<2,3>::InternalDataBase          &fedata,
@@ -311,12 +308,8 @@ FE_Poly<TensorProductPolynomials<2>,2,3>::fill_fe_values
     }
 
   if (flags & update_hessians && cell_similarity != CellSimilarity::translation)
-    {
-      AssertThrow(false, ExcNotImplemented());
-/*       this->compute_2nd (mapping, cell, QProjector<dim>::DataSetDescriptor::cell(), */
-/*                      mapping_data, fe_data, data); */
-    }
-
+    this->compute_2nd (mapping, cell, QProjector<2>::DataSetDescriptor::cell(),
+                      mapping_data, fe_data, data);
 }
 
 
@@ -325,7 +318,7 @@ inline
 void
 FE_Poly<PolynomialSpace<1>,1,2>::fill_fe_values (
   const Mapping<1,2>                      &mapping,
-  const Triangulation<1,2>::cell_iterator & /*cell*/,
+  const Triangulation<1,2>::cell_iterator &cell,
   const Quadrature<1>                     &quadrature,
   Mapping<1,2>::InternalDataBase          &mapping_data,
   Mapping<1,2>::InternalDataBase          &fedata,
@@ -347,20 +340,15 @@ FE_Poly<PolynomialSpace<1>,1,2>::fill_fe_values (
       if (flags & update_values)
        for (unsigned int i=0; i<quadrature.size(); ++i)
          data.shape_values(k,i) = fe_data.shape_values[k][i];
-
-// TODO: I would think this should work. Guido
-
+      
       if (flags & update_gradients && cell_similarity != CellSimilarity::translation)
        mapping.transform(fe_data.shape_gradients[k], data.shape_gradients[k],
                          mapping_data, mapping_covariant);
     }
 
   if (flags & update_hessians && cell_similarity != CellSimilarity::translation)
-    {
-      AssertThrow(false, ExcNotImplemented());
-/*       this->compute_2nd (mapping, cell, QProjector<dim>::DataSetDescriptor::cell(), */
-/*                      mapping_data, fe_data, data); */
-    }
+    this->compute_2nd (mapping, cell, QProjector<1>::DataSetDescriptor::cell(),
+                      mapping_data, fe_data, data);
 }
 
 
@@ -369,7 +357,7 @@ inline
 void
 FE_Poly<PolynomialSpace<2>,2,3>::fill_fe_values
   (const Mapping<2,3>                      &mapping,
-   const Triangulation<2,3>::cell_iterator &/* cell */,
+   const Triangulation<2,3>::cell_iterator &cell,
    const Quadrature<2>                     &quadrature,
    Mapping<2,3>::InternalDataBase          &mapping_data,
    Mapping<2,3>::InternalDataBase          &fedata,
@@ -394,12 +382,8 @@ FE_Poly<PolynomialSpace<2>,2,3>::fill_fe_values
     }
 
   if (flags & update_hessians && cell_similarity != CellSimilarity::translation)
-    {
-      AssertThrow(false, ExcNotImplemented());
-/*       this->compute_2nd (mapping, cell, QProjector<dim>::DataSetDescriptor::cell(), */
-/*                      mapping_data, fe_data, data); */
-    }
-
+    this->compute_2nd (mapping, cell, QProjector<2>::DataSetDescriptor::cell(),
+                      mapping_data, fe_data, data);
 }
 
 

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.