From 3c107f26371d4af803610f865b5e4201d0611799 Mon Sep 17 00:00:00 2001 From: allmaras Date: Thu, 8 Feb 2007 23:50:27 +0000 Subject: [PATCH] Fixed MatrixCreator::create_mass_matrix and create_laplace_matrix git-svn-id: https://svn.dealii.org/trunk@14428 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/numerics/matrices.cc | 24 ++++++++++----------- deal.II/doc/news/changes.html | 8 +++++++ 2 files changed, 20 insertions(+), 12 deletions(-) diff --git a/deal.II/deal.II/source/numerics/matrices.cc b/deal.II/deal.II/source/numerics/matrices.cc index 6959df271a..0cfa1ccbd4 100644 --- a/deal.II/deal.II/source/numerics/matrices.cc +++ b/deal.II/deal.II/source/numerics/matrices.cc @@ -821,9 +821,9 @@ MatrixCreator::create_mass_matrix_2 (const hp::MappingCollection &map fe_values.shape_value(j,point) * weight * coefficient_values[point]); - local_rhs(i) += fe_values.shape_value(i,point) * - rhs_values[point] * weight; } + local_rhs(i) += fe_values.shape_value(i,point) * + rhs_values[point] * weight; } } } @@ -846,9 +846,9 @@ MatrixCreator::create_mass_matrix_2 (const hp::MappingCollection &map fe_values.shape_value(j,point) * weight * coefficient_vector_values[point](component_i)); - local_rhs(i) += fe_values.shape_value(i,point) * - rhs_values[point] * weight; - } + } + local_rhs(i) += fe_values.shape_value(i,point) * + rhs_values[point] * weight; } } } @@ -868,8 +868,8 @@ MatrixCreator::create_mass_matrix_2 (const hp::MappingCollection &map cell_matrix(i,j) += (fe_values.shape_value(i,point) * fe_values.shape_value(j,point) * weight); - local_rhs(i) += v * rhs_values[point] * weight; } + local_rhs(i) += v * rhs_values[point] * weight; } } @@ -2115,8 +2115,8 @@ MatrixCreator::create_laplace_matrix_2 (const Mapping &mapping, fe.system_to_component_index(j).first)) cell_matrix(i,j) += (Du * Dv * weight * coefficient_values[point]); - local_rhs(i) += v * rhs_values[point] * weight; } + local_rhs(i) += v * rhs_values[point] * weight; } } } @@ -2140,8 +2140,8 @@ MatrixCreator::create_laplace_matrix_2 (const Mapping &mapping, (fe.system_to_component_index(j).first == component_i)) cell_matrix(i,j) += (Du * Dv * weight * coefficient_vector_values[point](component_i)); - local_rhs(i) += v * rhs_values[point] * weight; } + local_rhs(i) += v * rhs_values[point] * weight; } } } @@ -2161,8 +2161,8 @@ MatrixCreator::create_laplace_matrix_2 (const Mapping &mapping, (fe.system_to_component_index(i).first == fe.system_to_component_index(j).first)) cell_matrix(i,j) += (Du * Dv * weight); - local_rhs(i) += v * rhs_values[point] * weight; } + local_rhs(i) += v * rhs_values[point] * weight; } } @@ -2514,8 +2514,8 @@ MatrixCreator::create_laplace_matrix_2 (const hp::MappingCollection & fe.system_to_component_index(j).first)) cell_matrix(i,j) += (Du * Dv * weight * coefficient_values[point]); - local_rhs(i) += v * rhs_values[point] * weight; } + local_rhs(i) += v * rhs_values[point] * weight; } } } @@ -2539,8 +2539,8 @@ MatrixCreator::create_laplace_matrix_2 (const hp::MappingCollection & (fe.system_to_component_index(j).first == component_i)) cell_matrix(i,j) += (Du * Dv * weight * coefficient_vector_values[point](component_i)); - local_rhs(i) += v * rhs_values[point] * weight; } + local_rhs(i) += v * rhs_values[point] * weight; } } } @@ -2560,8 +2560,8 @@ MatrixCreator::create_laplace_matrix_2 (const hp::MappingCollection & (fe.system_to_component_index(i).first == fe.system_to_component_index(j).first)) cell_matrix(i,j) += (Du * Dv * weight); - local_rhs(i) += v * rhs_values[point] * weight; } + local_rhs(i) += v * rhs_values[point] * weight; } } diff --git a/deal.II/doc/news/changes.html b/deal.II/doc/news/changes.html index b9434e91c6..4dd40005e6 100644 --- a/deal.II/doc/news/changes.html +++ b/deal.II/doc/news/changes.html @@ -975,6 +975,14 @@ inconvenience this causes.

deal.II

    +
  1. Fixed: MatrixCreator::create_mass_matrix and + create_laplace_matrix computed wrong values for the right + hand sides. This has been fixed. +
    + (Moritz Allmaras 2007/02/08) +

    +
  2. Extended: DataOutBase::Patch has been extended by a new boolean flag points_are_available, which defaults to false. It is set to true if the -- 2.39.5