grid_in_* \
point_inside_? \
full_matrix_* \
+ complex_real_full_matrix_* \
block_matrix_array_* \
data_out_0* \
data_out_faces_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. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_01/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,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 running over const iterators. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_02/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+
+ for (typename FullMatrix<std::complex<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,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,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 running over const iterators starting at the second line. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_03/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+
+
+ for (typename FullMatrix<std::complex<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,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,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 using an array. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_04/output";
+
+
+template <typename number>
+void
+check ()
+{
+ const std::complex<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<std::complex<number> > m (5,5,array);
+
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (1.0,0.0)
+DEAL::0 1 (2.0,0.0)
+DEAL::0 2 (3.0,0.0)
+DEAL::0 3 (4.0,0.0)
+DEAL::0 4 (5.0,0.0)
+DEAL::1 0 (6.0,0.0)
+DEAL::1 1 (7.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (9.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::2 0 (1.0,0.0)
+DEAL::2 1 (2.0,0.0)
+DEAL::2 2 (3.0,0.0)
+DEAL::2 3 (4.0,0.0)
+DEAL::2 4 (5.0,0.0)
+DEAL::3 0 (6.0,0.0)
+DEAL::3 1 (7.0,0.0)
+DEAL::3 2 (8.0,0.0)
+DEAL::3 3 (9.0,0.0)
+DEAL::3 4 (0.0,0.0)
+DEAL::4 0 (1.0,0.0)
+DEAL::4 1 (2.0,0.0)
+DEAL::4 2 (3.0,0.0)
+DEAL::4 3 (4.0,0.0)
+DEAL::4 4 (5.0,0.0)
+DEAL::0 0 (1.0,0.0)
+DEAL::0 1 (2.0,0.0)
+DEAL::0 2 (3.0,0.0)
+DEAL::0 3 (4.0,0.0)
+DEAL::0 4 (5.0,0.0)
+DEAL::1 0 (6.0,0.0)
+DEAL::1 1 (7.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (9.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::2 0 (1.0,0.0)
+DEAL::2 1 (2.0,0.0)
+DEAL::2 2 (3.0,0.0)
+DEAL::2 3 (4.0,0.0)
+DEAL::2 4 (5.0,0.0)
+DEAL::3 0 (6.0,0.0)
+DEAL::3 1 (7.0,0.0)
+DEAL::3 2 (8.0,0.0)
+DEAL::3 3 (9.0,0.0)
+DEAL::3 4 (0.0,0.0)
+DEAL::4 0 (1.0,0.0)
+DEAL::4 1 (2.0,0.0)
+DEAL::4 2 (3.0,0.0)
+DEAL::4 3 (4.0,0.0)
+DEAL::4 4 (5.0,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 using IdentityMatrix. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_05/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m (IdentityMatrix(5));
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (1.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (1.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (1.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (1.0,0.0)
+DEAL::3 4 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (1.0,0.0)
+DEAL::0 0 (1.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (1.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (1.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (1.0,0.0)
+DEAL::3 4 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (1.0,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 matrix assignment. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_06/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+
+ FullMatrix<std::complex<number> > n;
+ n = m;
+
+ print_matrix (n);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,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 matrix=0. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_07/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m = std::complex<number>(0);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (0.0,0.0)
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (0.0,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. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_08/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+
+ FullMatrix<std::complex<number> > n;
+ n.copy_from (m);
+
+ print_matrix (n);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,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 assignment from IdenityMatrix. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_09/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ m = IdentityMatrix(5);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (1.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (1.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (1.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (1.0,0.0)
+DEAL::3 4 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (1.0,0.0)
+DEAL::0 0 (1.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (1.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (1.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (1.0,0.0)
+DEAL::3 4 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (1.0,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::fill. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_10/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+
+ FullMatrix<std::complex<number> > n (10, 10);
+ n.fill (m, 3, 4);
+ print_matrix (n);
+}
+
--- /dev/null
+
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (2.0,0.0)
+DEAL::3 5 (3.0,0.0)
+DEAL::3 6 (4.0,0.0)
+DEAL::3 7 (5.0,0.0)
+DEAL::3 8 (6.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (4.0,0.0)
+DEAL::4 5 (6.0,0.0)
+DEAL::4 6 (8.0,0.0)
+DEAL::4 7 (10.,0.0)
+DEAL::4 8 (12.,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (6.0,0.0)
+DEAL::5 5 (9.0,0.0)
+DEAL::5 6 (12.,0.0)
+DEAL::5 7 (15.,0.0)
+DEAL::5 8 (18.,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (8.0,0.0)
+DEAL::6 5 (12.,0.0)
+DEAL::6 6 (16.,0.0)
+DEAL::6 7 (20.,0.0)
+DEAL::6 8 (24.,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (10.,0.0)
+DEAL::7 5 (15.,0.0)
+DEAL::7 6 (20.,0.0)
+DEAL::7 7 (25.,0.0)
+DEAL::7 8 (30.,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,0.0)
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (2.0,0.0)
+DEAL::3 5 (3.0,0.0)
+DEAL::3 6 (4.0,0.0)
+DEAL::3 7 (5.0,0.0)
+DEAL::3 8 (6.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (4.0,0.0)
+DEAL::4 5 (6.0,0.0)
+DEAL::4 6 (8.0,0.0)
+DEAL::4 7 (10.,0.0)
+DEAL::4 8 (12.,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (6.0,0.0)
+DEAL::5 5 (9.0,0.0)
+DEAL::5 6 (12.,0.0)
+DEAL::5 7 (15.,0.0)
+DEAL::5 8 (18.,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (8.0,0.0)
+DEAL::6 5 (12.,0.0)
+DEAL::6 6 (16.,0.0)
+DEAL::6 7 (20.,0.0)
+DEAL::6 8 (24.,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (10.,0.0)
+DEAL::7 5 (15.,0.0)
+DEAL::7 6 (20.,0.0)
+DEAL::7 7 (25.,0.0)
+DEAL::7 8 (30.,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,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. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_11/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<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. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_12/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<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. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_13/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ Vector<std::complex<number> > v;
+ make_vector (v);
+
+ deallog << m.matrix_norm_square (v) << std::endl;
+}
+
--- /dev/null
+
+DEAL::(3.8e+03,0.0)
+DEAL::(3.8e+03,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::matrix_scalar_product. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_14/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ Vector<std::complex<number> > v,w;
+ make_vector (v);
+ make_vector (w);
+ for (unsigned int i=0; i<w.size(); ++i)
+ w(i) = w(i)+std::complex<number>(1.);
+
+ deallog << m.matrix_scalar_product (v,w) << std::endl;
+}
+
--- /dev/null
+
+DEAL::(5.0e+03,0.0)
+DEAL::(5.0e+03,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::l1_norm. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_15/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<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. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_16/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<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. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_17/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<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. like the full_matrix_* tests,
+// but use complex-valued matrices and vectors, even though we only store real
+// values in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_18/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<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. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_19/output";
+
+
+template <typename number>
+void
+check ()
+{
+ for (unsigned int n=1; n<=3; ++n)
+ {
+ const std::complex<number> array[] = { 1, 2, 3, 4, 5, 6, 7, 8, 9 };
+
+ FullMatrix<std::complex<number> > m(n,n, array);
+ print_matrix (m);
+ deallog << m.determinant () << std::endl;
+ }
+}
+
+
--- /dev/null
+
+DEAL::0 0 (1.0,0.0)
+DEAL::(1.0,0.0)
+DEAL::0 0 (1.0,0.0)
+DEAL::0 1 (2.0,0.0)
+DEAL::1 0 (3.0,0.0)
+DEAL::1 1 (4.0,0.0)
+DEAL::(-2.0,0.0)
+DEAL::0 0 (1.0,0.0)
+DEAL::0 1 (2.0,0.0)
+DEAL::0 2 (3.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (5.0,0.0)
+DEAL::1 2 (6.0,0.0)
+DEAL::2 0 (7.0,0.0)
+DEAL::2 1 (8.0,0.0)
+DEAL::2 2 (9.0,0.0)
+DEAL::(0.0,0.0)
+DEAL::0 0 (1.0,0.0)
+DEAL::(1.0,0.0)
+DEAL::0 0 (1.0,0.0)
+DEAL::0 1 (2.0,0.0)
+DEAL::1 0 (3.0,0.0)
+DEAL::1 1 (4.0,0.0)
+DEAL::(-2.0,0.0)
+DEAL::0 0 (1.0,0.0)
+DEAL::0 1 (2.0,0.0)
+DEAL::0 2 (3.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (5.0,0.0)
+DEAL::1 2 (6.0,0.0)
+DEAL::2 0 (7.0,0.0)
+DEAL::2 1 (8.0,0.0)
+DEAL::2 2 (9.0,0.0)
+DEAL::(0.0,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. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_20/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ deallog << m.trace() << std::endl;
+}
+
--- /dev/null
+
+DEAL::(70.,0.0)
+DEAL::(70.,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::print. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_21/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m.print (deallog.get_file_stream());
+}
+
--- /dev/null
+
+(2.0,0.0)(3.0,0.0)(4.0,0.0)(5.0,0.0)(6.0,0.0)
+(4.0,0.0)(6.0,0.0)(8.0,0.0)(10.,0.0)(12.,0.0)
+(6.0,0.0)(9.0,0.0)(12.,0.0)(15.,0.0)(18.,0.0)
+(8.0,0.0)(12.,0.0)(16.,0.0)(20.,0.0)(24.,0.0)
+(10.,0.0)(15.,0.0)(20.,0.0)(25.,0.0)(30.,0.0)
+(2.0,0.0)(3.0,0.0)(4.0,0.0)(5.0,0.0)(6.0,0.0)
+(4.0,0.0)(6.0,0.0)(8.0,0.0)(10.,0.0)(12.,0.0)
+(6.0,0.0)(9.0,0.0)(12.,0.0)(15.,0.0)(18.,0.0)
+(8.0,0.0)(12.,0.0)(16.,0.0)(20.,0.0)(24.,0.0)
+(10.,0.0)(15.,0.0)(20.,0.0)(25.,0.0)(30.,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::operator*=. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_22/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m *= 3.1415;
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (6.3,0.0)
+DEAL::0 1 (9.4,0.0)
+DEAL::0 2 (13.,0.0)
+DEAL::0 3 (16.,0.0)
+DEAL::0 4 (19.,0.0)
+DEAL::1 0 (13.,0.0)
+DEAL::1 1 (19.,0.0)
+DEAL::1 2 (25.,0.0)
+DEAL::1 3 (31.,0.0)
+DEAL::1 4 (38.,0.0)
+DEAL::2 0 (19.,0.0)
+DEAL::2 1 (28.,0.0)
+DEAL::2 2 (38.,0.0)
+DEAL::2 3 (47.,0.0)
+DEAL::2 4 (57.,0.0)
+DEAL::3 0 (25.,0.0)
+DEAL::3 1 (38.,0.0)
+DEAL::3 2 (50.,0.0)
+DEAL::3 3 (63.,0.0)
+DEAL::3 4 (75.,0.0)
+DEAL::4 0 (31.,0.0)
+DEAL::4 1 (47.,0.0)
+DEAL::4 2 (63.,0.0)
+DEAL::4 3 (79.,0.0)
+DEAL::4 4 (94.,0.0)
+DEAL::0 0 (6.3,0.0)
+DEAL::0 1 (9.4,0.0)
+DEAL::0 2 (13.,0.0)
+DEAL::0 3 (16.,0.0)
+DEAL::0 4 (19.,0.0)
+DEAL::1 0 (13.,0.0)
+DEAL::1 1 (19.,0.0)
+DEAL::1 2 (25.,0.0)
+DEAL::1 3 (31.,0.0)
+DEAL::1 4 (38.,0.0)
+DEAL::2 0 (19.,0.0)
+DEAL::2 1 (28.,0.0)
+DEAL::2 2 (38.,0.0)
+DEAL::2 3 (47.,0.0)
+DEAL::2 4 (57.,0.0)
+DEAL::3 0 (25.,0.0)
+DEAL::3 1 (38.,0.0)
+DEAL::3 2 (50.,0.0)
+DEAL::3 3 (63.,0.0)
+DEAL::3 4 (75.,0.0)
+DEAL::4 0 (31.,0.0)
+DEAL::4 1 (47.,0.0)
+DEAL::4 2 (63.,0.0)
+DEAL::4 3 (79.,0.0)
+DEAL::4 4 (94.,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::operator/=. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_23/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m /= 3.1415926;
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (0.64,0.0)
+DEAL::0 1 (0.95,0.0)
+DEAL::0 2 (1.3,0.0)
+DEAL::0 3 (1.6,0.0)
+DEAL::0 4 (1.9,0.0)
+DEAL::1 0 (1.3,0.0)
+DEAL::1 1 (1.9,0.0)
+DEAL::1 2 (2.5,0.0)
+DEAL::1 3 (3.2,0.0)
+DEAL::1 4 (3.8,0.0)
+DEAL::2 0 (1.9,0.0)
+DEAL::2 1 (2.9,0.0)
+DEAL::2 2 (3.8,0.0)
+DEAL::2 3 (4.8,0.0)
+DEAL::2 4 (5.7,0.0)
+DEAL::3 0 (2.5,0.0)
+DEAL::3 1 (3.8,0.0)
+DEAL::3 2 (5.1,0.0)
+DEAL::3 3 (6.4,0.0)
+DEAL::3 4 (7.6,0.0)
+DEAL::4 0 (3.2,0.0)
+DEAL::4 1 (4.8,0.0)
+DEAL::4 2 (6.4,0.0)
+DEAL::4 3 (8.0,0.0)
+DEAL::4 4 (9.5,0.0)
+DEAL::0 0 (0.64,0.0)
+DEAL::0 1 (0.95,0.0)
+DEAL::0 2 (1.3,0.0)
+DEAL::0 3 (1.6,0.0)
+DEAL::0 4 (1.9,0.0)
+DEAL::1 0 (1.3,0.0)
+DEAL::1 1 (1.9,0.0)
+DEAL::1 2 (2.5,0.0)
+DEAL::1 3 (3.2,0.0)
+DEAL::1 4 (3.8,0.0)
+DEAL::2 0 (1.9,0.0)
+DEAL::2 1 (2.9,0.0)
+DEAL::2 2 (3.8,0.0)
+DEAL::2 3 (4.8,0.0)
+DEAL::2 4 (5.7,0.0)
+DEAL::3 0 (2.5,0.0)
+DEAL::3 1 (3.8,0.0)
+DEAL::3 2 (5.1,0.0)
+DEAL::3 3 (6.4,0.0)
+DEAL::3 4 (7.6,0.0)
+DEAL::4 0 (3.2,0.0)
+DEAL::4 1 (4.8,0.0)
+DEAL::4 2 (6.4,0.0)
+DEAL::4 3 (8.0,0.0)
+DEAL::4 4 (9.5,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::add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_24/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n;
+ make_matrix (m);
+ make_matrix (n);
+ m.add (3.14159, n);
+
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (8.3,0.0)
+DEAL::0 1 (12.,0.0)
+DEAL::0 2 (17.,0.0)
+DEAL::0 3 (21.,0.0)
+DEAL::0 4 (25.,0.0)
+DEAL::1 0 (17.,0.0)
+DEAL::1 1 (25.,0.0)
+DEAL::1 2 (33.,0.0)
+DEAL::1 3 (41.,0.0)
+DEAL::1 4 (50.,0.0)
+DEAL::2 0 (25.,0.0)
+DEAL::2 1 (37.,0.0)
+DEAL::2 2 (50.,0.0)
+DEAL::2 3 (62.,0.0)
+DEAL::2 4 (75.,0.0)
+DEAL::3 0 (33.,0.0)
+DEAL::3 1 (50.,0.0)
+DEAL::3 2 (66.,0.0)
+DEAL::3 3 (83.,0.0)
+DEAL::3 4 (99.,0.0)
+DEAL::4 0 (41.,0.0)
+DEAL::4 1 (62.,0.0)
+DEAL::4 2 (83.,0.0)
+DEAL::4 3 (1.0e+02,0.0)
+DEAL::4 4 (1.2e+02,0.0)
+DEAL::0 0 (8.3,0.0)
+DEAL::0 1 (12.,0.0)
+DEAL::0 2 (17.,0.0)
+DEAL::0 3 (21.,0.0)
+DEAL::0 4 (25.,0.0)
+DEAL::1 0 (17.,0.0)
+DEAL::1 1 (25.,0.0)
+DEAL::1 2 (33.,0.0)
+DEAL::1 3 (41.,0.0)
+DEAL::1 4 (50.,0.0)
+DEAL::2 0 (25.,0.0)
+DEAL::2 1 (37.,0.0)
+DEAL::2 2 (50.,0.0)
+DEAL::2 3 (62.,0.0)
+DEAL::2 4 (75.,0.0)
+DEAL::3 0 (33.,0.0)
+DEAL::3 1 (50.,0.0)
+DEAL::3 2 (66.,0.0)
+DEAL::3 3 (83.,0.0)
+DEAL::3 4 (99.,0.0)
+DEAL::4 0 (41.,0.0)
+DEAL::4 1 (62.,0.0)
+DEAL::4 2 (83.,0.0)
+DEAL::4 3 (1.0e+02,0.0)
+DEAL::4 4 (1.2e+02,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::add(4). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_25/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n,o;
+ make_matrix (m);
+ make_matrix (n);
+ make_matrix (o);
+ m.add (3.1415, n, 2.718, o);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (14.,0.0)
+DEAL::0 1 (21.,0.0)
+DEAL::0 2 (27.,0.0)
+DEAL::0 3 (34.,0.0)
+DEAL::0 4 (41.,0.0)
+DEAL::1 0 (27.,0.0)
+DEAL::1 1 (41.,0.0)
+DEAL::1 2 (55.,0.0)
+DEAL::1 3 (69.,0.0)
+DEAL::1 4 (82.,0.0)
+DEAL::2 0 (41.,0.0)
+DEAL::2 1 (62.,0.0)
+DEAL::2 2 (82.,0.0)
+DEAL::2 3 (1.0e+02,0.0)
+DEAL::2 4 (1.2e+02,0.0)
+DEAL::3 0 (55.,0.0)
+DEAL::3 1 (82.,0.0)
+DEAL::3 2 (1.1e+02,0.0)
+DEAL::3 3 (1.4e+02,0.0)
+DEAL::3 4 (1.6e+02,0.0)
+DEAL::4 0 (69.,0.0)
+DEAL::4 1 (1.0e+02,0.0)
+DEAL::4 2 (1.4e+02,0.0)
+DEAL::4 3 (1.7e+02,0.0)
+DEAL::4 4 (2.1e+02,0.0)
+DEAL::0 0 (14.,0.0)
+DEAL::0 1 (21.,0.0)
+DEAL::0 2 (27.,0.0)
+DEAL::0 3 (34.,0.0)
+DEAL::0 4 (41.,0.0)
+DEAL::1 0 (27.,0.0)
+DEAL::1 1 (41.,0.0)
+DEAL::1 2 (55.,0.0)
+DEAL::1 3 (69.,0.0)
+DEAL::1 4 (82.,0.0)
+DEAL::2 0 (41.,0.0)
+DEAL::2 1 (62.,0.0)
+DEAL::2 2 (82.,0.0)
+DEAL::2 3 (1.0e+02,0.0)
+DEAL::2 4 (1.2e+02,0.0)
+DEAL::3 0 (55.,0.0)
+DEAL::3 1 (82.,0.0)
+DEAL::3 2 (1.1e+02,0.0)
+DEAL::3 3 (1.4e+02,0.0)
+DEAL::3 4 (1.6e+02,0.0)
+DEAL::4 0 (69.,0.0)
+DEAL::4 1 (1.0e+02,0.0)
+DEAL::4 2 (1.4e+02,0.0)
+DEAL::4 3 (1.7e+02,0.0)
+DEAL::4 4 (2.1e+02,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::add (6). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_26/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n,o,p;
+ make_matrix (m);
+ make_matrix (n);
+ make_matrix (o);
+ make_matrix (p);
+ m.add (3.1415, n, 2.718, o, 1.414, p);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (17.,0.0)
+DEAL::0 1 (25.,0.0)
+DEAL::0 2 (33.,0.0)
+DEAL::0 3 (41.,0.0)
+DEAL::0 4 (50.,0.0)
+DEAL::1 0 (33.,0.0)
+DEAL::1 1 (50.,0.0)
+DEAL::1 2 (66.,0.0)
+DEAL::1 3 (83.,0.0)
+DEAL::1 4 (99.,0.0)
+DEAL::2 0 (50.,0.0)
+DEAL::2 1 (74.,0.0)
+DEAL::2 2 (99.,0.0)
+DEAL::2 3 (1.2e+02,0.0)
+DEAL::2 4 (1.5e+02,0.0)
+DEAL::3 0 (66.,0.0)
+DEAL::3 1 (99.,0.0)
+DEAL::3 2 (1.3e+02,0.0)
+DEAL::3 3 (1.7e+02,0.0)
+DEAL::3 4 (2.0e+02,0.0)
+DEAL::4 0 (83.,0.0)
+DEAL::4 1 (1.2e+02,0.0)
+DEAL::4 2 (1.7e+02,0.0)
+DEAL::4 3 (2.1e+02,0.0)
+DEAL::4 4 (2.5e+02,0.0)
+DEAL::0 0 (17.,0.0)
+DEAL::0 1 (25.,0.0)
+DEAL::0 2 (33.,0.0)
+DEAL::0 3 (41.,0.0)
+DEAL::0 4 (50.,0.0)
+DEAL::1 0 (33.,0.0)
+DEAL::1 1 (50.,0.0)
+DEAL::1 2 (66.,0.0)
+DEAL::1 3 (83.,0.0)
+DEAL::1 4 (99.,0.0)
+DEAL::2 0 (50.,0.0)
+DEAL::2 1 (74.,0.0)
+DEAL::2 2 (99.,0.0)
+DEAL::2 3 (1.2e+02,0.0)
+DEAL::2 4 (1.5e+02,0.0)
+DEAL::3 0 (66.,0.0)
+DEAL::3 1 (99.,0.0)
+DEAL::3 2 (1.3e+02,0.0)
+DEAL::3 3 (1.7e+02,0.0)
+DEAL::3 4 (2.0e+02,0.0)
+DEAL::4 0 (83.,0.0)
+DEAL::4 1 (1.2e+02,0.0)
+DEAL::4 2 (1.7e+02,0.0)
+DEAL::4 3 (2.1e+02,0.0)
+DEAL::4 4 (2.5e+02,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::add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_27/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m(10,10),n;
+ make_matrix (n);
+ m.add (n, 2., 3, 4, 2, 3);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (30.,0.0)
+DEAL::3 5 (36.,0.0)
+DEAL::3 6 (0.0,0.0)
+DEAL::3 7 (0.0,0.0)
+DEAL::3 8 (0.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (40.,0.0)
+DEAL::4 5 (48.,0.0)
+DEAL::4 6 (0.0,0.0)
+DEAL::4 7 (0.0,0.0)
+DEAL::4 8 (0.0,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (50.,0.0)
+DEAL::5 5 (60.,0.0)
+DEAL::5 6 (0.0,0.0)
+DEAL::5 7 (0.0,0.0)
+DEAL::5 8 (0.0,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (0.0,0.0)
+DEAL::6 5 (0.0,0.0)
+DEAL::6 6 (0.0,0.0)
+DEAL::6 7 (0.0,0.0)
+DEAL::6 8 (0.0,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (0.0,0.0)
+DEAL::7 5 (0.0,0.0)
+DEAL::7 6 (0.0,0.0)
+DEAL::7 7 (0.0,0.0)
+DEAL::7 8 (0.0,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,0.0)
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (30.,0.0)
+DEAL::3 5 (36.,0.0)
+DEAL::3 6 (0.0,0.0)
+DEAL::3 7 (0.0,0.0)
+DEAL::3 8 (0.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (40.,0.0)
+DEAL::4 5 (48.,0.0)
+DEAL::4 6 (0.0,0.0)
+DEAL::4 7 (0.0,0.0)
+DEAL::4 8 (0.0,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (50.,0.0)
+DEAL::5 5 (60.,0.0)
+DEAL::5 6 (0.0,0.0)
+DEAL::5 7 (0.0,0.0)
+DEAL::5 8 (0.0,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (0.0,0.0)
+DEAL::6 5 (0.0,0.0)
+DEAL::6 6 (0.0,0.0)
+DEAL::6 7 (0.0,0.0)
+DEAL::6 8 (0.0,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (0.0,0.0)
+DEAL::7 5 (0.0,0.0)
+DEAL::7 6 (0.0,0.0)
+DEAL::7 7 (0.0,0.0)
+DEAL::7 8 (0.0,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,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::Tadd. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_28/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m(10,10),n;
+ make_matrix (n);
+ m.Tadd (n, 2., 3, 4, 2, 3);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (30.,0.0)
+DEAL::3 5 (40.,0.0)
+DEAL::3 6 (50.,0.0)
+DEAL::3 7 (0.0,0.0)
+DEAL::3 8 (0.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (36.,0.0)
+DEAL::4 5 (48.,0.0)
+DEAL::4 6 (60.,0.0)
+DEAL::4 7 (0.0,0.0)
+DEAL::4 8 (0.0,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (0.0,0.0)
+DEAL::5 5 (0.0,0.0)
+DEAL::5 6 (0.0,0.0)
+DEAL::5 7 (0.0,0.0)
+DEAL::5 8 (0.0,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (0.0,0.0)
+DEAL::6 5 (0.0,0.0)
+DEAL::6 6 (0.0,0.0)
+DEAL::6 7 (0.0,0.0)
+DEAL::6 8 (0.0,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (0.0,0.0)
+DEAL::7 5 (0.0,0.0)
+DEAL::7 6 (0.0,0.0)
+DEAL::7 7 (0.0,0.0)
+DEAL::7 8 (0.0,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,0.0)
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (30.,0.0)
+DEAL::3 5 (40.,0.0)
+DEAL::3 6 (50.,0.0)
+DEAL::3 7 (0.0,0.0)
+DEAL::3 8 (0.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (36.,0.0)
+DEAL::4 5 (48.,0.0)
+DEAL::4 6 (60.,0.0)
+DEAL::4 7 (0.0,0.0)
+DEAL::4 8 (0.0,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (0.0,0.0)
+DEAL::5 5 (0.0,0.0)
+DEAL::5 6 (0.0,0.0)
+DEAL::5 7 (0.0,0.0)
+DEAL::5 8 (0.0,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (0.0,0.0)
+DEAL::6 5 (0.0,0.0)
+DEAL::6 6 (0.0,0.0)
+DEAL::6 7 (0.0,0.0)
+DEAL::6 8 (0.0,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (0.0,0.0)
+DEAL::7 5 (0.0,0.0)
+DEAL::7 6 (0.0,0.0)
+DEAL::7 7 (0.0,0.0)
+DEAL::7 8 (0.0,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,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::add_row. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_29/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m.add_row (2, 3.141, 4);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (37.,0.0)
+DEAL::2 1 (56.,0.0)
+DEAL::2 2 (75.,0.0)
+DEAL::2 3 (94.,0.0)
+DEAL::2 4 (1.1e+02,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (37.,0.0)
+DEAL::2 1 (56.,0.0)
+DEAL::2 2 (75.,0.0)
+DEAL::2 3 (94.,0.0)
+DEAL::2 4 (1.1e+02,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,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::add_row. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_30/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m.add_row (2, 3.141, 4, 2.7182, 3);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (59.,0.0)
+DEAL::2 1 (89.,0.0)
+DEAL::2 2 (1.2e+02,0.0)
+DEAL::2 3 (1.5e+02,0.0)
+DEAL::2 4 (1.8e+02,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (59.,0.0)
+DEAL::2 1 (89.,0.0)
+DEAL::2 2 (1.2e+02,0.0)
+DEAL::2 3 (1.5e+02,0.0)
+DEAL::2 4 (1.8e+02,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,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::add_col. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_31/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m.add_col (2, 3.141, 4);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (23.,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (46.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (69.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (91.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (1.1e+02,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (23.,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (46.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (69.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (91.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (1.1e+02,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,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::add_col. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_32/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m.add_col (2, 3.141, 4, 2.718, 3);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (36.,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (73.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (1.1e+02,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (1.5e+02,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (1.8e+02,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (36.,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (73.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (1.1e+02,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (1.5e+02,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (1.8e+02,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,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::swap_row. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_33/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m.swap_row (2, 4);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (10.,0.0)
+DEAL::2 1 (15.,0.0)
+DEAL::2 2 (20.,0.0)
+DEAL::2 3 (25.,0.0)
+DEAL::2 4 (30.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (6.0,0.0)
+DEAL::4 1 (9.0,0.0)
+DEAL::4 2 (12.,0.0)
+DEAL::4 3 (15.,0.0)
+DEAL::4 4 (18.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (10.,0.0)
+DEAL::2 1 (15.,0.0)
+DEAL::2 2 (20.,0.0)
+DEAL::2 3 (25.,0.0)
+DEAL::2 4 (30.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (6.0,0.0)
+DEAL::4 1 (9.0,0.0)
+DEAL::4 2 (12.,0.0)
+DEAL::4 3 (15.,0.0)
+DEAL::4 4 (18.,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::swap_col. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_34/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m.swap_col (2, 4);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (6.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (4.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (12.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (8.0,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (18.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (12.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (24.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (16.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (30.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (20.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (6.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (4.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (12.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (8.0,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (18.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (12.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (24.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (16.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (30.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (20.,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::add_diag. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_35/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m, n;
+ make_matrix (m);
+ make_matrix (n);
+ m.add_diag (3.141, m);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (8.3,0.0)
+DEAL::0 1 (12.,0.0)
+DEAL::0 2 (17.,0.0)
+DEAL::0 3 (21.,0.0)
+DEAL::0 4 (25.,0.0)
+DEAL::1 0 (17.,0.0)
+DEAL::1 1 (25.,0.0)
+DEAL::1 2 (33.,0.0)
+DEAL::1 3 (41.,0.0)
+DEAL::1 4 (50.,0.0)
+DEAL::2 0 (25.,0.0)
+DEAL::2 1 (37.,0.0)
+DEAL::2 2 (50.,0.0)
+DEAL::2 3 (62.,0.0)
+DEAL::2 4 (75.,0.0)
+DEAL::3 0 (33.,0.0)
+DEAL::3 1 (50.,0.0)
+DEAL::3 2 (66.,0.0)
+DEAL::3 3 (83.,0.0)
+DEAL::3 4 (99.,0.0)
+DEAL::4 0 (41.,0.0)
+DEAL::4 1 (62.,0.0)
+DEAL::4 2 (83.,0.0)
+DEAL::4 3 (1.0e+02,0.0)
+DEAL::4 4 (1.2e+02,0.0)
+DEAL::0 0 (8.3,0.0)
+DEAL::0 1 (12.,0.0)
+DEAL::0 2 (17.,0.0)
+DEAL::0 3 (21.,0.0)
+DEAL::0 4 (25.,0.0)
+DEAL::1 0 (17.,0.0)
+DEAL::1 1 (25.,0.0)
+DEAL::1 2 (33.,0.0)
+DEAL::1 3 (41.,0.0)
+DEAL::1 4 (50.,0.0)
+DEAL::2 0 (25.,0.0)
+DEAL::2 1 (37.,0.0)
+DEAL::2 2 (50.,0.0)
+DEAL::2 3 (62.,0.0)
+DEAL::2 4 (75.,0.0)
+DEAL::3 0 (33.,0.0)
+DEAL::3 1 (50.,0.0)
+DEAL::3 2 (66.,0.0)
+DEAL::3 3 (83.,0.0)
+DEAL::3 4 (99.,0.0)
+DEAL::4 0 (41.,0.0)
+DEAL::4 1 (62.,0.0)
+DEAL::4 2 (83.,0.0)
+DEAL::4 3 (1.0e+02,0.0)
+DEAL::4 4 (1.2e+02,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::diag_add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_36/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m.diagadd (3.141);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (5.1,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (9.1,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (15.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (23.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (33.,0.0)
+DEAL::0 0 (5.1,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (9.1,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (15.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (23.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (33.,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::equ. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_37/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n;
+ make_matrix (m);
+ make_matrix (n);
+ m.equ (3.14159, n);
+
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (6.3,0.0)
+DEAL::0 1 (9.4,0.0)
+DEAL::0 2 (13.,0.0)
+DEAL::0 3 (16.,0.0)
+DEAL::0 4 (19.,0.0)
+DEAL::1 0 (13.,0.0)
+DEAL::1 1 (19.,0.0)
+DEAL::1 2 (25.,0.0)
+DEAL::1 3 (31.,0.0)
+DEAL::1 4 (38.,0.0)
+DEAL::2 0 (19.,0.0)
+DEAL::2 1 (28.,0.0)
+DEAL::2 2 (38.,0.0)
+DEAL::2 3 (47.,0.0)
+DEAL::2 4 (57.,0.0)
+DEAL::3 0 (25.,0.0)
+DEAL::3 1 (38.,0.0)
+DEAL::3 2 (50.,0.0)
+DEAL::3 3 (63.,0.0)
+DEAL::3 4 (75.,0.0)
+DEAL::4 0 (31.,0.0)
+DEAL::4 1 (47.,0.0)
+DEAL::4 2 (63.,0.0)
+DEAL::4 3 (79.,0.0)
+DEAL::4 4 (94.,0.0)
+DEAL::0 0 (6.3,0.0)
+DEAL::0 1 (9.4,0.0)
+DEAL::0 2 (13.,0.0)
+DEAL::0 3 (16.,0.0)
+DEAL::0 4 (19.,0.0)
+DEAL::1 0 (13.,0.0)
+DEAL::1 1 (19.,0.0)
+DEAL::1 2 (25.,0.0)
+DEAL::1 3 (31.,0.0)
+DEAL::1 4 (38.,0.0)
+DEAL::2 0 (19.,0.0)
+DEAL::2 1 (28.,0.0)
+DEAL::2 2 (38.,0.0)
+DEAL::2 3 (47.,0.0)
+DEAL::2 4 (57.,0.0)
+DEAL::3 0 (25.,0.0)
+DEAL::3 1 (38.,0.0)
+DEAL::3 2 (50.,0.0)
+DEAL::3 3 (63.,0.0)
+DEAL::3 4 (75.,0.0)
+DEAL::4 0 (31.,0.0)
+DEAL::4 1 (47.,0.0)
+DEAL::4 2 (63.,0.0)
+DEAL::4 3 (79.,0.0)
+DEAL::4 4 (94.,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::equ(4). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_38/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n,o;
+ make_matrix (m);
+ make_matrix (n);
+ make_matrix (o);
+ m.equ (3.1415, n, 2.718, o);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (12.,0.0)
+DEAL::0 1 (18.,0.0)
+DEAL::0 2 (23.,0.0)
+DEAL::0 3 (29.,0.0)
+DEAL::0 4 (35.,0.0)
+DEAL::1 0 (23.,0.0)
+DEAL::1 1 (35.,0.0)
+DEAL::1 2 (47.,0.0)
+DEAL::1 3 (59.,0.0)
+DEAL::1 4 (70.,0.0)
+DEAL::2 0 (35.,0.0)
+DEAL::2 1 (53.,0.0)
+DEAL::2 2 (70.,0.0)
+DEAL::2 3 (88.,0.0)
+DEAL::2 4 (1.1e+02,0.0)
+DEAL::3 0 (47.,0.0)
+DEAL::3 1 (70.,0.0)
+DEAL::3 2 (94.,0.0)
+DEAL::3 3 (1.2e+02,0.0)
+DEAL::3 4 (1.4e+02,0.0)
+DEAL::4 0 (59.,0.0)
+DEAL::4 1 (88.,0.0)
+DEAL::4 2 (1.2e+02,0.0)
+DEAL::4 3 (1.5e+02,0.0)
+DEAL::4 4 (1.8e+02,0.0)
+DEAL::0 0 (12.,0.0)
+DEAL::0 1 (18.,0.0)
+DEAL::0 2 (23.,0.0)
+DEAL::0 3 (29.,0.0)
+DEAL::0 4 (35.,0.0)
+DEAL::1 0 (23.,0.0)
+DEAL::1 1 (35.,0.0)
+DEAL::1 2 (47.,0.0)
+DEAL::1 3 (59.,0.0)
+DEAL::1 4 (70.,0.0)
+DEAL::2 0 (35.,0.0)
+DEAL::2 1 (53.,0.0)
+DEAL::2 2 (70.,0.0)
+DEAL::2 3 (88.,0.0)
+DEAL::2 4 (1.1e+02,0.0)
+DEAL::3 0 (47.,0.0)
+DEAL::3 1 (70.,0.0)
+DEAL::3 2 (94.,0.0)
+DEAL::3 3 (1.2e+02,0.0)
+DEAL::3 4 (1.4e+02,0.0)
+DEAL::4 0 (59.,0.0)
+DEAL::4 1 (88.,0.0)
+DEAL::4 2 (1.2e+02,0.0)
+DEAL::4 3 (1.5e+02,0.0)
+DEAL::4 4 (1.8e+02,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::equ (6). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_39/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n,o,p;
+ make_matrix (m);
+ make_matrix (n);
+ make_matrix (o);
+ make_matrix (p);
+ m.equ (3.1415, n, 2.718, o, 1.414, p);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (15.,0.0)
+DEAL::0 1 (22.,0.0)
+DEAL::0 2 (29.,0.0)
+DEAL::0 3 (36.,0.0)
+DEAL::0 4 (44.,0.0)
+DEAL::1 0 (29.,0.0)
+DEAL::1 1 (44.,0.0)
+DEAL::1 2 (58.,0.0)
+DEAL::1 3 (73.,0.0)
+DEAL::1 4 (87.,0.0)
+DEAL::2 0 (44.,0.0)
+DEAL::2 1 (65.,0.0)
+DEAL::2 2 (87.,0.0)
+DEAL::2 3 (1.1e+02,0.0)
+DEAL::2 4 (1.3e+02,0.0)
+DEAL::3 0 (58.,0.0)
+DEAL::3 1 (87.,0.0)
+DEAL::3 2 (1.2e+02,0.0)
+DEAL::3 3 (1.5e+02,0.0)
+DEAL::3 4 (1.7e+02,0.0)
+DEAL::4 0 (73.,0.0)
+DEAL::4 1 (1.1e+02,0.0)
+DEAL::4 2 (1.5e+02,0.0)
+DEAL::4 3 (1.8e+02,0.0)
+DEAL::4 4 (2.2e+02,0.0)
+DEAL::0 0 (15.,0.0)
+DEAL::0 1 (22.,0.0)
+DEAL::0 2 (29.,0.0)
+DEAL::0 3 (36.,0.0)
+DEAL::0 4 (44.,0.0)
+DEAL::1 0 (29.,0.0)
+DEAL::1 1 (44.,0.0)
+DEAL::1 2 (58.,0.0)
+DEAL::1 3 (73.,0.0)
+DEAL::1 4 (87.,0.0)
+DEAL::2 0 (44.,0.0)
+DEAL::2 1 (65.,0.0)
+DEAL::2 2 (87.,0.0)
+DEAL::2 3 (1.1e+02,0.0)
+DEAL::2 4 (1.3e+02,0.0)
+DEAL::3 0 (58.,0.0)
+DEAL::3 1 (87.,0.0)
+DEAL::3 2 (1.2e+02,0.0)
+DEAL::3 3 (1.5e+02,0.0)
+DEAL::3 4 (1.7e+02,0.0)
+DEAL::4 0 (73.,0.0)
+DEAL::4 1 (1.1e+02,0.0)
+DEAL::4 2 (1.5e+02,0.0)
+DEAL::4 3 (1.8e+02,0.0)
+DEAL::4 4 (2.2e+02,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::symmetrize. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_40/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m.symmetrize ();
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.5,0.0)
+DEAL::0 2 (5.0,0.0)
+DEAL::0 3 (6.5,0.0)
+DEAL::0 4 (8.0,0.0)
+DEAL::1 0 (3.5,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.5,0.0)
+DEAL::1 3 (11.,0.0)
+DEAL::1 4 (14.,0.0)
+DEAL::2 0 (5.0,0.0)
+DEAL::2 1 (8.5,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (16.,0.0)
+DEAL::2 4 (19.,0.0)
+DEAL::3 0 (6.5,0.0)
+DEAL::3 1 (11.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (8.0,0.0)
+DEAL::4 1 (14.,0.0)
+DEAL::4 2 (19.,0.0)
+DEAL::4 3 (24.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.5,0.0)
+DEAL::0 2 (5.0,0.0)
+DEAL::0 3 (6.5,0.0)
+DEAL::0 4 (8.0,0.0)
+DEAL::1 0 (3.5,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.5,0.0)
+DEAL::1 3 (11.,0.0)
+DEAL::1 4 (14.,0.0)
+DEAL::2 0 (5.0,0.0)
+DEAL::2 1 (8.5,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (16.,0.0)
+DEAL::2 4 (19.,0.0)
+DEAL::3 0 (6.5,0.0)
+DEAL::3 1 (11.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (8.0,0.0)
+DEAL::4 1 (14.,0.0)
+DEAL::4 2 (19.,0.0)
+DEAL::4 3 (24.,0.0)
+DEAL::4 4 (30.,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::gauss_jordan. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_41/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ m(0,0) = m(1,1) = m(2,2) = m(3,3) = m(4,4) = 50;
+ m.gauss_jordan ();
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (0.021,0.0)
+DEAL::0 1 (0.00039,-0.0)
+DEAL::0 2 (0.00060,0.0)
+DEAL::0 3 (0.00095,0.0)
+DEAL::0 4 (0.0017,0.0)
+DEAL::1 0 (0.00052,0.0)
+DEAL::1 1 (0.022,0.0)
+DEAL::1 2 (0.0013,-0.0)
+DEAL::1 3 (0.0021,0.0)
+DEAL::1 4 (0.0037,0.0)
+DEAL::2 0 (0.00090,0.0)
+DEAL::2 1 (0.0015,0.0)
+DEAL::2 2 (0.024,0.0)
+DEAL::2 3 (0.0036,-0.0)
+DEAL::2 4 (0.0065,0.0)
+DEAL::3 0 (0.0015,0.0)
+DEAL::3 1 (0.0025,0.0)
+DEAL::3 2 (0.0038,0.0)
+DEAL::3 3 (0.027,0.0)
+DEAL::3 4 (-0.011,-0.0)
+DEAL::4 0 (0.0028,-0.0)
+DEAL::4 1 (0.0047,-0.0)
+DEAL::4 2 (0.0072,-0.0)
+DEAL::4 3 (-0.011,-0.0)
+DEAL::4 4 (0.030,0.0)
+DEAL::0 0 (0.021,0.0)
+DEAL::0 1 (0.00039,-0.0)
+DEAL::0 2 (0.00060,0.0)
+DEAL::0 3 (0.00095,0.0)
+DEAL::0 4 (0.0017,0.0)
+DEAL::1 0 (0.00052,0.0)
+DEAL::1 1 (0.022,0.0)
+DEAL::1 2 (0.0013,-0.0)
+DEAL::1 3 (0.0021,0.0)
+DEAL::1 4 (0.0037,0.0)
+DEAL::2 0 (0.00090,0.0)
+DEAL::2 1 (0.0015,0.0)
+DEAL::2 2 (0.024,0.0)
+DEAL::2 3 (0.0036,-0.0)
+DEAL::2 4 (0.0065,0.0)
+DEAL::3 0 (0.0015,0.0)
+DEAL::3 1 (0.0025,0.0)
+DEAL::3 2 (0.0038,0.0)
+DEAL::3 3 (0.027,0.0)
+DEAL::3 4 (-0.011,-0.0)
+DEAL::4 0 (0.0028,-0.0)
+DEAL::4 1 (0.0047,-0.0)
+DEAL::4 2 (0.0072,-0.0)
+DEAL::4 3 (-0.011,-0.0)
+DEAL::4 4 (0.030,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::invert. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_42/output";
+
+
+template <typename number>
+void
+check ()
+{
+ for (unsigned int n=1; n<=5; ++n)
+ {
+ const std::complex<number> array[] = { 50, 2, 3, 4, 5,
+ 6, 50, 8, 9, 0,
+ 1, 2, 50, 4, 5,
+ 6, 7, 8, 50, 0,
+ 1, 2, 3, 4, 50 };
+
+ FullMatrix<std::complex<number> > m (n,n,array), p(n,n);
+ p.invert (m);
+ print_matrix (p);
+ }
+}
+
--- /dev/null
+
+DEAL::0 0 (0.020,0.0)
+DEAL::0 0 (0.021,0.0)
+DEAL::0 1 (-0.010,-0.0)
+DEAL::1 0 (-0.015,-0.0)
+DEAL::1 1 (0.26,0.0)
+DEAL::0 0 (0.011,0.0)
+DEAL::0 1 (-0.022,0.0)
+DEAL::0 2 (0.011,0.0)
+DEAL::1 0 (-0.96,0.0)
+DEAL::1 1 (-1.1,-0.0)
+DEAL::1 2 (1.0,0.0)
+DEAL::2 0 (0.79,0.0)
+DEAL::2 1 (1.1,0.0)
+DEAL::2 2 (-0.88,-0.0)
+DEAL::0 0 (0.049,0.0)
+DEAL::0 1 (0.0,-0.0)
+DEAL::0 2 (-0.024,0.0)
+DEAL::0 3 (-0.024,-0.0)
+DEAL::1 0 (-0.32,-0.0)
+DEAL::1 1 (-0.012,0.0)
+DEAL::1 2 (-0.60,-0.0)
+DEAL::1 3 (0.43,0.0)
+DEAL::2 0 (0.074,0.0)
+DEAL::2 1 (0.023,0.0)
+DEAL::2 2 (-0.026,0.0)
+DEAL::2 3 (-0.072,-0.0)
+DEAL::3 0 (-0.26,-0.0)
+DEAL::3 1 (-0.012,0.0)
+DEAL::3 2 (0.62,0.0)
+DEAL::3 3 (0.15,0.0)
+DEAL::0 0 (0.020,0.0)
+DEAL::0 1 (0.00052,0.0)
+DEAL::0 2 (0.00081,0.0)
+DEAL::0 3 (0.0013,0.0)
+DEAL::0 4 (0.0019,-0.0)
+DEAL::1 0 (0.0020,0.0)
+DEAL::1 1 (0.021,0.0)
+DEAL::1 2 (0.0027,-0.0)
+DEAL::1 3 (0.0034,0.0)
+DEAL::1 4 (0.00047,0.0)
+DEAL::2 0 (0.00014,0.0)
+DEAL::2 1 (0.00054,0.0)
+DEAL::2 2 (0.020,0.0)
+DEAL::2 3 (0.0014,0.0)
+DEAL::2 4 (0.0020,-0.0)
+DEAL::3 0 (0.0021,-0.0)
+DEAL::3 1 (0.0027,-0.0)
+DEAL::3 2 (0.0028,-0.0)
+DEAL::3 3 (0.021,0.0)
+DEAL::3 4 (0.00049,0.0)
+DEAL::4 0 (0.00015,-0.0)
+DEAL::4 1 (0.00056,-0.0)
+DEAL::4 2 (0.00088,-0.0)
+DEAL::4 3 (0.0014,-0.0)
+DEAL::4 4 (0.020,0.0)
+DEAL::0 0 (0.020,0.0)
+DEAL::0 0 (0.021,0.0)
+DEAL::0 1 (-0.010,-0.0)
+DEAL::1 0 (-0.015,-0.0)
+DEAL::1 1 (0.26,0.0)
+DEAL::0 0 (0.011,0.0)
+DEAL::0 1 (-0.022,0.0)
+DEAL::0 2 (0.011,0.0)
+DEAL::1 0 (-0.96,0.0)
+DEAL::1 1 (-1.1,-0.0)
+DEAL::1 2 (1.0,0.0)
+DEAL::2 0 (0.79,0.0)
+DEAL::2 1 (1.1,0.0)
+DEAL::2 2 (-0.88,-0.0)
+DEAL::0 0 (0.049,0.0)
+DEAL::0 1 (0.0,-0.0)
+DEAL::0 2 (-0.024,0.0)
+DEAL::0 3 (-0.024,-0.0)
+DEAL::1 0 (-0.32,-0.0)
+DEAL::1 1 (-0.012,0.0)
+DEAL::1 2 (-0.60,-0.0)
+DEAL::1 3 (0.43,0.0)
+DEAL::2 0 (0.074,0.0)
+DEAL::2 1 (0.023,0.0)
+DEAL::2 2 (-0.026,0.0)
+DEAL::2 3 (-0.072,-0.0)
+DEAL::3 0 (-0.26,-0.0)
+DEAL::3 1 (-0.012,0.0)
+DEAL::3 2 (0.62,0.0)
+DEAL::3 3 (0.15,0.0)
+DEAL::0 0 (0.020,0.0)
+DEAL::0 1 (0.00052,0.0)
+DEAL::0 2 (0.00081,0.0)
+DEAL::0 3 (0.0013,0.0)
+DEAL::0 4 (0.0019,-0.0)
+DEAL::1 0 (0.0020,0.0)
+DEAL::1 1 (0.021,0.0)
+DEAL::1 2 (0.0027,-0.0)
+DEAL::1 3 (0.0034,0.0)
+DEAL::1 4 (0.00047,0.0)
+DEAL::2 0 (0.00014,0.0)
+DEAL::2 1 (0.00054,0.0)
+DEAL::2 2 (0.020,0.0)
+DEAL::2 3 (0.0014,0.0)
+DEAL::2 4 (0.0020,-0.0)
+DEAL::3 0 (0.0021,-0.0)
+DEAL::3 1 (0.0027,-0.0)
+DEAL::3 2 (0.0028,-0.0)
+DEAL::3 3 (0.021,0.0)
+DEAL::3 4 (0.00049,0.0)
+DEAL::4 0 (0.00015,-0.0)
+DEAL::4 1 (0.00056,-0.0)
+DEAL::4 2 (0.00088,-0.0)
+DEAL::4 3 (0.0014,-0.0)
+DEAL::4 4 (0.020,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::mmult. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_43/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<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,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (1.4e+02,0.0)
+DEAL::0 1 (2.1e+02,0.0)
+DEAL::0 2 (2.8e+02,0.0)
+DEAL::0 3 (3.6e+02,0.0)
+DEAL::0 4 (4.3e+02,0.0)
+DEAL::1 0 (2.8e+02,0.0)
+DEAL::1 1 (4.3e+02,0.0)
+DEAL::1 2 (5.7e+02,0.0)
+DEAL::1 3 (7.1e+02,0.0)
+DEAL::1 4 (8.5e+02,0.0)
+DEAL::2 0 (4.3e+02,0.0)
+DEAL::2 1 (6.4e+02,0.0)
+DEAL::2 2 (8.5e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.7e+02,0.0)
+DEAL::3 1 (8.5e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (7.1e+02,0.0)
+DEAL::4 1 (1.1e+03,0.0)
+DEAL::4 2 (1.4e+03,0.0)
+DEAL::4 3 (1.8e+03,0.0)
+DEAL::4 4 (2.1e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (1.4e+02,0.0)
+DEAL::0 1 (2.1e+02,0.0)
+DEAL::0 2 (2.8e+02,0.0)
+DEAL::0 3 (3.5e+02,0.0)
+DEAL::0 4 (4.2e+02,0.0)
+DEAL::1 0 (2.8e+02,0.0)
+DEAL::1 1 (4.2e+02,0.0)
+DEAL::1 2 (5.6e+02,0.0)
+DEAL::1 3 (7.0e+02,0.0)
+DEAL::1 4 (8.4e+02,0.0)
+DEAL::2 0 (4.2e+02,0.0)
+DEAL::2 1 (6.3e+02,0.0)
+DEAL::2 2 (8.4e+02,0.0)
+DEAL::2 3 (1.0e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.6e+02,0.0)
+DEAL::3 1 (8.4e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (7.0e+02,0.0)
+DEAL::4 1 (1.0e+03,0.0)
+DEAL::4 2 (1.4e+03,0.0)
+DEAL::4 3 (1.8e+03,0.0)
+DEAL::4 4 (2.1e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (1.4e+02,0.0)
+DEAL::0 1 (2.1e+02,0.0)
+DEAL::0 2 (2.8e+02,0.0)
+DEAL::0 3 (3.6e+02,0.0)
+DEAL::0 4 (4.3e+02,0.0)
+DEAL::1 0 (2.8e+02,0.0)
+DEAL::1 1 (4.3e+02,0.0)
+DEAL::1 2 (5.7e+02,0.0)
+DEAL::1 3 (7.1e+02,0.0)
+DEAL::1 4 (8.5e+02,0.0)
+DEAL::2 0 (4.3e+02,0.0)
+DEAL::2 1 (6.4e+02,0.0)
+DEAL::2 2 (8.5e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.7e+02,0.0)
+DEAL::3 1 (8.5e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (7.1e+02,0.0)
+DEAL::4 1 (1.1e+03,0.0)
+DEAL::4 2 (1.4e+03,0.0)
+DEAL::4 3 (1.8e+03,0.0)
+DEAL::4 4 (2.1e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (1.4e+02,0.0)
+DEAL::0 1 (2.1e+02,0.0)
+DEAL::0 2 (2.8e+02,0.0)
+DEAL::0 3 (3.5e+02,0.0)
+DEAL::0 4 (4.2e+02,0.0)
+DEAL::1 0 (2.8e+02,0.0)
+DEAL::1 1 (4.2e+02,0.0)
+DEAL::1 2 (5.6e+02,0.0)
+DEAL::1 3 (7.0e+02,0.0)
+DEAL::1 4 (8.4e+02,0.0)
+DEAL::2 0 (4.2e+02,0.0)
+DEAL::2 1 (6.3e+02,0.0)
+DEAL::2 2 (8.4e+02,0.0)
+DEAL::2 3 (1.0e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.6e+02,0.0)
+DEAL::3 1 (8.4e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (7.0e+02,0.0)
+DEAL::4 1 (1.0e+03,0.0)
+DEAL::4 2 (1.4e+03,0.0)
+DEAL::4 3 (1.8e+03,0.0)
+DEAL::4 4 (2.1e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,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::Tmmult. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_44/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<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,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.2e+02,0.0)
+DEAL::0 1 (3.3e+02,0.0)
+DEAL::0 2 (4.4e+02,0.0)
+DEAL::0 3 (5.6e+02,0.0)
+DEAL::0 4 (6.7e+02,0.0)
+DEAL::1 0 (3.3e+02,0.0)
+DEAL::1 1 (5.0e+02,0.0)
+DEAL::1 2 (6.7e+02,0.0)
+DEAL::1 3 (8.4e+02,0.0)
+DEAL::1 4 (1.0e+03,0.0)
+DEAL::2 0 (4.5e+02,0.0)
+DEAL::2 1 (6.7e+02,0.0)
+DEAL::2 2 (8.9e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.6e+02,0.0)
+DEAL::3 1 (8.4e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (6.7e+02,0.0)
+DEAL::4 1 (1.0e+03,0.0)
+DEAL::4 2 (1.3e+03,0.0)
+DEAL::4 3 (1.7e+03,0.0)
+DEAL::4 4 (2.0e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.2e+02,0.0)
+DEAL::0 1 (3.3e+02,0.0)
+DEAL::0 2 (4.4e+02,0.0)
+DEAL::0 3 (5.5e+02,0.0)
+DEAL::0 4 (6.6e+02,0.0)
+DEAL::1 0 (3.3e+02,0.0)
+DEAL::1 1 (5.0e+02,0.0)
+DEAL::1 2 (6.6e+02,0.0)
+DEAL::1 3 (8.2e+02,0.0)
+DEAL::1 4 (9.9e+02,0.0)
+DEAL::2 0 (4.4e+02,0.0)
+DEAL::2 1 (6.6e+02,0.0)
+DEAL::2 2 (8.8e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.5e+02,0.0)
+DEAL::3 1 (8.2e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.6e+03,0.0)
+DEAL::4 0 (6.6e+02,0.0)
+DEAL::4 1 (9.9e+02,0.0)
+DEAL::4 2 (1.3e+03,0.0)
+DEAL::4 3 (1.6e+03,0.0)
+DEAL::4 4 (2.0e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.2e+02,0.0)
+DEAL::0 1 (3.3e+02,0.0)
+DEAL::0 2 (4.4e+02,0.0)
+DEAL::0 3 (5.6e+02,0.0)
+DEAL::0 4 (6.7e+02,0.0)
+DEAL::1 0 (3.3e+02,0.0)
+DEAL::1 1 (5.0e+02,0.0)
+DEAL::1 2 (6.7e+02,0.0)
+DEAL::1 3 (8.4e+02,0.0)
+DEAL::1 4 (1.0e+03,0.0)
+DEAL::2 0 (4.5e+02,0.0)
+DEAL::2 1 (6.7e+02,0.0)
+DEAL::2 2 (8.9e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.6e+02,0.0)
+DEAL::3 1 (8.4e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (6.7e+02,0.0)
+DEAL::4 1 (1.0e+03,0.0)
+DEAL::4 2 (1.3e+03,0.0)
+DEAL::4 3 (1.7e+03,0.0)
+DEAL::4 4 (2.0e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.2e+02,0.0)
+DEAL::0 1 (3.3e+02,0.0)
+DEAL::0 2 (4.4e+02,0.0)
+DEAL::0 3 (5.5e+02,0.0)
+DEAL::0 4 (6.6e+02,0.0)
+DEAL::1 0 (3.3e+02,0.0)
+DEAL::1 1 (5.0e+02,0.0)
+DEAL::1 2 (6.6e+02,0.0)
+DEAL::1 3 (8.2e+02,0.0)
+DEAL::1 4 (9.9e+02,0.0)
+DEAL::2 0 (4.4e+02,0.0)
+DEAL::2 1 (6.6e+02,0.0)
+DEAL::2 2 (8.8e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.5e+02,0.0)
+DEAL::3 1 (8.2e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.6e+03,0.0)
+DEAL::4 0 (6.6e+02,0.0)
+DEAL::4 1 (9.9e+02,0.0)
+DEAL::4 2 (1.3e+03,0.0)
+DEAL::4 3 (1.6e+03,0.0)
+DEAL::4 4 (2.0e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,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::vmult. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_45/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ Vector<std::complex<number> > 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.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.6e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (70.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.5e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (71.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.6e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (70.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.5e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,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::vmult_add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_46/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ Vector<std::complex<number> > v, w;
+ make_vector (v);
+ make_vector (w);
+
+ m.vmult_add (v, w);
+ print_vector (v);
+ print_vector (w);
+}
+
--- /dev/null
+
+DEAL::0 (71.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.6e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (71.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.6e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,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::Tvmult. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_47/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ Vector<std::complex<number> > 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,0.0)
+DEAL::1 (1.7e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.1e+02,0.0)
+DEAL::1 (1.6e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.3e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.1e+02,0.0)
+DEAL::1 (1.7e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.1e+02,0.0)
+DEAL::1 (1.6e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.3e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,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::Tvmult_add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_48/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ Vector<std::complex<number> > 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,0.0)
+DEAL::1 (1.7e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.1e+02,0.0)
+DEAL::1 (1.7e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,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::precondition_Jacobi. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_49/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ Vector<std::complex<number> > 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,0.0)
+DEAL::1 (1.0,0.0)
+DEAL::2 (0.79,0.0)
+DEAL::3 (0.63,0.0)
+DEAL::4 (0.52,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.6,0.0)
+DEAL::1 (1.0,0.0)
+DEAL::2 (0.79,0.0)
+DEAL::3 (0.63,0.0)
+DEAL::4 (0.52,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,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::residual. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_50/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ Vector<std::complex<number> > 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.,0.0)
+DEAL::1 (-1.4e+02,0.0)
+DEAL::2 (-2.1e+02,0.0)
+DEAL::3 (-2.8e+02,0.0)
+DEAL::4 (-3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (-69.,0.0)
+DEAL::1 (-1.4e+02,0.0)
+DEAL::2 (-2.1e+02,0.0)
+DEAL::3 (-2.8e+02,0.0)
+DEAL::4 (-3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,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::forward. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_51/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ Vector<std::complex<number> > x, y;
+ make_vector (x);
+ make_vector (y);
+
+ m.forward (x, y);
+ print_vector (x);
+ print_vector (y);
+}
+
--- /dev/null
+
+DEAL::0 (0.50,0.0)
+DEAL::1 (0.0,0.0)
+DEAL::2 (0.0,0.0)
+DEAL::3 (0.0,0.0)
+DEAL::4 (0.0,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (0.50,0.0)
+DEAL::1 (0.0,0.0)
+DEAL::2 (0.0,0.0)
+DEAL::3 (0.0,0.0)
+DEAL::4 (0.0,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,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::backward. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_52/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_matrix (m);
+ Vector<std::complex<number> > x, y;
+ make_vector (x);
+ make_vector (y);
+
+ m.backward (x, y);
+ print_vector (x);
+ print_vector (y);
+}
+
--- /dev/null
+
+DEAL::0 (0.0,0.0)
+DEAL::1 (0.0,0.0)
+DEAL::2 (0.0,0.0)
+DEAL::3 (0.0,0.0)
+DEAL::4 (0.17,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (0.0,0.0)
+DEAL::1 (0.0,0.0)
+DEAL::2 (0.0,0.0)
+DEAL::3 (0.0,0.0)
+DEAL::4 (0.17,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
{
FullMatrix<number> m;
make_matrix (m);
- Vector<double> v;
+ Vector<number> v;
make_vector (v);
deallog << m.matrix_norm_square (v) << std::endl;
{
FullMatrix<number> m;
make_matrix (m);
- Vector<double> v,w;
+ Vector<number> v,w;
make_vector (v);
make_vector (w);
for (unsigned int i=0; i<w.size(); ++i)
{
FullMatrix<number> m;
make_matrix (m);
- Vector<double> v, w;
+ Vector<number> v, w;
make_vector (v);
make_vector (w);
{
FullMatrix<number> m;
make_matrix (m);
- Vector<double> v, w;
+ Vector<number> v, w;
make_vector (v);
make_vector (w);
{
FullMatrix<number> m;
make_matrix (m);
- Vector<double> v, w;
+ Vector<number> v, w;
make_vector (v);
make_vector (w);
{
FullMatrix<number> m;
make_matrix (m);
- Vector<double> v, w;
+ Vector<number> v, w;
make_vector (v);
make_vector (w);
{
FullMatrix<number> m;
make_matrix (m);
- Vector<double> v, w;
+ Vector<number> v, w;
make_vector (v);
make_vector (w);
{
FullMatrix<number> m;
make_matrix (m);
- Vector<double> v, w, x;
+ Vector<number> v, w, x;
make_vector (v);
make_vector (w);
make_vector (x);
{
FullMatrix<number> m;
make_matrix (m);
- Vector<double> x, y;
+ Vector<number> x, y;
make_vector (x);
make_vector (y);
{
FullMatrix<number> m;
make_matrix (m);
- Vector<double> x, y;
+ Vector<number> x, y;
make_vector (x);
make_vector (y);
#include <vector>
-void test (Vector<double> &v,
- Vector<double> &w)
+void test (Vector<number> &v,
+ Vector<number> &w)
{
FullMatrix<double> m(v.size(), v.size());
for (unsigned int i=0; i<m.m(); ++i)
try
{
- Vector<double> v (100);
- Vector<double> w (100);
+ Vector<number> v (100);
+ Vector<number> w (100);
test (v,w);
}
catch (std::exception &exc)
#include <vector>
-void test (Vector<double> &v,
- Vector<double> &w)
+void test (Vector<number> &v,
+ Vector<number> &w)
{
FullMatrix<double> m(v.size(), v.size());
for (unsigned int i=0; i<m.m(); ++i)
try
{
- Vector<double> v (100);
- Vector<double> w (100);
+ Vector<number> v (100);
+ Vector<number> w (100);
test (v,w);
}
catch (std::exception &exc)
#include <vector>
-void test (Vector<double> &v,
- Vector<double> &w)
+void test (Vector<number> &v,
+ Vector<number> &w)
{
FullMatrix<double> m(v.size(), v.size());
for (unsigned int i=0; i<m.m(); ++i)
try
{
- Vector<double> v (100);
- Vector<double> w (100);
+ Vector<number> v (100);
+ Vector<number> w (100);
test (v,w);
}
catch (std::exception &exc)
#include <vector>
-void test (Vector<double> &v)
+void test (Vector<number> &v)
{
FullMatrix<double> m(v.size(), v.size());
for (unsigned int i=0; i<m.m(); ++i)
try
{
- Vector<double> v (100);
+ Vector<number> v (100);
test (v);
}
catch (std::exception &exc)
#include <vector>
-void test (Vector<double> &v,
- Vector<double> &w,
- Vector<double> &x)
+void test (Vector<number> &v,
+ Vector<number> &w,
+ Vector<number> &x)
{
FullMatrix<double> m(v.size(), v.size());
for (unsigned int i=0; i<m.m(); ++i)
try
{
- Vector<double> v (100);
- Vector<double> w (100);
- Vector<double> x (100);
+ Vector<number> v (100);
+ Vector<number> w (100);
+ Vector<number> x (100);
test (v,w,x);
}
catch (std::exception &exc)