From 8188be179ec4da1e7a8c4c80431250cdddd07693 Mon Sep 17 00:00:00 2001 From: Ralf Hartmann Date: Tue, 10 Apr 2001 12:54:12 +0000 Subject: [PATCH] Add mapping combatibility mode. git-svn-id: https://svn.dealii.org/trunk@4421 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/numerics/matrices.cc | 5 +++++ deal.II/deal.II/source/numerics/vectors.cc | 7 +++++++ 2 files changed, 12 insertions(+) diff --git a/deal.II/deal.II/source/numerics/matrices.cc b/deal.II/deal.II/source/numerics/matrices.cc index a6b2faf45c..2c1f2ca4e3 100644 --- a/deal.II/deal.II/source/numerics/matrices.cc +++ b/deal.II/deal.II/source/numerics/matrices.cc @@ -156,6 +156,7 @@ void MatrixCreator::create_mass_matrix (const DoFHandler &dof, SparseMatrix &matrix, const Function * const coefficient) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; create_mass_matrix(mapping, dof, q, matrix, coefficient); } @@ -303,6 +304,7 @@ void MatrixCreator::create_mass_matrix (const DoFHandler &dof, Vector &rhs_vector, const Function * const coefficient) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; create_mass_matrix(mapping, dof, q, matrix, rhs, rhs_vector, coefficient); } @@ -699,6 +701,7 @@ void MatrixCreator::create_boundary_mass_matrix (const DoFHandler std::vector &dof_to_boundary_mapping, const Function *a) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; create_boundary_mass_matrix(mapping, dof, q, matrix, rhs, rhs_vector, dof_to_boundary_mapping, a); @@ -824,6 +827,7 @@ void MatrixCreator::create_laplace_matrix (const DoFHandler &dof, SparseMatrix &matrix, const Function * const coefficient) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; create_laplace_matrix(mapping, dof, q, matrix, coefficient); } @@ -1012,6 +1016,7 @@ void MatrixCreator::create_laplace_matrix (const DoFHandler &dof, Vector &rhs_vector, const Function * const coefficient) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; create_laplace_matrix(mapping, dof, q, matrix, rhs, rhs_vector, coefficient); } diff --git a/deal.II/deal.II/source/numerics/vectors.cc b/deal.II/deal.II/source/numerics/vectors.cc index 6a2efed949..27aa6fa8e8 100644 --- a/deal.II/deal.II/source/numerics/vectors.cc +++ b/deal.II/deal.II/source/numerics/vectors.cc @@ -228,6 +228,7 @@ void VectorTools::interpolate (const DoFHandler &dof, const Function &function, Vector &vec) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; interpolate(mapping, dof, function, vec); } @@ -437,6 +438,7 @@ void VectorTools::project (const DoFHandler &dof, const Quadrature &q_boundary, const bool project_to_boundary_first) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; project(mapping, dof, constraints, quadrature, function, vec, enforce_zero_boundary, q_boundary, project_to_boundary_first); @@ -539,6 +541,7 @@ void VectorTools::create_right_hand_side (const DoFHandler &dof_handler, const Function &rhs_function, Vector &rhs_vector) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; create_right_hand_side(mapping, dof_handler, quadrature, rhs_function, rhs_vector); @@ -739,6 +742,7 @@ VectorTools::interpolate_boundary_values (const DoFHandler &dof, std::map &boundary_values, const std::vector &component_mask) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; interpolate_boundary_values(mapping, dof, boundary_component, boundary_function, boundary_values, component_mask); @@ -846,6 +850,7 @@ VectorTools::project_boundary_values (const DoFHandler &dof, const Quadrature &q, std::map &boundary_values) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; project_boundary_values(mapping, dof, boundary_functions, q, boundary_values); } @@ -1193,6 +1198,7 @@ VectorTools::integrate_difference (const DoFHandler &dof, const NormType &norm, const Function *weight) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; integrate_difference(mapping, dof, fe_function, exact_solution, difference, q, norm, weight); @@ -1245,6 +1251,7 @@ VectorTools::compute_mean_value (const DoFHandler &dof, Vector &v, const unsigned int component) { + Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping")); static const MappingQ1 mapping; return compute_mean_value(mapping, dof, quadrature, v, component); } -- 2.39.5