call to <code>interpolate_boundary_values</code> for boundary indicator one:
@code
VectorTools::interpolate_boundary_values(dof_handler,
- 1,
+ types::boundary_id(1),
Functions::ConstantFunction<2>(1.),
boundary_values);
@endcode
// boundary by indicators, and tell the interpolate_boundary_values
// function to only compute the boundary values on a certain part of the
// boundary (e.g. the clamped part, or the inflow boundary). By default,
- // all boundaries have a 0 boundary indicator, unless otherwise specified. If
- // sections of the boundary have different boundary conditions, you have to
+ // all boundaries have a 0 boundary indicator, unless otherwise specified.
+ // (For example, many functions in namespace GridGenerator specify otherwise.)
+ // If sections of the boundary have different boundary conditions, you have to
// number those parts with different boundary indicators. The function call
// below will then only determine boundary values for those parts of the
// boundary for which the boundary indicator is in fact the zero specified as
// class.
std::map<types::global_dof_index, double> boundary_values;
VectorTools::interpolate_boundary_values(dof_handler,
- 0,
+ types::boundary_id(0),
Functions::ZeroFunction<2>(),
boundary_values);
// Now that we got the list of boundary DoFs and their respective boundary
// on faces that have been marked with boundary indicator 0 (because that's
// what we say the function should work on with the second argument below).
// If there are faces with boundary id other than 0, then the function
- // interpolate_boundary_values will do nothing on these faces. For
+ // interpolate_boundary_values() will do nothing on these faces. For
// the Laplace equation doing nothing is equivalent to assuming that
// on those parts of the boundary a zero Neumann boundary condition holds.
std::map<types::global_dof_index, double> boundary_values;
VectorTools::interpolate_boundary_values(dof_handler,
- 0,
+ types::boundary_id(0),
BoundaryValues<dim>(),
boundary_values);
MatrixTools::apply_boundary_values(boundary_values,
// With the matrix so built, we use zero boundary values again:
std::map<types::global_dof_index, double> boundary_values;
VectorTools::interpolate_boundary_values(dof_handler,
- 0,
+ types::boundary_id(0),
Functions::ZeroFunction<dim>(),
boundary_values);
MatrixTools::apply_boundary_values(boundary_values,
// order: if two constraints conflict then the constraint matrix either abort
// or throw an exception via the Assert macro.
VectorTools::interpolate_boundary_values(dof_handler,
- 0,
+ types::boundary_id(0),
Functions::ZeroFunction<dim>(),
constraints);
std::map<types::global_dof_index, double> boundary_values;
VectorTools::interpolate_boundary_values(dof_handler,
- 0,
+ types::boundary_id(0),
Solution<dim>(),
boundary_values);
MatrixTools::apply_boundary_values(boundary_values,
constraints.clear();
DoFTools::make_hanging_node_constraints(dof_handler, constraints);
VectorTools::interpolate_boundary_values(dof_handler,
- 0,
+ types::boundary_id(0),
Functions::ZeroFunction<dim>(dim),
constraints);
constraints.close();