point_inside_? \
full_matrix_* \
complex_real_full_matrix_* \
+ complex_complex_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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_01/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_02/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_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,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_03/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_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,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_06/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+
+ FullMatrix<std::complex<number> > n;
+ n = m;
+
+ print_matrix (n);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_07/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_08/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+
+ FullMatrix<std::complex<number> > n;
+ n.copy_from (m);
+
+ print_matrix (n);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_10/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_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,12.)
+DEAL::3 5 (3.0,15.)
+DEAL::3 6 (4.0,18.)
+DEAL::3 7 (5.0,21.)
+DEAL::3 8 (6.0,24.)
+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,16.)
+DEAL::4 5 (6.0,20.)
+DEAL::4 6 (8.0,24.)
+DEAL::4 7 (10.,28.)
+DEAL::4 8 (12.,32.)
+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,20.)
+DEAL::5 5 (9.0,25.)
+DEAL::5 6 (12.,30.)
+DEAL::5 7 (15.,35.)
+DEAL::5 8 (18.,40.)
+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,24.)
+DEAL::6 5 (12.,30.)
+DEAL::6 6 (16.,36.)
+DEAL::6 7 (20.,42.)
+DEAL::6 8 (24.,48.)
+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.,28.)
+DEAL::7 5 (15.,35.)
+DEAL::7 6 (20.,42.)
+DEAL::7 7 (25.,49.)
+DEAL::7 8 (30.,56.)
+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,12.)
+DEAL::3 5 (3.0,15.)
+DEAL::3 6 (4.0,18.)
+DEAL::3 7 (5.0,21.)
+DEAL::3 8 (6.0,24.)
+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,16.)
+DEAL::4 5 (6.0,20.)
+DEAL::4 6 (8.0,24.)
+DEAL::4 7 (10.,28.)
+DEAL::4 8 (12.,32.)
+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,20.)
+DEAL::5 5 (9.0,25.)
+DEAL::5 6 (12.,30.)
+DEAL::5 7 (15.,35.)
+DEAL::5 8 (18.,40.)
+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,24.)
+DEAL::6 5 (12.,30.)
+DEAL::6 6 (16.,36.)
+DEAL::6 7 (20.,42.)
+DEAL::6 8 (24.,48.)
+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.,28.)
+DEAL::7 5 (15.,35.)
+DEAL::7 6 (20.,42.)
+DEAL::7 7 (25.,49.)
+DEAL::7 8 (30.,56.)
+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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_11/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m, n;
+ make_complex_matrix (m);
+ make_complex_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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_12/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m(3,3);
+ deallog << m.all_zero() << std::endl;
+
+ make_complex_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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_13/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ Vector<std::complex<number> > v;
+ make_complex_vector (v);
+
+ deallog << m.matrix_norm_square (v) << std::endl;
+}
+
--- /dev/null
+
+DEAL::(-3.3e+04,-1.1e+03)
+DEAL::(-3.3e+04,-1.1e+03)
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::matrix_scalar_product. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_14/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ Vector<std::complex<number> > v,w;
+ make_complex_vector (v);
+ make_complex_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::(-3.6e+04,3.2e+03)
+DEAL::(-3.6e+04,3.2e+03)
--- /dev/null
+//----------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2007 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::l1_norm. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_15/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ deallog << m.l1_norm() << std::endl;
+}
+
--- /dev/null
+
+DEAL::2.2e+02
+DEAL::2.2e+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::linfty_norm. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_16/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ deallog << m.linfty_norm() << std::endl;
+}
+
--- /dev/null
+
+DEAL::2.3e+02
+DEAL::2.3e+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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_17/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ deallog << m.frobenius_norm() << std::endl;
+}
+
--- /dev/null
+
+DEAL::1.7e+02
+DEAL::1.7e+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::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_complex_full_matrix_18/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ deallog << m.relative_symmetry_norm2() << std::endl;
+}
+
--- /dev/null
+
+DEAL::0.081
+DEAL::0.081
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_20/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ deallog << m.trace() << std::endl;
+}
+
--- /dev/null
+
+DEAL::(70.,1.6e+02)
+DEAL::(70.,1.6e+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::print. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_21/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m.print (deallog.get_file_stream());
+}
+
--- /dev/null
+
+(2.0,12.)(3.0,15.)(4.0,18.)(5.0,21.)(6.0,24.)
+(4.0,16.)(6.0,20.)(8.0,24.)(10.,28.)(12.,32.)
+(6.0,20.)(9.0,25.)(12.,30.)(15.,35.)(18.,40.)
+(8.0,24.)(12.,30.)(16.,36.)(20.,42.)(24.,48.)
+(10.,28.)(15.,35.)(20.,42.)(25.,49.)(30.,56.)
+(2.0,12.)(3.0,15.)(4.0,18.)(5.0,21.)(6.0,24.)
+(4.0,16.)(6.0,20.)(8.0,24.)(10.,28.)(12.,32.)
+(6.0,20.)(9.0,25.)(12.,30.)(15.,35.)(18.,40.)
+(8.0,24.)(12.,30.)(16.,36.)(20.,42.)(24.,48.)
+(10.,28.)(15.,35.)(20.,42.)(25.,49.)(30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_22/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m *= 3.1415;
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (6.3,38.)
+DEAL::0 1 (9.4,47.)
+DEAL::0 2 (13.,57.)
+DEAL::0 3 (16.,66.)
+DEAL::0 4 (19.,75.)
+DEAL::1 0 (13.,50.)
+DEAL::1 1 (19.,63.)
+DEAL::1 2 (25.,75.)
+DEAL::1 3 (31.,88.)
+DEAL::1 4 (38.,1.0e+02)
+DEAL::2 0 (19.,63.)
+DEAL::2 1 (28.,79.)
+DEAL::2 2 (38.,94.)
+DEAL::2 3 (47.,1.1e+02)
+DEAL::2 4 (57.,1.3e+02)
+DEAL::3 0 (25.,75.)
+DEAL::3 1 (38.,94.)
+DEAL::3 2 (50.,1.1e+02)
+DEAL::3 3 (63.,1.3e+02)
+DEAL::3 4 (75.,1.5e+02)
+DEAL::4 0 (31.,88.)
+DEAL::4 1 (47.,1.1e+02)
+DEAL::4 2 (63.,1.3e+02)
+DEAL::4 3 (79.,1.5e+02)
+DEAL::4 4 (94.,1.8e+02)
+DEAL::0 0 (6.3,38.)
+DEAL::0 1 (9.4,47.)
+DEAL::0 2 (13.,57.)
+DEAL::0 3 (16.,66.)
+DEAL::0 4 (19.,75.)
+DEAL::1 0 (13.,50.)
+DEAL::1 1 (19.,63.)
+DEAL::1 2 (25.,75.)
+DEAL::1 3 (31.,88.)
+DEAL::1 4 (38.,1.0e+02)
+DEAL::2 0 (19.,63.)
+DEAL::2 1 (28.,79.)
+DEAL::2 2 (38.,94.)
+DEAL::2 3 (47.,1.1e+02)
+DEAL::2 4 (57.,1.3e+02)
+DEAL::3 0 (25.,75.)
+DEAL::3 1 (38.,94.)
+DEAL::3 2 (50.,1.1e+02)
+DEAL::3 3 (63.,1.3e+02)
+DEAL::3 4 (75.,1.5e+02)
+DEAL::4 0 (31.,88.)
+DEAL::4 1 (47.,1.1e+02)
+DEAL::4 2 (63.,1.3e+02)
+DEAL::4 3 (79.,1.5e+02)
+DEAL::4 4 (94.,1.8e+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::operator/=. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_23/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m /= 3.1415926;
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (0.64,3.8)
+DEAL::0 1 (0.95,4.8)
+DEAL::0 2 (1.3,5.7)
+DEAL::0 3 (1.6,6.7)
+DEAL::0 4 (1.9,7.6)
+DEAL::1 0 (1.3,5.1)
+DEAL::1 1 (1.9,6.4)
+DEAL::1 2 (2.5,7.6)
+DEAL::1 3 (3.2,8.9)
+DEAL::1 4 (3.8,10.)
+DEAL::2 0 (1.9,6.4)
+DEAL::2 1 (2.9,8.0)
+DEAL::2 2 (3.8,9.5)
+DEAL::2 3 (4.8,11.)
+DEAL::2 4 (5.7,13.)
+DEAL::3 0 (2.5,7.6)
+DEAL::3 1 (3.8,9.5)
+DEAL::3 2 (5.1,11.)
+DEAL::3 3 (6.4,13.)
+DEAL::3 4 (7.6,15.)
+DEAL::4 0 (3.2,8.9)
+DEAL::4 1 (4.8,11.)
+DEAL::4 2 (6.4,13.)
+DEAL::4 3 (8.0,16.)
+DEAL::4 4 (9.5,18.)
+DEAL::0 0 (0.64,3.8)
+DEAL::0 1 (0.95,4.8)
+DEAL::0 2 (1.3,5.7)
+DEAL::0 3 (1.6,6.7)
+DEAL::0 4 (1.9,7.6)
+DEAL::1 0 (1.3,5.1)
+DEAL::1 1 (1.9,6.4)
+DEAL::1 2 (2.5,7.6)
+DEAL::1 3 (3.2,8.9)
+DEAL::1 4 (3.8,10.)
+DEAL::2 0 (1.9,6.4)
+DEAL::2 1 (2.9,8.0)
+DEAL::2 2 (3.8,9.5)
+DEAL::2 3 (4.8,11.)
+DEAL::2 4 (5.7,13.)
+DEAL::3 0 (2.5,7.6)
+DEAL::3 1 (3.8,9.5)
+DEAL::3 2 (5.1,11.)
+DEAL::3 3 (6.4,13.)
+DEAL::3 4 (7.6,15.)
+DEAL::4 0 (3.2,8.9)
+DEAL::4 1 (4.8,11.)
+DEAL::4 2 (6.4,13.)
+DEAL::4 3 (8.0,16.)
+DEAL::4 4 (9.5,18.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_24/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n;
+ make_complex_matrix (m);
+ make_complex_matrix (n);
+ m.add (3.14159, n);
+
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (8.3,50.)
+DEAL::0 1 (12.,62.)
+DEAL::0 2 (17.,75.)
+DEAL::0 3 (21.,87.)
+DEAL::0 4 (25.,99.)
+DEAL::1 0 (17.,66.)
+DEAL::1 1 (25.,83.)
+DEAL::1 2 (33.,99.)
+DEAL::1 3 (41.,1.2e+02)
+DEAL::1 4 (50.,1.3e+02)
+DEAL::2 0 (25.,83.)
+DEAL::2 1 (37.,1.0e+02)
+DEAL::2 2 (50.,1.2e+02)
+DEAL::2 3 (62.,1.4e+02)
+DEAL::2 4 (75.,1.7e+02)
+DEAL::3 0 (33.,99.)
+DEAL::3 1 (50.,1.2e+02)
+DEAL::3 2 (66.,1.5e+02)
+DEAL::3 3 (83.,1.7e+02)
+DEAL::3 4 (99.,2.0e+02)
+DEAL::4 0 (41.,1.2e+02)
+DEAL::4 1 (62.,1.4e+02)
+DEAL::4 2 (83.,1.7e+02)
+DEAL::4 3 (1.0e+02,2.0e+02)
+DEAL::4 4 (1.2e+02,2.3e+02)
+DEAL::0 0 (8.3,50.)
+DEAL::0 1 (12.,62.)
+DEAL::0 2 (17.,75.)
+DEAL::0 3 (21.,87.)
+DEAL::0 4 (25.,99.)
+DEAL::1 0 (17.,66.)
+DEAL::1 1 (25.,83.)
+DEAL::1 2 (33.,99.)
+DEAL::1 3 (41.,1.2e+02)
+DEAL::1 4 (50.,1.3e+02)
+DEAL::2 0 (25.,83.)
+DEAL::2 1 (37.,1.0e+02)
+DEAL::2 2 (50.,1.2e+02)
+DEAL::2 3 (62.,1.4e+02)
+DEAL::2 4 (75.,1.7e+02)
+DEAL::3 0 (33.,99.)
+DEAL::3 1 (50.,1.2e+02)
+DEAL::3 2 (66.,1.5e+02)
+DEAL::3 3 (83.,1.7e+02)
+DEAL::3 4 (99.,2.0e+02)
+DEAL::4 0 (41.,1.2e+02)
+DEAL::4 1 (62.,1.4e+02)
+DEAL::4 2 (83.,1.7e+02)
+DEAL::4 3 (1.0e+02,2.0e+02)
+DEAL::4 4 (1.2e+02,2.3e+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::add(4). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_25/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n,o;
+ make_complex_matrix (m);
+ make_complex_matrix (n);
+ make_complex_matrix (o);
+ m.add (3.1415, n, 2.718, o);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (14.,82.)
+DEAL::0 1 (21.,1.0e+02)
+DEAL::0 2 (27.,1.2e+02)
+DEAL::0 3 (34.,1.4e+02)
+DEAL::0 4 (41.,1.6e+02)
+DEAL::1 0 (27.,1.1e+02)
+DEAL::1 1 (41.,1.4e+02)
+DEAL::1 2 (55.,1.6e+02)
+DEAL::1 3 (69.,1.9e+02)
+DEAL::1 4 (82.,2.2e+02)
+DEAL::2 0 (41.,1.4e+02)
+DEAL::2 1 (62.,1.7e+02)
+DEAL::2 2 (82.,2.1e+02)
+DEAL::2 3 (1.0e+02,2.4e+02)
+DEAL::2 4 (1.2e+02,2.7e+02)
+DEAL::3 0 (55.,1.6e+02)
+DEAL::3 1 (82.,2.1e+02)
+DEAL::3 2 (1.1e+02,2.5e+02)
+DEAL::3 3 (1.4e+02,2.9e+02)
+DEAL::3 4 (1.6e+02,3.3e+02)
+DEAL::4 0 (69.,1.9e+02)
+DEAL::4 1 (1.0e+02,2.4e+02)
+DEAL::4 2 (1.4e+02,2.9e+02)
+DEAL::4 3 (1.7e+02,3.4e+02)
+DEAL::4 4 (2.1e+02,3.8e+02)
+DEAL::0 0 (14.,82.)
+DEAL::0 1 (21.,1.0e+02)
+DEAL::0 2 (27.,1.2e+02)
+DEAL::0 3 (34.,1.4e+02)
+DEAL::0 4 (41.,1.6e+02)
+DEAL::1 0 (27.,1.1e+02)
+DEAL::1 1 (41.,1.4e+02)
+DEAL::1 2 (55.,1.6e+02)
+DEAL::1 3 (69.,1.9e+02)
+DEAL::1 4 (82.,2.2e+02)
+DEAL::2 0 (41.,1.4e+02)
+DEAL::2 1 (62.,1.7e+02)
+DEAL::2 2 (82.,2.1e+02)
+DEAL::2 3 (1.0e+02,2.4e+02)
+DEAL::2 4 (1.2e+02,2.7e+02)
+DEAL::3 0 (55.,1.6e+02)
+DEAL::3 1 (82.,2.1e+02)
+DEAL::3 2 (1.1e+02,2.5e+02)
+DEAL::3 3 (1.4e+02,2.9e+02)
+DEAL::3 4 (1.6e+02,3.3e+02)
+DEAL::4 0 (69.,1.9e+02)
+DEAL::4 1 (1.0e+02,2.4e+02)
+DEAL::4 2 (1.4e+02,2.9e+02)
+DEAL::4 3 (1.7e+02,3.4e+02)
+DEAL::4 4 (2.1e+02,3.8e+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::add (6). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_26/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n,o,p;
+ make_complex_matrix (m);
+ make_complex_matrix (n);
+ make_complex_matrix (o);
+ make_complex_matrix (p);
+ m.add (3.1415, n, 2.718, o, 1.414, p);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (17.,99.)
+DEAL::0 1 (25.,1.2e+02)
+DEAL::0 2 (33.,1.5e+02)
+DEAL::0 3 (41.,1.7e+02)
+DEAL::0 4 (50.,2.0e+02)
+DEAL::1 0 (33.,1.3e+02)
+DEAL::1 1 (50.,1.7e+02)
+DEAL::1 2 (66.,2.0e+02)
+DEAL::1 3 (83.,2.3e+02)
+DEAL::1 4 (99.,2.6e+02)
+DEAL::2 0 (50.,1.7e+02)
+DEAL::2 1 (74.,2.1e+02)
+DEAL::2 2 (99.,2.5e+02)
+DEAL::2 3 (1.2e+02,2.9e+02)
+DEAL::2 4 (1.5e+02,3.3e+02)
+DEAL::3 0 (66.,2.0e+02)
+DEAL::3 1 (99.,2.5e+02)
+DEAL::3 2 (1.3e+02,3.0e+02)
+DEAL::3 3 (1.7e+02,3.5e+02)
+DEAL::3 4 (2.0e+02,4.0e+02)
+DEAL::4 0 (83.,2.3e+02)
+DEAL::4 1 (1.2e+02,2.9e+02)
+DEAL::4 2 (1.7e+02,3.5e+02)
+DEAL::4 3 (2.1e+02,4.1e+02)
+DEAL::4 4 (2.5e+02,4.6e+02)
+DEAL::0 0 (17.,99.)
+DEAL::0 1 (25.,1.2e+02)
+DEAL::0 2 (33.,1.5e+02)
+DEAL::0 3 (41.,1.7e+02)
+DEAL::0 4 (50.,2.0e+02)
+DEAL::1 0 (33.,1.3e+02)
+DEAL::1 1 (50.,1.7e+02)
+DEAL::1 2 (66.,2.0e+02)
+DEAL::1 3 (83.,2.3e+02)
+DEAL::1 4 (99.,2.6e+02)
+DEAL::2 0 (50.,1.7e+02)
+DEAL::2 1 (74.,2.1e+02)
+DEAL::2 2 (99.,2.5e+02)
+DEAL::2 3 (1.2e+02,2.9e+02)
+DEAL::2 4 (1.5e+02,3.3e+02)
+DEAL::3 0 (66.,2.0e+02)
+DEAL::3 1 (99.,2.5e+02)
+DEAL::3 2 (1.3e+02,3.0e+02)
+DEAL::3 3 (1.7e+02,3.5e+02)
+DEAL::3 4 (2.0e+02,4.0e+02)
+DEAL::4 0 (83.,2.3e+02)
+DEAL::4 1 (1.2e+02,2.9e+02)
+DEAL::4 2 (1.7e+02,3.5e+02)
+DEAL::4 3 (2.1e+02,4.1e+02)
+DEAL::4 4 (2.5e+02,4.6e+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::add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_27/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m(10,10),n;
+ make_complex_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.,70.)
+DEAL::3 5 (36.,80.)
+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.,84.)
+DEAL::4 5 (48.,96.)
+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.,98.)
+DEAL::5 5 (60.,1.1e+02)
+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.,70.)
+DEAL::3 5 (36.,80.)
+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.,84.)
+DEAL::4 5 (48.,96.)
+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.,98.)
+DEAL::5 5 (60.,1.1e+02)
+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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_28/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m(10,10),n;
+ make_complex_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.,70.)
+DEAL::3 5 (40.,84.)
+DEAL::3 6 (50.,98.)
+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.,80.)
+DEAL::4 5 (48.,96.)
+DEAL::4 6 (60.,1.1e+02)
+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.,70.)
+DEAL::3 5 (40.,84.)
+DEAL::3 6 (50.,98.)
+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.,80.)
+DEAL::4 5 (48.,96.)
+DEAL::4 6 (60.,1.1e+02)
+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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_29/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m.add_row (2, 3.141, 4);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (37.,1.1e+02)
+DEAL::2 1 (56.,1.3e+02)
+DEAL::2 2 (75.,1.6e+02)
+DEAL::2 3 (94.,1.9e+02)
+DEAL::2 4 (1.1e+02,2.2e+02)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (37.,1.1e+02)
+DEAL::2 1 (56.,1.3e+02)
+DEAL::2 2 (75.,1.6e+02)
+DEAL::2 3 (94.,1.9e+02)
+DEAL::2 4 (1.1e+02,2.2e+02)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_30/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m.add_row (2, 3.141, 4, 2.7182, 3);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (59.,1.7e+02)
+DEAL::2 1 (89.,2.2e+02)
+DEAL::2 2 (1.2e+02,2.6e+02)
+DEAL::2 3 (1.5e+02,3.0e+02)
+DEAL::2 4 (1.8e+02,3.5e+02)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (59.,1.7e+02)
+DEAL::2 1 (89.,2.2e+02)
+DEAL::2 2 (1.2e+02,2.6e+02)
+DEAL::2 3 (1.5e+02,3.0e+02)
+DEAL::2 4 (1.8e+02,3.5e+02)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_31/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m.add_col (2, 3.141, 4);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (23.,93.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (46.,1.2e+02)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (69.,1.6e+02)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (91.,1.9e+02)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (1.1e+02,2.2e+02)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (23.,93.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (46.,1.2e+02)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (69.,1.6e+02)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (91.,1.9e+02)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (1.1e+02,2.2e+02)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_32/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m.add_col (2, 3.141, 4, 2.718, 3);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (36.,1.5e+02)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (73.,2.0e+02)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (1.1e+02,2.5e+02)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (1.5e+02,3.0e+02)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (1.8e+02,3.5e+02)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (36.,1.5e+02)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (73.,2.0e+02)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (1.1e+02,2.5e+02)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (1.5e+02,3.0e+02)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (1.8e+02,3.5e+02)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_33/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m.swap_row (2, 4);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (10.,28.)
+DEAL::2 1 (15.,35.)
+DEAL::2 2 (20.,42.)
+DEAL::2 3 (25.,49.)
+DEAL::2 4 (30.,56.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (6.0,20.)
+DEAL::4 1 (9.0,25.)
+DEAL::4 2 (12.,30.)
+DEAL::4 3 (15.,35.)
+DEAL::4 4 (18.,40.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (10.,28.)
+DEAL::2 1 (15.,35.)
+DEAL::2 2 (20.,42.)
+DEAL::2 3 (25.,49.)
+DEAL::2 4 (30.,56.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (6.0,20.)
+DEAL::4 1 (9.0,25.)
+DEAL::4 2 (12.,30.)
+DEAL::4 3 (15.,35.)
+DEAL::4 4 (18.,40.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_34/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m.swap_col (2, 4);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (6.0,24.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (4.0,18.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (12.,32.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (8.0,24.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (18.,40.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (12.,30.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (24.,48.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (16.,36.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (30.,56.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (20.,42.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (6.0,24.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (4.0,18.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (12.,32.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (8.0,24.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (18.,40.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (12.,30.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (24.,48.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (16.,36.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (30.,56.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (20.,42.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_35/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m, n;
+ make_complex_matrix (m);
+ make_complex_matrix (n);
+ m.add_diag (3.141, m);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (8.3,50.)
+DEAL::0 1 (12.,62.)
+DEAL::0 2 (17.,75.)
+DEAL::0 3 (21.,87.)
+DEAL::0 4 (25.,99.)
+DEAL::1 0 (17.,66.)
+DEAL::1 1 (25.,83.)
+DEAL::1 2 (33.,99.)
+DEAL::1 3 (41.,1.2e+02)
+DEAL::1 4 (50.,1.3e+02)
+DEAL::2 0 (25.,83.)
+DEAL::2 1 (37.,1.0e+02)
+DEAL::2 2 (50.,1.2e+02)
+DEAL::2 3 (62.,1.4e+02)
+DEAL::2 4 (75.,1.7e+02)
+DEAL::3 0 (33.,99.)
+DEAL::3 1 (50.,1.2e+02)
+DEAL::3 2 (66.,1.5e+02)
+DEAL::3 3 (83.,1.7e+02)
+DEAL::3 4 (99.,2.0e+02)
+DEAL::4 0 (41.,1.2e+02)
+DEAL::4 1 (62.,1.4e+02)
+DEAL::4 2 (83.,1.7e+02)
+DEAL::4 3 (1.0e+02,2.0e+02)
+DEAL::4 4 (1.2e+02,2.3e+02)
+DEAL::0 0 (8.3,50.)
+DEAL::0 1 (12.,62.)
+DEAL::0 2 (17.,75.)
+DEAL::0 3 (21.,87.)
+DEAL::0 4 (25.,99.)
+DEAL::1 0 (17.,66.)
+DEAL::1 1 (25.,83.)
+DEAL::1 2 (33.,99.)
+DEAL::1 3 (41.,1.2e+02)
+DEAL::1 4 (50.,1.3e+02)
+DEAL::2 0 (25.,83.)
+DEAL::2 1 (37.,1.0e+02)
+DEAL::2 2 (50.,1.2e+02)
+DEAL::2 3 (62.,1.4e+02)
+DEAL::2 4 (75.,1.7e+02)
+DEAL::3 0 (33.,99.)
+DEAL::3 1 (50.,1.2e+02)
+DEAL::3 2 (66.,1.5e+02)
+DEAL::3 3 (83.,1.7e+02)
+DEAL::3 4 (99.,2.0e+02)
+DEAL::4 0 (41.,1.2e+02)
+DEAL::4 1 (62.,1.4e+02)
+DEAL::4 2 (83.,1.7e+02)
+DEAL::4 3 (1.0e+02,2.0e+02)
+DEAL::4 4 (1.2e+02,2.3e+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::diag_add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_36/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m.diagadd (3.141);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (5.1,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (9.1,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (15.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (23.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (33.,56.)
+DEAL::0 0 (5.1,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (9.1,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (15.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (23.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (33.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_37/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n;
+ make_complex_matrix (m);
+ make_complex_matrix (n);
+ m.equ (3.14159, n);
+
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (6.3,38.)
+DEAL::0 1 (9.4,47.)
+DEAL::0 2 (13.,57.)
+DEAL::0 3 (16.,66.)
+DEAL::0 4 (19.,75.)
+DEAL::1 0 (13.,50.)
+DEAL::1 1 (19.,63.)
+DEAL::1 2 (25.,75.)
+DEAL::1 3 (31.,88.)
+DEAL::1 4 (38.,1.0e+02)
+DEAL::2 0 (19.,63.)
+DEAL::2 1 (28.,79.)
+DEAL::2 2 (38.,94.)
+DEAL::2 3 (47.,1.1e+02)
+DEAL::2 4 (57.,1.3e+02)
+DEAL::3 0 (25.,75.)
+DEAL::3 1 (38.,94.)
+DEAL::3 2 (50.,1.1e+02)
+DEAL::3 3 (63.,1.3e+02)
+DEAL::3 4 (75.,1.5e+02)
+DEAL::4 0 (31.,88.)
+DEAL::4 1 (47.,1.1e+02)
+DEAL::4 2 (63.,1.3e+02)
+DEAL::4 3 (79.,1.5e+02)
+DEAL::4 4 (94.,1.8e+02)
+DEAL::0 0 (6.3,38.)
+DEAL::0 1 (9.4,47.)
+DEAL::0 2 (13.,57.)
+DEAL::0 3 (16.,66.)
+DEAL::0 4 (19.,75.)
+DEAL::1 0 (13.,50.)
+DEAL::1 1 (19.,63.)
+DEAL::1 2 (25.,75.)
+DEAL::1 3 (31.,88.)
+DEAL::1 4 (38.,1.0e+02)
+DEAL::2 0 (19.,63.)
+DEAL::2 1 (28.,79.)
+DEAL::2 2 (38.,94.)
+DEAL::2 3 (47.,1.1e+02)
+DEAL::2 4 (57.,1.3e+02)
+DEAL::3 0 (25.,75.)
+DEAL::3 1 (38.,94.)
+DEAL::3 2 (50.,1.1e+02)
+DEAL::3 3 (63.,1.3e+02)
+DEAL::3 4 (75.,1.5e+02)
+DEAL::4 0 (31.,88.)
+DEAL::4 1 (47.,1.1e+02)
+DEAL::4 2 (63.,1.3e+02)
+DEAL::4 3 (79.,1.5e+02)
+DEAL::4 4 (94.,1.8e+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::equ(4). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_38/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n,o;
+ make_complex_matrix (m);
+ make_complex_matrix (n);
+ make_complex_matrix (o);
+ m.equ (3.1415, n, 2.718, o);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (12.,70.)
+DEAL::0 1 (18.,88.)
+DEAL::0 2 (23.,1.1e+02)
+DEAL::0 3 (29.,1.2e+02)
+DEAL::0 4 (35.,1.4e+02)
+DEAL::1 0 (23.,94.)
+DEAL::1 1 (35.,1.2e+02)
+DEAL::1 2 (47.,1.4e+02)
+DEAL::1 3 (59.,1.6e+02)
+DEAL::1 4 (70.,1.9e+02)
+DEAL::2 0 (35.,1.2e+02)
+DEAL::2 1 (53.,1.5e+02)
+DEAL::2 2 (70.,1.8e+02)
+DEAL::2 3 (88.,2.1e+02)
+DEAL::2 4 (1.1e+02,2.3e+02)
+DEAL::3 0 (47.,1.4e+02)
+DEAL::3 1 (70.,1.8e+02)
+DEAL::3 2 (94.,2.1e+02)
+DEAL::3 3 (1.2e+02,2.5e+02)
+DEAL::3 4 (1.4e+02,2.8e+02)
+DEAL::4 0 (59.,1.6e+02)
+DEAL::4 1 (88.,2.1e+02)
+DEAL::4 2 (1.2e+02,2.5e+02)
+DEAL::4 3 (1.5e+02,2.9e+02)
+DEAL::4 4 (1.8e+02,3.3e+02)
+DEAL::0 0 (12.,70.)
+DEAL::0 1 (18.,88.)
+DEAL::0 2 (23.,1.1e+02)
+DEAL::0 3 (29.,1.2e+02)
+DEAL::0 4 (35.,1.4e+02)
+DEAL::1 0 (23.,94.)
+DEAL::1 1 (35.,1.2e+02)
+DEAL::1 2 (47.,1.4e+02)
+DEAL::1 3 (59.,1.6e+02)
+DEAL::1 4 (70.,1.9e+02)
+DEAL::2 0 (35.,1.2e+02)
+DEAL::2 1 (53.,1.5e+02)
+DEAL::2 2 (70.,1.8e+02)
+DEAL::2 3 (88.,2.1e+02)
+DEAL::2 4 (1.1e+02,2.3e+02)
+DEAL::3 0 (47.,1.4e+02)
+DEAL::3 1 (70.,1.8e+02)
+DEAL::3 2 (94.,2.1e+02)
+DEAL::3 3 (1.2e+02,2.5e+02)
+DEAL::3 4 (1.4e+02,2.8e+02)
+DEAL::4 0 (59.,1.6e+02)
+DEAL::4 1 (88.,2.1e+02)
+DEAL::4 2 (1.2e+02,2.5e+02)
+DEAL::4 3 (1.5e+02,2.9e+02)
+DEAL::4 4 (1.8e+02,3.3e+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::equ (6). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_39/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m,n,o,p;
+ make_complex_matrix (m);
+ make_complex_matrix (n);
+ make_complex_matrix (o);
+ make_complex_matrix (p);
+ m.equ (3.1415, n, 2.718, o, 1.414, p);
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (15.,87.)
+DEAL::0 1 (22.,1.1e+02)
+DEAL::0 2 (29.,1.3e+02)
+DEAL::0 3 (36.,1.5e+02)
+DEAL::0 4 (44.,1.7e+02)
+DEAL::1 0 (29.,1.2e+02)
+DEAL::1 1 (44.,1.5e+02)
+DEAL::1 2 (58.,1.7e+02)
+DEAL::1 3 (73.,2.0e+02)
+DEAL::1 4 (87.,2.3e+02)
+DEAL::2 0 (44.,1.5e+02)
+DEAL::2 1 (65.,1.8e+02)
+DEAL::2 2 (87.,2.2e+02)
+DEAL::2 3 (1.1e+02,2.5e+02)
+DEAL::2 4 (1.3e+02,2.9e+02)
+DEAL::3 0 (58.,1.7e+02)
+DEAL::3 1 (87.,2.2e+02)
+DEAL::3 2 (1.2e+02,2.6e+02)
+DEAL::3 3 (1.5e+02,3.1e+02)
+DEAL::3 4 (1.7e+02,3.5e+02)
+DEAL::4 0 (73.,2.0e+02)
+DEAL::4 1 (1.1e+02,2.5e+02)
+DEAL::4 2 (1.5e+02,3.1e+02)
+DEAL::4 3 (1.8e+02,3.6e+02)
+DEAL::4 4 (2.2e+02,4.1e+02)
+DEAL::0 0 (15.,87.)
+DEAL::0 1 (22.,1.1e+02)
+DEAL::0 2 (29.,1.3e+02)
+DEAL::0 3 (36.,1.5e+02)
+DEAL::0 4 (44.,1.7e+02)
+DEAL::1 0 (29.,1.2e+02)
+DEAL::1 1 (44.,1.5e+02)
+DEAL::1 2 (58.,1.7e+02)
+DEAL::1 3 (73.,2.0e+02)
+DEAL::1 4 (87.,2.3e+02)
+DEAL::2 0 (44.,1.5e+02)
+DEAL::2 1 (65.,1.8e+02)
+DEAL::2 2 (87.,2.2e+02)
+DEAL::2 3 (1.1e+02,2.5e+02)
+DEAL::2 4 (1.3e+02,2.9e+02)
+DEAL::3 0 (58.,1.7e+02)
+DEAL::3 1 (87.,2.2e+02)
+DEAL::3 2 (1.2e+02,2.6e+02)
+DEAL::3 3 (1.5e+02,3.1e+02)
+DEAL::3 4 (1.7e+02,3.5e+02)
+DEAL::4 0 (73.,2.0e+02)
+DEAL::4 1 (1.1e+02,2.5e+02)
+DEAL::4 2 (1.5e+02,3.1e+02)
+DEAL::4 3 (1.8e+02,3.6e+02)
+DEAL::4 4 (2.2e+02,4.1e+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::symmetrize. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_40/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ m.symmetrize ();
+ print_matrix (m);
+}
+
--- /dev/null
+
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.5,16.)
+DEAL::0 2 (5.0,19.)
+DEAL::0 3 (6.5,22.)
+DEAL::0 4 (8.0,26.)
+DEAL::1 0 (3.5,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.5,24.)
+DEAL::1 3 (11.,29.)
+DEAL::1 4 (14.,34.)
+DEAL::2 0 (5.0,19.)
+DEAL::2 1 (8.5,24.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (16.,36.)
+DEAL::2 4 (19.,41.)
+DEAL::3 0 (6.5,22.)
+DEAL::3 1 (11.,29.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (8.0,26.)
+DEAL::4 1 (14.,34.)
+DEAL::4 2 (19.,41.)
+DEAL::4 3 (24.,48.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.5,16.)
+DEAL::0 2 (5.0,19.)
+DEAL::0 3 (6.5,22.)
+DEAL::0 4 (8.0,26.)
+DEAL::1 0 (3.5,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.5,24.)
+DEAL::1 3 (11.,29.)
+DEAL::1 4 (14.,34.)
+DEAL::2 0 (5.0,19.)
+DEAL::2 1 (8.5,24.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (16.,36.)
+DEAL::2 4 (19.,41.)
+DEAL::3 0 (6.5,22.)
+DEAL::3 1 (11.,29.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (8.0,26.)
+DEAL::4 1 (14.,34.)
+DEAL::4 2 (19.,41.)
+DEAL::4 3 (24.,48.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_41/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_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.018,0.0043)
+DEAL::0 1 (0.0021,0.00086)
+DEAL::0 2 (0.0022,0.0014)
+DEAL::0 3 (0.0020,0.0020)
+DEAL::0 4 (0.0014,0.0024)
+DEAL::1 0 (0.0023,0.00080)
+DEAL::1 1 (0.016,0.0072)
+DEAL::1 2 (0.0030,0.0022)
+DEAL::1 3 (0.0026,0.0031)
+DEAL::1 4 (0.0017,0.0037)
+DEAL::2 0 (0.0026,0.0014)
+DEAL::2 1 (0.0032,0.0022)
+DEAL::2 2 (0.013,0.0094)
+DEAL::2 3 (0.0026,0.0046)
+DEAL::2 4 (0.0013,0.0051)
+DEAL::3 0 (0.0025,0.0022)
+DEAL::3 1 (0.0030,0.0033)
+DEAL::3 2 (0.0027,0.0047)
+DEAL::3 3 (0.0097,0.0099)
+DEAL::3 4 (-1.3e-05,0.0061)
+DEAL::4 0 (0.0020,0.0027)
+DEAL::4 1 (0.0021,0.0040)
+DEAL::4 2 (0.0015,0.0054)
+DEAL::4 3 (-6.6e-05,0.0062)
+DEAL::4 4 (0.0072,0.0098)
+DEAL::0 0 (0.018,0.0043)
+DEAL::0 1 (0.0021,0.00086)
+DEAL::0 2 (0.0022,0.0014)
+DEAL::0 3 (0.0020,0.0020)
+DEAL::0 4 (0.0014,0.0024)
+DEAL::1 0 (0.0023,0.00080)
+DEAL::1 1 (0.016,0.0072)
+DEAL::1 2 (0.0030,0.0022)
+DEAL::1 3 (0.0026,0.0031)
+DEAL::1 4 (0.0017,0.0037)
+DEAL::2 0 (0.0026,0.0014)
+DEAL::2 1 (0.0032,0.0022)
+DEAL::2 2 (0.013,0.0094)
+DEAL::2 3 (0.0026,0.0046)
+DEAL::2 4 (0.0013,0.0051)
+DEAL::3 0 (0.0025,0.0022)
+DEAL::3 1 (0.0030,0.0033)
+DEAL::3 2 (0.0027,0.0047)
+DEAL::3 3 (0.0097,0.0099)
+DEAL::3 4 (-1.3e-05,0.0061)
+DEAL::4 0 (0.0020,0.0027)
+DEAL::4 1 (0.0021,0.0040)
+DEAL::4 2 (0.0015,0.0054)
+DEAL::4 3 (-6.6e-05,0.0062)
+DEAL::4 4 (0.0072,0.0098)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_43/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m, n, o;
+ make_complex_matrix (m);
+ make_complex_matrix (n);
+ make_complex_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,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (-1.8e+03,1.1e+03)
+DEAL::0 1 (-2.2e+03,1.5e+03)
+DEAL::0 2 (-2.6e+03,1.9e+03)
+DEAL::0 3 (-3.0e+03,2.3e+03)
+DEAL::0 4 (-3.4e+03,2.7e+03)
+DEAL::1 0 (-2.3e+03,1.7e+03)
+DEAL::1 1 (-2.8e+03,2.3e+03)
+DEAL::1 2 (-3.3e+03,2.9e+03)
+DEAL::1 3 (-3.8e+03,3.6e+03)
+DEAL::1 4 (-4.3e+03,4.2e+03)
+DEAL::2 0 (-2.8e+03,2.3e+03)
+DEAL::2 1 (-3.4e+03,3.2e+03)
+DEAL::2 2 (-3.9e+03,4.0e+03)
+DEAL::2 3 (-4.5e+03,4.8e+03)
+DEAL::2 4 (-5.1e+03,5.7e+03)
+DEAL::3 0 (-3.3e+03,3.0e+03)
+DEAL::3 1 (-3.9e+03,4.0e+03)
+DEAL::3 2 (-4.6e+03,5.1e+03)
+DEAL::3 3 (-5.3e+03,6.1e+03)
+DEAL::3 4 (-6.0e+03,7.2e+03)
+DEAL::4 0 (-3.8e+03,3.6e+03)
+DEAL::4 1 (-4.5e+03,4.9e+03)
+DEAL::4 2 (-5.3e+03,6.1e+03)
+DEAL::4 3 (-6.1e+03,7.4e+03)
+DEAL::4 4 (-6.8e+03,8.7e+03)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (-1.8e+03,1.0e+03)
+DEAL::0 1 (-2.2e+03,1.4e+03)
+DEAL::0 2 (-2.6e+03,1.9e+03)
+DEAL::0 3 (-3.0e+03,2.3e+03)
+DEAL::0 4 (-3.4e+03,2.7e+03)
+DEAL::1 0 (-2.3e+03,1.7e+03)
+DEAL::1 1 (-2.8e+03,2.3e+03)
+DEAL::1 2 (-3.3e+03,2.9e+03)
+DEAL::1 3 (-3.8e+03,3.5e+03)
+DEAL::1 4 (-4.3e+03,4.2e+03)
+DEAL::2 0 (-2.8e+03,2.3e+03)
+DEAL::2 1 (-3.4e+03,3.2e+03)
+DEAL::2 2 (-4.0e+03,4.0e+03)
+DEAL::2 3 (-4.6e+03,4.8e+03)
+DEAL::2 4 (-5.1e+03,5.6e+03)
+DEAL::3 0 (-3.3e+03,3.0e+03)
+DEAL::3 1 (-4.0e+03,4.0e+03)
+DEAL::3 2 (-4.6e+03,5.0e+03)
+DEAL::3 3 (-5.3e+03,6.1e+03)
+DEAL::3 4 (-6.0e+03,7.1e+03)
+DEAL::4 0 (-3.8e+03,3.6e+03)
+DEAL::4 1 (-4.6e+03,4.8e+03)
+DEAL::4 2 (-5.3e+03,6.1e+03)
+DEAL::4 3 (-6.1e+03,7.4e+03)
+DEAL::4 4 (-6.9e+03,8.6e+03)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (-1.8e+03,1.1e+03)
+DEAL::0 1 (-2.2e+03,1.5e+03)
+DEAL::0 2 (-2.6e+03,1.9e+03)
+DEAL::0 3 (-3.0e+03,2.3e+03)
+DEAL::0 4 (-3.4e+03,2.7e+03)
+DEAL::1 0 (-2.3e+03,1.7e+03)
+DEAL::1 1 (-2.8e+03,2.3e+03)
+DEAL::1 2 (-3.3e+03,2.9e+03)
+DEAL::1 3 (-3.8e+03,3.6e+03)
+DEAL::1 4 (-4.3e+03,4.2e+03)
+DEAL::2 0 (-2.8e+03,2.3e+03)
+DEAL::2 1 (-3.4e+03,3.2e+03)
+DEAL::2 2 (-3.9e+03,4.0e+03)
+DEAL::2 3 (-4.5e+03,4.8e+03)
+DEAL::2 4 (-5.1e+03,5.7e+03)
+DEAL::3 0 (-3.3e+03,3.0e+03)
+DEAL::3 1 (-3.9e+03,4.0e+03)
+DEAL::3 2 (-4.6e+03,5.1e+03)
+DEAL::3 3 (-5.3e+03,6.1e+03)
+DEAL::3 4 (-6.0e+03,7.2e+03)
+DEAL::4 0 (-3.8e+03,3.6e+03)
+DEAL::4 1 (-4.5e+03,4.9e+03)
+DEAL::4 2 (-5.3e+03,6.1e+03)
+DEAL::4 3 (-6.1e+03,7.4e+03)
+DEAL::4 4 (-6.8e+03,8.7e+03)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (-1.8e+03,1.0e+03)
+DEAL::0 1 (-2.2e+03,1.4e+03)
+DEAL::0 2 (-2.6e+03,1.9e+03)
+DEAL::0 3 (-3.0e+03,2.3e+03)
+DEAL::0 4 (-3.4e+03,2.7e+03)
+DEAL::1 0 (-2.3e+03,1.7e+03)
+DEAL::1 1 (-2.8e+03,2.3e+03)
+DEAL::1 2 (-3.3e+03,2.9e+03)
+DEAL::1 3 (-3.8e+03,3.5e+03)
+DEAL::1 4 (-4.3e+03,4.2e+03)
+DEAL::2 0 (-2.8e+03,2.3e+03)
+DEAL::2 1 (-3.4e+03,3.2e+03)
+DEAL::2 2 (-4.0e+03,4.0e+03)
+DEAL::2 3 (-4.6e+03,4.8e+03)
+DEAL::2 4 (-5.1e+03,5.6e+03)
+DEAL::3 0 (-3.3e+03,3.0e+03)
+DEAL::3 1 (-4.0e+03,4.0e+03)
+DEAL::3 2 (-4.6e+03,5.0e+03)
+DEAL::3 3 (-5.3e+03,6.1e+03)
+DEAL::3 4 (-6.0e+03,7.1e+03)
+DEAL::4 0 (-3.8e+03,3.6e+03)
+DEAL::4 1 (-4.6e+03,4.8e+03)
+DEAL::4 2 (-5.3e+03,6.1e+03)
+DEAL::4 3 (-6.1e+03,7.4e+03)
+DEAL::4 4 (-6.9e+03,8.6e+03)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_44/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m, n, o;
+ make_complex_matrix (m);
+ make_complex_matrix (n);
+ make_complex_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,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (-1.9e+03,1.4e+03)
+DEAL::0 1 (-2.4e+03,1.9e+03)
+DEAL::0 2 (-2.8e+03,2.4e+03)
+DEAL::0 3 (-3.2e+03,2.9e+03)
+DEAL::0 4 (-3.7e+03,3.4e+03)
+DEAL::1 0 (-2.4e+03,1.9e+03)
+DEAL::1 1 (-2.9e+03,2.6e+03)
+DEAL::1 2 (-3.4e+03,3.3e+03)
+DEAL::1 3 (-3.9e+03,3.9e+03)
+DEAL::1 4 (-4.4e+03,4.6e+03)
+DEAL::2 0 (-2.8e+03,2.4e+03)
+DEAL::2 1 (-3.4e+03,3.3e+03)
+DEAL::2 2 (-4.0e+03,4.1e+03)
+DEAL::2 3 (-4.6e+03,5.0e+03)
+DEAL::2 4 (-5.1e+03,5.8e+03)
+DEAL::3 0 (-3.2e+03,2.9e+03)
+DEAL::3 1 (-3.9e+03,3.9e+03)
+DEAL::3 2 (-4.6e+03,5.0e+03)
+DEAL::3 3 (-5.2e+03,6.0e+03)
+DEAL::3 4 (-5.9e+03,7.0e+03)
+DEAL::4 0 (-3.6e+03,3.4e+03)
+DEAL::4 1 (-4.4e+03,4.6e+03)
+DEAL::4 2 (-5.1e+03,5.8e+03)
+DEAL::4 3 (-5.9e+03,7.0e+03)
+DEAL::4 4 (-6.6e+03,8.2e+03)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (-1.9e+03,1.4e+03)
+DEAL::0 1 (-2.4e+03,1.9e+03)
+DEAL::0 2 (-2.8e+03,2.4e+03)
+DEAL::0 3 (-3.2e+03,2.9e+03)
+DEAL::0 4 (-3.7e+03,3.4e+03)
+DEAL::1 0 (-2.4e+03,1.9e+03)
+DEAL::1 1 (-2.9e+03,2.6e+03)
+DEAL::1 2 (-3.4e+03,3.2e+03)
+DEAL::1 3 (-3.9e+03,3.9e+03)
+DEAL::1 4 (-4.4e+03,4.6e+03)
+DEAL::2 0 (-2.8e+03,2.4e+03)
+DEAL::2 1 (-3.4e+03,3.2e+03)
+DEAL::2 2 (-4.0e+03,4.1e+03)
+DEAL::2 3 (-4.6e+03,4.9e+03)
+DEAL::2 4 (-5.2e+03,5.8e+03)
+DEAL::3 0 (-3.2e+03,2.9e+03)
+DEAL::3 1 (-3.9e+03,3.9e+03)
+DEAL::3 2 (-4.6e+03,4.9e+03)
+DEAL::3 3 (-5.2e+03,6.0e+03)
+DEAL::3 4 (-5.9e+03,7.0e+03)
+DEAL::4 0 (-3.7e+03,3.4e+03)
+DEAL::4 1 (-4.4e+03,4.6e+03)
+DEAL::4 2 (-5.2e+03,5.8e+03)
+DEAL::4 3 (-5.9e+03,7.0e+03)
+DEAL::4 4 (-6.7e+03,8.2e+03)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (-1.9e+03,1.4e+03)
+DEAL::0 1 (-2.4e+03,1.9e+03)
+DEAL::0 2 (-2.8e+03,2.4e+03)
+DEAL::0 3 (-3.2e+03,2.9e+03)
+DEAL::0 4 (-3.7e+03,3.4e+03)
+DEAL::1 0 (-2.4e+03,1.9e+03)
+DEAL::1 1 (-2.9e+03,2.6e+03)
+DEAL::1 2 (-3.4e+03,3.3e+03)
+DEAL::1 3 (-3.9e+03,3.9e+03)
+DEAL::1 4 (-4.4e+03,4.6e+03)
+DEAL::2 0 (-2.8e+03,2.4e+03)
+DEAL::2 1 (-3.4e+03,3.3e+03)
+DEAL::2 2 (-4.0e+03,4.1e+03)
+DEAL::2 3 (-4.6e+03,5.0e+03)
+DEAL::2 4 (-5.1e+03,5.8e+03)
+DEAL::3 0 (-3.2e+03,2.9e+03)
+DEAL::3 1 (-3.9e+03,3.9e+03)
+DEAL::3 2 (-4.6e+03,5.0e+03)
+DEAL::3 3 (-5.2e+03,6.0e+03)
+DEAL::3 4 (-5.9e+03,7.0e+03)
+DEAL::4 0 (-3.6e+03,3.4e+03)
+DEAL::4 1 (-4.4e+03,4.6e+03)
+DEAL::4 2 (-5.1e+03,5.8e+03)
+DEAL::4 3 (-5.9e+03,7.0e+03)
+DEAL::4 4 (-6.6e+03,8.2e+03)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
+DEAL::0 0 (-1.9e+03,1.4e+03)
+DEAL::0 1 (-2.4e+03,1.9e+03)
+DEAL::0 2 (-2.8e+03,2.4e+03)
+DEAL::0 3 (-3.2e+03,2.9e+03)
+DEAL::0 4 (-3.7e+03,3.4e+03)
+DEAL::1 0 (-2.4e+03,1.9e+03)
+DEAL::1 1 (-2.9e+03,2.6e+03)
+DEAL::1 2 (-3.4e+03,3.2e+03)
+DEAL::1 3 (-3.9e+03,3.9e+03)
+DEAL::1 4 (-4.4e+03,4.6e+03)
+DEAL::2 0 (-2.8e+03,2.4e+03)
+DEAL::2 1 (-3.4e+03,3.2e+03)
+DEAL::2 2 (-4.0e+03,4.1e+03)
+DEAL::2 3 (-4.6e+03,4.9e+03)
+DEAL::2 4 (-5.2e+03,5.8e+03)
+DEAL::3 0 (-3.2e+03,2.9e+03)
+DEAL::3 1 (-3.9e+03,3.9e+03)
+DEAL::3 2 (-4.6e+03,4.9e+03)
+DEAL::3 3 (-5.2e+03,6.0e+03)
+DEAL::3 4 (-5.9e+03,7.0e+03)
+DEAL::4 0 (-3.7e+03,3.4e+03)
+DEAL::4 1 (-4.4e+03,4.6e+03)
+DEAL::4 2 (-5.2e+03,5.8e+03)
+DEAL::4 3 (-5.9e+03,7.0e+03)
+DEAL::4 4 (-6.7e+03,8.2e+03)
+DEAL::0 0 (2.0,12.)
+DEAL::0 1 (3.0,15.)
+DEAL::0 2 (4.0,18.)
+DEAL::0 3 (5.0,21.)
+DEAL::0 4 (6.0,24.)
+DEAL::1 0 (4.0,16.)
+DEAL::1 1 (6.0,20.)
+DEAL::1 2 (8.0,24.)
+DEAL::1 3 (10.,28.)
+DEAL::1 4 (12.,32.)
+DEAL::2 0 (6.0,20.)
+DEAL::2 1 (9.0,25.)
+DEAL::2 2 (12.,30.)
+DEAL::2 3 (15.,35.)
+DEAL::2 4 (18.,40.)
+DEAL::3 0 (8.0,24.)
+DEAL::3 1 (12.,30.)
+DEAL::3 2 (16.,36.)
+DEAL::3 3 (20.,42.)
+DEAL::3 4 (24.,48.)
+DEAL::4 0 (10.,28.)
+DEAL::4 1 (15.,35.)
+DEAL::4 2 (20.,42.)
+DEAL::4 3 (25.,49.)
+DEAL::4 4 (30.,56.)
--- /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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_45/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ Vector<std::complex<number> > v, w;
+ make_complex_vector (v);
+ make_complex_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 (-4.1e+02,4.1e+02)
+DEAL::1 (-5.0e+02,6.2e+02)
+DEAL::2 (-5.9e+02,8.4e+02)
+DEAL::3 (-6.8e+02,1.0e+03)
+DEAL::4 (-7.6e+02,1.3e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (-4.1e+02,4.1e+02)
+DEAL::1 (-5.0e+02,6.2e+02)
+DEAL::2 (-5.9e+02,8.3e+02)
+DEAL::3 (-6.8e+02,1.0e+03)
+DEAL::4 (-7.7e+02,1.2e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (-4.1e+02,4.1e+02)
+DEAL::1 (-5.0e+02,6.2e+02)
+DEAL::2 (-5.9e+02,8.4e+02)
+DEAL::3 (-6.8e+02,1.0e+03)
+DEAL::4 (-7.6e+02,1.3e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (-4.1e+02,4.1e+02)
+DEAL::1 (-5.0e+02,6.2e+02)
+DEAL::2 (-5.9e+02,8.3e+02)
+DEAL::3 (-6.8e+02,1.0e+03)
+DEAL::4 (-7.7e+02,1.2e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_46/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ Vector<std::complex<number> > v, w;
+ make_complex_vector (v);
+ make_complex_vector (w);
+
+ m.vmult_add (v, w);
+ print_vector (v);
+ print_vector (w);
+}
+
--- /dev/null
+
+DEAL::0 (-4.1e+02,4.1e+02)
+DEAL::1 (-5.0e+02,6.2e+02)
+DEAL::2 (-5.9e+02,8.4e+02)
+DEAL::3 (-6.8e+02,1.0e+03)
+DEAL::4 (-7.6e+02,1.3e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (-4.1e+02,4.1e+02)
+DEAL::1 (-5.0e+02,6.2e+02)
+DEAL::2 (-5.9e+02,8.4e+02)
+DEAL::3 (-6.8e+02,1.0e+03)
+DEAL::4 (-7.6e+02,1.3e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_47/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ Vector<std::complex<number> > v, w;
+ make_complex_vector (v);
+ make_complex_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 (-4.3e+02,5.1e+02)
+DEAL::1 (-5.1e+02,6.8e+02)
+DEAL::2 (-5.9e+02,8.6e+02)
+DEAL::3 (-6.7e+02,1.0e+03)
+DEAL::4 (-7.4e+02,1.2e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (-4.3e+02,5.1e+02)
+DEAL::1 (-5.1e+02,6.8e+02)
+DEAL::2 (-5.9e+02,8.5e+02)
+DEAL::3 (-6.7e+02,1.0e+03)
+DEAL::4 (-7.5e+02,1.2e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (-4.3e+02,5.1e+02)
+DEAL::1 (-5.1e+02,6.8e+02)
+DEAL::2 (-5.9e+02,8.6e+02)
+DEAL::3 (-6.7e+02,1.0e+03)
+DEAL::4 (-7.4e+02,1.2e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (-4.3e+02,5.1e+02)
+DEAL::1 (-5.1e+02,6.8e+02)
+DEAL::2 (-5.9e+02,8.5e+02)
+DEAL::3 (-6.7e+02,1.0e+03)
+DEAL::4 (-7.5e+02,1.2e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_48/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ Vector<std::complex<number> > v, w;
+ make_complex_vector (v);
+ make_complex_vector (w);
+
+ m.Tvmult_add (v, w);
+ print_vector (v);
+ print_vector (w);
+}
+
--- /dev/null
+
+DEAL::0 (-4.3e+02,5.1e+02)
+DEAL::1 (-5.1e+02,6.8e+02)
+DEAL::2 (-5.9e+02,8.6e+02)
+DEAL::3 (-6.7e+02,1.0e+03)
+DEAL::4 (-7.4e+02,1.2e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (-4.3e+02,5.1e+02)
+DEAL::1 (-5.1e+02,6.8e+02)
+DEAL::2 (-5.9e+02,8.6e+02)
+DEAL::3 (-6.7e+02,1.0e+03)
+DEAL::4 (-7.4e+02,1.2e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_49/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ Vector<std::complex<number> > v, w;
+ make_complex_vector (v);
+ make_complex_vector (w);
+
+ m.precondition_Jacobi (v, w, 3.141);
+ print_vector (v);
+ print_vector (w);
+}
+
--- /dev/null
+
+DEAL::0 (0.81,-0.13)
+DEAL::1 (0.66,-0.12)
+DEAL::2 (0.56,-0.090)
+DEAL::3 (0.48,-0.070)
+DEAL::4 (0.42,-0.054)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (0.81,-0.13)
+DEAL::1 (0.66,-0.12)
+DEAL::2 (0.56,-0.090)
+DEAL::3 (0.48,-0.070)
+DEAL::4 (0.42,-0.054)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_50/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ Vector<std::complex<number> > v, w, x;
+ make_complex_vector (v);
+ make_complex_vector (w);
+ make_complex_vector (x);
+
+ m.residual (v, w, x);
+ print_vector (v);
+ print_vector (w);
+ print_vector (x);
+}
+
--- /dev/null
+
+DEAL::0 (4.1e+02,-4.1e+02)
+DEAL::1 (5.0e+02,-6.2e+02)
+DEAL::2 (5.9e+02,-8.2e+02)
+DEAL::3 (6.8e+02,-1.0e+03)
+DEAL::4 (7.8e+02,-1.2e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (4.1e+02,-4.1e+02)
+DEAL::1 (5.0e+02,-6.2e+02)
+DEAL::2 (5.9e+02,-8.2e+02)
+DEAL::3 (6.8e+02,-1.0e+03)
+DEAL::4 (7.8e+02,-1.2e+03)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_51/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ Vector<std::complex<number> > x, y;
+ make_complex_vector (x);
+ make_complex_vector (y);
+
+ m.forward (x, y);
+ print_vector (x);
+ print_vector (y);
+}
+
--- /dev/null
+
+DEAL::0 (0.26,-0.041)
+DEAL::1 (0.0069,-0.014)
+DEAL::2 (0.0045,0.0060)
+DEAL::3 (0.0028,0.0029)
+DEAL::4 (0.0018,0.0016)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (0.26,-0.041)
+DEAL::1 (0.0069,-0.014)
+DEAL::2 (0.0045,0.0060)
+DEAL::3 (0.0028,0.0029)
+DEAL::4 (0.0018,0.0016)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.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; this time we actually store complex values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_complex_full_matrix_52/output";
+
+
+template <typename number>
+void
+check ()
+{
+ FullMatrix<std::complex<number> > m;
+ make_complex_matrix (m);
+ Vector<std::complex<number> > x, y;
+ make_complex_vector (x);
+ make_complex_vector (y);
+
+ m.backward (x, y);
+ print_vector (x);
+ print_vector (y);
+}
+
--- /dev/null
+
+DEAL::0 (0.0049,0.0082)
+DEAL::1 (0.0029,0.0033)
+DEAL::2 (0.0017,0.0015)
+DEAL::3 (0.0011,0.00080)
+DEAL::4 (0.13,-0.017)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)
+DEAL::0 (0.0049,0.0082)
+DEAL::1 (0.0029,0.0033)
+DEAL::2 (0.0017,0.0015)
+DEAL::3 (0.0011,0.00080)
+DEAL::4 (0.13,-0.017)
+DEAL::0 (1.0,3.0)
+DEAL::1 (2.0,4.0)
+DEAL::2 (3.0,5.0)
+DEAL::3 (4.0,6.0)
+DEAL::4 (5.0,7.0)