]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Throw an exception upon inversion of too large matrices in all cases, not only in...
authorhartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 13 Apr 1999 17:24:59 +0000 (17:24 +0000)
committerhartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 13 Apr 1999 17:24:59 +0000 (17:24 +0000)
git-svn-id: https://svn.dealii.org/trunk@1134 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/lac/include/lac/blocksparsematrix.templates.h

index fb737d4fb4511500fbc09b9953551b868ee9267e..7967fb9fd8ed821c489980ded006fdcf6f224cda 100644 (file)
@@ -169,9 +169,11 @@ void BlockSparseMatrix<number, inverse_type>::invert_diagblocks()
        for (unsigned int column_cell=0, column=cell*blocksize;
             column_cell<blocksize; ++column_cell, ++column)
          M_cell(row_cell,column_cell)=(*this)(row,column);
-      if (blocksize<=4)
-       inverse[cell].invert(M_cell);
-      else
+      try
+       {
+         inverse[cell].invert(M_cell);
+       }
+      catch (ExcNotImplemented &)
        {
          M_cell.gauss_jordan();
          inverse[cell]=M_cell;

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.