* Constructor for the P1 nonconforming element.
* It is only for 2D and codimension = 0.
*/
- FE_P1NC() ;
+ FE_P1NC();
- virtual std::string get_name () const ;
+ virtual std::string get_name () const;
- virtual UpdateFlags requires_update_flags (const UpdateFlags flags) const ;
+ virtual UpdateFlags requires_update_flags (const UpdateFlags flags) const;
- virtual FiniteElement<2,2> *clone () const ;
+ virtual FiniteElement<2,2> *clone () const;
/**
* Destructor.
get_data (const UpdateFlags update_flags,
const Mapping<2,2> &,
const Quadrature<2> &,
- dealii::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const ;
+ dealii::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const;
/**
* Compute the data on the current cell.
/**
* Create the constraints matrix for hanging edges.
*/
- void initialize_constraints () ;
+ void initialize_constraints ();
class InternalData : public FiniteElement<2,2>::InternalDataBase
{
// face support points: 2 end vertices
unit_face_support_points.resize(2);
- unit_face_support_points[0][0] = 0.0 ;
- unit_face_support_points[1][0] = 1.0 ;
+ unit_face_support_points[0][0] = 0.0;
+ unit_face_support_points[1][0] = 1.0;
// initialize constraints matrix
- initialize_constraints () ;
+ initialize_constraints ();
}
UpdateFlags out = update_default;
if (flags & update_values)
- out |= update_values | update_quadrature_points ;
+ out |= update_values | update_quadrature_points;
if (flags & update_gradients)
- out |= update_gradients ;
+ out |= update_gradients;
if (flags & update_cell_normal_vectors)
out |= update_cell_normal_vectors | update_JxW_values;
if (flags & update_hessians)
mpt[3](1) = (cell->vertex(2)(1) + cell->vertex(3)(1))*0.5;
// center point
- Point<2> cpt ;
+ Point<2> cpt;
cpt(0) = (mpt[0](0) + mpt[1](0) + mpt[2](0) + mpt[3](0))*0.25;
cpt(1) = (mpt[0](1) + mpt[1](1) + mpt[2](1) + mpt[3](1))*0.25;
Assert (dynamic_cast<const InternalData *> (&fe_internal) != 0, ExcInternalError());
const InternalData &fe_data = static_cast<const InternalData &> (fe_internal);
- const UpdateFlags flags(fe_internal.update_each) ;
+ const UpdateFlags flags(fe_internal.update_each);
const unsigned int n_q_points = mapping_data.quadrature_points.size();
{
if (flags & update_values)
{
- values[k] = coeffs[k][0]*mapping_data.quadrature_points[i](0) + coeffs[k][1]*mapping_data.quadrature_points[i](1) + coeffs[k][2] ;
+ values[k] = coeffs[k][0]*mapping_data.quadrature_points[i](0) + coeffs[k][1]*mapping_data.quadrature_points[i](1) + coeffs[k][2];
output_data.shape_values[k][i] = values[k];
}
if (flags & update_gradients)
{
- grads[k][0] = coeffs[k][0] ;
- grads[k][1] = coeffs[k][1] ;
+ grads[k][0] = coeffs[k][0];
+ grads[k][1] = coeffs[k][1];
output_data.shape_gradients[k][i] = grads[k];
}
}
dealii::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const
{
- const UpdateFlags flags(fe_internal.update_each) ;
+ const UpdateFlags flags(fe_internal.update_each);
const unsigned int n_q_points = mapping_data.quadrature_points.size();
{
if (flags & update_values)
{
- values[k] = coeffs[k][0]*mapping_data.quadrature_points[i](0) + coeffs[k][1]*mapping_data.quadrature_points[i](1) + coeffs[k][2] ;
+ values[k] = coeffs[k][0]*mapping_data.quadrature_points[i](0) + coeffs[k][1]*mapping_data.quadrature_points[i](1) + coeffs[k][2];
output_data.shape_values[k][i] = values[k];
}
if (flags & update_gradients)
{
- grads[k][0] = coeffs[k][0] ;
- grads[k][1] = coeffs[k][1] ;
+ grads[k][0] = coeffs[k][0];
+ grads[k][1] = coeffs[k][1];
output_data.shape_gradients[k][i] = grads[k];
}
}
dealii::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const
{
- const UpdateFlags flags(fe_internal.update_each) ;
+ const UpdateFlags flags(fe_internal.update_each);
const unsigned int n_q_points = mapping_data.quadrature_points.size();
{
if (flags & update_values)
{
- values[k] = coeffs[k][0]*mapping_data.quadrature_points[i](0) + coeffs[k][1]*mapping_data.quadrature_points[i](1) + coeffs[k][2] ;
+ values[k] = coeffs[k][0]*mapping_data.quadrature_points[i](0) + coeffs[k][1]*mapping_data.quadrature_points[i](1) + coeffs[k][2];
output_data.shape_values[k][i] = values[k];
}
if (flags & update_gradients)
{
- grads[k][0] = coeffs[k][0] ;
- grads[k][1] = coeffs[k][1] ;
+ grads[k][0] = coeffs[k][0];
+ grads[k][1] = coeffs[k][1];
output_data.shape_gradients[k][i] = grads[k];
}
}
interface_constraints
.TableBase<2,double>::reinit (interface_constraints_size());
- interface_constraints(0,0) = 0.5 ;
- interface_constraints(0,1) = 0.5 ;
+ interface_constraints(0,0) = 0.5;
+ interface_constraints(0,1) = 0.5;
}