* not allowed. In that case, use the initialization function with
* several DoFHandler arguments.
*/
- template <typename DoFHandlerType, typename QuadratureType>
+ template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
- reinit(const Mapping<dim> & mapping,
- const DoFHandlerType & dof_handler,
- const AffineConstraints<double> &constraint,
- const QuadratureType & quad,
- const AdditionalData additional_data = AdditionalData());
+ reinit(const Mapping<dim> & mapping,
+ const DoFHandlerType & dof_handler,
+ const AffineConstraints<number2> &constraint,
+ const QuadratureType & quad,
+ const AdditionalData additional_data = AdditionalData());
/**
* Initializes the data structures. Same as above, but using a $Q_1$
* mapping.
*/
- template <typename DoFHandlerType, typename QuadratureType>
+ template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
- reinit(const DoFHandlerType & dof_handler,
- const AffineConstraints<double> &constraint,
- const QuadratureType & quad,
- const AdditionalData additional_data = AdditionalData());
+ reinit(const DoFHandlerType & dof_handler,
+ const AffineConstraints<number2> &constraint,
+ const QuadratureType & quad,
+ const AdditionalData additional_data = AdditionalData());
/**
* Same as above.
* more. Use the reinit function without index set argument to choose the
* one provided by DoFHandler::locally_owned_dofs().
*/
- template <typename DoFHandlerType, typename QuadratureType>
+ template <typename DoFHandlerType, typename QuadratureType, typename number2>
DEAL_II_DEPRECATED void
- reinit(const Mapping<dim> & mapping,
- const DoFHandlerType & dof_handler,
- const AffineConstraints<double> &constraint,
- const IndexSet & locally_owned_dofs,
- const QuadratureType & quad,
- const AdditionalData additional_data = AdditionalData());
+ reinit(const Mapping<dim> & mapping,
+ const DoFHandlerType & dof_handler,
+ const AffineConstraints<number2> &constraint,
+ const IndexSet & locally_owned_dofs,
+ const QuadratureType & quad,
+ const AdditionalData additional_data = AdditionalData());
/**
* Extracts the information needed to perform loops over cells. The
* can be sets independently from the number of DoFHandlers, when several
* elements are always integrated with the same quadrature formula.
*/
- template <typename DoFHandlerType, typename QuadratureType>
+ template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
- reinit(const Mapping<dim> & mapping,
- const std::vector<const DoFHandlerType *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> &constraint,
- const std::vector<QuadratureType> & quad,
+ reinit(const Mapping<dim> & mapping,
+ const std::vector<const DoFHandlerType *> & dof_handler,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
+ const std::vector<QuadratureType> & quad,
const AdditionalData additional_data = AdditionalData());
/**
* Initializes the data structures. Same as above, but using a $Q_1$
* mapping.
*/
- template <typename DoFHandlerType, typename QuadratureType>
+ template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
- reinit(const std::vector<const DoFHandlerType *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> &constraint,
- const std::vector<QuadratureType> & quad,
+ reinit(const std::vector<const DoFHandlerType *> & dof_handler,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
+ const std::vector<QuadratureType> & quad,
const AdditionalData additional_data = AdditionalData());
/**
* more. Use the reinit function without index set argument to choose the
* one provided by DoFHandler::locally_owned_dofs().
*/
- template <typename DoFHandlerType, typename QuadratureType>
+ template <typename DoFHandlerType, typename QuadratureType, typename number2>
DEAL_II_DEPRECATED void
- reinit(const Mapping<dim> & mapping,
- const std::vector<const DoFHandlerType *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> &constraint,
+ reinit(const Mapping<dim> & mapping,
+ const std::vector<const DoFHandlerType *> & dof_handler,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
const std::vector<IndexSet> & locally_owned_set,
const std::vector<QuadratureType> &quad,
const AdditionalData additional_data = AdditionalData());
* as might be necessary when several components in a vector-valued problem
* are integrated together based on the same quadrature formula.
*/
- template <typename DoFHandlerType, typename QuadratureType>
+ template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
- reinit(const Mapping<dim> & mapping,
- const std::vector<const DoFHandlerType *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> &constraint,
- const QuadratureType & quad,
+ reinit(const Mapping<dim> & mapping,
+ const std::vector<const DoFHandlerType *> & dof_handler,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
+ const QuadratureType & quad,
const AdditionalData additional_data = AdditionalData());
/**
* Initializes the data structures. Same as above, but using a $Q_1$
* mapping.
*/
- template <typename DoFHandlerType, typename QuadratureType>
+ template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
- reinit(const std::vector<const DoFHandlerType *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> &constraint,
- const QuadratureType & quad,
+ reinit(const std::vector<const DoFHandlerType *> & dof_handler,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
+ const QuadratureType & quad,
const AdditionalData additional_data = AdditionalData());
/**
* This is the actual reinit function that sets up the indices for the
* DoFHandler case.
*/
+ template <typename number2>
void
internal_reinit(
- const Mapping<dim> & mapping,
- const std::vector<const DoFHandler<dim> *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> &constraint,
- const std::vector<IndexSet> & locally_owned_set,
- const std::vector<hp::QCollection<1>> & quad,
- const AdditionalData & additional_data);
+ const Mapping<dim> & mapping,
+ const std::vector<const DoFHandler<dim> *> & dof_handler,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
+ const std::vector<IndexSet> & locally_owned_set,
+ const std::vector<hp::QCollection<1>> & quad,
+ const AdditionalData & additional_data);
/**
* Same as before but for hp::DoFHandler instead of generic DoFHandler type.
*/
+ template <typename number2>
void
internal_reinit(
- const Mapping<dim> & mapping,
- const std::vector<const hp::DoFHandler<dim> *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> &constraint,
- const std::vector<IndexSet> & locally_owned_set,
- const std::vector<hp::QCollection<1>> & quad,
- const AdditionalData & additional_data);
+ const Mapping<dim> & mapping,
+ const std::vector<const hp::DoFHandler<dim> *> & dof_handler,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
+ const std::vector<IndexSet> & locally_owned_set,
+ const std::vector<hp::QCollection<1>> & quad,
+ const AdditionalData & additional_data);
/**
* Initializes the fields in DoFInfo together with the constraint pool that
* because several DoFInfo classes can have the same weights which
* consequently only need to be stored once).
*/
+ template <typename number2>
void
initialize_indices(
- const std::vector<const AffineConstraints<double> *> &constraint,
- const std::vector<IndexSet> & locally_owned_set,
- const AdditionalData & additional_data);
+ const std::vector<const AffineConstraints<number2> *> &constraint,
+ const std::vector<IndexSet> & locally_owned_set,
+ const AdditionalData & additional_data);
/**
* Initializes the DoFHandlers based on a DoFHandler<dim> argument.
template <int dim, typename Number>
-template <typename DoFHandlerType, typename QuadratureType>
+template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
MatrixFree<dim, Number>::reinit(
const DoFHandlerType & dof_handler,
- const AffineConstraints<double> & constraints_in,
+ const AffineConstraints<number2> & constraints_in,
const QuadratureType & quad,
const typename MatrixFree<dim, Number>::AdditionalData additional_data)
{
- std::vector<const DoFHandlerType *> dof_handlers;
- std::vector<const AffineConstraints<double> *> constraints;
- std::vector<QuadratureType> quads;
+ std::vector<const DoFHandlerType *> dof_handlers;
+ std::vector<const AffineConstraints<number2> *> constraints;
+ std::vector<QuadratureType> quads;
dof_handlers.push_back(&dof_handler);
constraints.push_back(&constraints_in);
template <int dim, typename Number>
-template <typename DoFHandlerType, typename QuadratureType>
+template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
MatrixFree<dim, Number>::reinit(
const Mapping<dim> & mapping,
const DoFHandlerType & dof_handler,
- const AffineConstraints<double> & constraints_in,
+ const AffineConstraints<number2> & constraints_in,
const QuadratureType & quad,
const typename MatrixFree<dim, Number>::AdditionalData additional_data)
{
- std::vector<const DoFHandlerType *> dof_handlers;
- std::vector<const AffineConstraints<double> *> constraints;
+ std::vector<const DoFHandlerType *> dof_handlers;
+ std::vector<const AffineConstraints<number2> *> constraints;
dof_handlers.push_back(&dof_handler);
constraints.push_back(&constraints_in);
template <int dim, typename Number>
-template <typename DoFHandlerType, typename QuadratureType>
+template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
MatrixFree<dim, Number>::reinit(
const std::vector<const DoFHandlerType *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> & constraint,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
const std::vector<QuadratureType> & quad,
const typename MatrixFree<dim, Number>::AdditionalData additional_data)
{
template <int dim, typename Number>
-template <typename DoFHandlerType, typename QuadratureType>
+template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
MatrixFree<dim, Number>::reinit(
const std::vector<const DoFHandlerType *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> & constraint,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
const QuadratureType & quad,
const typename MatrixFree<dim, Number>::AdditionalData additional_data)
{
template <int dim, typename Number>
-template <typename DoFHandlerType, typename QuadratureType>
+template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
MatrixFree<dim, Number>::reinit(
const Mapping<dim> & mapping,
const std::vector<const DoFHandlerType *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> & constraint,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
const QuadratureType & quad,
const typename MatrixFree<dim, Number>::AdditionalData additional_data)
{
template <int dim, typename Number>
-template <typename DoFHandlerType, typename QuadratureType>
+template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
MatrixFree<dim, Number>::reinit(
const Mapping<dim> & mapping,
const std::vector<const DoFHandlerType *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> & constraint,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
const std::vector<QuadratureType> & quad,
const typename MatrixFree<dim, Number>::AdditionalData additional_data)
{
template <int dim, typename Number>
-template <typename DoFHandlerType, typename QuadratureType>
+template <typename DoFHandlerType, typename QuadratureType, typename number2>
void
MatrixFree<dim, Number>::reinit(
const Mapping<dim> & mapping,
const std::vector<const DoFHandlerType *> & dof_handler,
- const std::vector<const AffineConstraints<double> *> & constraint,
+ const std::vector<const AffineConstraints<number2> *> &constraint,
const std::vector<IndexSet> & locally_owned_set,
const std::vector<QuadratureType> & quad,
const typename MatrixFree<dim, Number>::AdditionalData additional_data)