68 typedef dealii::Tensor<1, dim, dealii::VectorizedArray<Number>> vector;
70 typedef typename Base::VectorType VectorType;
71 typedef typename Base::IntegratorCell IntegratorCell;
72 typedef typename Base::IntegratorFace IntegratorFace;
78 typedef Number value_type;
85 initialize(dealii::MatrixFree<dim, Number>
const & matrix_free,
86 dealii::AffineConstraints<Number>
const & affine_constraints,
92 initialize(dealii::MatrixFree<dim, Number>
const & matrix_free,
93 dealii::AffineConstraints<Number>
const & affine_constraints,
95 std::shared_ptr<DivKernel> div_penalty_kernel,
96 std::shared_ptr<ContiKernel> conti_penalty_kernel);
102 get_divergence_kernel_data()
const;
105 get_continuity_kernel_data()
const;
108 get_time_step_size()
const;
110 dealii::LinearAlgebra::distributed::Vector<Number>
const &
111 get_velocity()
const;
114 update(VectorType
const & velocity,
double const & dt);
118 reinit_cell_derived(IntegratorCell & integrator,
unsigned int const cell)
const final;
121 reinit_face_derived(IntegratorFace & integrator_m,
122 IntegratorFace & integrator_p,
123 unsigned int const face)
const final;
126 reinit_boundary_face_derived(IntegratorFace & integrator_m,
unsigned int const face)
const final;
129 reinit_face_cell_based_derived(IntegratorFace & integrator_m,
130 IntegratorFace & integrator_p,
131 unsigned int const cell,
132 unsigned int const face,
133 dealii::types::boundary_id
const boundary_id)
const final;
136 do_cell_integral(IntegratorCell & integrator)
const final;
139 do_face_integral(IntegratorFace & integrator_m, IntegratorFace & integrator_p)
const final;
142 do_face_int_integral(IntegratorFace & integrator_m, IntegratorFace & integrator_p)
const final;
145 do_face_ext_integral(IntegratorFace & integrator_m, IntegratorFace & integrator_p)
const final;
148 do_boundary_integral(IntegratorFace & integrator_m,
149 OperatorType
const & operator_type,
150 dealii::types::boundary_id
const & boundary_id)
const final;
154 VectorType
const * velocity;
155 double time_step_size;
157 std::shared_ptr<Operators::DivergencePenaltyKernel<dim, Number>> div_kernel;
158 std::shared_ptr<Operators::ContinuityPenaltyKernel<dim, Number>> conti_kernel;