(const Mapping<dim, spacedim> &mapping,
const DoFHandlerType &dof,
const Quadrature<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask = ComponentMask(),
static void estimate
(const DoFHandlerType &dof,
const Quadrature<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask = ComponentMask(),
(const Mapping<dim, spacedim> &mapping,
const DoFHandlerType &dof,
const Quadrature<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask = ComponentMask(),
static void estimate
(const DoFHandlerType &dof,
const Quadrature<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask = ComponentMask(),
(const Mapping<dim, spacedim> &mapping,
const DoFHandlerType &dof,
const hp::QCollection<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask = ComponentMask(),
static void estimate
(const DoFHandlerType &dof,
const hp::QCollection<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask = ComponentMask(),
(const Mapping<dim, spacedim> &mapping,
const DoFHandlerType &dof,
const hp::QCollection<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask = ComponentMask(),
static void estimate
(const DoFHandlerType &dof,
const hp::QCollection<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask = ComponentMask(),
(const Mapping<1,spacedim> &mapping,
const DoFHandlerType &dof,
const Quadrature<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask = ComponentMask(),
static void estimate
(const DoFHandlerType &dof,
const Quadrature<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask = ComponentMask(),
(const Mapping<1,spacedim> &mapping,
const DoFHandlerType &dof,
const Quadrature<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask = ComponentMask(),
static void estimate
(const DoFHandlerType &dof,
const Quadrature<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask = ComponentMask(),
(const Mapping<1,spacedim> &mapping,
const DoFHandlerType &dof,
const hp::QCollection<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask = ComponentMask(),
static void estimate
(const DoFHandlerType &dof,
const hp::QCollection<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask = ComponentMask(),
(const Mapping<1,spacedim> &mapping,
const DoFHandlerType &dof,
const hp::QCollection<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask = ComponentMask(),
static void estimate
(const DoFHandlerType &dof,
const hp::QCollection<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask = ComponentMask(),
* Some more references to input data to the
* KellyErrorEstimator::estimate() function.
*/
- const typename FunctionMap<spacedim>::type *neumann_bc;
- const ComponentMask component_mask;
- const Function<spacedim> *coefficients;
+ const typename FunctionMap<spacedim,number>::type *neumann_bc;
+ const ComponentMask component_mask;
+ const Function<spacedim> *coefficients;
/**
* Constructor.
const unsigned int n_solution_vectors,
const types::subdomain_id subdomain_id,
const types::material_id material_id,
- const typename FunctionMap<spacedim>::type *neumann_bc,
+ const typename FunctionMap<spacedim,number>::type *neumann_bc,
const ComponentMask &component_mask,
const Function<spacedim> *coefficients);
const unsigned int n_solution_vectors,
const types::subdomain_id subdomain_id,
const types::material_id material_id,
- const typename FunctionMap<spacedim>::type *neumann_bc,
+ const typename FunctionMap<spacedim,number>::type *neumann_bc,
const ComponentMask &component_mask,
const Function<spacedim> *coefficients)
:
// get the values of the boundary function at the quadrature points
if (n_components == 1)
{
- std::vector<double> g(n_q_points);
+ std::vector<number> g(n_q_points);
parallel_data.neumann_bc->find(boundary_id)->second
->value_list (fe_face_values_cell.get_present_fe_values()
.get_quadrature_points(), g);
}
else
{
- std::vector<dealii::Vector<double> >
- g(n_q_points, dealii::Vector<double>(n_components));
+ std::vector<dealii::Vector<number> >
+ g(n_q_points, dealii::Vector<number>(n_components));
parallel_data.neumann_bc->find(boundary_id)->second
->vector_value_list (fe_face_values_cell.get_present_fe_values()
.get_quadrature_points(),
estimate (const Mapping<dim, spacedim> &mapping,
const DoFHandlerType &dof_handler,
const Quadrature<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask,
KellyErrorEstimator<dim,spacedim>::
estimate (const DoFHandlerType &dof_handler,
const Quadrature<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask,
estimate (const Mapping<dim, spacedim> &mapping,
const DoFHandlerType &dof_handler,
const hp::QCollection<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask,
KellyErrorEstimator<dim, spacedim>::
estimate (const DoFHandlerType &dof_handler,
const hp::QCollection<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask,
estimate (const Mapping<dim, spacedim> &mapping,
const DoFHandlerType &dof_handler,
const hp::QCollection<dim-1> &face_quadratures,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask,
Assert (solutions.size() == errors.size(),
ExcIncompatibleNumberOfElements(solutions.size(), errors.size()));
- for (typename FunctionMap<spacedim>::type::const_iterator i=neumann_bc.begin();
+ for (typename FunctionMap<spacedim,typename InputVector::value_type>::type::const_iterator i=neumann_bc.begin();
i!=neumann_bc.end(); ++i)
Assert (i->second->n_components == n_components,
ExcInvalidBoundaryFunction(i->first,
estimate (const Mapping<dim, spacedim> &mapping,
const DoFHandlerType &dof_handler,
const Quadrature<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask,
void KellyErrorEstimator<dim, spacedim>::estimate
(const DoFHandlerType &dof_handler,
const Quadrature<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask,
void KellyErrorEstimator<dim, spacedim>::estimate
(const DoFHandlerType &dof_handler,
const hp::QCollection<dim-1> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (const Mapping<deal_II_dimension, deal_II_space_dimension> &,
const DH<deal_II_dimension,deal_II_space_dimension> &,
const Quadrature<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const VEC &,
Vector<float> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (
const DH<deal_II_dimension,deal_II_space_dimension> &,
const Quadrature<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const VEC &,
Vector<float> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (const Mapping<deal_II_dimension, deal_II_space_dimension> &,
const DH<deal_II_dimension,deal_II_space_dimension> &,
const Quadrature<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const std::vector<const VEC *> &,
std::vector<Vector<float> *> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (
const DH<deal_II_dimension,deal_II_space_dimension> &,
const Quadrature<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const std::vector<const VEC *> &,
std::vector<Vector<float> *> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (const Mapping<deal_II_dimension, deal_II_space_dimension> &,
const DH<deal_II_dimension,deal_II_space_dimension> &,
const hp::QCollection<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const VEC &,
Vector<float> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (
const DH<deal_II_dimension,deal_II_space_dimension> &,
const hp::QCollection<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const VEC &,
Vector<float> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (const Mapping<deal_II_dimension, deal_II_space_dimension> &,
const DH<deal_II_dimension,deal_II_space_dimension> &,
const hp::QCollection<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const std::vector<const VEC *> &,
std::vector<Vector<float> *> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (
const DH<deal_II_dimension,deal_II_space_dimension> &,
const hp::QCollection<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const std::vector<const VEC *> &,
std::vector<Vector<float> *> &,
const ComponentMask &,
estimate (const Mapping<1,spacedim> &mapping,
const DoFHandlerType &dof_handler,
const Quadrature<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask,
KellyErrorEstimator<1,spacedim>::
estimate (const DoFHandlerType &dof_handler,
const Quadrature<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask,
KellyErrorEstimator<1,spacedim>::
estimate (const DoFHandlerType &dof_handler,
const Quadrature<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask,
estimate (const Mapping<1,spacedim> &mapping,
const DoFHandlerType &dof_handler,
const hp::QCollection<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask,
KellyErrorEstimator<1,spacedim>::
estimate (const DoFHandlerType &dof_handler,
const hp::QCollection<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const InputVector &solution,
Vector<float> &error,
const ComponentMask &component_mask,
KellyErrorEstimator<1,spacedim>::
estimate (const DoFHandlerType &dof_handler,
const hp::QCollection<0> &quadrature,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask,
estimate (const Mapping<1,spacedim> & /*mapping*/,
const DoFHandlerType & /*dof_handler*/,
const hp::QCollection<0> &,
- const typename FunctionMap<spacedim>::type & /*neumann_bc*/,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type & /*neumann_bc*/,
const std::vector<const InputVector *> & /*solutions*/,
std::vector<Vector<float>*> & /*errors*/,
const ComponentMask & /*component_mask_*/,
estimate (const Mapping<1,spacedim> &mapping,
const DoFHandlerType &dof_handler,
const Quadrature<0> &,
- const typename FunctionMap<spacedim>::type &neumann_bc,
+ const typename FunctionMap<spacedim,typename InputVector::value_type>::type &neumann_bc,
const std::vector<const InputVector *> &solutions,
std::vector<Vector<float>*> &errors,
const ComponentMask &component_mask,
"indicator for internal boundaries in your boundary "
"value map."));
- for (typename FunctionMap<spacedim>::type::const_iterator i=neumann_bc.begin();
+ for (typename FunctionMap<spacedim,typename InputVector::value_type>::type::const_iterator i=neumann_bc.begin();
i!=neumann_bc.end(); ++i)
Assert (i->second->n_components == n_components,
ExcInvalidBoundaryFunction(i->first,
(coefficient->n_components == 1),
ExcInvalidCoefficient());
- for (typename FunctionMap<spacedim>::type::const_iterator i=neumann_bc.begin();
+ for (typename FunctionMap<spacedim,typename InputVector::value_type>::type::const_iterator i=neumann_bc.begin();
i!=neumann_bc.end(); ++i)
Assert (i->second->n_components == n_components,
ExcInvalidBoundaryFunction(i->first,
{
if (n_components==1)
{
- const double
+ const typename InputVector::value_type
v = neumann_bc.find(n)->second->value(cell->vertex(n));
for (unsigned int s=0; s<n_solution_vectors; ++s)
}
else
{
- Vector<double> v(n_components);
+ Vector<typename InputVector::value_type> v(n_components);
neumann_bc.find(n)->second->vector_value(cell->vertex(n), v);
for (unsigned int s=0; s<n_solution_vectors; ++s)
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (const Mapping<deal_II_dimension, deal_II_space_dimension> &,
const DH<deal_II_dimension,deal_II_space_dimension> &,
const Quadrature<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const VEC &,
Vector<float> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (
const DH<deal_II_dimension,deal_II_space_dimension> &,
const Quadrature<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const VEC &,
Vector<float> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (const Mapping<deal_II_dimension, deal_II_space_dimension> &,
const DH<deal_II_dimension,deal_II_space_dimension> &,
const Quadrature<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const std::vector<const VEC *> &,
std::vector<Vector<float> *> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (
const DH<deal_II_dimension,deal_II_space_dimension> &,
const Quadrature<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const std::vector<const VEC *> &,
std::vector<Vector<float> *> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (const Mapping<deal_II_dimension, deal_II_space_dimension> &,
const DH<deal_II_dimension,deal_II_space_dimension> &,
const hp::QCollection<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const VEC &,
Vector<float> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (
const DH<deal_II_dimension,deal_II_space_dimension> &,
const hp::QCollection<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const VEC &,
Vector<float> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (const Mapping<deal_II_dimension, deal_II_space_dimension> &,
const DH<deal_II_dimension,deal_II_space_dimension> &,
const hp::QCollection<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const std::vector<const VEC *> &,
std::vector<Vector<float> *> &,
const ComponentMask &,
estimate<VEC,DH<deal_II_dimension,deal_II_space_dimension> > (
const DH<deal_II_dimension,deal_II_space_dimension> &,
const hp::QCollection<deal_II_dimension-1> &,
- const FunctionMap<deal_II_space_dimension>::type &,
+ const FunctionMap<deal_II_space_dimension,VEC::value_type>::type &,
const std::vector<const VEC *> &,
std::vector<Vector<float> *> &,
const ComponentMask &,