]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Say something on the singularity of the matrix.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 14 Apr 2009 21:48:20 +0000 (21:48 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 14 Apr 2009 21:48:20 +0000 (21:48 +0000)
git-svn-id: https://svn.dealii.org/trunk@18611 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-34/step-34.cc

index 7b44991a3ae5c4406b35631a32ddae26bc795313..07509c97c9e7c50fe7daa759abc246e823330f3b 100644 (file)
@@ -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());
 }
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.