*/
virtual void
fill_fe_values (const Triangulation<2,2>::cell_iterator &cell,
- const CellSimilarity::Similarity ,
+ const CellSimilarity::Similarity,
const Quadrature<2> &quadrature,
const Mapping<2,2> &mapping,
const Mapping<2,2>::InternalDataBase &,
* Compute the data on the face of the current cell.
*/
virtual void
- fill_fe_face_values (const Triangulation<2,2>::cell_iterator &cell,
- const unsigned int face_no,
- const Quadrature<1> &quadrature,
- const Mapping<2,2> &mapping,
- const Mapping<2,2>::InternalDataBase &mapping_internal,
+ fill_fe_face_values (const Triangulation<2,2>::cell_iterator &cell,
+ const unsigned int face_no,
+ const Quadrature<1> &quadrature,
+ const Mapping<2,2> &mapping,
+ const Mapping<2,2>::InternalDataBase &mapping_internal,
const dealii::internal::FEValues::MappingRelatedData<2,2> &mapping_data,
- const InternalDataBase &fe_internal,
+ const InternalDataBase &fe_internal,
dealii::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const;
/**
* Compute the data on the subface of the current cell.
*/
virtual void
- fill_fe_subface_values (const Triangulation<2,2>::cell_iterator &cell,
- const unsigned int face_no,
- const unsigned int sub_no,
- const Quadrature<1> &quadrature,
- const Mapping<2,2> &mapping,
- const Mapping<2,2>::InternalDataBase &mapping_internal,
+ fill_fe_subface_values (const Triangulation<2,2>::cell_iterator &cell,
+ const unsigned int face_no,
+ const unsigned int sub_no,
+ const Quadrature<1> &quadrature,
+ const Mapping<2,2> &mapping,
+ const Mapping<2,2>::InternalDataBase &mapping_internal,
const dealii::internal::FEValues::MappingRelatedData<2,2> &mapping_data,
- const InternalDataBase &fe_internal,
+ const InternalDataBase &fe_internal,
dealii::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const;
/**
std::vector<bool> (1, false),
std::vector<ComponentMask>(4, ComponentMask(1,true)))
{
-
// face support points: 2 end vertices
unit_face_support_points.resize(2);
unit_face_support_points[0][0] = 0.0;
void
FE_P1NC::fill_fe_values (const Triangulation<2,2>::cell_iterator &cell,
- const CellSimilarity::Similarity ,
+ const CellSimilarity::Similarity,
const Quadrature<2> &quadrature,
const Mapping<2,2> &mapping,
const Mapping<2,2>::InternalDataBase &,
// hessian
if (flags & update_hessians)
output_data.shape_hessians = fe_data.shape_hessians;
-
}
void
-FE_P1NC::fill_fe_face_values (const Triangulation<2,2>::cell_iterator &cell,
- const unsigned int face_no,
- const Quadrature<1> &quadrature,
- const Mapping<2,2> &mapping,
- const Mapping<2,2>::InternalDataBase &mapping_internal,
+FE_P1NC::fill_fe_face_values (const Triangulation<2,2>::cell_iterator &cell,
+ const unsigned int face_no,
+ const Quadrature<1> &quadrature,
+ const Mapping<2,2> &mapping,
+ const Mapping<2,2>::InternalDataBase &mapping_internal,
const dealii::internal::FEValues::MappingRelatedData<2,2> &mapping_data,
- const InternalDataBase &fe_internal,
+ const InternalDataBase &fe_internal,
dealii::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const
-
{
const UpdateFlags flags(fe_internal.update_each);
output_data.shape_gradients[k][i] = grads[k];
}
}
-
}
}
void
-FE_P1NC::fill_fe_subface_values (const Triangulation<2,2>::cell_iterator &cell,
- const unsigned int face_no,
- const unsigned int sub_no,
- const Quadrature<1> &quadrature,
- const Mapping<2,2> &mapping,
- const Mapping<2,2>::InternalDataBase &mapping_internal,
+FE_P1NC::fill_fe_subface_values (const Triangulation<2,2>::cell_iterator &cell,
+ const unsigned int face_no,
+ const unsigned int sub_no,
+ const Quadrature<1> &quadrature,
+ const Mapping<2,2> &mapping,
+ const Mapping<2,2>::InternalDataBase &mapping_internal,
const dealii::internal::FEValues::MappingRelatedData<2,2> &mapping_data,
- const InternalDataBase &fe_internal,
+ const InternalDataBase &fe_internal,
dealii::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const
-
{
const UpdateFlags flags(fe_internal.update_each);
void FE_P1NC::initialize_constraints ()
{
-
std::vector<Point<1> > constraint_points;
// Add midpoint
interface_constraints(0,0) = 0.5;
interface_constraints(0,1) = 0.5;
-
}