void set_refinement_cycle (const unsigned int refinement_cycle);
virtual void operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const = 0;
+ const Vector<double> &solution) const = 0;
protected:
unsigned int refinement_cycle;
};
TableHandler &results_table);
virtual void operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const;
+ const Vector<double> &solution) const;
DeclException1 (ExcEvaluationPointNotFound,
Point<dim>,
void
PointValueEvaluation<dim>::
operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const
+ const Vector<double> &solution) const
{
// First allocate a variable that will hold the point value. Initialize
// it with a value that is clearly bogus, so that if we fail to set it
const typename DataOut<dim>::OutputFormat output_format);
virtual void operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const;
+ const Vector<double> &solution) const;
private:
const std::string output_name_base;
const typename DataOut<dim>::OutputFormat output_format;
template <int dim>
void
SolutionOutput<dim>::operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const
+ const Vector<double> &solution) const
{
DataOut<dim> data_out;
data_out.attach_dof_handler (dof_handler);
void set_refinement_cycle (const unsigned int refinement_cycle);
virtual void operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const = 0;
+ const Vector<double> &solution) const = 0;
protected:
unsigned int refinement_cycle;
};
PointValueEvaluation (const Point<dim> &evaluation_point);
virtual void operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const;
+ const Vector<double> &solution) const;
DeclException1 (ExcEvaluationPointNotFound,
Point<dim>,
void
PointValueEvaluation<dim>::
operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const
+ const Vector<double> &solution) const
{
double point_value = 1e20;
PointXDerivativeEvaluation (const Point<dim> &evaluation_point);
virtual void operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const;
+ const Vector<double> &solution) const;
DeclException1 (ExcEvaluationPointNotFound,
Point<dim>,
void
PointXDerivativeEvaluation<dim>::
operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const
+ const Vector<double> &solution) const
{
// This time initialize the return value with something useful, since we
// will have to add up a number of contributions and take the mean value
GridOutput (const std::string &output_name_base);
virtual void operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution) const;
+ const Vector<double> &solution) const;
private:
const std::string output_name_base;
};
template <int dim>
void
GridOutput<dim>::operator () (const DoFHandler<dim> &dof_handler,
- const Vector<double> &/*solution*/) const
+ const Vector<double> & /*solution*/) const
{
std::ostringstream filename;
filename << output_name_base << "-"
template <int dim>
inline
void
- BodyForce<dim>::vector_value (const Point<dim> &/*p*/,
+ BodyForce<dim>::vector_value (const Point<dim> & /*p*/,
Vector<double> &values) const
{
Assert (values.size() == dim,
template <int dim>
void
IncrementalBoundaryValues<dim>::
- vector_value (const Point<dim> &/*p*/,
+ vector_value (const Point<dim> & /*p*/,
Vector<double> &values) const
{
Assert (values.size() == dim,
// course. Given our discussion in the introduction of how the solution
// should look like, the following computations should be straightforward:
template <int dim>
- double RightHandSide<dim>::value (const Point<dim> &/*p*/,
+ double RightHandSide<dim>::value (const Point<dim> & /*p*/,
const unsigned int /*component*/) const
{
return 0;
template <int dim>
- double PressureBoundaryValues<dim>::value (const Point<dim> &p,
+ double PressureBoundaryValues<dim>::value (const Point<dim> &p,
const unsigned int /*component*/) const
{
const double alpha = 0.3;
template <int dim>
double
- PressureRightHandSide<dim>::value (const Point<dim> &/*p*/,
+ PressureRightHandSide<dim>::value (const Point<dim> & /*p*/,
const unsigned int /*component*/) const
{
return 0;
template <int dim>
double
- PressureBoundaryValues<dim>::value (const Point<dim> &p,
+ PressureBoundaryValues<dim>::value (const Point<dim> &p,
const unsigned int /*component*/) const
{
return 1-p[0];
template <int dim>
double
- InitialValues<dim>::value (const Point<dim> &p,
+ InitialValues<dim>::value (const Point<dim> &p,
const unsigned int component) const
{
return ZeroFunction<dim>(dim+2).value (p, component);
template <int dim>
double
- BoundaryValues<dim>::value (const Point<dim> &p,
+ BoundaryValues<dim>::value (const Point<dim> &p,
const unsigned int component) const
{
Assert (component < this->n_components,
template <int dim>
double
- RightHandSide<dim>::value (const Point<dim> &/*p*/,
+ RightHandSide<dim>::value (const Point<dim> & /*p*/,
const unsigned int /*component*/) const
{
return 0;
template <int dim>
- double InitialValuesU<dim>::value (const Point<dim> &/*p*/,
+ double InitialValuesU<dim>::value (const Point<dim> & /*p*/,
const unsigned int component) const
{
Assert (component == 0, ExcInternalError());
template <int dim>
- double InitialValuesV<dim>::value (const Point<dim> &/*p*/,
+ double InitialValuesV<dim>::value (const Point<dim> & /*p*/,
const unsigned int component) const
{
Assert (component == 0, ExcInternalError());
template <int dim>
- double RightHandSide<dim>::value (const Point<dim> &/*p*/,
+ double RightHandSide<dim>::value (const Point<dim> & /*p*/,
const unsigned int component) const
{
Assert (component == 0, ExcInternalError());
// describing the material parameters for the number of energy groups
// requested in the input file, and (iii) the finite element to be used by
// all energy groups:
- const Parameters ¶meters;
+ const Parameters ¶meters;
const MaterialData material_data;
FE_Q<dim> fe;
template <int dim>
inline
- void DirichletBoundaryValues<dim>::vector_value (const Point<dim> &/*p*/,
+ void DirichletBoundaryValues<dim>::vector_value (const Point<dim> & /*p*/,
Vector<double> &values) const
{
Assert (values.size() == 2, ExcDimensionMismatch (values.size(), 2));
// system, which consists of two copies of the scalar Q1 field, one for $v$
// and one for $w$:
template <int dim>
- UltrasoundProblem<dim>::UltrasoundProblem (ParameterHandler ¶m)
+ UltrasoundProblem<dim>::UltrasoundProblem (ParameterHandler ¶m)
:
prm(param),
dof_handler(triangulation),
template <int dim>
double
- TemperatureRightHandSide<dim>::value (const Point<dim> &p,
+ TemperatureRightHandSide<dim>::value (const Point<dim> &p,
const unsigned int component) const
{
Assert (component == 0,
template <class PreconditionerA, class PreconditionerMp>
BlockSchurPreconditioner<PreconditionerA, PreconditionerMp>::
- BlockSchurPreconditioner(const TrilinosWrappers::BlockSparseMatrix &S,
+ BlockSchurPreconditioner(const TrilinosWrappers::BlockSparseMatrix &S,
const InverseMatrix<TrilinosWrappers::SparseMatrix,
PreconditionerMp> &Mpinv,
const PreconditionerA &Apreconditioner)
double
compute_viscosity(const std::vector<double> &old_temperature,
const std::vector<double> &old_old_temperature,
- const std::vector<Tensor<1,dim> > &old_temperature_grads,
- const std::vector<Tensor<1,dim> > &old_old_temperature_grads,
+ const std::vector<Tensor<1,dim> > &old_temperature_grads,
+ const std::vector<Tensor<1,dim> > &old_old_temperature_grads,
const std::vector<double> &old_temperature_laplacians,
const std::vector<double> &old_old_temperature_laplacians,
- const std::vector<Tensor<1,dim> > &old_velocity_values,
- const std::vector<Tensor<1,dim> > &old_old_velocity_values,
+ const std::vector<Tensor<1,dim> > &old_velocity_values,
+ const std::vector<Tensor<1,dim> > &old_old_velocity_values,
const std::vector<double> &gamma_values,
const double global_u_infty,
const double global_T_variation,
BoussinesqFlowProblem<dim>::
compute_viscosity (const std::vector<double> &old_temperature,
const std::vector<double> &old_old_temperature,
- const std::vector<Tensor<1,dim> > &old_temperature_grads,
- const std::vector<Tensor<1,dim> > &old_old_temperature_grads,
+ const std::vector<Tensor<1,dim> > &old_temperature_grads,
+ const std::vector<Tensor<1,dim> > &old_old_temperature_grads,
const std::vector<double> &old_temperature_laplacians,
const std::vector<double> &old_old_temperature_laplacians,
- const std::vector<Tensor<1,dim> > &old_velocity_values,
- const std::vector<Tensor<1,dim> > &old_old_velocity_values,
+ const std::vector<Tensor<1,dim> > &old_velocity_values,
+ const std::vector<Tensor<1,dim> > &old_old_velocity_values,
const std::vector<double> &gamma_values,
const double global_u_infty,
const double global_T_variation,
template <int dim>
double
- TemperatureInitialValues<dim>::value (const Point<dim> &p,
+ TemperatureInitialValues<dim>::value (const Point<dim> &p,
const unsigned int) const
{
const double r = p.norm();
class BlockSchurPreconditioner : public Subscriptor
{
public:
- BlockSchurPreconditioner (const TrilinosWrappers::BlockSparseMatrix &S,
- const TrilinosWrappers::BlockSparseMatrix &Spre,
+ BlockSchurPreconditioner (const TrilinosWrappers::BlockSparseMatrix &S,
+ const TrilinosWrappers::BlockSparseMatrix &Spre,
const PreconditionerMp &Mppreconditioner,
const PreconditionerA &Apreconditioner,
const bool do_solve_A)
const SmartPointer<const TrilinosWrappers::BlockSparseMatrix> stokes_matrix;
const SmartPointer<const TrilinosWrappers::BlockSparseMatrix> stokes_preconditioner_matrix;
const PreconditionerMp &mp_preconditioner;
- const PreconditionerA &a_preconditioner;
+ const PreconditionerA &a_preconditioner;
const bool do_solve_A;
};
}
double
compute_viscosity(const std::vector<double> &old_temperature,
const std::vector<double> &old_old_temperature,
- const std::vector<Tensor<1,dim> > &old_temperature_grads,
- const std::vector<Tensor<1,dim> > &old_old_temperature_grads,
+ const std::vector<Tensor<1,dim> > &old_temperature_grads,
+ const std::vector<Tensor<1,dim> > &old_old_temperature_grads,
const std::vector<double> &old_temperature_laplacians,
const std::vector<double> &old_old_temperature_laplacians,
- const std::vector<Tensor<1,dim> > &old_velocity_values,
- const std::vector<Tensor<1,dim> > &old_old_velocity_values,
- const std::vector<SymmetricTensor<2,dim> > &old_strain_rates,
- const std::vector<SymmetricTensor<2,dim> > &old_old_strain_rates,
+ const std::vector<Tensor<1,dim> > &old_velocity_values,
+ const std::vector<Tensor<1,dim> > &old_old_velocity_values,
+ const std::vector<SymmetricTensor<2,dim> > &old_strain_rates,
+ const std::vector<SymmetricTensor<2,dim> > &old_old_strain_rates,
const double global_u_infty,
const double global_T_variation,
const double average_temperature,
void
local_assemble_stokes_system (const typename DoFHandler<dim>::active_cell_iterator &cell,
- Assembly::Scratch::StokesSystem<dim> &scratch,
+ Assembly::Scratch::StokesSystem<dim> &scratch,
Assembly::CopyData::StokesSystem<dim> &data);
void
void
local_assemble_temperature_matrix (const typename DoFHandler<dim>::active_cell_iterator &cell,
- Assembly::Scratch::TemperatureMatrix<dim> &scratch,
+ Assembly::Scratch::TemperatureMatrix<dim> &scratch,
Assembly::CopyData::TemperatureMatrix<dim> &data);
void
BoussinesqFlowProblem<dim>::
compute_viscosity (const std::vector<double> &old_temperature,
const std::vector<double> &old_old_temperature,
- const std::vector<Tensor<1,dim> > &old_temperature_grads,
- const std::vector<Tensor<1,dim> > &old_old_temperature_grads,
+ const std::vector<Tensor<1,dim> > &old_temperature_grads,
+ const std::vector<Tensor<1,dim> > &old_old_temperature_grads,
const std::vector<double> &old_temperature_laplacians,
const std::vector<double> &old_old_temperature_laplacians,
- const std::vector<Tensor<1,dim> > &old_velocity_values,
- const std::vector<Tensor<1,dim> > &old_old_velocity_values,
- const std::vector<SymmetricTensor<2,dim> > &old_strain_rates,
- const std::vector<SymmetricTensor<2,dim> > &old_old_strain_rates,
+ const std::vector<Tensor<1,dim> > &old_velocity_values,
+ const std::vector<Tensor<1,dim> > &old_old_velocity_values,
+ const std::vector<SymmetricTensor<2,dim> > &old_strain_rates,
+ const std::vector<SymmetricTensor<2,dim> > &old_old_strain_rates,
const double global_u_infty,
const double global_T_variation,
const double average_temperature,
BoussinesqFlowProblem<dim>::Postprocessor::
compute_derived_quantities_vector (const std::vector<Vector<double> > &uh,
const std::vector<std::vector<Tensor<1,dim> > > &duh,
- const std::vector<std::vector<Tensor<2,dim> > > &/*dduh*/,
- const std::vector<Point<dim> > &/*normals*/,
- const std::vector<Point<dim> > &/*evaluation_points*/,
+ const std::vector<std::vector<Tensor<2,dim> > > & /*dduh*/,
+ const std::vector<Point<dim> > & /*normals*/,
+ const std::vector<Point<dim> > & /*evaluation_points*/,
std::vector<Vector<double> > &computed_quantities) const
{
const unsigned int n_quadrature_points = uh.size();
void
compute_refinement_indicators (const DoFHandler<dim> &dof_handler,
const Mapping<dim> &mapping,
- const Vector<double> &solution,
+ const Vector<double> &solution,
Vector<double> &refinement_indicators)
{
const unsigned int dofs_per_cell = dof_handler.get_fe().dofs_per_cell;
EulerEquations<dim>::Postprocessor::
compute_derived_quantities_vector (const std::vector<Vector<double> > &uh,
const std::vector<std::vector<Tensor<1,dim> > > &duh,
- const std::vector<std::vector<Tensor<2,dim> > > &/*dduh*/,
- const std::vector<Point<dim> > &/*normals*/,
- const std::vector<Point<dim> > &/*evaluation_points*/,
+ const std::vector<std::vector<Tensor<2,dim> > > & /*dduh*/,
+ const std::vector<Point<dim> > & /*normals*/,
+ const std::vector<Point<dim> > & /*evaluation_points*/,
std::vector<Vector<double> > &computed_quantities) const
{
// At the beginning of the function, let us make sure that all variables
{}
};
- void assemble_one_cell_of_gradient (const IteratorPair &SI,
+ void assemble_one_cell_of_gradient (const IteratorPair &SI,
InitGradScratchData &scratch,
InitGradPerTaskData &data);
template <int dim>
void
NavierStokesProjection<dim>::
- assemble_one_cell_of_gradient (const IteratorPair &SI,
+ assemble_one_cell_of_gradient (const IteratorPair &SI,
InitGradScratchData &scratch,
InitGradPerTaskData &data)
{
template <int dim>
- double Coefficient<dim>::value (const Point<dim> &p,
+ double Coefficient<dim>::value (const Point<dim> &p,
const unsigned int component) const
{
return value<double>(p,component);
void clear();
void reinit (const MGDoFHandler<dim> &dof_handler,
- const ConstraintMatrix &constraints,
+ const ConstraintMatrix &constraints,
const unsigned int level = numbers::invalid_unsigned_int);
unsigned int m () const;
template <int dim, int fe_degree, typename number>
void
LaplaceOperator<dim,fe_degree,number>::reinit (const MGDoFHandler<dim> &dof_handler,
- const ConstraintMatrix &constraints,
+ const ConstraintMatrix &constraints,
const unsigned int level)
{
typename MatrixFree<dim,number>::AdditionalData additional_data;
MPI_Comm _mpi_communicator,
ConditionalOStream _pcout);
- void plast_linear_hardening (SymmetricTensor<4,dim> &stress_strain_tensor,
- SymmetricTensor<2,dim> &strain_tensor,
+ void plast_linear_hardening (SymmetricTensor<4,dim> &stress_strain_tensor,
+ SymmetricTensor<2,dim> &strain_tensor,
unsigned int &elast_points,
unsigned int &plast_points,
double &yield);
- void linearized_plast_linear_hardening (SymmetricTensor<4,dim> &stress_strain_tensor_linearized,
- SymmetricTensor<4,dim> &stress_strain_tensor,
- SymmetricTensor<2,dim> &strain_tensor);
+ void linearized_plast_linear_hardening (SymmetricTensor<4,dim> &stress_strain_tensor_linearized,
+ SymmetricTensor<4,dim> &stress_strain_tensor,
+ SymmetricTensor<2,dim> &strain_tensor);
inline SymmetricTensor<2,dim> get_strain (const FEValues<dim> &fe_values,
const unsigned int shape_func,
const unsigned int q_point) const;
}
template <int dim>
- void ConstitutiveLaw<dim>::plast_linear_hardening (SymmetricTensor<4,dim> &stress_strain_tensor,
- SymmetricTensor<2,dim> &strain_tensor,
+ void ConstitutiveLaw<dim>::plast_linear_hardening (SymmetricTensor<4,dim> &stress_strain_tensor,
+ SymmetricTensor<2,dim> &strain_tensor,
unsigned int &elast_points,
unsigned int &plast_points,
double &yield)
}
template <int dim>
- void ConstitutiveLaw<dim>::linearized_plast_linear_hardening (SymmetricTensor<4,dim> &stress_strain_tensor_linearized,
- SymmetricTensor<4,dim> &stress_strain_tensor,
- SymmetricTensor<2,dim> &strain_tensor)
+ void ConstitutiveLaw<dim>::linearized_plast_linear_hardening (SymmetricTensor<4,dim> &stress_strain_tensor_linearized,
+ SymmetricTensor<4,dim> &stress_strain_tensor,
+ SymmetricTensor<2,dim> &strain_tensor)
{
if (dim == 3)
{
template <int dim>
double
- PressureRightHandSide<dim>::value (const Point<dim> &/*p*/,
+ PressureRightHandSide<dim>::value (const Point<dim> & /*p*/,
const unsigned int /*component*/) const
{
return 0;
template <int dim>
double
- PressureBoundaryValues<dim>::value (const Point<dim> &p,
+ PressureBoundaryValues<dim>::value (const Point<dim> &p,
const unsigned int /*component*/) const
{
return 1-p[0];
template <int dim>
double
- SaturationInitialValues<dim>::value (const Point<dim> &/*p*/,
+ SaturationInitialValues<dim>::value (const Point<dim> & /*p*/,
const unsigned int /*component*/) const
{
return 0.2;
template <class PreconditionerA, class PreconditionerMp>
BlockSchurPreconditioner<PreconditionerA, PreconditionerMp>::
- BlockSchurPreconditioner(const TrilinosWrappers::BlockSparseMatrix &S,
+ BlockSchurPreconditioner(const TrilinosWrappers::BlockSparseMatrix &S,
const InverseMatrix<TrilinosWrappers::SparseMatrix,
PreconditionerMp> &Mpinv,
const PreconditionerA &Apreconditioner)
void project_back_saturation ();
double compute_viscosity (const std::vector<double> &old_saturation,
const std::vector<double> &old_old_saturation,
- const std::vector<Tensor<1,dim> > &old_saturation_grads,
- const std::vector<Tensor<1,dim> > &old_old_saturation_grads,
+ const std::vector<Tensor<1,dim> > &old_saturation_grads,
+ const std::vector<Tensor<1,dim> > &old_old_saturation_grads,
const std::vector<Vector<double> > &present_darcy_values,
const double global_max_u_F_prime,
const double global_S_variation,
TwoPhaseFlowProblem<dim>::
compute_viscosity (const std::vector<double> &old_saturation,
const std::vector<double> &old_old_saturation,
- const std::vector<Tensor<1,dim> > &old_saturation_grads,
- const std::vector<Tensor<1,dim> > &old_old_saturation_grads,
+ const std::vector<Tensor<1,dim> > &old_saturation_grads,
+ const std::vector<Tensor<1,dim> > &old_old_saturation_grads,
const std::vector<Vector<double> > &present_darcy_values,
const double global_max_u_F_prime,
const double global_S_variation,
template <int dim>
double
- StokesBoundaryValues<dim>::value (const Point<dim> &p,
+ StokesBoundaryValues<dim>::value (const Point<dim> &p,
const unsigned int component) const
{
Assert (component < this->n_components,
template <int dim>
double
- RightHandSide<dim>::value (const Point<dim> &/*p*/,
+ RightHandSide<dim>::value (const Point<dim> & /*p*/,
const unsigned int /*component*/) const
{
return 0;
void
Postprocessor<dim>::
compute_derived_quantities_vector (const std::vector<Vector<double> > &uh,
- const std::vector<std::vector<Tensor<1,dim> > > &/*duh*/,
- const std::vector<std::vector<Tensor<2,dim> > > &/*dduh*/,
- const std::vector<Point<dim> > &/*normals*/,
+ const std::vector<std::vector<Tensor<1,dim> > > & /*duh*/,
+ const std::vector<std::vector<Tensor<2,dim> > > & /*dduh*/,
+ const std::vector<Point<dim> > & /*normals*/,
const std::vector<Point<dim> > &evaluation_points,
std::vector<Vector<double> > &computed_quantities) const
{
public:
template <int dim>
static void estimate (const DoFHandler<dim> &dof,
- const Vector<double> &solution,
+ const Vector<double> &solution,
Vector<float> &error_per_cell);
DeclException2 (ExcInvalidVectorLength,
template <int dim>
static void estimate_interval (const DoFHandler<dim> &dof,
- const Vector<double> &solution,
+ const Vector<double> &solution,
const IndexInterval &index_interval,
Vector<float> &error_per_cell);
};
template <int dim>
void
GradientEstimation::estimate (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution,
+ const Vector<double> &solution,
Vector<float> &error_per_cell)
{
// Before starting with the work, we check that the vector into which the
template <int dim>
void
GradientEstimation::estimate_interval (const DoFHandler<dim> &dof_handler,
- const Vector<double> &solution,
+ const Vector<double> &solution,
const IndexInterval &index_interval,
Vector<float> &error_per_cell)
{