From: kanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Mon, 4 Mar 2013 13:13:24 +0000 (+0000)
Subject: add more tests for NaN
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a4245c8e07e9125bd13fde943f7102c165644c74;p=dealii-svn.git

add more tests for NaN

git-svn-id: https://svn.dealii.org/trunk@28720 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/include/deal.II/lac/full_matrix.templates.h b/deal.II/include/deal.II/lac/full_matrix.templates.h
index d7adfc766f..8a9a61b5b6 100644
--- a/deal.II/include/deal.II/lac/full_matrix.templates.h
+++ b/deal.II/include/deal.II/lac/full_matrix.templates.h
@@ -1,7 +1,7 @@
 //---------------------------------------------------------------------------
 //    $Id$
 //
-//    Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012 by the deal.II authors
+//    Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -292,6 +292,7 @@ void FullMatrix<number>::forward (Vector<number2>       &dst,
       for (j=0; j<i; ++j)
         s -= number(dst(j)) * (*this)(i,j);
       dst(i) = s/(*this)(i,i);
+      Assert(numbers::is_finite(dst(i)), ExcNumberNotFinite());
     }
 }
 
@@ -312,6 +313,7 @@ void FullMatrix<number>::backward (Vector<number2>       &dst,
       for (j=i+1; j<nu; ++j)
         s -= dst(j) * number2((*this)(i,j));
       dst(i) = s/number2((*this)(i,i));
+      Assert(numbers::is_finite(dst(i)), ExcNumberNotFinite());
     }
 }
 
diff --git a/deal.II/include/deal.II/lac/relaxation_block.templates.h b/deal.II/include/deal.II/lac/relaxation_block.templates.h
index 0d4f2f903a..d326a6aac0 100644
--- a/deal.II/include/deal.II/lac/relaxation_block.templates.h
+++ b/deal.II/include/deal.II/lac/relaxation_block.templates.h
@@ -1,7 +1,7 @@
 //---------------------------------------------------------------------------
 //    $Id$
 //
-//    Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2009, 2010, 2011, 2012 by the deal.II authors
+//    Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2009, 2010, 2011, 2012, 2013 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -215,6 +215,12 @@ RelaxationBlock<MATRIX,inverse_type>::do_step (
             }
           // Apply inverse diagonal
           this->inverse_vmult(block, x_cell, b_cell);
+#ifdef DEBUG
+	  for (unsigned int i=0;i<x_cell.size();++i)
+	    {
+	      Assert(numbers::is_finite(x_cell(i)), ExcNumberNotFinite());
+	    }
+#endif
           // Store in result vector
           row=additional_data->block_list.begin(block);
           for (unsigned int row_cell=0; row_cell<bs; ++row_cell, ++row)