From: bangerth Date: Thu, 16 Sep 2010 16:02:56 +0000 (+0000) Subject: Add a TODO. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=abb236293ffbba5631d70b8a69cf942d3179d2cf;p=dealii-svn.git Add a TODO. git-svn-id: https://svn.dealii.org/trunk@21995 0785d39b-7218-0410-832d-ea1e28bc413d --- 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