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,
}
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);
}
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,
}
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);
}
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,
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);
}
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,
}
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);
}