--- /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::mmult
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_43/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m, n, o;
+ make_matrix (m);
+ make_matrix (n);
+ make_matrix (o);
+
+ m.mmult (n, o, true);
+ print_matrix (m);
+ print_matrix (n);
+ print_matrix (o);
+
+ m.mmult (n, o, false);
+ print_matrix (m);
+ print_matrix (n);
+ print_matrix (o);
+}
+
--- /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 1.4e+02
+DEAL::0 1 2.1e+02
+DEAL::0 2 2.8e+02
+DEAL::0 3 3.6e+02
+DEAL::0 4 4.3e+02
+DEAL::1 0 2.8e+02
+DEAL::1 1 4.3e+02
+DEAL::1 2 5.7e+02
+DEAL::1 3 7.1e+02
+DEAL::1 4 8.5e+02
+DEAL::2 0 4.3e+02
+DEAL::2 1 6.4e+02
+DEAL::2 2 8.5e+02
+DEAL::2 3 1.1e+03
+DEAL::2 4 1.3e+03
+DEAL::3 0 5.7e+02
+DEAL::3 1 8.5e+02
+DEAL::3 2 1.1e+03
+DEAL::3 3 1.4e+03
+DEAL::3 4 1.7e+03
+DEAL::4 0 7.1e+02
+DEAL::4 1 1.1e+03
+DEAL::4 2 1.4e+03
+DEAL::4 3 1.8e+03
+DEAL::4 4 2.1e+03
+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.
+DEAL::0 0 1.4e+02
+DEAL::0 1 2.1e+02
+DEAL::0 2 2.8e+02
+DEAL::0 3 3.5e+02
+DEAL::0 4 4.2e+02
+DEAL::1 0 2.8e+02
+DEAL::1 1 4.2e+02
+DEAL::1 2 5.6e+02
+DEAL::1 3 7.0e+02
+DEAL::1 4 8.4e+02
+DEAL::2 0 4.2e+02
+DEAL::2 1 6.3e+02
+DEAL::2 2 8.4e+02
+DEAL::2 3 1.0e+03
+DEAL::2 4 1.3e+03
+DEAL::3 0 5.6e+02
+DEAL::3 1 8.4e+02
+DEAL::3 2 1.1e+03
+DEAL::3 3 1.4e+03
+DEAL::3 4 1.7e+03
+DEAL::4 0 7.0e+02
+DEAL::4 1 1.0e+03
+DEAL::4 2 1.4e+03
+DEAL::4 3 1.8e+03
+DEAL::4 4 2.1e+03
+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.
+DEAL::0 0 1.4e+02
+DEAL::0 1 2.1e+02
+DEAL::0 2 2.8e+02
+DEAL::0 3 3.6e+02
+DEAL::0 4 4.3e+02
+DEAL::1 0 2.8e+02
+DEAL::1 1 4.3e+02
+DEAL::1 2 5.7e+02
+DEAL::1 3 7.1e+02
+DEAL::1 4 8.5e+02
+DEAL::2 0 4.3e+02
+DEAL::2 1 6.4e+02
+DEAL::2 2 8.5e+02
+DEAL::2 3 1.1e+03
+DEAL::2 4 1.3e+03
+DEAL::3 0 5.7e+02
+DEAL::3 1 8.5e+02
+DEAL::3 2 1.1e+03
+DEAL::3 3 1.4e+03
+DEAL::3 4 1.7e+03
+DEAL::4 0 7.1e+02
+DEAL::4 1 1.1e+03
+DEAL::4 2 1.4e+03
+DEAL::4 3 1.8e+03
+DEAL::4 4 2.1e+03
+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.
+DEAL::0 0 1.4e+02
+DEAL::0 1 2.1e+02
+DEAL::0 2 2.8e+02
+DEAL::0 3 3.5e+02
+DEAL::0 4 4.2e+02
+DEAL::1 0 2.8e+02
+DEAL::1 1 4.2e+02
+DEAL::1 2 5.6e+02
+DEAL::1 3 7.0e+02
+DEAL::1 4 8.4e+02
+DEAL::2 0 4.2e+02
+DEAL::2 1 6.3e+02
+DEAL::2 2 8.4e+02
+DEAL::2 3 1.0e+03
+DEAL::2 4 1.3e+03
+DEAL::3 0 5.6e+02
+DEAL::3 1 8.4e+02
+DEAL::3 2 1.1e+03
+DEAL::3 3 1.4e+03
+DEAL::3 4 1.7e+03
+DEAL::4 0 7.0e+02
+DEAL::4 1 1.0e+03
+DEAL::4 2 1.4e+03
+DEAL::4 3 1.8e+03
+DEAL::4 4 2.1e+03
+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 FullMatrix::Tmmult
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_44/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m, n, o;
+ make_matrix (m);
+ make_matrix (n);
+ make_matrix (o);
+
+ m.Tmmult (n, o, true);
+ print_matrix (m);
+ print_matrix (n);
+ print_matrix (o);
+
+ m.Tmmult (n, o, false);
+ print_matrix (m);
+ print_matrix (n);
+ print_matrix (o);
+}
+
--- /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.2e+02
+DEAL::0 1 3.3e+02
+DEAL::0 2 4.4e+02
+DEAL::0 3 5.6e+02
+DEAL::0 4 6.7e+02
+DEAL::1 0 3.3e+02
+DEAL::1 1 5.0e+02
+DEAL::1 2 6.7e+02
+DEAL::1 3 8.4e+02
+DEAL::1 4 1.0e+03
+DEAL::2 0 4.5e+02
+DEAL::2 1 6.7e+02
+DEAL::2 2 8.9e+02
+DEAL::2 3 1.1e+03
+DEAL::2 4 1.3e+03
+DEAL::3 0 5.6e+02
+DEAL::3 1 8.4e+02
+DEAL::3 2 1.1e+03
+DEAL::3 3 1.4e+03
+DEAL::3 4 1.7e+03
+DEAL::4 0 6.7e+02
+DEAL::4 1 1.0e+03
+DEAL::4 2 1.3e+03
+DEAL::4 3 1.7e+03
+DEAL::4 4 2.0e+03
+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.
+DEAL::0 0 2.2e+02
+DEAL::0 1 3.3e+02
+DEAL::0 2 4.4e+02
+DEAL::0 3 5.5e+02
+DEAL::0 4 6.6e+02
+DEAL::1 0 3.3e+02
+DEAL::1 1 5.0e+02
+DEAL::1 2 6.6e+02
+DEAL::1 3 8.2e+02
+DEAL::1 4 9.9e+02
+DEAL::2 0 4.4e+02
+DEAL::2 1 6.6e+02
+DEAL::2 2 8.8e+02
+DEAL::2 3 1.1e+03
+DEAL::2 4 1.3e+03
+DEAL::3 0 5.5e+02
+DEAL::3 1 8.2e+02
+DEAL::3 2 1.1e+03
+DEAL::3 3 1.4e+03
+DEAL::3 4 1.6e+03
+DEAL::4 0 6.6e+02
+DEAL::4 1 9.9e+02
+DEAL::4 2 1.3e+03
+DEAL::4 3 1.6e+03
+DEAL::4 4 2.0e+03
+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.
+DEAL::0 0 2.2e+02
+DEAL::0 1 3.3e+02
+DEAL::0 2 4.4e+02
+DEAL::0 3 5.6e+02
+DEAL::0 4 6.7e+02
+DEAL::1 0 3.3e+02
+DEAL::1 1 5.0e+02
+DEAL::1 2 6.7e+02
+DEAL::1 3 8.4e+02
+DEAL::1 4 1.0e+03
+DEAL::2 0 4.5e+02
+DEAL::2 1 6.7e+02
+DEAL::2 2 8.9e+02
+DEAL::2 3 1.1e+03
+DEAL::2 4 1.3e+03
+DEAL::3 0 5.6e+02
+DEAL::3 1 8.4e+02
+DEAL::3 2 1.1e+03
+DEAL::3 3 1.4e+03
+DEAL::3 4 1.7e+03
+DEAL::4 0 6.7e+02
+DEAL::4 1 1.0e+03
+DEAL::4 2 1.3e+03
+DEAL::4 3 1.7e+03
+DEAL::4 4 2.0e+03
+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.
+DEAL::0 0 2.2e+02
+DEAL::0 1 3.3e+02
+DEAL::0 2 4.4e+02
+DEAL::0 3 5.5e+02
+DEAL::0 4 6.6e+02
+DEAL::1 0 3.3e+02
+DEAL::1 1 5.0e+02
+DEAL::1 2 6.6e+02
+DEAL::1 3 8.2e+02
+DEAL::1 4 9.9e+02
+DEAL::2 0 4.4e+02
+DEAL::2 1 6.6e+02
+DEAL::2 2 8.8e+02
+DEAL::2 3 1.1e+03
+DEAL::2 4 1.3e+03
+DEAL::3 0 5.5e+02
+DEAL::3 1 8.2e+02
+DEAL::3 2 1.1e+03
+DEAL::3 3 1.4e+03
+DEAL::3 4 1.6e+03
+DEAL::4 0 6.6e+02
+DEAL::4 1 9.9e+02
+DEAL::4 2 1.3e+03
+DEAL::4 3 1.6e+03
+DEAL::4 4 2.0e+03
+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 FullMatrix::vmult
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_45/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ Vector<double> v, w;
+ make_vector (v);
+ make_vector (w);
+
+ m.vmult (v, w, true);
+ print_vector (v);
+ print_vector (w);
+
+ m.vmult (v, w, false);
+ print_vector (v);
+ print_vector (w);
+}
+
--- /dev/null
+
+DEAL::0 71.
+DEAL::1 1.4e+02
+DEAL::2 2.1e+02
+DEAL::3 2.8e+02
+DEAL::4 3.6e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 70.
+DEAL::1 1.4e+02
+DEAL::2 2.1e+02
+DEAL::3 2.8e+02
+DEAL::4 3.5e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 71.
+DEAL::1 1.4e+02
+DEAL::2 2.1e+02
+DEAL::3 2.8e+02
+DEAL::4 3.6e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 70.
+DEAL::1 1.4e+02
+DEAL::2 2.1e+02
+DEAL::3 2.8e+02
+DEAL::4 3.5e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::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 FullMatrix::vmult_add
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_46/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ Vector<double> v, w;
+ make_vector (v);
+ make_vector (w);
+
+ m.vmult_add (v, w);
+ print_vector (v);
+ print_vector (w);
+}
+
--- /dev/null
+
+DEAL::0 71.
+DEAL::1 1.4e+02
+DEAL::2 2.1e+02
+DEAL::3 2.8e+02
+DEAL::4 3.6e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 71.
+DEAL::1 1.4e+02
+DEAL::2 2.1e+02
+DEAL::3 2.8e+02
+DEAL::4 3.6e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::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 FullMatrix::Tvmult
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_47/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ Vector<double> v, w;
+ make_vector (v);
+ make_vector (w);
+
+ m.Tvmult (v, w, true);
+ print_vector (v);
+ print_vector (w);
+
+ m.Tvmult (v, w, false);
+ print_vector (v);
+ print_vector (w);
+}
+
--- /dev/null
+
+DEAL::0 1.1e+02
+DEAL::1 1.7e+02
+DEAL::2 2.2e+02
+DEAL::3 2.8e+02
+DEAL::4 3.4e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 1.1e+02
+DEAL::1 1.6e+02
+DEAL::2 2.2e+02
+DEAL::3 2.8e+02
+DEAL::4 3.3e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 1.1e+02
+DEAL::1 1.7e+02
+DEAL::2 2.2e+02
+DEAL::3 2.8e+02
+DEAL::4 3.4e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 1.1e+02
+DEAL::1 1.6e+02
+DEAL::2 2.2e+02
+DEAL::3 2.8e+02
+DEAL::4 3.3e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::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 FullMatrix::Tvmult_add
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_48/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ Vector<double> v, w;
+ make_vector (v);
+ make_vector (w);
+
+ m.Tvmult_add (v, w);
+ print_vector (v);
+ print_vector (w);
+}
+
--- /dev/null
+
+DEAL::0 1.1e+02
+DEAL::1 1.7e+02
+DEAL::2 2.2e+02
+DEAL::3 2.8e+02
+DEAL::4 3.4e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 1.1e+02
+DEAL::1 1.7e+02
+DEAL::2 2.2e+02
+DEAL::3 2.8e+02
+DEAL::4 3.4e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::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 FullMatrix::precondition_Jacobi
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_49/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ Vector<double> v, w;
+ make_vector (v);
+ make_vector (w);
+
+ m.precondition_Jacobi (v, w, 3.141);
+ print_vector (v);
+ print_vector (w);
+}
+
--- /dev/null
+
+DEAL::0 1.6
+DEAL::1 1.0
+DEAL::2 0.79
+DEAL::3 0.63
+DEAL::4 0.52
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 1.6
+DEAL::1 1.0
+DEAL::2 0.79
+DEAL::3 0.63
+DEAL::4 0.52
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::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 FullMatrix::residual
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_50/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ Vector<double> v, w, x;
+ make_vector (v);
+ make_vector (w);
+ make_vector (x);
+
+ m.residual (v, w, x);
+ print_vector (v);
+ print_vector (w);
+ print_vector (x);
+}
+
--- /dev/null
+
+DEAL::0 -69.
+DEAL::1 -1.4e+02
+DEAL::2 -2.1e+02
+DEAL::3 -2.8e+02
+DEAL::4 -3.4e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 -69.
+DEAL::1 -1.4e+02
+DEAL::2 -2.1e+02
+DEAL::3 -2.8e+02
+DEAL::4 -3.4e+02
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::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 FullMatrix::forward
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_51/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ Vector<double> x, y;
+ make_vector (x);
+ make_vector (y);
+
+ m.forward (x, y);
+ print_vector (x);
+ print_vector (y);
+}
+
--- /dev/null
+
+DEAL::0 0.50
+DEAL::1 0
+DEAL::2 0
+DEAL::3 0
+DEAL::4 0
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 0.50
+DEAL::1 0
+DEAL::2 0
+DEAL::3 0
+DEAL::4 0
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::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 FullMatrix::backward
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "full_matrix_52/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<number> m;
+ make_matrix (m);
+ Vector<double> x, y;
+ make_vector (x);
+ make_vector (y);
+
+ m.backward (x, y);
+ print_vector (x);
+ print_vector (y);
+}
+
--- /dev/null
+
+DEAL::0 0
+DEAL::1 0
+DEAL::2 0
+DEAL::3 0
+DEAL::4 0.17
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+DEAL::0 0
+DEAL::1 0
+DEAL::2 0
+DEAL::3 0
+DEAL::4 0.17
+DEAL::0 1.0
+DEAL::1 2.0
+DEAL::2 3.0
+DEAL::3 4.0
+DEAL::4 5.0
+template <typename number>
+void
+print_vector (const Vector<number> &v)
+{
+ for (unsigned int i=0; i<v.size(); ++i)
+ deallog << i << ' ' << v(i)
+ << std::endl;
+}
+
+
+
int
main()
{