From: guido Date: Thu, 21 Nov 2002 17:43:33 +0000 (+0000) Subject: point_difference instantiated X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=24e83d04ade4dcdc2efde6229c417ae9d0c3e244;p=dealii-svn.git point_difference instantiated git-svn-id: https://svn.dealii.org/trunk@6779 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/numerics/vectors.cc b/deal.II/deal.II/source/numerics/vectors.cc index 63b0ad0090..d4c4ae8a0d 100644 --- a/deal.II/deal.II/source/numerics/vectors.cc +++ b/deal.II/deal.II/source/numerics/vectors.cc @@ -1373,9 +1373,9 @@ VectorTools::point_difference (const DoFHandler& dof, ExcDimensionMismatch(difference.size(), fe.n_components())); typename DoFHandler::active_cell_iterator - cell = dof_handler->begin_active(); + cell = dof.begin_active(); const typename DoFHandler::active_cell_iterator - endc = dof_handler->end(); + endc = dof.end(); std::vector > u_value(1, Vector (fe.n_components())); @@ -1383,8 +1383,8 @@ VectorTools::point_difference (const DoFHandler& dof, { if (cell->point_inside(point)) { - Point unit_point = mapping->transform_real_to_unit_cell( - cell, evaluation_point); + Point unit_point = mapping.transform_real_to_unit_cell( + cell, point); for (unsigned int i=0; i > quadrature_point(1, unit_point); @@ -1393,7 +1393,7 @@ VectorTools::point_difference (const DoFHandler& dof, FEValues fe_values( mapping, fe, quadrature, UpdateFlags(update_values)); fe_values.reinit(cell); - fe_values.get_function_values(u, u_value); + fe_values.get_function_values(fe_function, u_value); break; } } @@ -1405,7 +1405,9 @@ VectorTools::point_difference (const DoFHandler& dof, exact_function.vector_value(point, difference); } for (unsigned int i=0;i Vector &, const std::set &); template -void VectorTools::interpolate_boundary_values (const DoFHandler &, - const unsigned char, - const Function &, - std::map &, - const std::vector &); +void VectorTools::interpolate_boundary_values ( + const DoFHandler &, + const unsigned char, + const Function &, + std::map &, + const std::vector &); template -void VectorTools::integrate_difference - -(const Mapping &, - const DoFHandler &, - const Vector &, - const Function &, - Vector &, - const Quadrature &, - const NormType &, - const Function *, - const double); +void VectorTools::integrate_difference ( + const Mapping &, + const DoFHandler &, + const Vector &, + const Function &, + Vector &, + const Quadrature &, + const NormType &, + const Function *, + const double); template void VectorTools::integrate_difference (const DoFHandler &, @@ -1740,6 +1742,37 @@ void VectorTools::integrate_difference const NormType &, const Function *, const double); + +template +void VectorTools::point_difference ( + const DoFHandler&, + const Vector&, + const Function&, + Vector&, + const Point&); +template +void VectorTools::point_difference ( + const DoFHandler&, + const Vector&, + const Function&, + Vector&, + const Point&); +template +void VectorTools::point_difference ( + const DoFHandler&, + const BlockVector&, + const Function&, + Vector&, + const Point&); +template +void VectorTools::point_difference ( + const DoFHandler&, + const BlockVector&, + const Function&, + Vector&, + const Point&); + + #if deal_II_dimension != 1 template void VectorTools::project_boundary_values