]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add a set of tests where we store real values in complex storage.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Sun, 18 Nov 2007 23:47:02 +0000 (23:47 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Sun, 18 Nov 2007 23:47:02 +0000 (23:47 +0000)
git-svn-id: https://svn.dealii.org/trunk@15519 0785d39b-7218-0410-832d-ea1e28bc413d

120 files changed:
tests/bits/Makefile
tests/bits/complex_real_full_matrix_01.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_01/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_02.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_02/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_03.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_03/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_04.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_04/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_05.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_05/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_06.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_06/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_07.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_07/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_08.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_08/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_09.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_09/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_10.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_10/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_11.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_11/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_12.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_12/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_13.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_13/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_14.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_14/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_15.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_15/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_16.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_16/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_17.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_17/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_18.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_18/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_19.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_19/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_20.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_20/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_21.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_21/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_22.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_22/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_23.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_23/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_24.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_24/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_25.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_25/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_26.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_26/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_27.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_27/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_28.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_28/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_29.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_29/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_30.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_30/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_31.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_31/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_32.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_32/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_33.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_33/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_34.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_34/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_35.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_35/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_36.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_36/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_37.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_37/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_38.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_38/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_39.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_39/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_40.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_40/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_41.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_41/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_42.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_42/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_43.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_43/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_44.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_44/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_45.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_45/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_46.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_46/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_47.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_47/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_48.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_48/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_49.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_49/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_50.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_50/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_51.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_51/cmp/generic [new file with mode: 0644]
tests/bits/complex_real_full_matrix_52.cc [new file with mode: 0644]
tests/bits/complex_real_full_matrix_52/cmp/generic [new file with mode: 0644]
tests/bits/full_matrix_13.cc
tests/bits/full_matrix_14.cc
tests/bits/full_matrix_45.cc
tests/bits/full_matrix_46.cc
tests/bits/full_matrix_47.cc
tests/bits/full_matrix_48.cc
tests/bits/full_matrix_49.cc
tests/bits/full_matrix_50.cc
tests/bits/full_matrix_51.cc
tests/bits/full_matrix_52.cc
tests/bits/full_matrix_vector_01.cc
tests/bits/full_matrix_vector_02.cc
tests/bits/full_matrix_vector_05.cc
tests/bits/full_matrix_vector_06.cc
tests/bits/full_matrix_vector_07.cc

index 12eb6720e2c841d269ec5b36185f9593c1d17b5d..6e04691d358c9b7362587785ef4f0b3ac067be6b 100644 (file)
@@ -25,6 +25,7 @@ tests_x = grid_generator_?? \
          grid_in_* \
          point_inside_? \
          full_matrix_* \
+         complex_real_full_matrix_* \
          block_matrix_array_* \
          data_out_0* \
           data_out_faces_0* \
diff --git a/tests/bits/complex_real_full_matrix_01.cc b/tests/bits/complex_real_full_matrix_01.cc
new file mode 100644 (file)
index 0000000..e3bf176
--- /dev/null
@@ -0,0 +1,35 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check creation and output of a matrix. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_01/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_01/cmp/generic b/tests/bits/complex_real_full_matrix_01/cmp/generic
new file mode 100644 (file)
index 0000000..9decec3
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_02.cc b/tests/bits/complex_real_full_matrix_02.cc
new file mode 100644 (file)
index 0000000..a93966c
--- /dev/null
@@ -0,0 +1,40 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check running over const iterators. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_02/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+
+  for (typename FullMatrix<std::complex<number> >::const_iterator
+        p = m.begin(); p!=m.end(); ++p)
+    deallog << p->row() << ' ' << p->column() << ' '
+           << p->value()
+           << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_02/cmp/generic b/tests/bits/complex_real_full_matrix_02/cmp/generic
new file mode 100644 (file)
index 0000000..9decec3
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_03.cc b/tests/bits/complex_real_full_matrix_03.cc
new file mode 100644 (file)
index 0000000..f0be1ff
--- /dev/null
@@ -0,0 +1,41 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check running over const iterators starting at the second line. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_03/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+
+  
+  for (typename FullMatrix<std::complex<number> >::const_iterator
+        p = m.begin(1); p!=m.end(1); ++p)
+    deallog << p->row() << ' ' << p->column() << ' '
+           << p->value()
+           << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_03/cmp/generic b/tests/bits/complex_real_full_matrix_03/cmp/generic
new file mode 100644 (file)
index 0000000..1234139
--- /dev/null
@@ -0,0 +1,11 @@
+
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_04.cc b/tests/bits/complex_real_full_matrix_04.cc
new file mode 100644 (file)
index 0000000..230adb8
--- /dev/null
@@ -0,0 +1,41 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check creation and output of a matrix using an array. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_04/output";
+
+
+template <typename number>
+void
+check ()
+{
+  const std::complex<number> array[] = { 1, 2, 3, 4, 5,
+                          6, 7, 8, 9, 0,
+                          1, 2, 3, 4, 5,
+                          6, 7, 8, 9, 0,
+                          1, 2, 3, 4, 5 };
+  
+  FullMatrix<std::complex<number> > m (5,5,array);
+
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_04/cmp/generic b/tests/bits/complex_real_full_matrix_04/cmp/generic
new file mode 100644 (file)
index 0000000..a0e80cc
--- /dev/null
@@ -0,0 +1,51 @@
+
+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)
diff --git a/tests/bits/complex_real_full_matrix_05.cc b/tests/bits/complex_real_full_matrix_05.cc
new file mode 100644 (file)
index 0000000..2d19332
--- /dev/null
@@ -0,0 +1,34 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check creation using IdentityMatrix. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_05/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m (IdentityMatrix(5));
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_05/cmp/generic b/tests/bits/complex_real_full_matrix_05/cmp/generic
new file mode 100644 (file)
index 0000000..e804e89
--- /dev/null
@@ -0,0 +1,51 @@
+
+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)
diff --git a/tests/bits/complex_real_full_matrix_06.cc b/tests/bits/complex_real_full_matrix_06.cc
new file mode 100644 (file)
index 0000000..d82a0e0
--- /dev/null
@@ -0,0 +1,39 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check matrix assignment. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_06/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+
+  FullMatrix<std::complex<number> > n;
+  n = m;
+  
+  print_matrix (n);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_06/cmp/generic b/tests/bits/complex_real_full_matrix_06/cmp/generic
new file mode 100644 (file)
index 0000000..9decec3
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_07.cc b/tests/bits/complex_real_full_matrix_07.cc
new file mode 100644 (file)
index 0000000..8db949a
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check matrix=0. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_07/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m = std::complex<number>(0);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_07/cmp/generic b/tests/bits/complex_real_full_matrix_07/cmp/generic
new file mode 100644 (file)
index 0000000..5b6383b
--- /dev/null
@@ -0,0 +1,51 @@
+
+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)
diff --git a/tests/bits/complex_real_full_matrix_08.cc b/tests/bits/complex_real_full_matrix_08.cc
new file mode 100644 (file)
index 0000000..3b5cd4c
--- /dev/null
@@ -0,0 +1,39 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::copy_from. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_08/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+
+  FullMatrix<std::complex<number> > n;
+  n.copy_from (m);
+  
+  print_matrix (n);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_08/cmp/generic b/tests/bits/complex_real_full_matrix_08/cmp/generic
new file mode 100644 (file)
index 0000000..9decec3
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_09.cc b/tests/bits/complex_real_full_matrix_09.cc
new file mode 100644 (file)
index 0000000..fb1a7bb
--- /dev/null
@@ -0,0 +1,35 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check assignment from IdenityMatrix. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_09/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  m = IdentityMatrix(5);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_09/cmp/generic b/tests/bits/complex_real_full_matrix_09/cmp/generic
new file mode 100644 (file)
index 0000000..e804e89
--- /dev/null
@@ -0,0 +1,51 @@
+
+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)
diff --git a/tests/bits/complex_real_full_matrix_10.cc b/tests/bits/complex_real_full_matrix_10.cc
new file mode 100644 (file)
index 0000000..1e2df74
--- /dev/null
@@ -0,0 +1,38 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::fill. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_10/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+
+  FullMatrix<std::complex<number> > n (10, 10);
+  n.fill (m, 3, 4);
+  print_matrix (n);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_10/cmp/generic b/tests/bits/complex_real_full_matrix_10/cmp/generic
new file mode 100644 (file)
index 0000000..ec8fa95
--- /dev/null
@@ -0,0 +1,201 @@
+
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (2.0,0.0)
+DEAL::3 5 (3.0,0.0)
+DEAL::3 6 (4.0,0.0)
+DEAL::3 7 (5.0,0.0)
+DEAL::3 8 (6.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (4.0,0.0)
+DEAL::4 5 (6.0,0.0)
+DEAL::4 6 (8.0,0.0)
+DEAL::4 7 (10.,0.0)
+DEAL::4 8 (12.,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (6.0,0.0)
+DEAL::5 5 (9.0,0.0)
+DEAL::5 6 (12.,0.0)
+DEAL::5 7 (15.,0.0)
+DEAL::5 8 (18.,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (8.0,0.0)
+DEAL::6 5 (12.,0.0)
+DEAL::6 6 (16.,0.0)
+DEAL::6 7 (20.,0.0)
+DEAL::6 8 (24.,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (10.,0.0)
+DEAL::7 5 (15.,0.0)
+DEAL::7 6 (20.,0.0)
+DEAL::7 7 (25.,0.0)
+DEAL::7 8 (30.,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,0.0)
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (2.0,0.0)
+DEAL::3 5 (3.0,0.0)
+DEAL::3 6 (4.0,0.0)
+DEAL::3 7 (5.0,0.0)
+DEAL::3 8 (6.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (4.0,0.0)
+DEAL::4 5 (6.0,0.0)
+DEAL::4 6 (8.0,0.0)
+DEAL::4 7 (10.,0.0)
+DEAL::4 8 (12.,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (6.0,0.0)
+DEAL::5 5 (9.0,0.0)
+DEAL::5 6 (12.,0.0)
+DEAL::5 7 (15.,0.0)
+DEAL::5 8 (18.,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (8.0,0.0)
+DEAL::6 5 (12.,0.0)
+DEAL::6 6 (16.,0.0)
+DEAL::6 7 (20.,0.0)
+DEAL::6 8 (24.,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (10.,0.0)
+DEAL::7 5 (15.,0.0)
+DEAL::7 6 (20.,0.0)
+DEAL::7 7 (25.,0.0)
+DEAL::7 8 (30.,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,0.0)
diff --git a/tests/bits/complex_real_full_matrix_11.cc b/tests/bits/complex_real_full_matrix_11.cc
new file mode 100644 (file)
index 0000000..bff2bbc
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check creation and output of a matrix. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_11/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m, n;
+  make_matrix (m);
+  make_matrix (n);
+  deallog << (m==n) << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_11/cmp/generic b/tests/bits/complex_real_full_matrix_11/cmp/generic
new file mode 100644 (file)
index 0000000..e46c0ac
--- /dev/null
@@ -0,0 +1,3 @@
+
+DEAL::1
+DEAL::1
diff --git a/tests/bits/complex_real_full_matrix_12.cc b/tests/bits/complex_real_full_matrix_12.cc
new file mode 100644 (file)
index 0000000..f2489e5
--- /dev/null
@@ -0,0 +1,37 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::all_zero. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_12/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m(3,3);
+  deallog << m.all_zero() << std::endl;
+  
+  make_matrix (m);
+  deallog << m.all_zero() << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_12/cmp/generic b/tests/bits/complex_real_full_matrix_12/cmp/generic
new file mode 100644 (file)
index 0000000..de01f2c
--- /dev/null
@@ -0,0 +1,5 @@
+
+DEAL::1
+DEAL::0
+DEAL::1
+DEAL::0
diff --git a/tests/bits/complex_real_full_matrix_13.cc b/tests/bits/complex_real_full_matrix_13.cc
new file mode 100644 (file)
index 0000000..186c792
--- /dev/null
@@ -0,0 +1,38 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::matrix_norm_square. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_13/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  Vector<std::complex<number> > v;
+  make_vector (v);
+  
+  deallog << m.matrix_norm_square (v) << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_13/cmp/generic b/tests/bits/complex_real_full_matrix_13/cmp/generic
new file mode 100644 (file)
index 0000000..4a96c33
--- /dev/null
@@ -0,0 +1,3 @@
+
+DEAL::(3.8e+03,0.0)
+DEAL::(3.8e+03,0.0)
diff --git a/tests/bits/complex_real_full_matrix_14.cc b/tests/bits/complex_real_full_matrix_14.cc
new file mode 100644 (file)
index 0000000..3ba6af4
--- /dev/null
@@ -0,0 +1,41 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::matrix_scalar_product. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_14/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  Vector<std::complex<number> > v,w;
+  make_vector (v);
+  make_vector (w);
+  for (unsigned int i=0; i<w.size(); ++i)
+    w(i) = w(i)+std::complex<number>(1.);
+  
+  deallog << m.matrix_scalar_product (v,w) << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_14/cmp/generic b/tests/bits/complex_real_full_matrix_14/cmp/generic
new file mode 100644 (file)
index 0000000..ba1c954
--- /dev/null
@@ -0,0 +1,3 @@
+
+DEAL::(5.0e+03,0.0)
+DEAL::(5.0e+03,0.0)
diff --git a/tests/bits/complex_real_full_matrix_15.cc b/tests/bits/complex_real_full_matrix_15.cc
new file mode 100644 (file)
index 0000000..29311f8
--- /dev/null
@@ -0,0 +1,35 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::l1_norm. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_15/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  deallog << m.l1_norm() << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_15/cmp/generic b/tests/bits/complex_real_full_matrix_15/cmp/generic
new file mode 100644 (file)
index 0000000..5219ffb
--- /dev/null
@@ -0,0 +1,3 @@
+
+DEAL::90.
+DEAL::90.
diff --git a/tests/bits/complex_real_full_matrix_16.cc b/tests/bits/complex_real_full_matrix_16.cc
new file mode 100644 (file)
index 0000000..7f17d77
--- /dev/null
@@ -0,0 +1,35 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::linfty_norm. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_16/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  deallog << m.linfty_norm() << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_16/cmp/generic b/tests/bits/complex_real_full_matrix_16/cmp/generic
new file mode 100644 (file)
index 0000000..bf7e83e
--- /dev/null
@@ -0,0 +1,3 @@
+
+DEAL::1.0e+02
+DEAL::1.0e+02
diff --git a/tests/bits/complex_real_full_matrix_17.cc b/tests/bits/complex_real_full_matrix_17.cc
new file mode 100644 (file)
index 0000000..eea4fe0
--- /dev/null
@@ -0,0 +1,35 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::frobenius_norm. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_17/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  deallog << m.frobenius_norm() << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_17/cmp/generic b/tests/bits/complex_real_full_matrix_17/cmp/generic
new file mode 100644 (file)
index 0000000..d106e77
--- /dev/null
@@ -0,0 +1,3 @@
+
+DEAL::70.
+DEAL::70.
diff --git a/tests/bits/complex_real_full_matrix_18.cc b/tests/bits/complex_real_full_matrix_18.cc
new file mode 100644 (file)
index 0000000..d97b492
--- /dev/null
@@ -0,0 +1,35 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::relative_symmetry_norm2. like the full_matrix_* tests,
+// but use complex-valued matrices and vectors, even though we only store real
+// values in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_18/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  deallog << m.relative_symmetry_norm2() << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_18/cmp/generic b/tests/bits/complex_real_full_matrix_18/cmp/generic
new file mode 100644 (file)
index 0000000..7a5d74c
--- /dev/null
@@ -0,0 +1,3 @@
+
+DEAL::0.14
+DEAL::0.14
diff --git a/tests/bits/complex_real_full_matrix_19.cc b/tests/bits/complex_real_full_matrix_19.cc
new file mode 100644 (file)
index 0000000..7aa86a7
--- /dev/null
@@ -0,0 +1,41 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::determinant. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_19/output";
+
+
+template <typename number>
+void
+check ()
+{
+  for (unsigned int n=1; n<=3; ++n)
+    {
+      const std::complex<number> array[] = { 1, 2, 3, 4, 5, 6, 7, 8, 9 };
+    
+      FullMatrix<std::complex<number> > m(n,n, array);
+      print_matrix (m);
+      deallog << m.determinant () << std::endl;
+    }
+}
+
+
diff --git a/tests/bits/complex_real_full_matrix_19/cmp/generic b/tests/bits/complex_real_full_matrix_19/cmp/generic
new file mode 100644 (file)
index 0000000..39f150c
--- /dev/null
@@ -0,0 +1,35 @@
+
+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)
diff --git a/tests/bits/complex_real_full_matrix_20.cc b/tests/bits/complex_real_full_matrix_20.cc
new file mode 100644 (file)
index 0000000..20c85cd
--- /dev/null
@@ -0,0 +1,35 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::trace. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_20/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  deallog << m.trace() << std::endl;
+}
+
diff --git a/tests/bits/complex_real_full_matrix_20/cmp/generic b/tests/bits/complex_real_full_matrix_20/cmp/generic
new file mode 100644 (file)
index 0000000..3ca88bc
--- /dev/null
@@ -0,0 +1,3 @@
+
+DEAL::(70.,0.0)
+DEAL::(70.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_21.cc b/tests/bits/complex_real_full_matrix_21.cc
new file mode 100644 (file)
index 0000000..b175fb1
--- /dev/null
@@ -0,0 +1,35 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::print. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_21/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m.print (deallog.get_file_stream());
+}
+
diff --git a/tests/bits/complex_real_full_matrix_21/cmp/generic b/tests/bits/complex_real_full_matrix_21/cmp/generic
new file mode 100644 (file)
index 0000000..8c4940c
--- /dev/null
@@ -0,0 +1,11 @@
+
+(2.0,0.0)(3.0,0.0)(4.0,0.0)(5.0,0.0)(6.0,0.0)
+(4.0,0.0)(6.0,0.0)(8.0,0.0)(10.,0.0)(12.,0.0)
+(6.0,0.0)(9.0,0.0)(12.,0.0)(15.,0.0)(18.,0.0)
+(8.0,0.0)(12.,0.0)(16.,0.0)(20.,0.0)(24.,0.0)
+(10.,0.0)(15.,0.0)(20.,0.0)(25.,0.0)(30.,0.0)
+(2.0,0.0)(3.0,0.0)(4.0,0.0)(5.0,0.0)(6.0,0.0)
+(4.0,0.0)(6.0,0.0)(8.0,0.0)(10.,0.0)(12.,0.0)
+(6.0,0.0)(9.0,0.0)(12.,0.0)(15.,0.0)(18.,0.0)
+(8.0,0.0)(12.,0.0)(16.,0.0)(20.,0.0)(24.,0.0)
+(10.,0.0)(15.,0.0)(20.,0.0)(25.,0.0)(30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_22.cc b/tests/bits/complex_real_full_matrix_22.cc
new file mode 100644 (file)
index 0000000..1c33d29
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::operator*=. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_22/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m *= 3.1415;
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_22/cmp/generic b/tests/bits/complex_real_full_matrix_22/cmp/generic
new file mode 100644 (file)
index 0000000..570d6dd
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (6.3,0.0)
+DEAL::0 1 (9.4,0.0)
+DEAL::0 2 (13.,0.0)
+DEAL::0 3 (16.,0.0)
+DEAL::0 4 (19.,0.0)
+DEAL::1 0 (13.,0.0)
+DEAL::1 1 (19.,0.0)
+DEAL::1 2 (25.,0.0)
+DEAL::1 3 (31.,0.0)
+DEAL::1 4 (38.,0.0)
+DEAL::2 0 (19.,0.0)
+DEAL::2 1 (28.,0.0)
+DEAL::2 2 (38.,0.0)
+DEAL::2 3 (47.,0.0)
+DEAL::2 4 (57.,0.0)
+DEAL::3 0 (25.,0.0)
+DEAL::3 1 (38.,0.0)
+DEAL::3 2 (50.,0.0)
+DEAL::3 3 (63.,0.0)
+DEAL::3 4 (75.,0.0)
+DEAL::4 0 (31.,0.0)
+DEAL::4 1 (47.,0.0)
+DEAL::4 2 (63.,0.0)
+DEAL::4 3 (79.,0.0)
+DEAL::4 4 (94.,0.0)
+DEAL::0 0 (6.3,0.0)
+DEAL::0 1 (9.4,0.0)
+DEAL::0 2 (13.,0.0)
+DEAL::0 3 (16.,0.0)
+DEAL::0 4 (19.,0.0)
+DEAL::1 0 (13.,0.0)
+DEAL::1 1 (19.,0.0)
+DEAL::1 2 (25.,0.0)
+DEAL::1 3 (31.,0.0)
+DEAL::1 4 (38.,0.0)
+DEAL::2 0 (19.,0.0)
+DEAL::2 1 (28.,0.0)
+DEAL::2 2 (38.,0.0)
+DEAL::2 3 (47.,0.0)
+DEAL::2 4 (57.,0.0)
+DEAL::3 0 (25.,0.0)
+DEAL::3 1 (38.,0.0)
+DEAL::3 2 (50.,0.0)
+DEAL::3 3 (63.,0.0)
+DEAL::3 4 (75.,0.0)
+DEAL::4 0 (31.,0.0)
+DEAL::4 1 (47.,0.0)
+DEAL::4 2 (63.,0.0)
+DEAL::4 3 (79.,0.0)
+DEAL::4 4 (94.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_23.cc b/tests/bits/complex_real_full_matrix_23.cc
new file mode 100644 (file)
index 0000000..1a1b26f
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::operator/=. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_23/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m /= 3.1415926;
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_23/cmp/generic b/tests/bits/complex_real_full_matrix_23/cmp/generic
new file mode 100644 (file)
index 0000000..f5c61ce
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (0.64,0.0)
+DEAL::0 1 (0.95,0.0)
+DEAL::0 2 (1.3,0.0)
+DEAL::0 3 (1.6,0.0)
+DEAL::0 4 (1.9,0.0)
+DEAL::1 0 (1.3,0.0)
+DEAL::1 1 (1.9,0.0)
+DEAL::1 2 (2.5,0.0)
+DEAL::1 3 (3.2,0.0)
+DEAL::1 4 (3.8,0.0)
+DEAL::2 0 (1.9,0.0)
+DEAL::2 1 (2.9,0.0)
+DEAL::2 2 (3.8,0.0)
+DEAL::2 3 (4.8,0.0)
+DEAL::2 4 (5.7,0.0)
+DEAL::3 0 (2.5,0.0)
+DEAL::3 1 (3.8,0.0)
+DEAL::3 2 (5.1,0.0)
+DEAL::3 3 (6.4,0.0)
+DEAL::3 4 (7.6,0.0)
+DEAL::4 0 (3.2,0.0)
+DEAL::4 1 (4.8,0.0)
+DEAL::4 2 (6.4,0.0)
+DEAL::4 3 (8.0,0.0)
+DEAL::4 4 (9.5,0.0)
+DEAL::0 0 (0.64,0.0)
+DEAL::0 1 (0.95,0.0)
+DEAL::0 2 (1.3,0.0)
+DEAL::0 3 (1.6,0.0)
+DEAL::0 4 (1.9,0.0)
+DEAL::1 0 (1.3,0.0)
+DEAL::1 1 (1.9,0.0)
+DEAL::1 2 (2.5,0.0)
+DEAL::1 3 (3.2,0.0)
+DEAL::1 4 (3.8,0.0)
+DEAL::2 0 (1.9,0.0)
+DEAL::2 1 (2.9,0.0)
+DEAL::2 2 (3.8,0.0)
+DEAL::2 3 (4.8,0.0)
+DEAL::2 4 (5.7,0.0)
+DEAL::3 0 (2.5,0.0)
+DEAL::3 1 (3.8,0.0)
+DEAL::3 2 (5.1,0.0)
+DEAL::3 3 (6.4,0.0)
+DEAL::3 4 (7.6,0.0)
+DEAL::4 0 (3.2,0.0)
+DEAL::4 1 (4.8,0.0)
+DEAL::4 2 (6.4,0.0)
+DEAL::4 3 (8.0,0.0)
+DEAL::4 4 (9.5,0.0)
diff --git a/tests/bits/complex_real_full_matrix_24.cc b/tests/bits/complex_real_full_matrix_24.cc
new file mode 100644 (file)
index 0000000..6e6b545
--- /dev/null
@@ -0,0 +1,38 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_24/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m,n;
+  make_matrix (m);
+  make_matrix (n);
+  m.add (3.14159, n);
+  
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_24/cmp/generic b/tests/bits/complex_real_full_matrix_24/cmp/generic
new file mode 100644 (file)
index 0000000..c538f0d
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (8.3,0.0)
+DEAL::0 1 (12.,0.0)
+DEAL::0 2 (17.,0.0)
+DEAL::0 3 (21.,0.0)
+DEAL::0 4 (25.,0.0)
+DEAL::1 0 (17.,0.0)
+DEAL::1 1 (25.,0.0)
+DEAL::1 2 (33.,0.0)
+DEAL::1 3 (41.,0.0)
+DEAL::1 4 (50.,0.0)
+DEAL::2 0 (25.,0.0)
+DEAL::2 1 (37.,0.0)
+DEAL::2 2 (50.,0.0)
+DEAL::2 3 (62.,0.0)
+DEAL::2 4 (75.,0.0)
+DEAL::3 0 (33.,0.0)
+DEAL::3 1 (50.,0.0)
+DEAL::3 2 (66.,0.0)
+DEAL::3 3 (83.,0.0)
+DEAL::3 4 (99.,0.0)
+DEAL::4 0 (41.,0.0)
+DEAL::4 1 (62.,0.0)
+DEAL::4 2 (83.,0.0)
+DEAL::4 3 (1.0e+02,0.0)
+DEAL::4 4 (1.2e+02,0.0)
+DEAL::0 0 (8.3,0.0)
+DEAL::0 1 (12.,0.0)
+DEAL::0 2 (17.,0.0)
+DEAL::0 3 (21.,0.0)
+DEAL::0 4 (25.,0.0)
+DEAL::1 0 (17.,0.0)
+DEAL::1 1 (25.,0.0)
+DEAL::1 2 (33.,0.0)
+DEAL::1 3 (41.,0.0)
+DEAL::1 4 (50.,0.0)
+DEAL::2 0 (25.,0.0)
+DEAL::2 1 (37.,0.0)
+DEAL::2 2 (50.,0.0)
+DEAL::2 3 (62.,0.0)
+DEAL::2 4 (75.,0.0)
+DEAL::3 0 (33.,0.0)
+DEAL::3 1 (50.,0.0)
+DEAL::3 2 (66.,0.0)
+DEAL::3 3 (83.,0.0)
+DEAL::3 4 (99.,0.0)
+DEAL::4 0 (41.,0.0)
+DEAL::4 1 (62.,0.0)
+DEAL::4 2 (83.,0.0)
+DEAL::4 3 (1.0e+02,0.0)
+DEAL::4 4 (1.2e+02,0.0)
diff --git a/tests/bits/complex_real_full_matrix_25.cc b/tests/bits/complex_real_full_matrix_25.cc
new file mode 100644 (file)
index 0000000..cf06c59
--- /dev/null
@@ -0,0 +1,38 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::add(4). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_25/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m,n,o;
+  make_matrix (m);
+  make_matrix (n);
+  make_matrix (o);
+  m.add (3.1415, n, 2.718, o);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_25/cmp/generic b/tests/bits/complex_real_full_matrix_25/cmp/generic
new file mode 100644 (file)
index 0000000..e14efd0
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (14.,0.0)
+DEAL::0 1 (21.,0.0)
+DEAL::0 2 (27.,0.0)
+DEAL::0 3 (34.,0.0)
+DEAL::0 4 (41.,0.0)
+DEAL::1 0 (27.,0.0)
+DEAL::1 1 (41.,0.0)
+DEAL::1 2 (55.,0.0)
+DEAL::1 3 (69.,0.0)
+DEAL::1 4 (82.,0.0)
+DEAL::2 0 (41.,0.0)
+DEAL::2 1 (62.,0.0)
+DEAL::2 2 (82.,0.0)
+DEAL::2 3 (1.0e+02,0.0)
+DEAL::2 4 (1.2e+02,0.0)
+DEAL::3 0 (55.,0.0)
+DEAL::3 1 (82.,0.0)
+DEAL::3 2 (1.1e+02,0.0)
+DEAL::3 3 (1.4e+02,0.0)
+DEAL::3 4 (1.6e+02,0.0)
+DEAL::4 0 (69.,0.0)
+DEAL::4 1 (1.0e+02,0.0)
+DEAL::4 2 (1.4e+02,0.0)
+DEAL::4 3 (1.7e+02,0.0)
+DEAL::4 4 (2.1e+02,0.0)
+DEAL::0 0 (14.,0.0)
+DEAL::0 1 (21.,0.0)
+DEAL::0 2 (27.,0.0)
+DEAL::0 3 (34.,0.0)
+DEAL::0 4 (41.,0.0)
+DEAL::1 0 (27.,0.0)
+DEAL::1 1 (41.,0.0)
+DEAL::1 2 (55.,0.0)
+DEAL::1 3 (69.,0.0)
+DEAL::1 4 (82.,0.0)
+DEAL::2 0 (41.,0.0)
+DEAL::2 1 (62.,0.0)
+DEAL::2 2 (82.,0.0)
+DEAL::2 3 (1.0e+02,0.0)
+DEAL::2 4 (1.2e+02,0.0)
+DEAL::3 0 (55.,0.0)
+DEAL::3 1 (82.,0.0)
+DEAL::3 2 (1.1e+02,0.0)
+DEAL::3 3 (1.4e+02,0.0)
+DEAL::3 4 (1.6e+02,0.0)
+DEAL::4 0 (69.,0.0)
+DEAL::4 1 (1.0e+02,0.0)
+DEAL::4 2 (1.4e+02,0.0)
+DEAL::4 3 (1.7e+02,0.0)
+DEAL::4 4 (2.1e+02,0.0)
diff --git a/tests/bits/complex_real_full_matrix_26.cc b/tests/bits/complex_real_full_matrix_26.cc
new file mode 100644 (file)
index 0000000..e560637
--- /dev/null
@@ -0,0 +1,39 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::add (6). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_26/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m,n,o,p;
+  make_matrix (m);
+  make_matrix (n);
+  make_matrix (o);
+  make_matrix (p);
+  m.add (3.1415, n, 2.718, o, 1.414, p);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_26/cmp/generic b/tests/bits/complex_real_full_matrix_26/cmp/generic
new file mode 100644 (file)
index 0000000..767c310
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (17.,0.0)
+DEAL::0 1 (25.,0.0)
+DEAL::0 2 (33.,0.0)
+DEAL::0 3 (41.,0.0)
+DEAL::0 4 (50.,0.0)
+DEAL::1 0 (33.,0.0)
+DEAL::1 1 (50.,0.0)
+DEAL::1 2 (66.,0.0)
+DEAL::1 3 (83.,0.0)
+DEAL::1 4 (99.,0.0)
+DEAL::2 0 (50.,0.0)
+DEAL::2 1 (74.,0.0)
+DEAL::2 2 (99.,0.0)
+DEAL::2 3 (1.2e+02,0.0)
+DEAL::2 4 (1.5e+02,0.0)
+DEAL::3 0 (66.,0.0)
+DEAL::3 1 (99.,0.0)
+DEAL::3 2 (1.3e+02,0.0)
+DEAL::3 3 (1.7e+02,0.0)
+DEAL::3 4 (2.0e+02,0.0)
+DEAL::4 0 (83.,0.0)
+DEAL::4 1 (1.2e+02,0.0)
+DEAL::4 2 (1.7e+02,0.0)
+DEAL::4 3 (2.1e+02,0.0)
+DEAL::4 4 (2.5e+02,0.0)
+DEAL::0 0 (17.,0.0)
+DEAL::0 1 (25.,0.0)
+DEAL::0 2 (33.,0.0)
+DEAL::0 3 (41.,0.0)
+DEAL::0 4 (50.,0.0)
+DEAL::1 0 (33.,0.0)
+DEAL::1 1 (50.,0.0)
+DEAL::1 2 (66.,0.0)
+DEAL::1 3 (83.,0.0)
+DEAL::1 4 (99.,0.0)
+DEAL::2 0 (50.,0.0)
+DEAL::2 1 (74.,0.0)
+DEAL::2 2 (99.,0.0)
+DEAL::2 3 (1.2e+02,0.0)
+DEAL::2 4 (1.5e+02,0.0)
+DEAL::3 0 (66.,0.0)
+DEAL::3 1 (99.,0.0)
+DEAL::3 2 (1.3e+02,0.0)
+DEAL::3 3 (1.7e+02,0.0)
+DEAL::3 4 (2.0e+02,0.0)
+DEAL::4 0 (83.,0.0)
+DEAL::4 1 (1.2e+02,0.0)
+DEAL::4 2 (1.7e+02,0.0)
+DEAL::4 3 (2.1e+02,0.0)
+DEAL::4 4 (2.5e+02,0.0)
diff --git a/tests/bits/complex_real_full_matrix_27.cc b/tests/bits/complex_real_full_matrix_27.cc
new file mode 100644 (file)
index 0000000..3203b45
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_27/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m(10,10),n;
+  make_matrix (n);
+  m.add (n, 2., 3, 4, 2, 3);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_27/cmp/generic b/tests/bits/complex_real_full_matrix_27/cmp/generic
new file mode 100644 (file)
index 0000000..af1e79e
--- /dev/null
@@ -0,0 +1,201 @@
+
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (30.,0.0)
+DEAL::3 5 (36.,0.0)
+DEAL::3 6 (0.0,0.0)
+DEAL::3 7 (0.0,0.0)
+DEAL::3 8 (0.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (40.,0.0)
+DEAL::4 5 (48.,0.0)
+DEAL::4 6 (0.0,0.0)
+DEAL::4 7 (0.0,0.0)
+DEAL::4 8 (0.0,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (50.,0.0)
+DEAL::5 5 (60.,0.0)
+DEAL::5 6 (0.0,0.0)
+DEAL::5 7 (0.0,0.0)
+DEAL::5 8 (0.0,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (0.0,0.0)
+DEAL::6 5 (0.0,0.0)
+DEAL::6 6 (0.0,0.0)
+DEAL::6 7 (0.0,0.0)
+DEAL::6 8 (0.0,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (0.0,0.0)
+DEAL::7 5 (0.0,0.0)
+DEAL::7 6 (0.0,0.0)
+DEAL::7 7 (0.0,0.0)
+DEAL::7 8 (0.0,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,0.0)
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (30.,0.0)
+DEAL::3 5 (36.,0.0)
+DEAL::3 6 (0.0,0.0)
+DEAL::3 7 (0.0,0.0)
+DEAL::3 8 (0.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (40.,0.0)
+DEAL::4 5 (48.,0.0)
+DEAL::4 6 (0.0,0.0)
+DEAL::4 7 (0.0,0.0)
+DEAL::4 8 (0.0,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (50.,0.0)
+DEAL::5 5 (60.,0.0)
+DEAL::5 6 (0.0,0.0)
+DEAL::5 7 (0.0,0.0)
+DEAL::5 8 (0.0,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (0.0,0.0)
+DEAL::6 5 (0.0,0.0)
+DEAL::6 6 (0.0,0.0)
+DEAL::6 7 (0.0,0.0)
+DEAL::6 8 (0.0,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (0.0,0.0)
+DEAL::7 5 (0.0,0.0)
+DEAL::7 6 (0.0,0.0)
+DEAL::7 7 (0.0,0.0)
+DEAL::7 8 (0.0,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,0.0)
diff --git a/tests/bits/complex_real_full_matrix_28.cc b/tests/bits/complex_real_full_matrix_28.cc
new file mode 100644 (file)
index 0000000..7be4a76
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::Tadd. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_28/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m(10,10),n;
+  make_matrix (n);
+  m.Tadd (n, 2., 3, 4, 2, 3);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_28/cmp/generic b/tests/bits/complex_real_full_matrix_28/cmp/generic
new file mode 100644 (file)
index 0000000..c964424
--- /dev/null
@@ -0,0 +1,201 @@
+
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (30.,0.0)
+DEAL::3 5 (40.,0.0)
+DEAL::3 6 (50.,0.0)
+DEAL::3 7 (0.0,0.0)
+DEAL::3 8 (0.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (36.,0.0)
+DEAL::4 5 (48.,0.0)
+DEAL::4 6 (60.,0.0)
+DEAL::4 7 (0.0,0.0)
+DEAL::4 8 (0.0,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (0.0,0.0)
+DEAL::5 5 (0.0,0.0)
+DEAL::5 6 (0.0,0.0)
+DEAL::5 7 (0.0,0.0)
+DEAL::5 8 (0.0,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (0.0,0.0)
+DEAL::6 5 (0.0,0.0)
+DEAL::6 6 (0.0,0.0)
+DEAL::6 7 (0.0,0.0)
+DEAL::6 8 (0.0,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (0.0,0.0)
+DEAL::7 5 (0.0,0.0)
+DEAL::7 6 (0.0,0.0)
+DEAL::7 7 (0.0,0.0)
+DEAL::7 8 (0.0,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,0.0)
+DEAL::0 0 (0.0,0.0)
+DEAL::0 1 (0.0,0.0)
+DEAL::0 2 (0.0,0.0)
+DEAL::0 3 (0.0,0.0)
+DEAL::0 4 (0.0,0.0)
+DEAL::0 5 (0.0,0.0)
+DEAL::0 6 (0.0,0.0)
+DEAL::0 7 (0.0,0.0)
+DEAL::0 8 (0.0,0.0)
+DEAL::0 9 (0.0,0.0)
+DEAL::1 0 (0.0,0.0)
+DEAL::1 1 (0.0,0.0)
+DEAL::1 2 (0.0,0.0)
+DEAL::1 3 (0.0,0.0)
+DEAL::1 4 (0.0,0.0)
+DEAL::1 5 (0.0,0.0)
+DEAL::1 6 (0.0,0.0)
+DEAL::1 7 (0.0,0.0)
+DEAL::1 8 (0.0,0.0)
+DEAL::1 9 (0.0,0.0)
+DEAL::2 0 (0.0,0.0)
+DEAL::2 1 (0.0,0.0)
+DEAL::2 2 (0.0,0.0)
+DEAL::2 3 (0.0,0.0)
+DEAL::2 4 (0.0,0.0)
+DEAL::2 5 (0.0,0.0)
+DEAL::2 6 (0.0,0.0)
+DEAL::2 7 (0.0,0.0)
+DEAL::2 8 (0.0,0.0)
+DEAL::2 9 (0.0,0.0)
+DEAL::3 0 (0.0,0.0)
+DEAL::3 1 (0.0,0.0)
+DEAL::3 2 (0.0,0.0)
+DEAL::3 3 (0.0,0.0)
+DEAL::3 4 (30.,0.0)
+DEAL::3 5 (40.,0.0)
+DEAL::3 6 (50.,0.0)
+DEAL::3 7 (0.0,0.0)
+DEAL::3 8 (0.0,0.0)
+DEAL::3 9 (0.0,0.0)
+DEAL::4 0 (0.0,0.0)
+DEAL::4 1 (0.0,0.0)
+DEAL::4 2 (0.0,0.0)
+DEAL::4 3 (0.0,0.0)
+DEAL::4 4 (36.,0.0)
+DEAL::4 5 (48.,0.0)
+DEAL::4 6 (60.,0.0)
+DEAL::4 7 (0.0,0.0)
+DEAL::4 8 (0.0,0.0)
+DEAL::4 9 (0.0,0.0)
+DEAL::5 0 (0.0,0.0)
+DEAL::5 1 (0.0,0.0)
+DEAL::5 2 (0.0,0.0)
+DEAL::5 3 (0.0,0.0)
+DEAL::5 4 (0.0,0.0)
+DEAL::5 5 (0.0,0.0)
+DEAL::5 6 (0.0,0.0)
+DEAL::5 7 (0.0,0.0)
+DEAL::5 8 (0.0,0.0)
+DEAL::5 9 (0.0,0.0)
+DEAL::6 0 (0.0,0.0)
+DEAL::6 1 (0.0,0.0)
+DEAL::6 2 (0.0,0.0)
+DEAL::6 3 (0.0,0.0)
+DEAL::6 4 (0.0,0.0)
+DEAL::6 5 (0.0,0.0)
+DEAL::6 6 (0.0,0.0)
+DEAL::6 7 (0.0,0.0)
+DEAL::6 8 (0.0,0.0)
+DEAL::6 9 (0.0,0.0)
+DEAL::7 0 (0.0,0.0)
+DEAL::7 1 (0.0,0.0)
+DEAL::7 2 (0.0,0.0)
+DEAL::7 3 (0.0,0.0)
+DEAL::7 4 (0.0,0.0)
+DEAL::7 5 (0.0,0.0)
+DEAL::7 6 (0.0,0.0)
+DEAL::7 7 (0.0,0.0)
+DEAL::7 8 (0.0,0.0)
+DEAL::7 9 (0.0,0.0)
+DEAL::8 0 (0.0,0.0)
+DEAL::8 1 (0.0,0.0)
+DEAL::8 2 (0.0,0.0)
+DEAL::8 3 (0.0,0.0)
+DEAL::8 4 (0.0,0.0)
+DEAL::8 5 (0.0,0.0)
+DEAL::8 6 (0.0,0.0)
+DEAL::8 7 (0.0,0.0)
+DEAL::8 8 (0.0,0.0)
+DEAL::8 9 (0.0,0.0)
+DEAL::9 0 (0.0,0.0)
+DEAL::9 1 (0.0,0.0)
+DEAL::9 2 (0.0,0.0)
+DEAL::9 3 (0.0,0.0)
+DEAL::9 4 (0.0,0.0)
+DEAL::9 5 (0.0,0.0)
+DEAL::9 6 (0.0,0.0)
+DEAL::9 7 (0.0,0.0)
+DEAL::9 8 (0.0,0.0)
+DEAL::9 9 (0.0,0.0)
diff --git a/tests/bits/complex_real_full_matrix_29.cc b/tests/bits/complex_real_full_matrix_29.cc
new file mode 100644 (file)
index 0000000..b6ebf2b
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::add_row. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_29/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m.add_row (2, 3.141, 4);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_29/cmp/generic b/tests/bits/complex_real_full_matrix_29/cmp/generic
new file mode 100644 (file)
index 0000000..c1fb921
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (37.,0.0)
+DEAL::2 1 (56.,0.0)
+DEAL::2 2 (75.,0.0)
+DEAL::2 3 (94.,0.0)
+DEAL::2 4 (1.1e+02,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (37.,0.0)
+DEAL::2 1 (56.,0.0)
+DEAL::2 2 (75.,0.0)
+DEAL::2 3 (94.,0.0)
+DEAL::2 4 (1.1e+02,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_30.cc b/tests/bits/complex_real_full_matrix_30.cc
new file mode 100644 (file)
index 0000000..67ae96e
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::add_row. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_30/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m.add_row (2, 3.141, 4, 2.7182, 3);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_30/cmp/generic b/tests/bits/complex_real_full_matrix_30/cmp/generic
new file mode 100644 (file)
index 0000000..87889b0
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (59.,0.0)
+DEAL::2 1 (89.,0.0)
+DEAL::2 2 (1.2e+02,0.0)
+DEAL::2 3 (1.5e+02,0.0)
+DEAL::2 4 (1.8e+02,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (59.,0.0)
+DEAL::2 1 (89.,0.0)
+DEAL::2 2 (1.2e+02,0.0)
+DEAL::2 3 (1.5e+02,0.0)
+DEAL::2 4 (1.8e+02,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_31.cc b/tests/bits/complex_real_full_matrix_31.cc
new file mode 100644 (file)
index 0000000..0b35d30
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::add_col. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_31/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m.add_col (2, 3.141, 4);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_31/cmp/generic b/tests/bits/complex_real_full_matrix_31/cmp/generic
new file mode 100644 (file)
index 0000000..cace4be
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (23.,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (46.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (69.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (91.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (1.1e+02,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (23.,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (46.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (69.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (91.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (1.1e+02,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_32.cc b/tests/bits/complex_real_full_matrix_32.cc
new file mode 100644 (file)
index 0000000..ecbf69d
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::add_col. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_32/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m.add_col (2, 3.141, 4, 2.718, 3);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_32/cmp/generic b/tests/bits/complex_real_full_matrix_32/cmp/generic
new file mode 100644 (file)
index 0000000..4015cee
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (36.,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (73.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (1.1e+02,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (1.5e+02,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (1.8e+02,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (36.,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (73.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (1.1e+02,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (1.5e+02,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (1.8e+02,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_33.cc b/tests/bits/complex_real_full_matrix_33.cc
new file mode 100644 (file)
index 0000000..e2cc745
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::swap_row. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_33/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m.swap_row (2, 4);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_33/cmp/generic b/tests/bits/complex_real_full_matrix_33/cmp/generic
new file mode 100644 (file)
index 0000000..feb3645
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (10.,0.0)
+DEAL::2 1 (15.,0.0)
+DEAL::2 2 (20.,0.0)
+DEAL::2 3 (25.,0.0)
+DEAL::2 4 (30.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (6.0,0.0)
+DEAL::4 1 (9.0,0.0)
+DEAL::4 2 (12.,0.0)
+DEAL::4 3 (15.,0.0)
+DEAL::4 4 (18.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (10.,0.0)
+DEAL::2 1 (15.,0.0)
+DEAL::2 2 (20.,0.0)
+DEAL::2 3 (25.,0.0)
+DEAL::2 4 (30.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (6.0,0.0)
+DEAL::4 1 (9.0,0.0)
+DEAL::4 2 (12.,0.0)
+DEAL::4 3 (15.,0.0)
+DEAL::4 4 (18.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_34.cc b/tests/bits/complex_real_full_matrix_34.cc
new file mode 100644 (file)
index 0000000..7f06486
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::swap_col. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_34/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m.swap_col (2, 4);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_34/cmp/generic b/tests/bits/complex_real_full_matrix_34/cmp/generic
new file mode 100644 (file)
index 0000000..a60126c
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (6.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (4.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (12.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (8.0,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (18.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (12.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (24.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (16.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (30.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (20.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (6.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (4.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (12.,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (8.0,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (18.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (12.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (24.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (16.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (30.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (20.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_35.cc b/tests/bits/complex_real_full_matrix_35.cc
new file mode 100644 (file)
index 0000000..3571889
--- /dev/null
@@ -0,0 +1,37 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::add_diag. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_35/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m, n;
+  make_matrix (m);
+  make_matrix (n);
+  m.add_diag (3.141, m);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_35/cmp/generic b/tests/bits/complex_real_full_matrix_35/cmp/generic
new file mode 100644 (file)
index 0000000..c538f0d
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (8.3,0.0)
+DEAL::0 1 (12.,0.0)
+DEAL::0 2 (17.,0.0)
+DEAL::0 3 (21.,0.0)
+DEAL::0 4 (25.,0.0)
+DEAL::1 0 (17.,0.0)
+DEAL::1 1 (25.,0.0)
+DEAL::1 2 (33.,0.0)
+DEAL::1 3 (41.,0.0)
+DEAL::1 4 (50.,0.0)
+DEAL::2 0 (25.,0.0)
+DEAL::2 1 (37.,0.0)
+DEAL::2 2 (50.,0.0)
+DEAL::2 3 (62.,0.0)
+DEAL::2 4 (75.,0.0)
+DEAL::3 0 (33.,0.0)
+DEAL::3 1 (50.,0.0)
+DEAL::3 2 (66.,0.0)
+DEAL::3 3 (83.,0.0)
+DEAL::3 4 (99.,0.0)
+DEAL::4 0 (41.,0.0)
+DEAL::4 1 (62.,0.0)
+DEAL::4 2 (83.,0.0)
+DEAL::4 3 (1.0e+02,0.0)
+DEAL::4 4 (1.2e+02,0.0)
+DEAL::0 0 (8.3,0.0)
+DEAL::0 1 (12.,0.0)
+DEAL::0 2 (17.,0.0)
+DEAL::0 3 (21.,0.0)
+DEAL::0 4 (25.,0.0)
+DEAL::1 0 (17.,0.0)
+DEAL::1 1 (25.,0.0)
+DEAL::1 2 (33.,0.0)
+DEAL::1 3 (41.,0.0)
+DEAL::1 4 (50.,0.0)
+DEAL::2 0 (25.,0.0)
+DEAL::2 1 (37.,0.0)
+DEAL::2 2 (50.,0.0)
+DEAL::2 3 (62.,0.0)
+DEAL::2 4 (75.,0.0)
+DEAL::3 0 (33.,0.0)
+DEAL::3 1 (50.,0.0)
+DEAL::3 2 (66.,0.0)
+DEAL::3 3 (83.,0.0)
+DEAL::3 4 (99.,0.0)
+DEAL::4 0 (41.,0.0)
+DEAL::4 1 (62.,0.0)
+DEAL::4 2 (83.,0.0)
+DEAL::4 3 (1.0e+02,0.0)
+DEAL::4 4 (1.2e+02,0.0)
diff --git a/tests/bits/complex_real_full_matrix_36.cc b/tests/bits/complex_real_full_matrix_36.cc
new file mode 100644 (file)
index 0000000..58e071c
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::diag_add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_36/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m.diagadd (3.141);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_36/cmp/generic b/tests/bits/complex_real_full_matrix_36/cmp/generic
new file mode 100644 (file)
index 0000000..19f141d
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (5.1,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (9.1,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (15.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (23.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (33.,0.0)
+DEAL::0 0 (5.1,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (9.1,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (15.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (23.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (33.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_37.cc b/tests/bits/complex_real_full_matrix_37.cc
new file mode 100644 (file)
index 0000000..3c7e3c7
--- /dev/null
@@ -0,0 +1,38 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::equ. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_37/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m,n;
+  make_matrix (m);
+  make_matrix (n);
+  m.equ (3.14159, n);
+  
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_37/cmp/generic b/tests/bits/complex_real_full_matrix_37/cmp/generic
new file mode 100644 (file)
index 0000000..570d6dd
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (6.3,0.0)
+DEAL::0 1 (9.4,0.0)
+DEAL::0 2 (13.,0.0)
+DEAL::0 3 (16.,0.0)
+DEAL::0 4 (19.,0.0)
+DEAL::1 0 (13.,0.0)
+DEAL::1 1 (19.,0.0)
+DEAL::1 2 (25.,0.0)
+DEAL::1 3 (31.,0.0)
+DEAL::1 4 (38.,0.0)
+DEAL::2 0 (19.,0.0)
+DEAL::2 1 (28.,0.0)
+DEAL::2 2 (38.,0.0)
+DEAL::2 3 (47.,0.0)
+DEAL::2 4 (57.,0.0)
+DEAL::3 0 (25.,0.0)
+DEAL::3 1 (38.,0.0)
+DEAL::3 2 (50.,0.0)
+DEAL::3 3 (63.,0.0)
+DEAL::3 4 (75.,0.0)
+DEAL::4 0 (31.,0.0)
+DEAL::4 1 (47.,0.0)
+DEAL::4 2 (63.,0.0)
+DEAL::4 3 (79.,0.0)
+DEAL::4 4 (94.,0.0)
+DEAL::0 0 (6.3,0.0)
+DEAL::0 1 (9.4,0.0)
+DEAL::0 2 (13.,0.0)
+DEAL::0 3 (16.,0.0)
+DEAL::0 4 (19.,0.0)
+DEAL::1 0 (13.,0.0)
+DEAL::1 1 (19.,0.0)
+DEAL::1 2 (25.,0.0)
+DEAL::1 3 (31.,0.0)
+DEAL::1 4 (38.,0.0)
+DEAL::2 0 (19.,0.0)
+DEAL::2 1 (28.,0.0)
+DEAL::2 2 (38.,0.0)
+DEAL::2 3 (47.,0.0)
+DEAL::2 4 (57.,0.0)
+DEAL::3 0 (25.,0.0)
+DEAL::3 1 (38.,0.0)
+DEAL::3 2 (50.,0.0)
+DEAL::3 3 (63.,0.0)
+DEAL::3 4 (75.,0.0)
+DEAL::4 0 (31.,0.0)
+DEAL::4 1 (47.,0.0)
+DEAL::4 2 (63.,0.0)
+DEAL::4 3 (79.,0.0)
+DEAL::4 4 (94.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_38.cc b/tests/bits/complex_real_full_matrix_38.cc
new file mode 100644 (file)
index 0000000..52d9895
--- /dev/null
@@ -0,0 +1,38 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::equ(4). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_38/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m,n,o;
+  make_matrix (m);
+  make_matrix (n);
+  make_matrix (o);
+  m.equ (3.1415, n, 2.718, o);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_38/cmp/generic b/tests/bits/complex_real_full_matrix_38/cmp/generic
new file mode 100644 (file)
index 0000000..9449cfb
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (12.,0.0)
+DEAL::0 1 (18.,0.0)
+DEAL::0 2 (23.,0.0)
+DEAL::0 3 (29.,0.0)
+DEAL::0 4 (35.,0.0)
+DEAL::1 0 (23.,0.0)
+DEAL::1 1 (35.,0.0)
+DEAL::1 2 (47.,0.0)
+DEAL::1 3 (59.,0.0)
+DEAL::1 4 (70.,0.0)
+DEAL::2 0 (35.,0.0)
+DEAL::2 1 (53.,0.0)
+DEAL::2 2 (70.,0.0)
+DEAL::2 3 (88.,0.0)
+DEAL::2 4 (1.1e+02,0.0)
+DEAL::3 0 (47.,0.0)
+DEAL::3 1 (70.,0.0)
+DEAL::3 2 (94.,0.0)
+DEAL::3 3 (1.2e+02,0.0)
+DEAL::3 4 (1.4e+02,0.0)
+DEAL::4 0 (59.,0.0)
+DEAL::4 1 (88.,0.0)
+DEAL::4 2 (1.2e+02,0.0)
+DEAL::4 3 (1.5e+02,0.0)
+DEAL::4 4 (1.8e+02,0.0)
+DEAL::0 0 (12.,0.0)
+DEAL::0 1 (18.,0.0)
+DEAL::0 2 (23.,0.0)
+DEAL::0 3 (29.,0.0)
+DEAL::0 4 (35.,0.0)
+DEAL::1 0 (23.,0.0)
+DEAL::1 1 (35.,0.0)
+DEAL::1 2 (47.,0.0)
+DEAL::1 3 (59.,0.0)
+DEAL::1 4 (70.,0.0)
+DEAL::2 0 (35.,0.0)
+DEAL::2 1 (53.,0.0)
+DEAL::2 2 (70.,0.0)
+DEAL::2 3 (88.,0.0)
+DEAL::2 4 (1.1e+02,0.0)
+DEAL::3 0 (47.,0.0)
+DEAL::3 1 (70.,0.0)
+DEAL::3 2 (94.,0.0)
+DEAL::3 3 (1.2e+02,0.0)
+DEAL::3 4 (1.4e+02,0.0)
+DEAL::4 0 (59.,0.0)
+DEAL::4 1 (88.,0.0)
+DEAL::4 2 (1.2e+02,0.0)
+DEAL::4 3 (1.5e+02,0.0)
+DEAL::4 4 (1.8e+02,0.0)
diff --git a/tests/bits/complex_real_full_matrix_39.cc b/tests/bits/complex_real_full_matrix_39.cc
new file mode 100644 (file)
index 0000000..c11f212
--- /dev/null
@@ -0,0 +1,39 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::equ (6). like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_39/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m,n,o,p;
+  make_matrix (m);
+  make_matrix (n);
+  make_matrix (o);
+  make_matrix (p);
+  m.equ (3.1415, n, 2.718, o, 1.414, p);
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_39/cmp/generic b/tests/bits/complex_real_full_matrix_39/cmp/generic
new file mode 100644 (file)
index 0000000..f4c9095
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (15.,0.0)
+DEAL::0 1 (22.,0.0)
+DEAL::0 2 (29.,0.0)
+DEAL::0 3 (36.,0.0)
+DEAL::0 4 (44.,0.0)
+DEAL::1 0 (29.,0.0)
+DEAL::1 1 (44.,0.0)
+DEAL::1 2 (58.,0.0)
+DEAL::1 3 (73.,0.0)
+DEAL::1 4 (87.,0.0)
+DEAL::2 0 (44.,0.0)
+DEAL::2 1 (65.,0.0)
+DEAL::2 2 (87.,0.0)
+DEAL::2 3 (1.1e+02,0.0)
+DEAL::2 4 (1.3e+02,0.0)
+DEAL::3 0 (58.,0.0)
+DEAL::3 1 (87.,0.0)
+DEAL::3 2 (1.2e+02,0.0)
+DEAL::3 3 (1.5e+02,0.0)
+DEAL::3 4 (1.7e+02,0.0)
+DEAL::4 0 (73.,0.0)
+DEAL::4 1 (1.1e+02,0.0)
+DEAL::4 2 (1.5e+02,0.0)
+DEAL::4 3 (1.8e+02,0.0)
+DEAL::4 4 (2.2e+02,0.0)
+DEAL::0 0 (15.,0.0)
+DEAL::0 1 (22.,0.0)
+DEAL::0 2 (29.,0.0)
+DEAL::0 3 (36.,0.0)
+DEAL::0 4 (44.,0.0)
+DEAL::1 0 (29.,0.0)
+DEAL::1 1 (44.,0.0)
+DEAL::1 2 (58.,0.0)
+DEAL::1 3 (73.,0.0)
+DEAL::1 4 (87.,0.0)
+DEAL::2 0 (44.,0.0)
+DEAL::2 1 (65.,0.0)
+DEAL::2 2 (87.,0.0)
+DEAL::2 3 (1.1e+02,0.0)
+DEAL::2 4 (1.3e+02,0.0)
+DEAL::3 0 (58.,0.0)
+DEAL::3 1 (87.,0.0)
+DEAL::3 2 (1.2e+02,0.0)
+DEAL::3 3 (1.5e+02,0.0)
+DEAL::3 4 (1.7e+02,0.0)
+DEAL::4 0 (73.,0.0)
+DEAL::4 1 (1.1e+02,0.0)
+DEAL::4 2 (1.5e+02,0.0)
+DEAL::4 3 (1.8e+02,0.0)
+DEAL::4 4 (2.2e+02,0.0)
diff --git a/tests/bits/complex_real_full_matrix_40.cc b/tests/bits/complex_real_full_matrix_40.cc
new file mode 100644 (file)
index 0000000..e4a83d4
--- /dev/null
@@ -0,0 +1,36 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::symmetrize. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_40/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m.symmetrize ();
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_40/cmp/generic b/tests/bits/complex_real_full_matrix_40/cmp/generic
new file mode 100644 (file)
index 0000000..fce83ff
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.5,0.0)
+DEAL::0 2 (5.0,0.0)
+DEAL::0 3 (6.5,0.0)
+DEAL::0 4 (8.0,0.0)
+DEAL::1 0 (3.5,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.5,0.0)
+DEAL::1 3 (11.,0.0)
+DEAL::1 4 (14.,0.0)
+DEAL::2 0 (5.0,0.0)
+DEAL::2 1 (8.5,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (16.,0.0)
+DEAL::2 4 (19.,0.0)
+DEAL::3 0 (6.5,0.0)
+DEAL::3 1 (11.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (8.0,0.0)
+DEAL::4 1 (14.,0.0)
+DEAL::4 2 (19.,0.0)
+DEAL::4 3 (24.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.5,0.0)
+DEAL::0 2 (5.0,0.0)
+DEAL::0 3 (6.5,0.0)
+DEAL::0 4 (8.0,0.0)
+DEAL::1 0 (3.5,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.5,0.0)
+DEAL::1 3 (11.,0.0)
+DEAL::1 4 (14.,0.0)
+DEAL::2 0 (5.0,0.0)
+DEAL::2 1 (8.5,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (16.,0.0)
+DEAL::2 4 (19.,0.0)
+DEAL::3 0 (6.5,0.0)
+DEAL::3 1 (11.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (8.0,0.0)
+DEAL::4 1 (14.,0.0)
+DEAL::4 2 (19.,0.0)
+DEAL::4 3 (24.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_41.cc b/tests/bits/complex_real_full_matrix_41.cc
new file mode 100644 (file)
index 0000000..e122bf6
--- /dev/null
@@ -0,0 +1,37 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::gauss_jordan. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_41/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  m(0,0) = m(1,1) = m(2,2) = m(3,3) = m(4,4) = 50;
+  m.gauss_jordan ();
+  print_matrix (m);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_41/cmp/generic b/tests/bits/complex_real_full_matrix_41/cmp/generic
new file mode 100644 (file)
index 0000000..649f30d
--- /dev/null
@@ -0,0 +1,51 @@
+
+DEAL::0 0 (0.021,0.0)
+DEAL::0 1 (0.00039,-0.0)
+DEAL::0 2 (0.00060,0.0)
+DEAL::0 3 (0.00095,0.0)
+DEAL::0 4 (0.0017,0.0)
+DEAL::1 0 (0.00052,0.0)
+DEAL::1 1 (0.022,0.0)
+DEAL::1 2 (0.0013,-0.0)
+DEAL::1 3 (0.0021,0.0)
+DEAL::1 4 (0.0037,0.0)
+DEAL::2 0 (0.00090,0.0)
+DEAL::2 1 (0.0015,0.0)
+DEAL::2 2 (0.024,0.0)
+DEAL::2 3 (0.0036,-0.0)
+DEAL::2 4 (0.0065,0.0)
+DEAL::3 0 (0.0015,0.0)
+DEAL::3 1 (0.0025,0.0)
+DEAL::3 2 (0.0038,0.0)
+DEAL::3 3 (0.027,0.0)
+DEAL::3 4 (-0.011,-0.0)
+DEAL::4 0 (0.0028,-0.0)
+DEAL::4 1 (0.0047,-0.0)
+DEAL::4 2 (0.0072,-0.0)
+DEAL::4 3 (-0.011,-0.0)
+DEAL::4 4 (0.030,0.0)
+DEAL::0 0 (0.021,0.0)
+DEAL::0 1 (0.00039,-0.0)
+DEAL::0 2 (0.00060,0.0)
+DEAL::0 3 (0.00095,0.0)
+DEAL::0 4 (0.0017,0.0)
+DEAL::1 0 (0.00052,0.0)
+DEAL::1 1 (0.022,0.0)
+DEAL::1 2 (0.0013,-0.0)
+DEAL::1 3 (0.0021,0.0)
+DEAL::1 4 (0.0037,0.0)
+DEAL::2 0 (0.00090,0.0)
+DEAL::2 1 (0.0015,0.0)
+DEAL::2 2 (0.024,0.0)
+DEAL::2 3 (0.0036,-0.0)
+DEAL::2 4 (0.0065,0.0)
+DEAL::3 0 (0.0015,0.0)
+DEAL::3 1 (0.0025,0.0)
+DEAL::3 2 (0.0038,0.0)
+DEAL::3 3 (0.027,0.0)
+DEAL::3 4 (-0.011,-0.0)
+DEAL::4 0 (0.0028,-0.0)
+DEAL::4 1 (0.0047,-0.0)
+DEAL::4 2 (0.0072,-0.0)
+DEAL::4 3 (-0.011,-0.0)
+DEAL::4 4 (0.030,0.0)
diff --git a/tests/bits/complex_real_full_matrix_42.cc b/tests/bits/complex_real_full_matrix_42.cc
new file mode 100644 (file)
index 0000000..8868062
--- /dev/null
@@ -0,0 +1,44 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::invert. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_42/output";
+
+
+template <typename number>
+void
+check ()
+{
+  for (unsigned int n=1; n<=5; ++n)
+    {
+      const std::complex<number> array[] = { 50, 2, 3, 4, 5,
+                              6, 50, 8, 9, 0,
+                              1, 2, 50, 4, 5,
+                              6, 7, 8, 50, 0,
+                              1, 2, 3, 4, 50 };
+  
+      FullMatrix<std::complex<number> > m (n,n,array), p(n,n);
+      p.invert (m);
+      print_matrix (p);
+    }
+}
+
diff --git a/tests/bits/complex_real_full_matrix_42/cmp/generic b/tests/bits/complex_real_full_matrix_42/cmp/generic
new file mode 100644 (file)
index 0000000..47435b6
--- /dev/null
@@ -0,0 +1,111 @@
+
+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)
diff --git a/tests/bits/complex_real_full_matrix_43.cc b/tests/bits/complex_real_full_matrix_43.cc
new file mode 100644 (file)
index 0000000..80c40a0
--- /dev/null
@@ -0,0 +1,46 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::mmult. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_43/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m, n, o;
+  make_matrix (m);
+  make_matrix (n);
+  make_matrix (o);
+
+  m.mmult (n, o, true);
+  print_matrix (m);
+  print_matrix (n);
+  print_matrix (o);
+
+  m.mmult (n, o, false);
+  print_matrix (m);
+  print_matrix (n);
+  print_matrix (o);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_43/cmp/generic b/tests/bits/complex_real_full_matrix_43/cmp/generic
new file mode 100644 (file)
index 0000000..c957cb8
--- /dev/null
@@ -0,0 +1,301 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (1.4e+02,0.0)
+DEAL::0 1 (2.1e+02,0.0)
+DEAL::0 2 (2.8e+02,0.0)
+DEAL::0 3 (3.6e+02,0.0)
+DEAL::0 4 (4.3e+02,0.0)
+DEAL::1 0 (2.8e+02,0.0)
+DEAL::1 1 (4.3e+02,0.0)
+DEAL::1 2 (5.7e+02,0.0)
+DEAL::1 3 (7.1e+02,0.0)
+DEAL::1 4 (8.5e+02,0.0)
+DEAL::2 0 (4.3e+02,0.0)
+DEAL::2 1 (6.4e+02,0.0)
+DEAL::2 2 (8.5e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.7e+02,0.0)
+DEAL::3 1 (8.5e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (7.1e+02,0.0)
+DEAL::4 1 (1.1e+03,0.0)
+DEAL::4 2 (1.4e+03,0.0)
+DEAL::4 3 (1.8e+03,0.0)
+DEAL::4 4 (2.1e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (1.4e+02,0.0)
+DEAL::0 1 (2.1e+02,0.0)
+DEAL::0 2 (2.8e+02,0.0)
+DEAL::0 3 (3.5e+02,0.0)
+DEAL::0 4 (4.2e+02,0.0)
+DEAL::1 0 (2.8e+02,0.0)
+DEAL::1 1 (4.2e+02,0.0)
+DEAL::1 2 (5.6e+02,0.0)
+DEAL::1 3 (7.0e+02,0.0)
+DEAL::1 4 (8.4e+02,0.0)
+DEAL::2 0 (4.2e+02,0.0)
+DEAL::2 1 (6.3e+02,0.0)
+DEAL::2 2 (8.4e+02,0.0)
+DEAL::2 3 (1.0e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.6e+02,0.0)
+DEAL::3 1 (8.4e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (7.0e+02,0.0)
+DEAL::4 1 (1.0e+03,0.0)
+DEAL::4 2 (1.4e+03,0.0)
+DEAL::4 3 (1.8e+03,0.0)
+DEAL::4 4 (2.1e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (1.4e+02,0.0)
+DEAL::0 1 (2.1e+02,0.0)
+DEAL::0 2 (2.8e+02,0.0)
+DEAL::0 3 (3.6e+02,0.0)
+DEAL::0 4 (4.3e+02,0.0)
+DEAL::1 0 (2.8e+02,0.0)
+DEAL::1 1 (4.3e+02,0.0)
+DEAL::1 2 (5.7e+02,0.0)
+DEAL::1 3 (7.1e+02,0.0)
+DEAL::1 4 (8.5e+02,0.0)
+DEAL::2 0 (4.3e+02,0.0)
+DEAL::2 1 (6.4e+02,0.0)
+DEAL::2 2 (8.5e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.7e+02,0.0)
+DEAL::3 1 (8.5e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (7.1e+02,0.0)
+DEAL::4 1 (1.1e+03,0.0)
+DEAL::4 2 (1.4e+03,0.0)
+DEAL::4 3 (1.8e+03,0.0)
+DEAL::4 4 (2.1e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (1.4e+02,0.0)
+DEAL::0 1 (2.1e+02,0.0)
+DEAL::0 2 (2.8e+02,0.0)
+DEAL::0 3 (3.5e+02,0.0)
+DEAL::0 4 (4.2e+02,0.0)
+DEAL::1 0 (2.8e+02,0.0)
+DEAL::1 1 (4.2e+02,0.0)
+DEAL::1 2 (5.6e+02,0.0)
+DEAL::1 3 (7.0e+02,0.0)
+DEAL::1 4 (8.4e+02,0.0)
+DEAL::2 0 (4.2e+02,0.0)
+DEAL::2 1 (6.3e+02,0.0)
+DEAL::2 2 (8.4e+02,0.0)
+DEAL::2 3 (1.0e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.6e+02,0.0)
+DEAL::3 1 (8.4e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (7.0e+02,0.0)
+DEAL::4 1 (1.0e+03,0.0)
+DEAL::4 2 (1.4e+03,0.0)
+DEAL::4 3 (1.8e+03,0.0)
+DEAL::4 4 (2.1e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_44.cc b/tests/bits/complex_real_full_matrix_44.cc
new file mode 100644 (file)
index 0000000..6c1dc9c
--- /dev/null
@@ -0,0 +1,46 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::Tmmult. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_44/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m, n, o;
+  make_matrix (m);
+  make_matrix (n);
+  make_matrix (o);
+
+  m.Tmmult (n, o, true);
+  print_matrix (m);
+  print_matrix (n);
+  print_matrix (o);
+
+  m.Tmmult (n, o, false);
+  print_matrix (m);
+  print_matrix (n);
+  print_matrix (o);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_44/cmp/generic b/tests/bits/complex_real_full_matrix_44/cmp/generic
new file mode 100644 (file)
index 0000000..4e09c9d
--- /dev/null
@@ -0,0 +1,301 @@
+
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.2e+02,0.0)
+DEAL::0 1 (3.3e+02,0.0)
+DEAL::0 2 (4.4e+02,0.0)
+DEAL::0 3 (5.6e+02,0.0)
+DEAL::0 4 (6.7e+02,0.0)
+DEAL::1 0 (3.3e+02,0.0)
+DEAL::1 1 (5.0e+02,0.0)
+DEAL::1 2 (6.7e+02,0.0)
+DEAL::1 3 (8.4e+02,0.0)
+DEAL::1 4 (1.0e+03,0.0)
+DEAL::2 0 (4.5e+02,0.0)
+DEAL::2 1 (6.7e+02,0.0)
+DEAL::2 2 (8.9e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.6e+02,0.0)
+DEAL::3 1 (8.4e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (6.7e+02,0.0)
+DEAL::4 1 (1.0e+03,0.0)
+DEAL::4 2 (1.3e+03,0.0)
+DEAL::4 3 (1.7e+03,0.0)
+DEAL::4 4 (2.0e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.2e+02,0.0)
+DEAL::0 1 (3.3e+02,0.0)
+DEAL::0 2 (4.4e+02,0.0)
+DEAL::0 3 (5.5e+02,0.0)
+DEAL::0 4 (6.6e+02,0.0)
+DEAL::1 0 (3.3e+02,0.0)
+DEAL::1 1 (5.0e+02,0.0)
+DEAL::1 2 (6.6e+02,0.0)
+DEAL::1 3 (8.2e+02,0.0)
+DEAL::1 4 (9.9e+02,0.0)
+DEAL::2 0 (4.4e+02,0.0)
+DEAL::2 1 (6.6e+02,0.0)
+DEAL::2 2 (8.8e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.5e+02,0.0)
+DEAL::3 1 (8.2e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.6e+03,0.0)
+DEAL::4 0 (6.6e+02,0.0)
+DEAL::4 1 (9.9e+02,0.0)
+DEAL::4 2 (1.3e+03,0.0)
+DEAL::4 3 (1.6e+03,0.0)
+DEAL::4 4 (2.0e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.2e+02,0.0)
+DEAL::0 1 (3.3e+02,0.0)
+DEAL::0 2 (4.4e+02,0.0)
+DEAL::0 3 (5.6e+02,0.0)
+DEAL::0 4 (6.7e+02,0.0)
+DEAL::1 0 (3.3e+02,0.0)
+DEAL::1 1 (5.0e+02,0.0)
+DEAL::1 2 (6.7e+02,0.0)
+DEAL::1 3 (8.4e+02,0.0)
+DEAL::1 4 (1.0e+03,0.0)
+DEAL::2 0 (4.5e+02,0.0)
+DEAL::2 1 (6.7e+02,0.0)
+DEAL::2 2 (8.9e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.6e+02,0.0)
+DEAL::3 1 (8.4e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.7e+03,0.0)
+DEAL::4 0 (6.7e+02,0.0)
+DEAL::4 1 (1.0e+03,0.0)
+DEAL::4 2 (1.3e+03,0.0)
+DEAL::4 3 (1.7e+03,0.0)
+DEAL::4 4 (2.0e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
+DEAL::0 0 (2.2e+02,0.0)
+DEAL::0 1 (3.3e+02,0.0)
+DEAL::0 2 (4.4e+02,0.0)
+DEAL::0 3 (5.5e+02,0.0)
+DEAL::0 4 (6.6e+02,0.0)
+DEAL::1 0 (3.3e+02,0.0)
+DEAL::1 1 (5.0e+02,0.0)
+DEAL::1 2 (6.6e+02,0.0)
+DEAL::1 3 (8.2e+02,0.0)
+DEAL::1 4 (9.9e+02,0.0)
+DEAL::2 0 (4.4e+02,0.0)
+DEAL::2 1 (6.6e+02,0.0)
+DEAL::2 2 (8.8e+02,0.0)
+DEAL::2 3 (1.1e+03,0.0)
+DEAL::2 4 (1.3e+03,0.0)
+DEAL::3 0 (5.5e+02,0.0)
+DEAL::3 1 (8.2e+02,0.0)
+DEAL::3 2 (1.1e+03,0.0)
+DEAL::3 3 (1.4e+03,0.0)
+DEAL::3 4 (1.6e+03,0.0)
+DEAL::4 0 (6.6e+02,0.0)
+DEAL::4 1 (9.9e+02,0.0)
+DEAL::4 2 (1.3e+03,0.0)
+DEAL::4 3 (1.6e+03,0.0)
+DEAL::4 4 (2.0e+03,0.0)
+DEAL::0 0 (2.0,0.0)
+DEAL::0 1 (3.0,0.0)
+DEAL::0 2 (4.0,0.0)
+DEAL::0 3 (5.0,0.0)
+DEAL::0 4 (6.0,0.0)
+DEAL::1 0 (4.0,0.0)
+DEAL::1 1 (6.0,0.0)
+DEAL::1 2 (8.0,0.0)
+DEAL::1 3 (10.,0.0)
+DEAL::1 4 (12.,0.0)
+DEAL::2 0 (6.0,0.0)
+DEAL::2 1 (9.0,0.0)
+DEAL::2 2 (12.,0.0)
+DEAL::2 3 (15.,0.0)
+DEAL::2 4 (18.,0.0)
+DEAL::3 0 (8.0,0.0)
+DEAL::3 1 (12.,0.0)
+DEAL::3 2 (16.,0.0)
+DEAL::3 3 (20.,0.0)
+DEAL::3 4 (24.,0.0)
+DEAL::4 0 (10.,0.0)
+DEAL::4 1 (15.,0.0)
+DEAL::4 2 (20.,0.0)
+DEAL::4 3 (25.,0.0)
+DEAL::4 4 (30.,0.0)
diff --git a/tests/bits/complex_real_full_matrix_45.cc b/tests/bits/complex_real_full_matrix_45.cc
new file mode 100644 (file)
index 0000000..27b45cf
--- /dev/null
@@ -0,0 +1,45 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::vmult. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_45/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  Vector<std::complex<number> > v, w;
+  make_vector (v);
+  make_vector (w);
+
+  m.vmult (v, w, true);
+  print_vector (v);
+  print_vector (w);
+
+  m.vmult (v, w, false);
+  print_vector (v);
+  print_vector (w);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_45/cmp/generic b/tests/bits/complex_real_full_matrix_45/cmp/generic
new file mode 100644 (file)
index 0000000..ac31e79
--- /dev/null
@@ -0,0 +1,41 @@
+
+DEAL::0 (71.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.6e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (70.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.5e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (71.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.6e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (70.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.5e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
diff --git a/tests/bits/complex_real_full_matrix_46.cc b/tests/bits/complex_real_full_matrix_46.cc
new file mode 100644 (file)
index 0000000..e850b96
--- /dev/null
@@ -0,0 +1,41 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::vmult_add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_46/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  Vector<std::complex<number> > v, w;
+  make_vector (v);
+  make_vector (w);
+
+  m.vmult_add (v, w);
+  print_vector (v);
+  print_vector (w);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_46/cmp/generic b/tests/bits/complex_real_full_matrix_46/cmp/generic
new file mode 100644 (file)
index 0000000..3b9b54e
--- /dev/null
@@ -0,0 +1,21 @@
+
+DEAL::0 (71.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.6e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (71.,0.0)
+DEAL::1 (1.4e+02,0.0)
+DEAL::2 (2.1e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.6e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
diff --git a/tests/bits/complex_real_full_matrix_47.cc b/tests/bits/complex_real_full_matrix_47.cc
new file mode 100644 (file)
index 0000000..983f3f6
--- /dev/null
@@ -0,0 +1,45 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::Tvmult. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_47/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  Vector<std::complex<number> > v, w;
+  make_vector (v);
+  make_vector (w);
+
+  m.Tvmult (v, w, true);
+  print_vector (v);
+  print_vector (w);
+
+  m.Tvmult (v, w, false);
+  print_vector (v);
+  print_vector (w);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_47/cmp/generic b/tests/bits/complex_real_full_matrix_47/cmp/generic
new file mode 100644 (file)
index 0000000..5fd7d3f
--- /dev/null
@@ -0,0 +1,41 @@
+
+DEAL::0 (1.1e+02,0.0)
+DEAL::1 (1.7e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.1e+02,0.0)
+DEAL::1 (1.6e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.3e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.1e+02,0.0)
+DEAL::1 (1.7e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.1e+02,0.0)
+DEAL::1 (1.6e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.3e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
diff --git a/tests/bits/complex_real_full_matrix_48.cc b/tests/bits/complex_real_full_matrix_48.cc
new file mode 100644 (file)
index 0000000..8401962
--- /dev/null
@@ -0,0 +1,41 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::Tvmult_add. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_48/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  Vector<std::complex<number> > v, w;
+  make_vector (v);
+  make_vector (w);
+
+  m.Tvmult_add (v, w);
+  print_vector (v);
+  print_vector (w);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_48/cmp/generic b/tests/bits/complex_real_full_matrix_48/cmp/generic
new file mode 100644 (file)
index 0000000..4a337a7
--- /dev/null
@@ -0,0 +1,21 @@
+
+DEAL::0 (1.1e+02,0.0)
+DEAL::1 (1.7e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.1e+02,0.0)
+DEAL::1 (1.7e+02,0.0)
+DEAL::2 (2.2e+02,0.0)
+DEAL::3 (2.8e+02,0.0)
+DEAL::4 (3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
diff --git a/tests/bits/complex_real_full_matrix_49.cc b/tests/bits/complex_real_full_matrix_49.cc
new file mode 100644 (file)
index 0000000..85dbdf4
--- /dev/null
@@ -0,0 +1,41 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::precondition_Jacobi. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_49/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  Vector<std::complex<number> > v, w;
+  make_vector (v);
+  make_vector (w);
+
+  m.precondition_Jacobi (v, w, 3.141);
+  print_vector (v);
+  print_vector (w);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_49/cmp/generic b/tests/bits/complex_real_full_matrix_49/cmp/generic
new file mode 100644 (file)
index 0000000..7d8270c
--- /dev/null
@@ -0,0 +1,21 @@
+
+DEAL::0 (1.6,0.0)
+DEAL::1 (1.0,0.0)
+DEAL::2 (0.79,0.0)
+DEAL::3 (0.63,0.0)
+DEAL::4 (0.52,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.6,0.0)
+DEAL::1 (1.0,0.0)
+DEAL::2 (0.79,0.0)
+DEAL::3 (0.63,0.0)
+DEAL::4 (0.52,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
diff --git a/tests/bits/complex_real_full_matrix_50.cc b/tests/bits/complex_real_full_matrix_50.cc
new file mode 100644 (file)
index 0000000..7ca464e
--- /dev/null
@@ -0,0 +1,43 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::residual. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_50/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  Vector<std::complex<number> > v, w, x;
+  make_vector (v);
+  make_vector (w);
+  make_vector (x);
+
+  m.residual (v, w, x);
+  print_vector (v);
+  print_vector (w);
+  print_vector (x);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_50/cmp/generic b/tests/bits/complex_real_full_matrix_50/cmp/generic
new file mode 100644 (file)
index 0000000..215ee59
--- /dev/null
@@ -0,0 +1,31 @@
+
+DEAL::0 (-69.,0.0)
+DEAL::1 (-1.4e+02,0.0)
+DEAL::2 (-2.1e+02,0.0)
+DEAL::3 (-2.8e+02,0.0)
+DEAL::4 (-3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (-69.,0.0)
+DEAL::1 (-1.4e+02,0.0)
+DEAL::2 (-2.1e+02,0.0)
+DEAL::3 (-2.8e+02,0.0)
+DEAL::4 (-3.4e+02,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
diff --git a/tests/bits/complex_real_full_matrix_51.cc b/tests/bits/complex_real_full_matrix_51.cc
new file mode 100644 (file)
index 0000000..bc6e3b1
--- /dev/null
@@ -0,0 +1,41 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::forward. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_51/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  Vector<std::complex<number> > x, y;
+  make_vector (x);
+  make_vector (y);
+
+  m.forward (x, y);
+  print_vector (x);
+  print_vector (y);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_51/cmp/generic b/tests/bits/complex_real_full_matrix_51/cmp/generic
new file mode 100644 (file)
index 0000000..4ee3233
--- /dev/null
@@ -0,0 +1,21 @@
+
+DEAL::0 (0.50,0.0)
+DEAL::1 (0.0,0.0)
+DEAL::2 (0.0,0.0)
+DEAL::3 (0.0,0.0)
+DEAL::4 (0.0,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (0.50,0.0)
+DEAL::1 (0.0,0.0)
+DEAL::2 (0.0,0.0)
+DEAL::3 (0.0,0.0)
+DEAL::4 (0.0,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
diff --git a/tests/bits/complex_real_full_matrix_52.cc b/tests/bits/complex_real_full_matrix_52.cc
new file mode 100644 (file)
index 0000000..32914bc
--- /dev/null
@@ -0,0 +1,41 @@
+//----------------------------------------------------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2007 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------------------------------------------------
+
+
+// check FullMatrix::backward. like the full_matrix_* tests, but use
+// complex-valued matrices and vectors, even though we only store real values
+// in them
+
+
+#include "../tests.h"
+#include "full_matrix_common.h"
+
+
+std::string output_file_name = "complex_real_full_matrix_52/output";
+
+
+template <typename number>
+void
+check ()
+{
+  FullMatrix<std::complex<number> > m;
+  make_matrix (m);
+  Vector<std::complex<number> > x, y;
+  make_vector (x);
+  make_vector (y);
+
+  m.backward (x, y);
+  print_vector (x);
+  print_vector (y);
+}
+
diff --git a/tests/bits/complex_real_full_matrix_52/cmp/generic b/tests/bits/complex_real_full_matrix_52/cmp/generic
new file mode 100644 (file)
index 0000000..7864d6a
--- /dev/null
@@ -0,0 +1,21 @@
+
+DEAL::0 (0.0,0.0)
+DEAL::1 (0.0,0.0)
+DEAL::2 (0.0,0.0)
+DEAL::3 (0.0,0.0)
+DEAL::4 (0.17,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
+DEAL::0 (0.0,0.0)
+DEAL::1 (0.0,0.0)
+DEAL::2 (0.0,0.0)
+DEAL::3 (0.0,0.0)
+DEAL::4 (0.17,0.0)
+DEAL::0 (1.0,0.0)
+DEAL::1 (2.0,0.0)
+DEAL::2 (3.0,0.0)
+DEAL::3 (4.0,0.0)
+DEAL::4 (5.0,0.0)
index b6a7b59fb449e0ae9ad1e542fb0357722bd5a789..ae2f986eec7aac5a1f27e593d72c95746e69fbfc 100644 (file)
@@ -28,7 +28,7 @@ check ()
 {
   FullMatrix<number> m;
   make_matrix (m);
-  Vector<double> v;
+  Vector<number> v;
   make_vector (v);
   
   deallog << m.matrix_norm_square (v) << std::endl;
index 24569a620ad09a305aaf2cd7a458d170f6eefeef..a18128771929502f402661e084c708b2a2d564d5 100644 (file)
@@ -28,7 +28,7 @@ check ()
 {
   FullMatrix<number> m;
   make_matrix (m);
-  Vector<double> v,w;
+  Vector<number> v,w;
   make_vector (v);
   make_vector (w);
   for (unsigned int i=0; i<w.size(); ++i)
index 899cf002190a00e7fd79b90a680e47e0b6193bee..6dcd2e8110598b14e8fe68ce3f43aee1704d6ffd 100644 (file)
@@ -28,7 +28,7 @@ check ()
 {
   FullMatrix<number> m;
   make_matrix (m);
-  Vector<double> v, w;
+  Vector<number> v, w;
   make_vector (v);
   make_vector (w);
 
index b943470ad0712dfc3f51e131c77f36659f977ecc..e9365546354b2c64c0a1d4927192cee73f579da0 100644 (file)
@@ -28,7 +28,7 @@ check ()
 {
   FullMatrix<number> m;
   make_matrix (m);
-  Vector<double> v, w;
+  Vector<number> v, w;
   make_vector (v);
   make_vector (w);
 
index 2add1958e4148347cfd07d2df4081c24fcc13cb6..45d2158ccf73ddbe8e9a6c79873701b78762b674 100644 (file)
@@ -28,7 +28,7 @@ check ()
 {
   FullMatrix<number> m;
   make_matrix (m);
-  Vector<double> v, w;
+  Vector<number> v, w;
   make_vector (v);
   make_vector (w);
 
index 1f9bed89ab30ccfd90da587dd16a719fdadcb438..2cc4af5b8ab6760a1f5aafe3b576739a0a4f277e 100644 (file)
@@ -28,7 +28,7 @@ check ()
 {
   FullMatrix<number> m;
   make_matrix (m);
-  Vector<double> v, w;
+  Vector<number> v, w;
   make_vector (v);
   make_vector (w);
 
index 73020ba723450604d35405b450d62ba360c524e6..9bf58c62c3e1bfcd59537c56ef3b9a9e4fe578dc 100644 (file)
@@ -28,7 +28,7 @@ check ()
 {
   FullMatrix<number> m;
   make_matrix (m);
-  Vector<double> v, w;
+  Vector<number> v, w;
   make_vector (v);
   make_vector (w);
 
index f44a4a8375b00acbac6d9867e7d48d0556e2e2ce..8bdbee02ace1c96935fcf04f107e12dc8c31d4af 100644 (file)
@@ -28,7 +28,7 @@ check ()
 {
   FullMatrix<number> m;
   make_matrix (m);
-  Vector<double> v, w, x;
+  Vector<number> v, w, x;
   make_vector (v);
   make_vector (w);
   make_vector (x);
index 39121d657a9650a8c2bcc3fda6876625e34487d6..f134d7d1bb0c2676548f31bd61d044ababce2417 100644 (file)
@@ -28,7 +28,7 @@ check ()
 {
   FullMatrix<number> m;
   make_matrix (m);
-  Vector<double> x, y;
+  Vector<number> x, y;
   make_vector (x);
   make_vector (y);
 
index 230ba3fa22f46e615d7a5033b2fd29b4816d8232..eaf1439fa73ccf0e5e55f1777963613d308aa6a8 100644 (file)
@@ -28,7 +28,7 @@ check ()
 {
   FullMatrix<number> m;
   make_matrix (m);
-  Vector<double> x, y;
+  Vector<number> x, y;
   make_vector (x);
   make_vector (y);
 
index a0278f51f05e6b7c828da55f2331a439871f7dab..d4268b1fbee3fc34fa685909d7d8b1ed1e8e1e7a 100644 (file)
@@ -22,8 +22,8 @@
 #include <vector>
 
 
-void test (Vector<double> &v,
-           Vector<double> &w)
+void test (Vector<number> &v,
+           Vector<number> &w)
 {
   FullMatrix<double> m(v.size(), v.size());
   for (unsigned int i=0; i<m.m(); ++i)
@@ -64,8 +64,8 @@ int main ()
 
   try
     {
-      Vector<double> v (100);
-      Vector<double> w (100);
+      Vector<number> v (100);
+      Vector<number> w (100);
       test (v,w);
     }
   catch (std::exception &exc)
index 3a493e427b5cdf105de6315dff70c5e5c513651c..1b11b75abf1f49ba427728c39aee0ca1a2bb06cd 100644 (file)
@@ -22,8 +22,8 @@
 #include <vector>
 
 
-void test (Vector<double> &v,
-           Vector<double> &w)
+void test (Vector<number> &v,
+           Vector<number> &w)
 {
   FullMatrix<double> m(v.size(), v.size());
   for (unsigned int i=0; i<m.m(); ++i)
@@ -64,8 +64,8 @@ int main ()
 
   try
     {
-      Vector<double> v (100);
-      Vector<double> w (100);
+      Vector<number> v (100);
+      Vector<number> w (100);
       test (v,w);
     }
   catch (std::exception &exc)
index 0619f826f467d5518c6ab3c56de0271d173185cb..dd813503789c50ae98b75358ad79357fad908f29 100644 (file)
@@ -22,8 +22,8 @@
 #include <vector>
 
 
-void test (Vector<double> &v,
-           Vector<double> &w)
+void test (Vector<number> &v,
+           Vector<number> &w)
 {
   FullMatrix<double> m(v.size(), v.size());
   for (unsigned int i=0; i<m.m(); ++i)
@@ -70,8 +70,8 @@ int main ()
 
   try
     {
-      Vector<double> v (100);
-      Vector<double> w (100);
+      Vector<number> v (100);
+      Vector<number> w (100);
       test (v,w);
     }
   catch (std::exception &exc)
index 2544ee002005e04473fb1c139ffbebe99fe25450..46bbb3ebed480a4861479f70f169f9c0da1b71ef 100644 (file)
@@ -22,7 +22,7 @@
 #include <vector>
 
 
-void test (Vector<double> &v)
+void test (Vector<number> &v)
 {
   FullMatrix<double> m(v.size(), v.size());
   for (unsigned int i=0; i<m.m(); ++i)
@@ -62,7 +62,7 @@ int main ()
 
   try
     {
-      Vector<double> v (100);
+      Vector<number> v (100);
       test (v);
     }
   catch (std::exception &exc)
index 8115c35491b4d00e4ac489cfea3c21e3140c15c8..1e6b8db69eff1e6495c5f67b3524812c744c2f07 100644 (file)
@@ -22,9 +22,9 @@
 #include <vector>
 
 
-void test (Vector<double> &v,
-           Vector<double> &w,
-           Vector<double> &x)
+void test (Vector<number> &v,
+           Vector<number> &w,
+           Vector<number> &x)
 {
   FullMatrix<double> m(v.size(), v.size());
   for (unsigned int i=0; i<m.m(); ++i)
@@ -72,9 +72,9 @@ int main ()
 
   try
     {
-      Vector<double> v (100);
-      Vector<double> w (100);
-      Vector<double> x (100);
+      Vector<number> v (100);
+      Vector<number> w (100);
+      Vector<number> x (100);
       test (v,w,x);
     }
   catch (std::exception &exc)

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.