From 43f9c3a88af29b7d4d1c768e280524b4d476d75d Mon Sep 17 00:00:00 2001
From: bangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Tue, 14 Apr 2009 21:48:20 +0000
Subject: [PATCH] Say something on the singularity of the matrix.

git-svn-id: https://svn.dealii.org/trunk@18611 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/examples/step-34/step-34.cc | 14 ++++++++------
 1 file changed, 8 insertions(+), 6 deletions(-)

diff --git a/deal.II/examples/step-34/step-34.cc b/deal.II/examples/step-34/step-34.cc
index 7b44991a3a..07509c97c9 100644
--- a/deal.II/examples/step-34/step-34.cc
+++ b/deal.II/examples/step-34/step-34.cc
@@ -880,18 +880,20 @@ void BEMProblem<dim>::assemble_system() {
 				 // reality it only produces an LU
 				 // decomposition) and then apply this inverse
 				 // to the right hand side to yield the
-				 // solution:
+				 // solution.
+				 //
+				 // As mentioned in the introduction,
+				 // the solution is only known up to a
+				 // constant potential. We solve this
+				 // issue by subtracting the mean
+				 // value of the vector from each
+				 // vector entry to normalize it.
 template <int dim>
 void BEMProblem<dim>::solve_system() {
     SparseDirectUMFPACK inverse_matrix;
     inverse_matrix.initialize (system_matrix);
     inverse_matrix.vmult (phi, system_rhs);
 
-//TODO: is this true? it seems to me that the BIE is definite...    
-    // Since we are solving a purely Neumann problem, the solution is
-    // only known up to a constant potential. We solve this issue by
-    // subtracting the mean value of the vector from each vector
-    // entry.
     phi.add(-phi.mean_value());
 }
 
-- 
2.39.5