]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
invert_diagblocks not necessary anymore
authorguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 17 Jan 2003 08:00:30 +0000 (08:00 +0000)
committerguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 17 Jan 2003 08:00:30 +0000 (08:00 +0000)
git-svn-id: https://svn.dealii.org/trunk@6928 0785d39b-7218-0410-832d-ea1e28bc413d

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

index d51c232a6bb6fc4da2d2272bffcf661d74800647..85fd3cc06d7b91ecaba05d048e3d770c02bfa949 100644 (file)
@@ -1449,19 +1449,6 @@ void DGMethod<dim>::solve (Vector<double> &solution)
                                   // set the right block size.
   preconditioner.initialize(system_matrix, fe.dofs_per_cell);
 
-                                  // As the inverses of the diagonal
-                                  // blocks are needed in each
-                                  // preconditioner step, it is wise
-                                  // to invert the diagonal blocks of
-                                  // the matrix before starting the
-                                  // solver. Otherwise, it takes less
-                                  // memory, but the diagonal blocks
-                                  // are inverted in each
-                                  // preconditioner step,
-                                  // significantly slowing down the
-                                  // linear solving process.
-  preconditioner.invert_diagblocks();
-
                                   // After these preparations we are
                                   // ready to start the linear solver.
   solver.solve (system_matrix, solution, right_hand_side,

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.