From 9fbf5c1009839a33c1c8853448d722b62ae77fa7 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Thu, 16 Sep 2010 16:02:56 +0000 Subject: [PATCH] Add a TODO. git-svn-id: https://svn.dealii.org/trunk@21995 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/lac/include/lac/constraint_matrix.h | 11 ++++++++++- 1 file changed, 10 insertions(+), 1 deletion(-) diff --git a/deal.II/lac/include/lac/constraint_matrix.h b/deal.II/lac/include/lac/constraint_matrix.h index 49e14719ef..a30e41d48a 100644 --- a/deal.II/lac/include/lac/constraint_matrix.h +++ b/deal.II/lac/include/lac/constraint_matrix.h @@ -54,6 +54,15 @@ namespace internals } +//TODO[WB]: We should have a function of the kind +// ConstraintMatrix::add_constraint (const unsigned int constrained_dof, +// const std::vector > &entries, +// const double inhomogeneity = 0); +// rather than building up constraints piecemeal through add_line/add_entry +// etc. This would also eliminate the possibility of accidentally changing +// existing constraints into something pointless, see the discussion on the +// mailing list on "Tiny bug in interpolate_boundary_values" in Sept. 2010. + /** * This class implements dealing with linear (possibly inhomogeneous) * constraints on degrees of freedom. In particular, it handles constraints of @@ -100,7 +109,7 @@ namespace internals * boundary conditions (the standard way that is extensively discussed in the * tutorial programs is to use MatrixTools::apply_boundary_values). The * general principle of Dirichlet conditions are algebraic constraints of the - * form $x_{i} = b_i$, which fits into the form described above. + * form $x_{i} = b_i$, which fits into the form described above. * * Please note that when combining adaptively refined meshes with hanging node * constraints and inhomogeneous Dirichlet boundary conditions within one -- 2.39.5