From 38bc69399441f15f1d4c1b9dfd4e6f03335322a6 Mon Sep 17 00:00:00 2001 From: Matthias Maier Date: Sun, 25 Oct 2015 21:40:08 -0500 Subject: [PATCH] ... --- include/deal.II/lac/constrained_linear_operator.h | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/include/deal.II/lac/constrained_linear_operator.h b/include/deal.II/lac/constrained_linear_operator.h index 0597ae2baa..be63b086bc 100644 --- a/include/deal.II/lac/constrained_linear_operator.h +++ b/include/deal.II/lac/constrained_linear_operator.h @@ -41,7 +41,7 @@ DEAL_II_NAMESPACE_OPEN * Applying the LinearOperator object on a vector u results in * a vector v that stores the result of calling * ConstraintMatrix::distribute() on u - with one important - * difference: inhomogeneities are note applied, but always treated as 0 + * difference: inhomogeneities are not applied, but always treated as 0 * instead. * * The LinearOperator object created by this function is primarily used @@ -154,7 +154,7 @@ LinearOperator project_to_constrained_linear_operator( v = 0.; for (auto i : v.locally_owned_elements()) if (constraint_matrix.is_constrained(i)) - v(i) += u(i); + v(i) = u(i); }; return_op.Tvmult = [&constraint_matrix](Domain &v, const Range &u) @@ -162,7 +162,7 @@ LinearOperator project_to_constrained_linear_operator( v = 0.; for (auto i : v.locally_owned_elements()) if (constraint_matrix.is_constrained(i)) - v(i) += u(i); + v(i) = u(i); }; return return_op; -- 2.39.5