--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check running over const iterators
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_02/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+
+ for (typename FullMatrix<number>::const_iterator
+ p = m.begin(); p!=m.end(); ++p)
+ deallog << p->row() << ' ' << p->column() << ' '
+ << p->value()
+ << std::endl;
+}
+
--- /dev/null
+
+DEAL::0 0 2.0
+DEAL::0 1 3.0
+DEAL::0 2 4.0
+DEAL::0 3 5.0
+DEAL::0 4 6.0
+DEAL::1 0 4.0
+DEAL::1 1 6.0
+DEAL::1 2 8.0
+DEAL::1 3 10.
+DEAL::1 4 12.
+DEAL::2 0 6.0
+DEAL::2 1 9.0
+DEAL::2 2 12.
+DEAL::2 3 15.
+DEAL::2 4 18.
+DEAL::3 0 8.0
+DEAL::3 1 12.
+DEAL::3 2 16.
+DEAL::3 3 20.
+DEAL::3 4 24.
+DEAL::4 0 10.
+DEAL::4 1 15.
+DEAL::4 2 20.
+DEAL::4 3 25.
+DEAL::4 4 30.
+DEAL::0 0 2.0
+DEAL::0 1 3.0
+DEAL::0 2 4.0
+DEAL::0 3 5.0
+DEAL::0 4 6.0
+DEAL::1 0 4.0
+DEAL::1 1 6.0
+DEAL::1 2 8.0
+DEAL::1 3 10.
+DEAL::1 4 12.
+DEAL::2 0 6.0
+DEAL::2 1 9.0
+DEAL::2 2 12.
+DEAL::2 3 15.
+DEAL::2 4 18.
+DEAL::3 0 8.0
+DEAL::3 1 12.
+DEAL::3 2 16.
+DEAL::3 3 20.
+DEAL::3 4 24.
+DEAL::4 0 10.
+DEAL::4 1 15.
+DEAL::4 2 20.
+DEAL::4 3 25.
+DEAL::4 4 30.
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check running over const iterators starting at the second line
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_03/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+
+
+ for (typename FullMatrix<number>::const_iterator
+ p = m.begin(1); p!=m.end(1); ++p)
+ deallog << p->row() << ' ' << p->column() << ' '
+ << p->value()
+ << std::endl;
+}
+
--- /dev/null
+
+DEAL::1 0 4.0
+DEAL::1 1 6.0
+DEAL::1 2 8.0
+DEAL::1 3 10.
+DEAL::1 4 12.
+DEAL::1 0 4.0
+DEAL::1 1 6.0
+DEAL::1 2 8.0
+DEAL::1 3 10.
+DEAL::1 4 12.
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check creation and output of a matrix using an array
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_04/output";
+
+
+template <typename number>
+void
+check ()
+{
+ const number array[] = { 1, 2, 3, 4, 5,
+ 6, 7, 8, 9, 0,
+ 1, 2, 3, 4, 5,
+ 6, 7, 8, 9, 0,
+ 1, 2, 3, 4, 5 };
+
+ FullMatrix<number> m (5,5,array);
+
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 1.0
+DEAL::0 1 2.0
+DEAL::0 2 3.0
+DEAL::0 3 4.0
+DEAL::0 4 5.0
+DEAL::1 0 6.0
+DEAL::1 1 7.0
+DEAL::1 2 8.0
+DEAL::1 3 9.0
+DEAL::1 4 0
+DEAL::2 0 1.0
+DEAL::2 1 2.0
+DEAL::2 2 3.0
+DEAL::2 3 4.0
+DEAL::2 4 5.0
+DEAL::3 0 6.0
+DEAL::3 1 7.0
+DEAL::3 2 8.0
+DEAL::3 3 9.0
+DEAL::3 4 0
+DEAL::4 0 1.0
+DEAL::4 1 2.0
+DEAL::4 2 3.0
+DEAL::4 3 4.0
+DEAL::4 4 5.0
+DEAL::0 0 1.0
+DEAL::0 1 2.0
+DEAL::0 2 3.0
+DEAL::0 3 4.0
+DEAL::0 4 5.0
+DEAL::1 0 6.0
+DEAL::1 1 7.0
+DEAL::1 2 8.0
+DEAL::1 3 9.0
+DEAL::1 4 0.0
+DEAL::2 0 1.0
+DEAL::2 1 2.0
+DEAL::2 2 3.0
+DEAL::2 3 4.0
+DEAL::2 4 5.0
+DEAL::3 0 6.0
+DEAL::3 1 7.0
+DEAL::3 2 8.0
+DEAL::3 3 9.0
+DEAL::3 4 0.0
+DEAL::4 0 1.0
+DEAL::4 1 2.0
+DEAL::4 2 3.0
+DEAL::4 3 4.0
+DEAL::4 4 5.0
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check creation using IdentityMatrix
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_05/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m (IdentityMatrix(5));
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 1.0
+DEAL::0 1 0
+DEAL::0 2 0
+DEAL::0 3 0
+DEAL::0 4 0
+DEAL::1 0 0
+DEAL::1 1 1.0
+DEAL::1 2 0
+DEAL::1 3 0
+DEAL::1 4 0
+DEAL::2 0 0
+DEAL::2 1 0
+DEAL::2 2 1.0
+DEAL::2 3 0
+DEAL::2 4 0
+DEAL::3 0 0
+DEAL::3 1 0
+DEAL::3 2 0
+DEAL::3 3 1.0
+DEAL::3 4 0
+DEAL::4 0 0
+DEAL::4 1 0
+DEAL::4 2 0
+DEAL::4 3 0
+DEAL::4 4 1.0
+DEAL::0 0 1.0
+DEAL::0 1 0.0
+DEAL::0 2 0.0
+DEAL::0 3 0.0
+DEAL::0 4 0.0
+DEAL::1 0 0.0
+DEAL::1 1 1.0
+DEAL::1 2 0.0
+DEAL::1 3 0.0
+DEAL::1 4 0.0
+DEAL::2 0 0.0
+DEAL::2 1 0.0
+DEAL::2 2 1.0
+DEAL::2 3 0.0
+DEAL::2 4 0.0
+DEAL::3 0 0.0
+DEAL::3 1 0.0
+DEAL::3 2 0.0
+DEAL::3 3 1.0
+DEAL::3 4 0.0
+DEAL::4 0 0.0
+DEAL::4 1 0.0
+DEAL::4 2 0.0
+DEAL::4 3 0.0
+DEAL::4 4 1.0
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check matrix assignment
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_06/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+
+ FullMatrix<number> n;
+ n = m;
+
+ print_matrix (n);
+}
+
--- /dev/null
+
+DEAL::0 0 2.0
+DEAL::0 1 3.0
+DEAL::0 2 4.0
+DEAL::0 3 5.0
+DEAL::0 4 6.0
+DEAL::1 0 4.0
+DEAL::1 1 6.0
+DEAL::1 2 8.0
+DEAL::1 3 10.
+DEAL::1 4 12.
+DEAL::2 0 6.0
+DEAL::2 1 9.0
+DEAL::2 2 12.
+DEAL::2 3 15.
+DEAL::2 4 18.
+DEAL::3 0 8.0
+DEAL::3 1 12.
+DEAL::3 2 16.
+DEAL::3 3 20.
+DEAL::3 4 24.
+DEAL::4 0 10.
+DEAL::4 1 15.
+DEAL::4 2 20.
+DEAL::4 3 25.
+DEAL::4 4 30.
+DEAL::0 0 2.0
+DEAL::0 1 3.0
+DEAL::0 2 4.0
+DEAL::0 3 5.0
+DEAL::0 4 6.0
+DEAL::1 0 4.0
+DEAL::1 1 6.0
+DEAL::1 2 8.0
+DEAL::1 3 10.
+DEAL::1 4 12.
+DEAL::2 0 6.0
+DEAL::2 1 9.0
+DEAL::2 2 12.
+DEAL::2 3 15.
+DEAL::2 4 18.
+DEAL::3 0 8.0
+DEAL::3 1 12.
+DEAL::3 2 16.
+DEAL::3 3 20.
+DEAL::3 4 24.
+DEAL::4 0 10.
+DEAL::4 1 15.
+DEAL::4 2 20.
+DEAL::4 3 25.
+DEAL::4 4 30.
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check matrix=0
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_07/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ m = 0;
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 0
+DEAL::0 1 0
+DEAL::0 2 0
+DEAL::0 3 0
+DEAL::0 4 0
+DEAL::1 0 0
+DEAL::1 1 0
+DEAL::1 2 0
+DEAL::1 3 0
+DEAL::1 4 0
+DEAL::2 0 0
+DEAL::2 1 0
+DEAL::2 2 0
+DEAL::2 3 0
+DEAL::2 4 0
+DEAL::3 0 0
+DEAL::3 1 0
+DEAL::3 2 0
+DEAL::3 3 0
+DEAL::3 4 0
+DEAL::4 0 0
+DEAL::4 1 0
+DEAL::4 2 0
+DEAL::4 3 0
+DEAL::4 4 0
+DEAL::0 0 0.0
+DEAL::0 1 0.0
+DEAL::0 2 0.0
+DEAL::0 3 0.0
+DEAL::0 4 0.0
+DEAL::1 0 0.0
+DEAL::1 1 0.0
+DEAL::1 2 0.0
+DEAL::1 3 0.0
+DEAL::1 4 0.0
+DEAL::2 0 0.0
+DEAL::2 1 0.0
+DEAL::2 2 0.0
+DEAL::2 3 0.0
+DEAL::2 4 0.0
+DEAL::3 0 0.0
+DEAL::3 1 0.0
+DEAL::3 2 0.0
+DEAL::3 3 0.0
+DEAL::3 4 0.0
+DEAL::4 0 0.0
+DEAL::4 1 0.0
+DEAL::4 2 0.0
+DEAL::4 3 0.0
+DEAL::4 4 0.0
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::copy_from
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_08/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+
+ FullMatrix<number> n;
+ n.copy_from (m);
+
+ print_matrix (n);
+}
+
--- /dev/null
+
+DEAL::0 0 2.0
+DEAL::0 1 3.0
+DEAL::0 2 4.0
+DEAL::0 3 5.0
+DEAL::0 4 6.0
+DEAL::1 0 4.0
+DEAL::1 1 6.0
+DEAL::1 2 8.0
+DEAL::1 3 10.
+DEAL::1 4 12.
+DEAL::2 0 6.0
+DEAL::2 1 9.0
+DEAL::2 2 12.
+DEAL::2 3 15.
+DEAL::2 4 18.
+DEAL::3 0 8.0
+DEAL::3 1 12.
+DEAL::3 2 16.
+DEAL::3 3 20.
+DEAL::3 4 24.
+DEAL::4 0 10.
+DEAL::4 1 15.
+DEAL::4 2 20.
+DEAL::4 3 25.
+DEAL::4 4 30.
+DEAL::0 0 2.0
+DEAL::0 1 3.0
+DEAL::0 2 4.0
+DEAL::0 3 5.0
+DEAL::0 4 6.0
+DEAL::1 0 4.0
+DEAL::1 1 6.0
+DEAL::1 2 8.0
+DEAL::1 3 10.
+DEAL::1 4 12.
+DEAL::2 0 6.0
+DEAL::2 1 9.0
+DEAL::2 2 12.
+DEAL::2 3 15.
+DEAL::2 4 18.
+DEAL::3 0 8.0
+DEAL::3 1 12.
+DEAL::3 2 16.
+DEAL::3 3 20.
+DEAL::3 4 24.
+DEAL::4 0 10.
+DEAL::4 1 15.
+DEAL::4 2 20.
+DEAL::4 3 25.
+DEAL::4 4 30.
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check assignment from IdenityMatrix
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_09/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ m = IdentityMatrix(5);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 1.0
+DEAL::0 1 0
+DEAL::0 2 0
+DEAL::0 3 0
+DEAL::0 4 0
+DEAL::1 0 0
+DEAL::1 1 1.0
+DEAL::1 2 0
+DEAL::1 3 0
+DEAL::1 4 0
+DEAL::2 0 0
+DEAL::2 1 0
+DEAL::2 2 1.0
+DEAL::2 3 0
+DEAL::2 4 0
+DEAL::3 0 0
+DEAL::3 1 0
+DEAL::3 2 0
+DEAL::3 3 1.0
+DEAL::3 4 0
+DEAL::4 0 0
+DEAL::4 1 0
+DEAL::4 2 0
+DEAL::4 3 0
+DEAL::4 4 1.0
+DEAL::0 0 1.0
+DEAL::0 1 0.0
+DEAL::0 2 0.0
+DEAL::0 3 0.0
+DEAL::0 4 0.0
+DEAL::1 0 0.0
+DEAL::1 1 1.0
+DEAL::1 2 0.0
+DEAL::1 3 0.0
+DEAL::1 4 0.0
+DEAL::2 0 0.0
+DEAL::2 1 0.0
+DEAL::2 2 1.0
+DEAL::2 3 0.0
+DEAL::2 4 0.0
+DEAL::3 0 0.0
+DEAL::3 1 0.0
+DEAL::3 2 0.0
+DEAL::3 3 1.0
+DEAL::3 4 0.0
+DEAL::4 0 0.0
+DEAL::4 1 0.0
+DEAL::4 2 0.0
+DEAL::4 3 0.0
+DEAL::4 4 1.0
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::fill
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_10/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+
+ FullMatrix<number> n (10, 10);
+ n.fill (m, 3, 4);
+ print_matrix (n);
+}
+
--- /dev/null
+
+DEAL::0 0 0
+DEAL::0 1 0
+DEAL::0 2 0
+DEAL::0 3 0
+DEAL::0 4 0
+DEAL::0 5 0
+DEAL::0 6 0
+DEAL::0 7 0
+DEAL::0 8 0
+DEAL::0 9 0
+DEAL::1 0 0
+DEAL::1 1 0
+DEAL::1 2 0
+DEAL::1 3 0
+DEAL::1 4 0
+DEAL::1 5 0
+DEAL::1 6 0
+DEAL::1 7 0
+DEAL::1 8 0
+DEAL::1 9 0
+DEAL::2 0 0
+DEAL::2 1 0
+DEAL::2 2 0
+DEAL::2 3 0
+DEAL::2 4 0
+DEAL::2 5 0
+DEAL::2 6 0
+DEAL::2 7 0
+DEAL::2 8 0
+DEAL::2 9 0
+DEAL::3 0 0
+DEAL::3 1 0
+DEAL::3 2 0
+DEAL::3 3 0
+DEAL::3 4 2.0
+DEAL::3 5 3.0
+DEAL::3 6 4.0
+DEAL::3 7 5.0
+DEAL::3 8 6.0
+DEAL::3 9 0
+DEAL::4 0 0
+DEAL::4 1 0
+DEAL::4 2 0
+DEAL::4 3 0
+DEAL::4 4 4.0
+DEAL::4 5 6.0
+DEAL::4 6 8.0
+DEAL::4 7 10.
+DEAL::4 8 12.
+DEAL::4 9 0
+DEAL::5 0 0
+DEAL::5 1 0
+DEAL::5 2 0
+DEAL::5 3 0
+DEAL::5 4 6.0
+DEAL::5 5 9.0
+DEAL::5 6 12.
+DEAL::5 7 15.
+DEAL::5 8 18.
+DEAL::5 9 0
+DEAL::6 0 0
+DEAL::6 1 0
+DEAL::6 2 0
+DEAL::6 3 0
+DEAL::6 4 8.0
+DEAL::6 5 12.
+DEAL::6 6 16.
+DEAL::6 7 20.
+DEAL::6 8 24.
+DEAL::6 9 0
+DEAL::7 0 0
+DEAL::7 1 0
+DEAL::7 2 0
+DEAL::7 3 0
+DEAL::7 4 10.
+DEAL::7 5 15.
+DEAL::7 6 20.
+DEAL::7 7 25.
+DEAL::7 8 30.
+DEAL::7 9 0
+DEAL::8 0 0
+DEAL::8 1 0
+DEAL::8 2 0
+DEAL::8 3 0
+DEAL::8 4 0
+DEAL::8 5 0
+DEAL::8 6 0
+DEAL::8 7 0
+DEAL::8 8 0
+DEAL::8 9 0
+DEAL::9 0 0
+DEAL::9 1 0
+DEAL::9 2 0
+DEAL::9 3 0
+DEAL::9 4 0
+DEAL::9 5 0
+DEAL::9 6 0
+DEAL::9 7 0
+DEAL::9 8 0
+DEAL::9 9 0
+DEAL::0 0 0.0
+DEAL::0 1 0.0
+DEAL::0 2 0.0
+DEAL::0 3 0.0
+DEAL::0 4 0.0
+DEAL::0 5 0.0
+DEAL::0 6 0.0
+DEAL::0 7 0.0
+DEAL::0 8 0.0
+DEAL::0 9 0.0
+DEAL::1 0 0.0
+DEAL::1 1 0.0
+DEAL::1 2 0.0
+DEAL::1 3 0.0
+DEAL::1 4 0.0
+DEAL::1 5 0.0
+DEAL::1 6 0.0
+DEAL::1 7 0.0
+DEAL::1 8 0.0
+DEAL::1 9 0.0
+DEAL::2 0 0.0
+DEAL::2 1 0.0
+DEAL::2 2 0.0
+DEAL::2 3 0.0
+DEAL::2 4 0.0
+DEAL::2 5 0.0
+DEAL::2 6 0.0
+DEAL::2 7 0.0
+DEAL::2 8 0.0
+DEAL::2 9 0.0
+DEAL::3 0 0.0
+DEAL::3 1 0.0
+DEAL::3 2 0.0
+DEAL::3 3 0.0
+DEAL::3 4 2.0
+DEAL::3 5 3.0
+DEAL::3 6 4.0
+DEAL::3 7 5.0
+DEAL::3 8 6.0
+DEAL::3 9 0.0
+DEAL::4 0 0.0
+DEAL::4 1 0.0
+DEAL::4 2 0.0
+DEAL::4 3 0.0
+DEAL::4 4 4.0
+DEAL::4 5 6.0
+DEAL::4 6 8.0
+DEAL::4 7 10.
+DEAL::4 8 12.
+DEAL::4 9 0.0
+DEAL::5 0 0.0
+DEAL::5 1 0.0
+DEAL::5 2 0.0
+DEAL::5 3 0.0
+DEAL::5 4 6.0
+DEAL::5 5 9.0
+DEAL::5 6 12.
+DEAL::5 7 15.
+DEAL::5 8 18.
+DEAL::5 9 0.0
+DEAL::6 0 0.0
+DEAL::6 1 0.0
+DEAL::6 2 0.0
+DEAL::6 3 0.0
+DEAL::6 4 8.0
+DEAL::6 5 12.
+DEAL::6 6 16.
+DEAL::6 7 20.
+DEAL::6 8 24.
+DEAL::6 9 0.0
+DEAL::7 0 0.0
+DEAL::7 1 0.0
+DEAL::7 2 0.0
+DEAL::7 3 0.0
+DEAL::7 4 10.
+DEAL::7 5 15.
+DEAL::7 6 20.
+DEAL::7 7 25.
+DEAL::7 8 30.
+DEAL::7 9 0.0
+DEAL::8 0 0.0
+DEAL::8 1 0.0
+DEAL::8 2 0.0
+DEAL::8 3 0.0
+DEAL::8 4 0.0
+DEAL::8 5 0.0
+DEAL::8 6 0.0
+DEAL::8 7 0.0
+DEAL::8 8 0.0
+DEAL::8 9 0.0
+DEAL::9 0 0.0
+DEAL::9 1 0.0
+DEAL::9 2 0.0
+DEAL::9 3 0.0
+DEAL::9 4 0.0
+DEAL::9 5 0.0
+DEAL::9 6 0.0
+DEAL::9 7 0.0
+DEAL::9 8 0.0
+DEAL::9 9 0.0
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check creation and output of a matrix
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_11/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m, n;
+ make_matrix (m);
+ make_matrix (n);
+ deallog << (m==n) << std::endl;
+}
+
--- /dev/null
+
+DEAL::1
+DEAL::1
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::all_zero
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_12/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m(3,3);
+ deallog << m.all_zero() << std::endl;
+
+ make_matrix (m);
+ deallog << m.all_zero() << std::endl;
+}
+
--- /dev/null
+
+DEAL::1
+DEAL::0
+DEAL::1
+DEAL::0
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::matrix_norm_square
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_13/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ Vector<double> v;
+ make_vector (v);
+
+ deallog << m.matrix_norm_square (v) << std::endl;
+}
+
--- /dev/null
+
+DEAL::3.8e+03
+DEAL::3.8e+03
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::matrix_scalar_product
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_14/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ Vector<double> v,w;
+ make_vector (v);
+ make_vector (w);
+ for (unsigned int i=0; i<w.size(); ++i)
+ w(i) = w(i)+1.;
+
+ deallog << m.matrix_scalar_product (v,w) << std::endl;
+}
+
--- /dev/null
+
+DEAL::5.0e+03
+DEAL::5.0e+03
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::l1_norm
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_15/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ deallog << m.l1_norm() << std::endl;
+}
+
--- /dev/null
+
+DEAL::90.
+DEAL::90.
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::linfty_norm
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_16/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ deallog << m.linfty_norm() << std::endl;
+}
+
--- /dev/null
+
+DEAL::1.0e+02
+DEAL::1.0e+02
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::frobenius_norm
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_17/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ deallog << m.frobenius_norm() << std::endl;
+}
+
--- /dev/null
+
+DEAL::70.
+DEAL::70.
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::relative_symmetry_norm2
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_18/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ deallog << m.relative_symmetry_norm2() << std::endl;
+}
+
--- /dev/null
+
+DEAL::0.14
+DEAL::0.14
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::determinant
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_19/output";
+
+
+template <typename number>
+void
+check ()
+{
+ for (unsigned int n=1; n<=3; ++n)
+ {
+ const number array[] = { 1, 2, 3, 4, 5, 6, 7, 8, 9 };
+
+ FullMatrix<number> m(n,n, array);
+ print_matrix (m);
+ deallog << m.determinant () << std::endl;
+ }
+}
+
+
--- /dev/null
+
+DEAL::0 0 1.0
+DEAL::1.0
+DEAL::0 0 1.0
+DEAL::0 1 2.0
+DEAL::1 0 3.0
+DEAL::1 1 4.0
+DEAL::-2.0
+DEAL::0 0 1.0
+DEAL::0 1 2.0
+DEAL::0 2 3.0
+DEAL::1 0 4.0
+DEAL::1 1 5.0
+DEAL::1 2 6.0
+DEAL::2 0 7.0
+DEAL::2 1 8.0
+DEAL::2 2 9.0
+DEAL::0
+DEAL::0 0 1.0
+DEAL::1.0
+DEAL::0 0 1.0
+DEAL::0 1 2.0
+DEAL::1 0 3.0
+DEAL::1 1 4.0
+DEAL::-2.0
+DEAL::0 0 1.0
+DEAL::0 1 2.0
+DEAL::0 2 3.0
+DEAL::1 0 4.0
+DEAL::1 1 5.0
+DEAL::1 2 6.0
+DEAL::2 0 7.0
+DEAL::2 1 8.0
+DEAL::2 2 9.0
+DEAL::0.0
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::trace
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_20/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ deallog << m.trace() << std::endl;
+}
+
--- /dev/null
+
+DEAL::70.
+DEAL::70.
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::print
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_21/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ m.print (deallog.get_file_stream());
+}
+
--- /dev/null
+
+2.0 3.0 4.0 5.0 6.0
+4.0 6.0 8.0 10. 12.
+6.0 9.0 12. 15. 18.
+8.0 12. 16. 20. 24.
+10. 15. 20. 25. 30.
+2.0 3.0 4.0 5.0 6.0
+4.0 6.0 8.0 10. 12.
+6.0 9.0 12. 15. 18.
+8.0 12. 16. 20. 24.
+10. 15. 20. 25. 30.
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::operator*=
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_22/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ m *= 3.1415;
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 6.3
+DEAL::0 1 9.4
+DEAL::0 2 13.
+DEAL::0 3 16.
+DEAL::0 4 19.
+DEAL::1 0 13.
+DEAL::1 1 19.
+DEAL::1 2 25.
+DEAL::1 3 31.
+DEAL::1 4 38.
+DEAL::2 0 19.
+DEAL::2 1 28.
+DEAL::2 2 38.
+DEAL::2 3 47.
+DEAL::2 4 57.
+DEAL::3 0 25.
+DEAL::3 1 38.
+DEAL::3 2 50.
+DEAL::3 3 63.
+DEAL::3 4 75.
+DEAL::4 0 31.
+DEAL::4 1 47.
+DEAL::4 2 63.
+DEAL::4 3 79.
+DEAL::4 4 94.
+DEAL::0 0 6.3
+DEAL::0 1 9.4
+DEAL::0 2 13.
+DEAL::0 3 16.
+DEAL::0 4 19.
+DEAL::1 0 13.
+DEAL::1 1 19.
+DEAL::1 2 25.
+DEAL::1 3 31.
+DEAL::1 4 38.
+DEAL::2 0 19.
+DEAL::2 1 28.
+DEAL::2 2 38.
+DEAL::2 3 47.
+DEAL::2 4 57.
+DEAL::3 0 25.
+DEAL::3 1 38.
+DEAL::3 2 50.
+DEAL::3 3 63.
+DEAL::3 4 75.
+DEAL::4 0 31.
+DEAL::4 1 47.
+DEAL::4 2 63.
+DEAL::4 3 79.
+DEAL::4 4 94.
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::operator/=
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_23/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ m /= 3.1415926;
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 0.64
+DEAL::0 1 0.95
+DEAL::0 2 1.3
+DEAL::0 3 1.6
+DEAL::0 4 1.9
+DEAL::1 0 1.3
+DEAL::1 1 1.9
+DEAL::1 2 2.5
+DEAL::1 3 3.2
+DEAL::1 4 3.8
+DEAL::2 0 1.9
+DEAL::2 1 2.9
+DEAL::2 2 3.8
+DEAL::2 3 4.8
+DEAL::2 4 5.7
+DEAL::3 0 2.5
+DEAL::3 1 3.8
+DEAL::3 2 5.1
+DEAL::3 3 6.4
+DEAL::3 4 7.6
+DEAL::4 0 3.2
+DEAL::4 1 4.8
+DEAL::4 2 6.4
+DEAL::4 3 8.0
+DEAL::4 4 9.5
+DEAL::0 0 0.64
+DEAL::0 1 0.95
+DEAL::0 2 1.3
+DEAL::0 3 1.6
+DEAL::0 4 1.9
+DEAL::1 0 1.3
+DEAL::1 1 1.9
+DEAL::1 2 2.5
+DEAL::1 3 3.2
+DEAL::1 4 3.8
+DEAL::2 0 1.9
+DEAL::2 1 2.9
+DEAL::2 2 3.8
+DEAL::2 3 4.8
+DEAL::2 4 5.7
+DEAL::3 0 2.5
+DEAL::3 1 3.8
+DEAL::3 2 5.1
+DEAL::3 3 6.4
+DEAL::3 4 7.6
+DEAL::4 0 3.2
+DEAL::4 1 4.8
+DEAL::4 2 6.4
+DEAL::4 3 8.0
+DEAL::4 4 9.5