From: Matthias Maier Date: Mon, 26 Oct 2015 02:40:08 +0000 (-0500) Subject: ... X-Git-Tag: v8.4.0-rc2~126^2~2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=38bc69399441f15f1d4c1b9dfd4e6f03335322a6;p=dealii.git ... --- 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;