From: Wolfgang Bangerth Date: Sun, 18 Nov 2007 23:55:17 +0000 (+0000) Subject: Another set of tests that now actually use complex values in complex storage. X-Git-Tag: v8.0.0~9584 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c48e665f1552b64813ddcedba390ee21d2fca2ad;p=dealii.git Another set of tests that now actually use complex values in complex storage. git-svn-id: https://svn.dealii.org/trunk@15520 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/bits/Makefile b/tests/bits/Makefile index 6e04691d35..7aa33fc36b 100644 --- a/tests/bits/Makefile +++ b/tests/bits/Makefile @@ -26,6 +26,7 @@ tests_x = grid_generator_?? \ point_inside_? \ full_matrix_* \ complex_real_full_matrix_* \ + complex_complex_full_matrix_* \ block_matrix_array_* \ data_out_0* \ data_out_faces_0* \ diff --git a/tests/bits/complex_complex_full_matrix_01.cc b/tests/bits/complex_complex_full_matrix_01.cc new file mode 100644 index 0000000000..fd3a4f3a4a --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_01.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_01/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_01/cmp/generic b/tests/bits/complex_complex_full_matrix_01/cmp/generic new file mode 100644 index 0000000000..9d9b54f64c --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_01/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_02.cc b/tests/bits/complex_complex_full_matrix_02.cc new file mode 100644 index 0000000000..ff49667e5b --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_02.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_02/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + + for (typename FullMatrix >::const_iterator + p = m.begin(); p!=m.end(); ++p) + deallog << p->row() << ' ' << p->column() << ' ' + << p->value() + << std::endl; +} + diff --git a/tests/bits/complex_complex_full_matrix_02/cmp/generic b/tests/bits/complex_complex_full_matrix_02/cmp/generic new file mode 100644 index 0000000000..9d9b54f64c --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_02/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_03.cc b/tests/bits/complex_complex_full_matrix_03.cc new file mode 100644 index 0000000000..dacb83db53 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_03.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_03/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + + + for (typename FullMatrix >::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_complex_full_matrix_03/cmp/generic b/tests/bits/complex_complex_full_matrix_03/cmp/generic new file mode 100644 index 0000000000..e1876979bd --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_03/cmp/generic @@ -0,0 +1,11 @@ + +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) diff --git a/tests/bits/complex_complex_full_matrix_04.cc b/tests/bits/complex_complex_full_matrix_04.cc new file mode 100644 index 0000000000..bc5240e110 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_04.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_04/output"; + + +template +void +check () +{ + const std::complex 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 > m (5,5,array); + + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_04/cmp/generic b/tests/bits/complex_complex_full_matrix_04/cmp/generic new file mode 100644 index 0000000000..a0e80cc63f --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_04/cmp/generic @@ -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_complex_full_matrix_05.cc b/tests/bits/complex_complex_full_matrix_05.cc new file mode 100644 index 0000000000..d49f30d839 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_05.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_05/output"; + + +template +void +check () +{ + FullMatrix > m (IdentityMatrix(5)); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_05/cmp/generic b/tests/bits/complex_complex_full_matrix_05/cmp/generic new file mode 100644 index 0000000000..e804e89061 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_05/cmp/generic @@ -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_complex_full_matrix_06.cc b/tests/bits/complex_complex_full_matrix_06.cc new file mode 100644 index 0000000000..8145cf8d89 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_06.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_06/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + + FullMatrix > n; + n = m; + + print_matrix (n); +} + diff --git a/tests/bits/complex_complex_full_matrix_06/cmp/generic b/tests/bits/complex_complex_full_matrix_06/cmp/generic new file mode 100644 index 0000000000..9d9b54f64c --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_06/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_07.cc b/tests/bits/complex_complex_full_matrix_07.cc new file mode 100644 index 0000000000..070c431e2b --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_07.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_07/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m = std::complex(0); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_07/cmp/generic b/tests/bits/complex_complex_full_matrix_07/cmp/generic new file mode 100644 index 0000000000..5b6383b352 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_07/cmp/generic @@ -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_complex_full_matrix_08.cc b/tests/bits/complex_complex_full_matrix_08.cc new file mode 100644 index 0000000000..9e0ffda965 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_08.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_08/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + + FullMatrix > n; + n.copy_from (m); + + print_matrix (n); +} + diff --git a/tests/bits/complex_complex_full_matrix_08/cmp/generic b/tests/bits/complex_complex_full_matrix_08/cmp/generic new file mode 100644 index 0000000000..9d9b54f64c --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_08/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_09.cc b/tests/bits/complex_complex_full_matrix_09.cc new file mode 100644 index 0000000000..4f63c95792 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_09.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_09/output"; + + +template +void +check () +{ + FullMatrix > m; + m = IdentityMatrix(5); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_09/cmp/generic b/tests/bits/complex_complex_full_matrix_09/cmp/generic new file mode 100644 index 0000000000..e804e89061 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_09/cmp/generic @@ -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_complex_full_matrix_10.cc b/tests/bits/complex_complex_full_matrix_10.cc new file mode 100644 index 0000000000..84508ef5e6 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_10.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_10/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + + FullMatrix > n (10, 10); + n.fill (m, 3, 4); + print_matrix (n); +} + diff --git a/tests/bits/complex_complex_full_matrix_10/cmp/generic b/tests/bits/complex_complex_full_matrix_10/cmp/generic new file mode 100644 index 0000000000..184403999a --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_10/cmp/generic @@ -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,12.) +DEAL::3 5 (3.0,15.) +DEAL::3 6 (4.0,18.) +DEAL::3 7 (5.0,21.) +DEAL::3 8 (6.0,24.) +DEAL::3 9 (0.0,0.0) +DEAL::4 0 (0.0,0.0) +DEAL::4 1 (0.0,0.0) +DEAL::4 2 (0.0,0.0) +DEAL::4 3 (0.0,0.0) +DEAL::4 4 (4.0,16.) +DEAL::4 5 (6.0,20.) +DEAL::4 6 (8.0,24.) +DEAL::4 7 (10.,28.) +DEAL::4 8 (12.,32.) +DEAL::4 9 (0.0,0.0) +DEAL::5 0 (0.0,0.0) +DEAL::5 1 (0.0,0.0) +DEAL::5 2 (0.0,0.0) +DEAL::5 3 (0.0,0.0) +DEAL::5 4 (6.0,20.) +DEAL::5 5 (9.0,25.) +DEAL::5 6 (12.,30.) +DEAL::5 7 (15.,35.) +DEAL::5 8 (18.,40.) +DEAL::5 9 (0.0,0.0) +DEAL::6 0 (0.0,0.0) +DEAL::6 1 (0.0,0.0) +DEAL::6 2 (0.0,0.0) +DEAL::6 3 (0.0,0.0) +DEAL::6 4 (8.0,24.) +DEAL::6 5 (12.,30.) +DEAL::6 6 (16.,36.) +DEAL::6 7 (20.,42.) +DEAL::6 8 (24.,48.) +DEAL::6 9 (0.0,0.0) +DEAL::7 0 (0.0,0.0) +DEAL::7 1 (0.0,0.0) +DEAL::7 2 (0.0,0.0) +DEAL::7 3 (0.0,0.0) +DEAL::7 4 (10.,28.) +DEAL::7 5 (15.,35.) +DEAL::7 6 (20.,42.) +DEAL::7 7 (25.,49.) +DEAL::7 8 (30.,56.) +DEAL::7 9 (0.0,0.0) +DEAL::8 0 (0.0,0.0) +DEAL::8 1 (0.0,0.0) +DEAL::8 2 (0.0,0.0) +DEAL::8 3 (0.0,0.0) +DEAL::8 4 (0.0,0.0) +DEAL::8 5 (0.0,0.0) +DEAL::8 6 (0.0,0.0) +DEAL::8 7 (0.0,0.0) +DEAL::8 8 (0.0,0.0) +DEAL::8 9 (0.0,0.0) +DEAL::9 0 (0.0,0.0) +DEAL::9 1 (0.0,0.0) +DEAL::9 2 (0.0,0.0) +DEAL::9 3 (0.0,0.0) +DEAL::9 4 (0.0,0.0) +DEAL::9 5 (0.0,0.0) +DEAL::9 6 (0.0,0.0) +DEAL::9 7 (0.0,0.0) +DEAL::9 8 (0.0,0.0) +DEAL::9 9 (0.0,0.0) +DEAL::0 0 (0.0,0.0) +DEAL::0 1 (0.0,0.0) +DEAL::0 2 (0.0,0.0) +DEAL::0 3 (0.0,0.0) +DEAL::0 4 (0.0,0.0) +DEAL::0 5 (0.0,0.0) +DEAL::0 6 (0.0,0.0) +DEAL::0 7 (0.0,0.0) +DEAL::0 8 (0.0,0.0) +DEAL::0 9 (0.0,0.0) +DEAL::1 0 (0.0,0.0) +DEAL::1 1 (0.0,0.0) +DEAL::1 2 (0.0,0.0) +DEAL::1 3 (0.0,0.0) +DEAL::1 4 (0.0,0.0) +DEAL::1 5 (0.0,0.0) +DEAL::1 6 (0.0,0.0) +DEAL::1 7 (0.0,0.0) +DEAL::1 8 (0.0,0.0) +DEAL::1 9 (0.0,0.0) +DEAL::2 0 (0.0,0.0) +DEAL::2 1 (0.0,0.0) +DEAL::2 2 (0.0,0.0) +DEAL::2 3 (0.0,0.0) +DEAL::2 4 (0.0,0.0) +DEAL::2 5 (0.0,0.0) +DEAL::2 6 (0.0,0.0) +DEAL::2 7 (0.0,0.0) +DEAL::2 8 (0.0,0.0) +DEAL::2 9 (0.0,0.0) +DEAL::3 0 (0.0,0.0) +DEAL::3 1 (0.0,0.0) +DEAL::3 2 (0.0,0.0) +DEAL::3 3 (0.0,0.0) +DEAL::3 4 (2.0,12.) +DEAL::3 5 (3.0,15.) +DEAL::3 6 (4.0,18.) +DEAL::3 7 (5.0,21.) +DEAL::3 8 (6.0,24.) +DEAL::3 9 (0.0,0.0) +DEAL::4 0 (0.0,0.0) +DEAL::4 1 (0.0,0.0) +DEAL::4 2 (0.0,0.0) +DEAL::4 3 (0.0,0.0) +DEAL::4 4 (4.0,16.) +DEAL::4 5 (6.0,20.) +DEAL::4 6 (8.0,24.) +DEAL::4 7 (10.,28.) +DEAL::4 8 (12.,32.) +DEAL::4 9 (0.0,0.0) +DEAL::5 0 (0.0,0.0) +DEAL::5 1 (0.0,0.0) +DEAL::5 2 (0.0,0.0) +DEAL::5 3 (0.0,0.0) +DEAL::5 4 (6.0,20.) +DEAL::5 5 (9.0,25.) +DEAL::5 6 (12.,30.) +DEAL::5 7 (15.,35.) +DEAL::5 8 (18.,40.) +DEAL::5 9 (0.0,0.0) +DEAL::6 0 (0.0,0.0) +DEAL::6 1 (0.0,0.0) +DEAL::6 2 (0.0,0.0) +DEAL::6 3 (0.0,0.0) +DEAL::6 4 (8.0,24.) +DEAL::6 5 (12.,30.) +DEAL::6 6 (16.,36.) +DEAL::6 7 (20.,42.) +DEAL::6 8 (24.,48.) +DEAL::6 9 (0.0,0.0) +DEAL::7 0 (0.0,0.0) +DEAL::7 1 (0.0,0.0) +DEAL::7 2 (0.0,0.0) +DEAL::7 3 (0.0,0.0) +DEAL::7 4 (10.,28.) +DEAL::7 5 (15.,35.) +DEAL::7 6 (20.,42.) +DEAL::7 7 (25.,49.) +DEAL::7 8 (30.,56.) +DEAL::7 9 (0.0,0.0) +DEAL::8 0 (0.0,0.0) +DEAL::8 1 (0.0,0.0) +DEAL::8 2 (0.0,0.0) +DEAL::8 3 (0.0,0.0) +DEAL::8 4 (0.0,0.0) +DEAL::8 5 (0.0,0.0) +DEAL::8 6 (0.0,0.0) +DEAL::8 7 (0.0,0.0) +DEAL::8 8 (0.0,0.0) +DEAL::8 9 (0.0,0.0) +DEAL::9 0 (0.0,0.0) +DEAL::9 1 (0.0,0.0) +DEAL::9 2 (0.0,0.0) +DEAL::9 3 (0.0,0.0) +DEAL::9 4 (0.0,0.0) +DEAL::9 5 (0.0,0.0) +DEAL::9 6 (0.0,0.0) +DEAL::9 7 (0.0,0.0) +DEAL::9 8 (0.0,0.0) +DEAL::9 9 (0.0,0.0) diff --git a/tests/bits/complex_complex_full_matrix_11.cc b/tests/bits/complex_complex_full_matrix_11.cc new file mode 100644 index 0000000000..630cad921f --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_11.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_11/output"; + + +template +void +check () +{ + FullMatrix > m, n; + make_complex_matrix (m); + make_complex_matrix (n); + deallog << (m==n) << std::endl; +} + diff --git a/tests/bits/complex_complex_full_matrix_11/cmp/generic b/tests/bits/complex_complex_full_matrix_11/cmp/generic new file mode 100644 index 0000000000..e46c0ac61e --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_11/cmp/generic @@ -0,0 +1,3 @@ + +DEAL::1 +DEAL::1 diff --git a/tests/bits/complex_complex_full_matrix_12.cc b/tests/bits/complex_complex_full_matrix_12.cc new file mode 100644 index 0000000000..8b507f1a07 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_12.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_12/output"; + + +template +void +check () +{ + FullMatrix > m(3,3); + deallog << m.all_zero() << std::endl; + + make_complex_matrix (m); + deallog << m.all_zero() << std::endl; +} + diff --git a/tests/bits/complex_complex_full_matrix_12/cmp/generic b/tests/bits/complex_complex_full_matrix_12/cmp/generic new file mode 100644 index 0000000000..de01f2c730 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_12/cmp/generic @@ -0,0 +1,5 @@ + +DEAL::1 +DEAL::0 +DEAL::1 +DEAL::0 diff --git a/tests/bits/complex_complex_full_matrix_13.cc b/tests/bits/complex_complex_full_matrix_13.cc new file mode 100644 index 0000000000..17aac560bf --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_13.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_13/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + Vector > v; + make_complex_vector (v); + + deallog << m.matrix_norm_square (v) << std::endl; +} + diff --git a/tests/bits/complex_complex_full_matrix_13/cmp/generic b/tests/bits/complex_complex_full_matrix_13/cmp/generic new file mode 100644 index 0000000000..e01042f17c --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_13/cmp/generic @@ -0,0 +1,3 @@ + +DEAL::(-3.3e+04,-1.1e+03) +DEAL::(-3.3e+04,-1.1e+03) diff --git a/tests/bits/complex_complex_full_matrix_14.cc b/tests/bits/complex_complex_full_matrix_14.cc new file mode 100644 index 0000000000..6219eda099 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_14.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_14/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + Vector > v,w; + make_complex_vector (v); + make_complex_vector (w); + for (unsigned int i=0; i(1.); + + deallog << m.matrix_scalar_product (v,w) << std::endl; +} + diff --git a/tests/bits/complex_complex_full_matrix_14/cmp/generic b/tests/bits/complex_complex_full_matrix_14/cmp/generic new file mode 100644 index 0000000000..4289f0caa7 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_14/cmp/generic @@ -0,0 +1,3 @@ + +DEAL::(-3.6e+04,3.2e+03) +DEAL::(-3.6e+04,3.2e+03) diff --git a/tests/bits/complex_complex_full_matrix_15.cc b/tests/bits/complex_complex_full_matrix_15.cc new file mode 100644 index 0000000000..dd7ce3bb97 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_15.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_15/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + deallog << m.l1_norm() << std::endl; +} + diff --git a/tests/bits/complex_complex_full_matrix_15/cmp/generic b/tests/bits/complex_complex_full_matrix_15/cmp/generic new file mode 100644 index 0000000000..37e590a8e0 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_15/cmp/generic @@ -0,0 +1,3 @@ + +DEAL::2.2e+02 +DEAL::2.2e+02 diff --git a/tests/bits/complex_complex_full_matrix_16.cc b/tests/bits/complex_complex_full_matrix_16.cc new file mode 100644 index 0000000000..c3535ce011 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_16.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_16/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + deallog << m.linfty_norm() << std::endl; +} + diff --git a/tests/bits/complex_complex_full_matrix_16/cmp/generic b/tests/bits/complex_complex_full_matrix_16/cmp/generic new file mode 100644 index 0000000000..a70842070c --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_16/cmp/generic @@ -0,0 +1,3 @@ + +DEAL::2.3e+02 +DEAL::2.3e+02 diff --git a/tests/bits/complex_complex_full_matrix_17.cc b/tests/bits/complex_complex_full_matrix_17.cc new file mode 100644 index 0000000000..0738721f41 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_17.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_17/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + deallog << m.frobenius_norm() << std::endl; +} + diff --git a/tests/bits/complex_complex_full_matrix_17/cmp/generic b/tests/bits/complex_complex_full_matrix_17/cmp/generic new file mode 100644 index 0000000000..6e2a22e7b5 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_17/cmp/generic @@ -0,0 +1,3 @@ + +DEAL::1.7e+02 +DEAL::1.7e+02 diff --git a/tests/bits/complex_complex_full_matrix_18.cc b/tests/bits/complex_complex_full_matrix_18.cc new file mode 100644 index 0000000000..4d4601d11a --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_18.cc @@ -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_complex_full_matrix_18/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + deallog << m.relative_symmetry_norm2() << std::endl; +} + diff --git a/tests/bits/complex_complex_full_matrix_18/cmp/generic b/tests/bits/complex_complex_full_matrix_18/cmp/generic new file mode 100644 index 0000000000..95e574c7d0 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_18/cmp/generic @@ -0,0 +1,3 @@ + +DEAL::0.081 +DEAL::0.081 diff --git a/tests/bits/complex_complex_full_matrix_19.cc b/tests/bits/complex_complex_full_matrix_19.cc new file mode 100644 index 0000000000..ab40aa4ba5 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_19.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_19/output"; + + +template +void +check () +{ + for (unsigned int n=1; n<=3; ++n) + { + const std::complex array[] = { 1, 2, 3, 4, 5, 6, 7, 8, 9 }; + + FullMatrix > m(n,n, array); + print_matrix (m); + deallog << m.determinant () << std::endl; + } +} + + diff --git a/tests/bits/complex_complex_full_matrix_19/cmp/generic b/tests/bits/complex_complex_full_matrix_19/cmp/generic new file mode 100644 index 0000000000..39f150c940 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_19/cmp/generic @@ -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_complex_full_matrix_20.cc b/tests/bits/complex_complex_full_matrix_20.cc new file mode 100644 index 0000000000..9c30ea371f --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_20.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_20/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + deallog << m.trace() << std::endl; +} + diff --git a/tests/bits/complex_complex_full_matrix_20/cmp/generic b/tests/bits/complex_complex_full_matrix_20/cmp/generic new file mode 100644 index 0000000000..3658e79966 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_20/cmp/generic @@ -0,0 +1,3 @@ + +DEAL::(70.,1.6e+02) +DEAL::(70.,1.6e+02) diff --git a/tests/bits/complex_complex_full_matrix_21.cc b/tests/bits/complex_complex_full_matrix_21.cc new file mode 100644 index 0000000000..070def092f --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_21.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_21/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m.print (deallog.get_file_stream()); +} + diff --git a/tests/bits/complex_complex_full_matrix_21/cmp/generic b/tests/bits/complex_complex_full_matrix_21/cmp/generic new file mode 100644 index 0000000000..16bdf08f56 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_21/cmp/generic @@ -0,0 +1,11 @@ + +(2.0,12.)(3.0,15.)(4.0,18.)(5.0,21.)(6.0,24.) +(4.0,16.)(6.0,20.)(8.0,24.)(10.,28.)(12.,32.) +(6.0,20.)(9.0,25.)(12.,30.)(15.,35.)(18.,40.) +(8.0,24.)(12.,30.)(16.,36.)(20.,42.)(24.,48.) +(10.,28.)(15.,35.)(20.,42.)(25.,49.)(30.,56.) +(2.0,12.)(3.0,15.)(4.0,18.)(5.0,21.)(6.0,24.) +(4.0,16.)(6.0,20.)(8.0,24.)(10.,28.)(12.,32.) +(6.0,20.)(9.0,25.)(12.,30.)(15.,35.)(18.,40.) +(8.0,24.)(12.,30.)(16.,36.)(20.,42.)(24.,48.) +(10.,28.)(15.,35.)(20.,42.)(25.,49.)(30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_22.cc b/tests/bits/complex_complex_full_matrix_22.cc new file mode 100644 index 0000000000..d3f1f65ad6 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_22.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_22/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m *= 3.1415; + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_22/cmp/generic b/tests/bits/complex_complex_full_matrix_22/cmp/generic new file mode 100644 index 0000000000..97bce9567b --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_22/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (6.3,38.) +DEAL::0 1 (9.4,47.) +DEAL::0 2 (13.,57.) +DEAL::0 3 (16.,66.) +DEAL::0 4 (19.,75.) +DEAL::1 0 (13.,50.) +DEAL::1 1 (19.,63.) +DEAL::1 2 (25.,75.) +DEAL::1 3 (31.,88.) +DEAL::1 4 (38.,1.0e+02) +DEAL::2 0 (19.,63.) +DEAL::2 1 (28.,79.) +DEAL::2 2 (38.,94.) +DEAL::2 3 (47.,1.1e+02) +DEAL::2 4 (57.,1.3e+02) +DEAL::3 0 (25.,75.) +DEAL::3 1 (38.,94.) +DEAL::3 2 (50.,1.1e+02) +DEAL::3 3 (63.,1.3e+02) +DEAL::3 4 (75.,1.5e+02) +DEAL::4 0 (31.,88.) +DEAL::4 1 (47.,1.1e+02) +DEAL::4 2 (63.,1.3e+02) +DEAL::4 3 (79.,1.5e+02) +DEAL::4 4 (94.,1.8e+02) +DEAL::0 0 (6.3,38.) +DEAL::0 1 (9.4,47.) +DEAL::0 2 (13.,57.) +DEAL::0 3 (16.,66.) +DEAL::0 4 (19.,75.) +DEAL::1 0 (13.,50.) +DEAL::1 1 (19.,63.) +DEAL::1 2 (25.,75.) +DEAL::1 3 (31.,88.) +DEAL::1 4 (38.,1.0e+02) +DEAL::2 0 (19.,63.) +DEAL::2 1 (28.,79.) +DEAL::2 2 (38.,94.) +DEAL::2 3 (47.,1.1e+02) +DEAL::2 4 (57.,1.3e+02) +DEAL::3 0 (25.,75.) +DEAL::3 1 (38.,94.) +DEAL::3 2 (50.,1.1e+02) +DEAL::3 3 (63.,1.3e+02) +DEAL::3 4 (75.,1.5e+02) +DEAL::4 0 (31.,88.) +DEAL::4 1 (47.,1.1e+02) +DEAL::4 2 (63.,1.3e+02) +DEAL::4 3 (79.,1.5e+02) +DEAL::4 4 (94.,1.8e+02) diff --git a/tests/bits/complex_complex_full_matrix_23.cc b/tests/bits/complex_complex_full_matrix_23.cc new file mode 100644 index 0000000000..b3bb4fe38a --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_23.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_23/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m /= 3.1415926; + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_23/cmp/generic b/tests/bits/complex_complex_full_matrix_23/cmp/generic new file mode 100644 index 0000000000..acfda45045 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_23/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (0.64,3.8) +DEAL::0 1 (0.95,4.8) +DEAL::0 2 (1.3,5.7) +DEAL::0 3 (1.6,6.7) +DEAL::0 4 (1.9,7.6) +DEAL::1 0 (1.3,5.1) +DEAL::1 1 (1.9,6.4) +DEAL::1 2 (2.5,7.6) +DEAL::1 3 (3.2,8.9) +DEAL::1 4 (3.8,10.) +DEAL::2 0 (1.9,6.4) +DEAL::2 1 (2.9,8.0) +DEAL::2 2 (3.8,9.5) +DEAL::2 3 (4.8,11.) +DEAL::2 4 (5.7,13.) +DEAL::3 0 (2.5,7.6) +DEAL::3 1 (3.8,9.5) +DEAL::3 2 (5.1,11.) +DEAL::3 3 (6.4,13.) +DEAL::3 4 (7.6,15.) +DEAL::4 0 (3.2,8.9) +DEAL::4 1 (4.8,11.) +DEAL::4 2 (6.4,13.) +DEAL::4 3 (8.0,16.) +DEAL::4 4 (9.5,18.) +DEAL::0 0 (0.64,3.8) +DEAL::0 1 (0.95,4.8) +DEAL::0 2 (1.3,5.7) +DEAL::0 3 (1.6,6.7) +DEAL::0 4 (1.9,7.6) +DEAL::1 0 (1.3,5.1) +DEAL::1 1 (1.9,6.4) +DEAL::1 2 (2.5,7.6) +DEAL::1 3 (3.2,8.9) +DEAL::1 4 (3.8,10.) +DEAL::2 0 (1.9,6.4) +DEAL::2 1 (2.9,8.0) +DEAL::2 2 (3.8,9.5) +DEAL::2 3 (4.8,11.) +DEAL::2 4 (5.7,13.) +DEAL::3 0 (2.5,7.6) +DEAL::3 1 (3.8,9.5) +DEAL::3 2 (5.1,11.) +DEAL::3 3 (6.4,13.) +DEAL::3 4 (7.6,15.) +DEAL::4 0 (3.2,8.9) +DEAL::4 1 (4.8,11.) +DEAL::4 2 (6.4,13.) +DEAL::4 3 (8.0,16.) +DEAL::4 4 (9.5,18.) diff --git a/tests/bits/complex_complex_full_matrix_24.cc b/tests/bits/complex_complex_full_matrix_24.cc new file mode 100644 index 0000000000..814681d67a --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_24.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_24/output"; + + +template +void +check () +{ + FullMatrix > m,n; + make_complex_matrix (m); + make_complex_matrix (n); + m.add (3.14159, n); + + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_24/cmp/generic b/tests/bits/complex_complex_full_matrix_24/cmp/generic new file mode 100644 index 0000000000..f3848e7b88 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_24/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (8.3,50.) +DEAL::0 1 (12.,62.) +DEAL::0 2 (17.,75.) +DEAL::0 3 (21.,87.) +DEAL::0 4 (25.,99.) +DEAL::1 0 (17.,66.) +DEAL::1 1 (25.,83.) +DEAL::1 2 (33.,99.) +DEAL::1 3 (41.,1.2e+02) +DEAL::1 4 (50.,1.3e+02) +DEAL::2 0 (25.,83.) +DEAL::2 1 (37.,1.0e+02) +DEAL::2 2 (50.,1.2e+02) +DEAL::2 3 (62.,1.4e+02) +DEAL::2 4 (75.,1.7e+02) +DEAL::3 0 (33.,99.) +DEAL::3 1 (50.,1.2e+02) +DEAL::3 2 (66.,1.5e+02) +DEAL::3 3 (83.,1.7e+02) +DEAL::3 4 (99.,2.0e+02) +DEAL::4 0 (41.,1.2e+02) +DEAL::4 1 (62.,1.4e+02) +DEAL::4 2 (83.,1.7e+02) +DEAL::4 3 (1.0e+02,2.0e+02) +DEAL::4 4 (1.2e+02,2.3e+02) +DEAL::0 0 (8.3,50.) +DEAL::0 1 (12.,62.) +DEAL::0 2 (17.,75.) +DEAL::0 3 (21.,87.) +DEAL::0 4 (25.,99.) +DEAL::1 0 (17.,66.) +DEAL::1 1 (25.,83.) +DEAL::1 2 (33.,99.) +DEAL::1 3 (41.,1.2e+02) +DEAL::1 4 (50.,1.3e+02) +DEAL::2 0 (25.,83.) +DEAL::2 1 (37.,1.0e+02) +DEAL::2 2 (50.,1.2e+02) +DEAL::2 3 (62.,1.4e+02) +DEAL::2 4 (75.,1.7e+02) +DEAL::3 0 (33.,99.) +DEAL::3 1 (50.,1.2e+02) +DEAL::3 2 (66.,1.5e+02) +DEAL::3 3 (83.,1.7e+02) +DEAL::3 4 (99.,2.0e+02) +DEAL::4 0 (41.,1.2e+02) +DEAL::4 1 (62.,1.4e+02) +DEAL::4 2 (83.,1.7e+02) +DEAL::4 3 (1.0e+02,2.0e+02) +DEAL::4 4 (1.2e+02,2.3e+02) diff --git a/tests/bits/complex_complex_full_matrix_25.cc b/tests/bits/complex_complex_full_matrix_25.cc new file mode 100644 index 0000000000..6ab1c8f3e2 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_25.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_25/output"; + + +template +void +check () +{ + FullMatrix > m,n,o; + make_complex_matrix (m); + make_complex_matrix (n); + make_complex_matrix (o); + m.add (3.1415, n, 2.718, o); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_25/cmp/generic b/tests/bits/complex_complex_full_matrix_25/cmp/generic new file mode 100644 index 0000000000..46ce3ab9f9 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_25/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (14.,82.) +DEAL::0 1 (21.,1.0e+02) +DEAL::0 2 (27.,1.2e+02) +DEAL::0 3 (34.,1.4e+02) +DEAL::0 4 (41.,1.6e+02) +DEAL::1 0 (27.,1.1e+02) +DEAL::1 1 (41.,1.4e+02) +DEAL::1 2 (55.,1.6e+02) +DEAL::1 3 (69.,1.9e+02) +DEAL::1 4 (82.,2.2e+02) +DEAL::2 0 (41.,1.4e+02) +DEAL::2 1 (62.,1.7e+02) +DEAL::2 2 (82.,2.1e+02) +DEAL::2 3 (1.0e+02,2.4e+02) +DEAL::2 4 (1.2e+02,2.7e+02) +DEAL::3 0 (55.,1.6e+02) +DEAL::3 1 (82.,2.1e+02) +DEAL::3 2 (1.1e+02,2.5e+02) +DEAL::3 3 (1.4e+02,2.9e+02) +DEAL::3 4 (1.6e+02,3.3e+02) +DEAL::4 0 (69.,1.9e+02) +DEAL::4 1 (1.0e+02,2.4e+02) +DEAL::4 2 (1.4e+02,2.9e+02) +DEAL::4 3 (1.7e+02,3.4e+02) +DEAL::4 4 (2.1e+02,3.8e+02) +DEAL::0 0 (14.,82.) +DEAL::0 1 (21.,1.0e+02) +DEAL::0 2 (27.,1.2e+02) +DEAL::0 3 (34.,1.4e+02) +DEAL::0 4 (41.,1.6e+02) +DEAL::1 0 (27.,1.1e+02) +DEAL::1 1 (41.,1.4e+02) +DEAL::1 2 (55.,1.6e+02) +DEAL::1 3 (69.,1.9e+02) +DEAL::1 4 (82.,2.2e+02) +DEAL::2 0 (41.,1.4e+02) +DEAL::2 1 (62.,1.7e+02) +DEAL::2 2 (82.,2.1e+02) +DEAL::2 3 (1.0e+02,2.4e+02) +DEAL::2 4 (1.2e+02,2.7e+02) +DEAL::3 0 (55.,1.6e+02) +DEAL::3 1 (82.,2.1e+02) +DEAL::3 2 (1.1e+02,2.5e+02) +DEAL::3 3 (1.4e+02,2.9e+02) +DEAL::3 4 (1.6e+02,3.3e+02) +DEAL::4 0 (69.,1.9e+02) +DEAL::4 1 (1.0e+02,2.4e+02) +DEAL::4 2 (1.4e+02,2.9e+02) +DEAL::4 3 (1.7e+02,3.4e+02) +DEAL::4 4 (2.1e+02,3.8e+02) diff --git a/tests/bits/complex_complex_full_matrix_26.cc b/tests/bits/complex_complex_full_matrix_26.cc new file mode 100644 index 0000000000..4890a8a270 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_26.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_26/output"; + + +template +void +check () +{ + FullMatrix > m,n,o,p; + make_complex_matrix (m); + make_complex_matrix (n); + make_complex_matrix (o); + make_complex_matrix (p); + m.add (3.1415, n, 2.718, o, 1.414, p); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_26/cmp/generic b/tests/bits/complex_complex_full_matrix_26/cmp/generic new file mode 100644 index 0000000000..7a069f7dde --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_26/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (17.,99.) +DEAL::0 1 (25.,1.2e+02) +DEAL::0 2 (33.,1.5e+02) +DEAL::0 3 (41.,1.7e+02) +DEAL::0 4 (50.,2.0e+02) +DEAL::1 0 (33.,1.3e+02) +DEAL::1 1 (50.,1.7e+02) +DEAL::1 2 (66.,2.0e+02) +DEAL::1 3 (83.,2.3e+02) +DEAL::1 4 (99.,2.6e+02) +DEAL::2 0 (50.,1.7e+02) +DEAL::2 1 (74.,2.1e+02) +DEAL::2 2 (99.,2.5e+02) +DEAL::2 3 (1.2e+02,2.9e+02) +DEAL::2 4 (1.5e+02,3.3e+02) +DEAL::3 0 (66.,2.0e+02) +DEAL::3 1 (99.,2.5e+02) +DEAL::3 2 (1.3e+02,3.0e+02) +DEAL::3 3 (1.7e+02,3.5e+02) +DEAL::3 4 (2.0e+02,4.0e+02) +DEAL::4 0 (83.,2.3e+02) +DEAL::4 1 (1.2e+02,2.9e+02) +DEAL::4 2 (1.7e+02,3.5e+02) +DEAL::4 3 (2.1e+02,4.1e+02) +DEAL::4 4 (2.5e+02,4.6e+02) +DEAL::0 0 (17.,99.) +DEAL::0 1 (25.,1.2e+02) +DEAL::0 2 (33.,1.5e+02) +DEAL::0 3 (41.,1.7e+02) +DEAL::0 4 (50.,2.0e+02) +DEAL::1 0 (33.,1.3e+02) +DEAL::1 1 (50.,1.7e+02) +DEAL::1 2 (66.,2.0e+02) +DEAL::1 3 (83.,2.3e+02) +DEAL::1 4 (99.,2.6e+02) +DEAL::2 0 (50.,1.7e+02) +DEAL::2 1 (74.,2.1e+02) +DEAL::2 2 (99.,2.5e+02) +DEAL::2 3 (1.2e+02,2.9e+02) +DEAL::2 4 (1.5e+02,3.3e+02) +DEAL::3 0 (66.,2.0e+02) +DEAL::3 1 (99.,2.5e+02) +DEAL::3 2 (1.3e+02,3.0e+02) +DEAL::3 3 (1.7e+02,3.5e+02) +DEAL::3 4 (2.0e+02,4.0e+02) +DEAL::4 0 (83.,2.3e+02) +DEAL::4 1 (1.2e+02,2.9e+02) +DEAL::4 2 (1.7e+02,3.5e+02) +DEAL::4 3 (2.1e+02,4.1e+02) +DEAL::4 4 (2.5e+02,4.6e+02) diff --git a/tests/bits/complex_complex_full_matrix_27.cc b/tests/bits/complex_complex_full_matrix_27.cc new file mode 100644 index 0000000000..28961aee57 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_27.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_27/output"; + + +template +void +check () +{ + FullMatrix > m(10,10),n; + make_complex_matrix (n); + m.add (n, 2., 3, 4, 2, 3); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_27/cmp/generic b/tests/bits/complex_complex_full_matrix_27/cmp/generic new file mode 100644 index 0000000000..af6eac0d0e --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_27/cmp/generic @@ -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.,70.) +DEAL::3 5 (36.,80.) +DEAL::3 6 (0.0,0.0) +DEAL::3 7 (0.0,0.0) +DEAL::3 8 (0.0,0.0) +DEAL::3 9 (0.0,0.0) +DEAL::4 0 (0.0,0.0) +DEAL::4 1 (0.0,0.0) +DEAL::4 2 (0.0,0.0) +DEAL::4 3 (0.0,0.0) +DEAL::4 4 (40.,84.) +DEAL::4 5 (48.,96.) +DEAL::4 6 (0.0,0.0) +DEAL::4 7 (0.0,0.0) +DEAL::4 8 (0.0,0.0) +DEAL::4 9 (0.0,0.0) +DEAL::5 0 (0.0,0.0) +DEAL::5 1 (0.0,0.0) +DEAL::5 2 (0.0,0.0) +DEAL::5 3 (0.0,0.0) +DEAL::5 4 (50.,98.) +DEAL::5 5 (60.,1.1e+02) +DEAL::5 6 (0.0,0.0) +DEAL::5 7 (0.0,0.0) +DEAL::5 8 (0.0,0.0) +DEAL::5 9 (0.0,0.0) +DEAL::6 0 (0.0,0.0) +DEAL::6 1 (0.0,0.0) +DEAL::6 2 (0.0,0.0) +DEAL::6 3 (0.0,0.0) +DEAL::6 4 (0.0,0.0) +DEAL::6 5 (0.0,0.0) +DEAL::6 6 (0.0,0.0) +DEAL::6 7 (0.0,0.0) +DEAL::6 8 (0.0,0.0) +DEAL::6 9 (0.0,0.0) +DEAL::7 0 (0.0,0.0) +DEAL::7 1 (0.0,0.0) +DEAL::7 2 (0.0,0.0) +DEAL::7 3 (0.0,0.0) +DEAL::7 4 (0.0,0.0) +DEAL::7 5 (0.0,0.0) +DEAL::7 6 (0.0,0.0) +DEAL::7 7 (0.0,0.0) +DEAL::7 8 (0.0,0.0) +DEAL::7 9 (0.0,0.0) +DEAL::8 0 (0.0,0.0) +DEAL::8 1 (0.0,0.0) +DEAL::8 2 (0.0,0.0) +DEAL::8 3 (0.0,0.0) +DEAL::8 4 (0.0,0.0) +DEAL::8 5 (0.0,0.0) +DEAL::8 6 (0.0,0.0) +DEAL::8 7 (0.0,0.0) +DEAL::8 8 (0.0,0.0) +DEAL::8 9 (0.0,0.0) +DEAL::9 0 (0.0,0.0) +DEAL::9 1 (0.0,0.0) +DEAL::9 2 (0.0,0.0) +DEAL::9 3 (0.0,0.0) +DEAL::9 4 (0.0,0.0) +DEAL::9 5 (0.0,0.0) +DEAL::9 6 (0.0,0.0) +DEAL::9 7 (0.0,0.0) +DEAL::9 8 (0.0,0.0) +DEAL::9 9 (0.0,0.0) +DEAL::0 0 (0.0,0.0) +DEAL::0 1 (0.0,0.0) +DEAL::0 2 (0.0,0.0) +DEAL::0 3 (0.0,0.0) +DEAL::0 4 (0.0,0.0) +DEAL::0 5 (0.0,0.0) +DEAL::0 6 (0.0,0.0) +DEAL::0 7 (0.0,0.0) +DEAL::0 8 (0.0,0.0) +DEAL::0 9 (0.0,0.0) +DEAL::1 0 (0.0,0.0) +DEAL::1 1 (0.0,0.0) +DEAL::1 2 (0.0,0.0) +DEAL::1 3 (0.0,0.0) +DEAL::1 4 (0.0,0.0) +DEAL::1 5 (0.0,0.0) +DEAL::1 6 (0.0,0.0) +DEAL::1 7 (0.0,0.0) +DEAL::1 8 (0.0,0.0) +DEAL::1 9 (0.0,0.0) +DEAL::2 0 (0.0,0.0) +DEAL::2 1 (0.0,0.0) +DEAL::2 2 (0.0,0.0) +DEAL::2 3 (0.0,0.0) +DEAL::2 4 (0.0,0.0) +DEAL::2 5 (0.0,0.0) +DEAL::2 6 (0.0,0.0) +DEAL::2 7 (0.0,0.0) +DEAL::2 8 (0.0,0.0) +DEAL::2 9 (0.0,0.0) +DEAL::3 0 (0.0,0.0) +DEAL::3 1 (0.0,0.0) +DEAL::3 2 (0.0,0.0) +DEAL::3 3 (0.0,0.0) +DEAL::3 4 (30.,70.) +DEAL::3 5 (36.,80.) +DEAL::3 6 (0.0,0.0) +DEAL::3 7 (0.0,0.0) +DEAL::3 8 (0.0,0.0) +DEAL::3 9 (0.0,0.0) +DEAL::4 0 (0.0,0.0) +DEAL::4 1 (0.0,0.0) +DEAL::4 2 (0.0,0.0) +DEAL::4 3 (0.0,0.0) +DEAL::4 4 (40.,84.) +DEAL::4 5 (48.,96.) +DEAL::4 6 (0.0,0.0) +DEAL::4 7 (0.0,0.0) +DEAL::4 8 (0.0,0.0) +DEAL::4 9 (0.0,0.0) +DEAL::5 0 (0.0,0.0) +DEAL::5 1 (0.0,0.0) +DEAL::5 2 (0.0,0.0) +DEAL::5 3 (0.0,0.0) +DEAL::5 4 (50.,98.) +DEAL::5 5 (60.,1.1e+02) +DEAL::5 6 (0.0,0.0) +DEAL::5 7 (0.0,0.0) +DEAL::5 8 (0.0,0.0) +DEAL::5 9 (0.0,0.0) +DEAL::6 0 (0.0,0.0) +DEAL::6 1 (0.0,0.0) +DEAL::6 2 (0.0,0.0) +DEAL::6 3 (0.0,0.0) +DEAL::6 4 (0.0,0.0) +DEAL::6 5 (0.0,0.0) +DEAL::6 6 (0.0,0.0) +DEAL::6 7 (0.0,0.0) +DEAL::6 8 (0.0,0.0) +DEAL::6 9 (0.0,0.0) +DEAL::7 0 (0.0,0.0) +DEAL::7 1 (0.0,0.0) +DEAL::7 2 (0.0,0.0) +DEAL::7 3 (0.0,0.0) +DEAL::7 4 (0.0,0.0) +DEAL::7 5 (0.0,0.0) +DEAL::7 6 (0.0,0.0) +DEAL::7 7 (0.0,0.0) +DEAL::7 8 (0.0,0.0) +DEAL::7 9 (0.0,0.0) +DEAL::8 0 (0.0,0.0) +DEAL::8 1 (0.0,0.0) +DEAL::8 2 (0.0,0.0) +DEAL::8 3 (0.0,0.0) +DEAL::8 4 (0.0,0.0) +DEAL::8 5 (0.0,0.0) +DEAL::8 6 (0.0,0.0) +DEAL::8 7 (0.0,0.0) +DEAL::8 8 (0.0,0.0) +DEAL::8 9 (0.0,0.0) +DEAL::9 0 (0.0,0.0) +DEAL::9 1 (0.0,0.0) +DEAL::9 2 (0.0,0.0) +DEAL::9 3 (0.0,0.0) +DEAL::9 4 (0.0,0.0) +DEAL::9 5 (0.0,0.0) +DEAL::9 6 (0.0,0.0) +DEAL::9 7 (0.0,0.0) +DEAL::9 8 (0.0,0.0) +DEAL::9 9 (0.0,0.0) diff --git a/tests/bits/complex_complex_full_matrix_28.cc b/tests/bits/complex_complex_full_matrix_28.cc new file mode 100644 index 0000000000..4037afbc9b --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_28.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_28/output"; + + +template +void +check () +{ + FullMatrix > m(10,10),n; + make_complex_matrix (n); + m.Tadd (n, 2., 3, 4, 2, 3); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_28/cmp/generic b/tests/bits/complex_complex_full_matrix_28/cmp/generic new file mode 100644 index 0000000000..83e9a512e8 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_28/cmp/generic @@ -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.,70.) +DEAL::3 5 (40.,84.) +DEAL::3 6 (50.,98.) +DEAL::3 7 (0.0,0.0) +DEAL::3 8 (0.0,0.0) +DEAL::3 9 (0.0,0.0) +DEAL::4 0 (0.0,0.0) +DEAL::4 1 (0.0,0.0) +DEAL::4 2 (0.0,0.0) +DEAL::4 3 (0.0,0.0) +DEAL::4 4 (36.,80.) +DEAL::4 5 (48.,96.) +DEAL::4 6 (60.,1.1e+02) +DEAL::4 7 (0.0,0.0) +DEAL::4 8 (0.0,0.0) +DEAL::4 9 (0.0,0.0) +DEAL::5 0 (0.0,0.0) +DEAL::5 1 (0.0,0.0) +DEAL::5 2 (0.0,0.0) +DEAL::5 3 (0.0,0.0) +DEAL::5 4 (0.0,0.0) +DEAL::5 5 (0.0,0.0) +DEAL::5 6 (0.0,0.0) +DEAL::5 7 (0.0,0.0) +DEAL::5 8 (0.0,0.0) +DEAL::5 9 (0.0,0.0) +DEAL::6 0 (0.0,0.0) +DEAL::6 1 (0.0,0.0) +DEAL::6 2 (0.0,0.0) +DEAL::6 3 (0.0,0.0) +DEAL::6 4 (0.0,0.0) +DEAL::6 5 (0.0,0.0) +DEAL::6 6 (0.0,0.0) +DEAL::6 7 (0.0,0.0) +DEAL::6 8 (0.0,0.0) +DEAL::6 9 (0.0,0.0) +DEAL::7 0 (0.0,0.0) +DEAL::7 1 (0.0,0.0) +DEAL::7 2 (0.0,0.0) +DEAL::7 3 (0.0,0.0) +DEAL::7 4 (0.0,0.0) +DEAL::7 5 (0.0,0.0) +DEAL::7 6 (0.0,0.0) +DEAL::7 7 (0.0,0.0) +DEAL::7 8 (0.0,0.0) +DEAL::7 9 (0.0,0.0) +DEAL::8 0 (0.0,0.0) +DEAL::8 1 (0.0,0.0) +DEAL::8 2 (0.0,0.0) +DEAL::8 3 (0.0,0.0) +DEAL::8 4 (0.0,0.0) +DEAL::8 5 (0.0,0.0) +DEAL::8 6 (0.0,0.0) +DEAL::8 7 (0.0,0.0) +DEAL::8 8 (0.0,0.0) +DEAL::8 9 (0.0,0.0) +DEAL::9 0 (0.0,0.0) +DEAL::9 1 (0.0,0.0) +DEAL::9 2 (0.0,0.0) +DEAL::9 3 (0.0,0.0) +DEAL::9 4 (0.0,0.0) +DEAL::9 5 (0.0,0.0) +DEAL::9 6 (0.0,0.0) +DEAL::9 7 (0.0,0.0) +DEAL::9 8 (0.0,0.0) +DEAL::9 9 (0.0,0.0) +DEAL::0 0 (0.0,0.0) +DEAL::0 1 (0.0,0.0) +DEAL::0 2 (0.0,0.0) +DEAL::0 3 (0.0,0.0) +DEAL::0 4 (0.0,0.0) +DEAL::0 5 (0.0,0.0) +DEAL::0 6 (0.0,0.0) +DEAL::0 7 (0.0,0.0) +DEAL::0 8 (0.0,0.0) +DEAL::0 9 (0.0,0.0) +DEAL::1 0 (0.0,0.0) +DEAL::1 1 (0.0,0.0) +DEAL::1 2 (0.0,0.0) +DEAL::1 3 (0.0,0.0) +DEAL::1 4 (0.0,0.0) +DEAL::1 5 (0.0,0.0) +DEAL::1 6 (0.0,0.0) +DEAL::1 7 (0.0,0.0) +DEAL::1 8 (0.0,0.0) +DEAL::1 9 (0.0,0.0) +DEAL::2 0 (0.0,0.0) +DEAL::2 1 (0.0,0.0) +DEAL::2 2 (0.0,0.0) +DEAL::2 3 (0.0,0.0) +DEAL::2 4 (0.0,0.0) +DEAL::2 5 (0.0,0.0) +DEAL::2 6 (0.0,0.0) +DEAL::2 7 (0.0,0.0) +DEAL::2 8 (0.0,0.0) +DEAL::2 9 (0.0,0.0) +DEAL::3 0 (0.0,0.0) +DEAL::3 1 (0.0,0.0) +DEAL::3 2 (0.0,0.0) +DEAL::3 3 (0.0,0.0) +DEAL::3 4 (30.,70.) +DEAL::3 5 (40.,84.) +DEAL::3 6 (50.,98.) +DEAL::3 7 (0.0,0.0) +DEAL::3 8 (0.0,0.0) +DEAL::3 9 (0.0,0.0) +DEAL::4 0 (0.0,0.0) +DEAL::4 1 (0.0,0.0) +DEAL::4 2 (0.0,0.0) +DEAL::4 3 (0.0,0.0) +DEAL::4 4 (36.,80.) +DEAL::4 5 (48.,96.) +DEAL::4 6 (60.,1.1e+02) +DEAL::4 7 (0.0,0.0) +DEAL::4 8 (0.0,0.0) +DEAL::4 9 (0.0,0.0) +DEAL::5 0 (0.0,0.0) +DEAL::5 1 (0.0,0.0) +DEAL::5 2 (0.0,0.0) +DEAL::5 3 (0.0,0.0) +DEAL::5 4 (0.0,0.0) +DEAL::5 5 (0.0,0.0) +DEAL::5 6 (0.0,0.0) +DEAL::5 7 (0.0,0.0) +DEAL::5 8 (0.0,0.0) +DEAL::5 9 (0.0,0.0) +DEAL::6 0 (0.0,0.0) +DEAL::6 1 (0.0,0.0) +DEAL::6 2 (0.0,0.0) +DEAL::6 3 (0.0,0.0) +DEAL::6 4 (0.0,0.0) +DEAL::6 5 (0.0,0.0) +DEAL::6 6 (0.0,0.0) +DEAL::6 7 (0.0,0.0) +DEAL::6 8 (0.0,0.0) +DEAL::6 9 (0.0,0.0) +DEAL::7 0 (0.0,0.0) +DEAL::7 1 (0.0,0.0) +DEAL::7 2 (0.0,0.0) +DEAL::7 3 (0.0,0.0) +DEAL::7 4 (0.0,0.0) +DEAL::7 5 (0.0,0.0) +DEAL::7 6 (0.0,0.0) +DEAL::7 7 (0.0,0.0) +DEAL::7 8 (0.0,0.0) +DEAL::7 9 (0.0,0.0) +DEAL::8 0 (0.0,0.0) +DEAL::8 1 (0.0,0.0) +DEAL::8 2 (0.0,0.0) +DEAL::8 3 (0.0,0.0) +DEAL::8 4 (0.0,0.0) +DEAL::8 5 (0.0,0.0) +DEAL::8 6 (0.0,0.0) +DEAL::8 7 (0.0,0.0) +DEAL::8 8 (0.0,0.0) +DEAL::8 9 (0.0,0.0) +DEAL::9 0 (0.0,0.0) +DEAL::9 1 (0.0,0.0) +DEAL::9 2 (0.0,0.0) +DEAL::9 3 (0.0,0.0) +DEAL::9 4 (0.0,0.0) +DEAL::9 5 (0.0,0.0) +DEAL::9 6 (0.0,0.0) +DEAL::9 7 (0.0,0.0) +DEAL::9 8 (0.0,0.0) +DEAL::9 9 (0.0,0.0) diff --git a/tests/bits/complex_complex_full_matrix_29.cc b/tests/bits/complex_complex_full_matrix_29.cc new file mode 100644 index 0000000000..2fefe0c776 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_29.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_29/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m.add_row (2, 3.141, 4); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_29/cmp/generic b/tests/bits/complex_complex_full_matrix_29/cmp/generic new file mode 100644 index 0000000000..a805b1b0f7 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_29/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (37.,1.1e+02) +DEAL::2 1 (56.,1.3e+02) +DEAL::2 2 (75.,1.6e+02) +DEAL::2 3 (94.,1.9e+02) +DEAL::2 4 (1.1e+02,2.2e+02) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (37.,1.1e+02) +DEAL::2 1 (56.,1.3e+02) +DEAL::2 2 (75.,1.6e+02) +DEAL::2 3 (94.,1.9e+02) +DEAL::2 4 (1.1e+02,2.2e+02) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_30.cc b/tests/bits/complex_complex_full_matrix_30.cc new file mode 100644 index 0000000000..ea44c0f453 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_30.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_30/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m.add_row (2, 3.141, 4, 2.7182, 3); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_30/cmp/generic b/tests/bits/complex_complex_full_matrix_30/cmp/generic new file mode 100644 index 0000000000..64ec81cb70 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_30/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (59.,1.7e+02) +DEAL::2 1 (89.,2.2e+02) +DEAL::2 2 (1.2e+02,2.6e+02) +DEAL::2 3 (1.5e+02,3.0e+02) +DEAL::2 4 (1.8e+02,3.5e+02) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (59.,1.7e+02) +DEAL::2 1 (89.,2.2e+02) +DEAL::2 2 (1.2e+02,2.6e+02) +DEAL::2 3 (1.5e+02,3.0e+02) +DEAL::2 4 (1.8e+02,3.5e+02) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_31.cc b/tests/bits/complex_complex_full_matrix_31.cc new file mode 100644 index 0000000000..f69c0fd6c1 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_31.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_31/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m.add_col (2, 3.141, 4); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_31/cmp/generic b/tests/bits/complex_complex_full_matrix_31/cmp/generic new file mode 100644 index 0000000000..e631f5bbf9 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_31/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (23.,93.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (46.,1.2e+02) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (69.,1.6e+02) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (91.,1.9e+02) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (1.1e+02,2.2e+02) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (23.,93.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (46.,1.2e+02) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (69.,1.6e+02) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (91.,1.9e+02) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (1.1e+02,2.2e+02) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_32.cc b/tests/bits/complex_complex_full_matrix_32.cc new file mode 100644 index 0000000000..a8e479a4a6 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_32.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_32/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m.add_col (2, 3.141, 4, 2.718, 3); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_32/cmp/generic b/tests/bits/complex_complex_full_matrix_32/cmp/generic new file mode 100644 index 0000000000..3bb4ab8838 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_32/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (36.,1.5e+02) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (73.,2.0e+02) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (1.1e+02,2.5e+02) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (1.5e+02,3.0e+02) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (1.8e+02,3.5e+02) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (36.,1.5e+02) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (73.,2.0e+02) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (1.1e+02,2.5e+02) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (1.5e+02,3.0e+02) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (1.8e+02,3.5e+02) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_33.cc b/tests/bits/complex_complex_full_matrix_33.cc new file mode 100644 index 0000000000..a06cfa5cd1 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_33.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_33/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m.swap_row (2, 4); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_33/cmp/generic b/tests/bits/complex_complex_full_matrix_33/cmp/generic new file mode 100644 index 0000000000..03dfd43261 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_33/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (10.,28.) +DEAL::2 1 (15.,35.) +DEAL::2 2 (20.,42.) +DEAL::2 3 (25.,49.) +DEAL::2 4 (30.,56.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (6.0,20.) +DEAL::4 1 (9.0,25.) +DEAL::4 2 (12.,30.) +DEAL::4 3 (15.,35.) +DEAL::4 4 (18.,40.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (10.,28.) +DEAL::2 1 (15.,35.) +DEAL::2 2 (20.,42.) +DEAL::2 3 (25.,49.) +DEAL::2 4 (30.,56.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (6.0,20.) +DEAL::4 1 (9.0,25.) +DEAL::4 2 (12.,30.) +DEAL::4 3 (15.,35.) +DEAL::4 4 (18.,40.) diff --git a/tests/bits/complex_complex_full_matrix_34.cc b/tests/bits/complex_complex_full_matrix_34.cc new file mode 100644 index 0000000000..5499c7e1b1 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_34.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_34/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m.swap_col (2, 4); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_34/cmp/generic b/tests/bits/complex_complex_full_matrix_34/cmp/generic new file mode 100644 index 0000000000..1ab9dfaa34 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_34/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (6.0,24.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (4.0,18.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (12.,32.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (8.0,24.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (18.,40.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (12.,30.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (24.,48.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (16.,36.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (30.,56.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (20.,42.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (6.0,24.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (4.0,18.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (12.,32.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (8.0,24.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (18.,40.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (12.,30.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (24.,48.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (16.,36.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (30.,56.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (20.,42.) diff --git a/tests/bits/complex_complex_full_matrix_35.cc b/tests/bits/complex_complex_full_matrix_35.cc new file mode 100644 index 0000000000..399b958a67 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_35.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_35/output"; + + +template +void +check () +{ + FullMatrix > m, n; + make_complex_matrix (m); + make_complex_matrix (n); + m.add_diag (3.141, m); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_35/cmp/generic b/tests/bits/complex_complex_full_matrix_35/cmp/generic new file mode 100644 index 0000000000..f3848e7b88 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_35/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (8.3,50.) +DEAL::0 1 (12.,62.) +DEAL::0 2 (17.,75.) +DEAL::0 3 (21.,87.) +DEAL::0 4 (25.,99.) +DEAL::1 0 (17.,66.) +DEAL::1 1 (25.,83.) +DEAL::1 2 (33.,99.) +DEAL::1 3 (41.,1.2e+02) +DEAL::1 4 (50.,1.3e+02) +DEAL::2 0 (25.,83.) +DEAL::2 1 (37.,1.0e+02) +DEAL::2 2 (50.,1.2e+02) +DEAL::2 3 (62.,1.4e+02) +DEAL::2 4 (75.,1.7e+02) +DEAL::3 0 (33.,99.) +DEAL::3 1 (50.,1.2e+02) +DEAL::3 2 (66.,1.5e+02) +DEAL::3 3 (83.,1.7e+02) +DEAL::3 4 (99.,2.0e+02) +DEAL::4 0 (41.,1.2e+02) +DEAL::4 1 (62.,1.4e+02) +DEAL::4 2 (83.,1.7e+02) +DEAL::4 3 (1.0e+02,2.0e+02) +DEAL::4 4 (1.2e+02,2.3e+02) +DEAL::0 0 (8.3,50.) +DEAL::0 1 (12.,62.) +DEAL::0 2 (17.,75.) +DEAL::0 3 (21.,87.) +DEAL::0 4 (25.,99.) +DEAL::1 0 (17.,66.) +DEAL::1 1 (25.,83.) +DEAL::1 2 (33.,99.) +DEAL::1 3 (41.,1.2e+02) +DEAL::1 4 (50.,1.3e+02) +DEAL::2 0 (25.,83.) +DEAL::2 1 (37.,1.0e+02) +DEAL::2 2 (50.,1.2e+02) +DEAL::2 3 (62.,1.4e+02) +DEAL::2 4 (75.,1.7e+02) +DEAL::3 0 (33.,99.) +DEAL::3 1 (50.,1.2e+02) +DEAL::3 2 (66.,1.5e+02) +DEAL::3 3 (83.,1.7e+02) +DEAL::3 4 (99.,2.0e+02) +DEAL::4 0 (41.,1.2e+02) +DEAL::4 1 (62.,1.4e+02) +DEAL::4 2 (83.,1.7e+02) +DEAL::4 3 (1.0e+02,2.0e+02) +DEAL::4 4 (1.2e+02,2.3e+02) diff --git a/tests/bits/complex_complex_full_matrix_36.cc b/tests/bits/complex_complex_full_matrix_36.cc new file mode 100644 index 0000000000..867cb71ed8 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_36.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_36/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m.diagadd (3.141); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_36/cmp/generic b/tests/bits/complex_complex_full_matrix_36/cmp/generic new file mode 100644 index 0000000000..e16df449ca --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_36/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (5.1,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (9.1,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (15.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (23.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (33.,56.) +DEAL::0 0 (5.1,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (9.1,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (15.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (23.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (33.,56.) diff --git a/tests/bits/complex_complex_full_matrix_37.cc b/tests/bits/complex_complex_full_matrix_37.cc new file mode 100644 index 0000000000..0bd600e301 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_37.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_37/output"; + + +template +void +check () +{ + FullMatrix > m,n; + make_complex_matrix (m); + make_complex_matrix (n); + m.equ (3.14159, n); + + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_37/cmp/generic b/tests/bits/complex_complex_full_matrix_37/cmp/generic new file mode 100644 index 0000000000..97bce9567b --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_37/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (6.3,38.) +DEAL::0 1 (9.4,47.) +DEAL::0 2 (13.,57.) +DEAL::0 3 (16.,66.) +DEAL::0 4 (19.,75.) +DEAL::1 0 (13.,50.) +DEAL::1 1 (19.,63.) +DEAL::1 2 (25.,75.) +DEAL::1 3 (31.,88.) +DEAL::1 4 (38.,1.0e+02) +DEAL::2 0 (19.,63.) +DEAL::2 1 (28.,79.) +DEAL::2 2 (38.,94.) +DEAL::2 3 (47.,1.1e+02) +DEAL::2 4 (57.,1.3e+02) +DEAL::3 0 (25.,75.) +DEAL::3 1 (38.,94.) +DEAL::3 2 (50.,1.1e+02) +DEAL::3 3 (63.,1.3e+02) +DEAL::3 4 (75.,1.5e+02) +DEAL::4 0 (31.,88.) +DEAL::4 1 (47.,1.1e+02) +DEAL::4 2 (63.,1.3e+02) +DEAL::4 3 (79.,1.5e+02) +DEAL::4 4 (94.,1.8e+02) +DEAL::0 0 (6.3,38.) +DEAL::0 1 (9.4,47.) +DEAL::0 2 (13.,57.) +DEAL::0 3 (16.,66.) +DEAL::0 4 (19.,75.) +DEAL::1 0 (13.,50.) +DEAL::1 1 (19.,63.) +DEAL::1 2 (25.,75.) +DEAL::1 3 (31.,88.) +DEAL::1 4 (38.,1.0e+02) +DEAL::2 0 (19.,63.) +DEAL::2 1 (28.,79.) +DEAL::2 2 (38.,94.) +DEAL::2 3 (47.,1.1e+02) +DEAL::2 4 (57.,1.3e+02) +DEAL::3 0 (25.,75.) +DEAL::3 1 (38.,94.) +DEAL::3 2 (50.,1.1e+02) +DEAL::3 3 (63.,1.3e+02) +DEAL::3 4 (75.,1.5e+02) +DEAL::4 0 (31.,88.) +DEAL::4 1 (47.,1.1e+02) +DEAL::4 2 (63.,1.3e+02) +DEAL::4 3 (79.,1.5e+02) +DEAL::4 4 (94.,1.8e+02) diff --git a/tests/bits/complex_complex_full_matrix_38.cc b/tests/bits/complex_complex_full_matrix_38.cc new file mode 100644 index 0000000000..6bd0255b98 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_38.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_38/output"; + + +template +void +check () +{ + FullMatrix > m,n,o; + make_complex_matrix (m); + make_complex_matrix (n); + make_complex_matrix (o); + m.equ (3.1415, n, 2.718, o); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_38/cmp/generic b/tests/bits/complex_complex_full_matrix_38/cmp/generic new file mode 100644 index 0000000000..7516425c7d --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_38/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (12.,70.) +DEAL::0 1 (18.,88.) +DEAL::0 2 (23.,1.1e+02) +DEAL::0 3 (29.,1.2e+02) +DEAL::0 4 (35.,1.4e+02) +DEAL::1 0 (23.,94.) +DEAL::1 1 (35.,1.2e+02) +DEAL::1 2 (47.,1.4e+02) +DEAL::1 3 (59.,1.6e+02) +DEAL::1 4 (70.,1.9e+02) +DEAL::2 0 (35.,1.2e+02) +DEAL::2 1 (53.,1.5e+02) +DEAL::2 2 (70.,1.8e+02) +DEAL::2 3 (88.,2.1e+02) +DEAL::2 4 (1.1e+02,2.3e+02) +DEAL::3 0 (47.,1.4e+02) +DEAL::3 1 (70.,1.8e+02) +DEAL::3 2 (94.,2.1e+02) +DEAL::3 3 (1.2e+02,2.5e+02) +DEAL::3 4 (1.4e+02,2.8e+02) +DEAL::4 0 (59.,1.6e+02) +DEAL::4 1 (88.,2.1e+02) +DEAL::4 2 (1.2e+02,2.5e+02) +DEAL::4 3 (1.5e+02,2.9e+02) +DEAL::4 4 (1.8e+02,3.3e+02) +DEAL::0 0 (12.,70.) +DEAL::0 1 (18.,88.) +DEAL::0 2 (23.,1.1e+02) +DEAL::0 3 (29.,1.2e+02) +DEAL::0 4 (35.,1.4e+02) +DEAL::1 0 (23.,94.) +DEAL::1 1 (35.,1.2e+02) +DEAL::1 2 (47.,1.4e+02) +DEAL::1 3 (59.,1.6e+02) +DEAL::1 4 (70.,1.9e+02) +DEAL::2 0 (35.,1.2e+02) +DEAL::2 1 (53.,1.5e+02) +DEAL::2 2 (70.,1.8e+02) +DEAL::2 3 (88.,2.1e+02) +DEAL::2 4 (1.1e+02,2.3e+02) +DEAL::3 0 (47.,1.4e+02) +DEAL::3 1 (70.,1.8e+02) +DEAL::3 2 (94.,2.1e+02) +DEAL::3 3 (1.2e+02,2.5e+02) +DEAL::3 4 (1.4e+02,2.8e+02) +DEAL::4 0 (59.,1.6e+02) +DEAL::4 1 (88.,2.1e+02) +DEAL::4 2 (1.2e+02,2.5e+02) +DEAL::4 3 (1.5e+02,2.9e+02) +DEAL::4 4 (1.8e+02,3.3e+02) diff --git a/tests/bits/complex_complex_full_matrix_39.cc b/tests/bits/complex_complex_full_matrix_39.cc new file mode 100644 index 0000000000..d853368271 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_39.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_39/output"; + + +template +void +check () +{ + FullMatrix > m,n,o,p; + make_complex_matrix (m); + make_complex_matrix (n); + make_complex_matrix (o); + make_complex_matrix (p); + m.equ (3.1415, n, 2.718, o, 1.414, p); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_39/cmp/generic b/tests/bits/complex_complex_full_matrix_39/cmp/generic new file mode 100644 index 0000000000..5c5b1fb22e --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_39/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (15.,87.) +DEAL::0 1 (22.,1.1e+02) +DEAL::0 2 (29.,1.3e+02) +DEAL::0 3 (36.,1.5e+02) +DEAL::0 4 (44.,1.7e+02) +DEAL::1 0 (29.,1.2e+02) +DEAL::1 1 (44.,1.5e+02) +DEAL::1 2 (58.,1.7e+02) +DEAL::1 3 (73.,2.0e+02) +DEAL::1 4 (87.,2.3e+02) +DEAL::2 0 (44.,1.5e+02) +DEAL::2 1 (65.,1.8e+02) +DEAL::2 2 (87.,2.2e+02) +DEAL::2 3 (1.1e+02,2.5e+02) +DEAL::2 4 (1.3e+02,2.9e+02) +DEAL::3 0 (58.,1.7e+02) +DEAL::3 1 (87.,2.2e+02) +DEAL::3 2 (1.2e+02,2.6e+02) +DEAL::3 3 (1.5e+02,3.1e+02) +DEAL::3 4 (1.7e+02,3.5e+02) +DEAL::4 0 (73.,2.0e+02) +DEAL::4 1 (1.1e+02,2.5e+02) +DEAL::4 2 (1.5e+02,3.1e+02) +DEAL::4 3 (1.8e+02,3.6e+02) +DEAL::4 4 (2.2e+02,4.1e+02) +DEAL::0 0 (15.,87.) +DEAL::0 1 (22.,1.1e+02) +DEAL::0 2 (29.,1.3e+02) +DEAL::0 3 (36.,1.5e+02) +DEAL::0 4 (44.,1.7e+02) +DEAL::1 0 (29.,1.2e+02) +DEAL::1 1 (44.,1.5e+02) +DEAL::1 2 (58.,1.7e+02) +DEAL::1 3 (73.,2.0e+02) +DEAL::1 4 (87.,2.3e+02) +DEAL::2 0 (44.,1.5e+02) +DEAL::2 1 (65.,1.8e+02) +DEAL::2 2 (87.,2.2e+02) +DEAL::2 3 (1.1e+02,2.5e+02) +DEAL::2 4 (1.3e+02,2.9e+02) +DEAL::3 0 (58.,1.7e+02) +DEAL::3 1 (87.,2.2e+02) +DEAL::3 2 (1.2e+02,2.6e+02) +DEAL::3 3 (1.5e+02,3.1e+02) +DEAL::3 4 (1.7e+02,3.5e+02) +DEAL::4 0 (73.,2.0e+02) +DEAL::4 1 (1.1e+02,2.5e+02) +DEAL::4 2 (1.5e+02,3.1e+02) +DEAL::4 3 (1.8e+02,3.6e+02) +DEAL::4 4 (2.2e+02,4.1e+02) diff --git a/tests/bits/complex_complex_full_matrix_40.cc b/tests/bits/complex_complex_full_matrix_40.cc new file mode 100644 index 0000000000..5cdabe888e --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_40.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_40/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m.symmetrize (); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_40/cmp/generic b/tests/bits/complex_complex_full_matrix_40/cmp/generic new file mode 100644 index 0000000000..e8cf1827f4 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_40/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.5,16.) +DEAL::0 2 (5.0,19.) +DEAL::0 3 (6.5,22.) +DEAL::0 4 (8.0,26.) +DEAL::1 0 (3.5,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.5,24.) +DEAL::1 3 (11.,29.) +DEAL::1 4 (14.,34.) +DEAL::2 0 (5.0,19.) +DEAL::2 1 (8.5,24.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (16.,36.) +DEAL::2 4 (19.,41.) +DEAL::3 0 (6.5,22.) +DEAL::3 1 (11.,29.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (8.0,26.) +DEAL::4 1 (14.,34.) +DEAL::4 2 (19.,41.) +DEAL::4 3 (24.,48.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.5,16.) +DEAL::0 2 (5.0,19.) +DEAL::0 3 (6.5,22.) +DEAL::0 4 (8.0,26.) +DEAL::1 0 (3.5,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.5,24.) +DEAL::1 3 (11.,29.) +DEAL::1 4 (14.,34.) +DEAL::2 0 (5.0,19.) +DEAL::2 1 (8.5,24.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (16.,36.) +DEAL::2 4 (19.,41.) +DEAL::3 0 (6.5,22.) +DEAL::3 1 (11.,29.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (8.0,26.) +DEAL::4 1 (14.,34.) +DEAL::4 2 (19.,41.) +DEAL::4 3 (24.,48.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_41.cc b/tests/bits/complex_complex_full_matrix_41.cc new file mode 100644 index 0000000000..9dbed44f5f --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_41.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_41/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + m(0,0) = m(1,1) = m(2,2) = m(3,3) = m(4,4) = 50; + m.gauss_jordan (); + print_matrix (m); +} + diff --git a/tests/bits/complex_complex_full_matrix_41/cmp/generic b/tests/bits/complex_complex_full_matrix_41/cmp/generic new file mode 100644 index 0000000000..f0700a8e3c --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_41/cmp/generic @@ -0,0 +1,51 @@ + +DEAL::0 0 (0.018,0.0043) +DEAL::0 1 (0.0021,0.00086) +DEAL::0 2 (0.0022,0.0014) +DEAL::0 3 (0.0020,0.0020) +DEAL::0 4 (0.0014,0.0024) +DEAL::1 0 (0.0023,0.00080) +DEAL::1 1 (0.016,0.0072) +DEAL::1 2 (0.0030,0.0022) +DEAL::1 3 (0.0026,0.0031) +DEAL::1 4 (0.0017,0.0037) +DEAL::2 0 (0.0026,0.0014) +DEAL::2 1 (0.0032,0.0022) +DEAL::2 2 (0.013,0.0094) +DEAL::2 3 (0.0026,0.0046) +DEAL::2 4 (0.0013,0.0051) +DEAL::3 0 (0.0025,0.0022) +DEAL::3 1 (0.0030,0.0033) +DEAL::3 2 (0.0027,0.0047) +DEAL::3 3 (0.0097,0.0099) +DEAL::3 4 (-1.3e-05,0.0061) +DEAL::4 0 (0.0020,0.0027) +DEAL::4 1 (0.0021,0.0040) +DEAL::4 2 (0.0015,0.0054) +DEAL::4 3 (-6.6e-05,0.0062) +DEAL::4 4 (0.0072,0.0098) +DEAL::0 0 (0.018,0.0043) +DEAL::0 1 (0.0021,0.00086) +DEAL::0 2 (0.0022,0.0014) +DEAL::0 3 (0.0020,0.0020) +DEAL::0 4 (0.0014,0.0024) +DEAL::1 0 (0.0023,0.00080) +DEAL::1 1 (0.016,0.0072) +DEAL::1 2 (0.0030,0.0022) +DEAL::1 3 (0.0026,0.0031) +DEAL::1 4 (0.0017,0.0037) +DEAL::2 0 (0.0026,0.0014) +DEAL::2 1 (0.0032,0.0022) +DEAL::2 2 (0.013,0.0094) +DEAL::2 3 (0.0026,0.0046) +DEAL::2 4 (0.0013,0.0051) +DEAL::3 0 (0.0025,0.0022) +DEAL::3 1 (0.0030,0.0033) +DEAL::3 2 (0.0027,0.0047) +DEAL::3 3 (0.0097,0.0099) +DEAL::3 4 (-1.3e-05,0.0061) +DEAL::4 0 (0.0020,0.0027) +DEAL::4 1 (0.0021,0.0040) +DEAL::4 2 (0.0015,0.0054) +DEAL::4 3 (-6.6e-05,0.0062) +DEAL::4 4 (0.0072,0.0098) diff --git a/tests/bits/complex_complex_full_matrix_42.cc b/tests/bits/complex_complex_full_matrix_42.cc new file mode 100644 index 0000000000..fdd12a0f49 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_42.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_42/output"; + + +template +void +check () +{ + for (unsigned int n=1; n<=5; ++n) + { + const std::complex 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 > m (n,n,array), p(n,n); + p.invert (m); + print_matrix (p); + } +} + diff --git a/tests/bits/complex_complex_full_matrix_42/cmp/generic b/tests/bits/complex_complex_full_matrix_42/cmp/generic new file mode 100644 index 0000000000..47435b6391 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_42/cmp/generic @@ -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_complex_full_matrix_43.cc b/tests/bits/complex_complex_full_matrix_43.cc new file mode 100644 index 0000000000..a61ba3f44b --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_43.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_43/output"; + + +template +void +check () +{ + FullMatrix > m, n, o; + make_complex_matrix (m); + make_complex_matrix (n); + make_complex_matrix (o); + + m.mmult (n, o, true); + print_matrix (m); + print_matrix (n); + print_matrix (o); + + m.mmult (n, o, false); + print_matrix (m); + print_matrix (n); + print_matrix (o); +} + diff --git a/tests/bits/complex_complex_full_matrix_43/cmp/generic b/tests/bits/complex_complex_full_matrix_43/cmp/generic new file mode 100644 index 0000000000..88656fba42 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_43/cmp/generic @@ -0,0 +1,301 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (-1.8e+03,1.1e+03) +DEAL::0 1 (-2.2e+03,1.5e+03) +DEAL::0 2 (-2.6e+03,1.9e+03) +DEAL::0 3 (-3.0e+03,2.3e+03) +DEAL::0 4 (-3.4e+03,2.7e+03) +DEAL::1 0 (-2.3e+03,1.7e+03) +DEAL::1 1 (-2.8e+03,2.3e+03) +DEAL::1 2 (-3.3e+03,2.9e+03) +DEAL::1 3 (-3.8e+03,3.6e+03) +DEAL::1 4 (-4.3e+03,4.2e+03) +DEAL::2 0 (-2.8e+03,2.3e+03) +DEAL::2 1 (-3.4e+03,3.2e+03) +DEAL::2 2 (-3.9e+03,4.0e+03) +DEAL::2 3 (-4.5e+03,4.8e+03) +DEAL::2 4 (-5.1e+03,5.7e+03) +DEAL::3 0 (-3.3e+03,3.0e+03) +DEAL::3 1 (-3.9e+03,4.0e+03) +DEAL::3 2 (-4.6e+03,5.1e+03) +DEAL::3 3 (-5.3e+03,6.1e+03) +DEAL::3 4 (-6.0e+03,7.2e+03) +DEAL::4 0 (-3.8e+03,3.6e+03) +DEAL::4 1 (-4.5e+03,4.9e+03) +DEAL::4 2 (-5.3e+03,6.1e+03) +DEAL::4 3 (-6.1e+03,7.4e+03) +DEAL::4 4 (-6.8e+03,8.7e+03) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (-1.8e+03,1.0e+03) +DEAL::0 1 (-2.2e+03,1.4e+03) +DEAL::0 2 (-2.6e+03,1.9e+03) +DEAL::0 3 (-3.0e+03,2.3e+03) +DEAL::0 4 (-3.4e+03,2.7e+03) +DEAL::1 0 (-2.3e+03,1.7e+03) +DEAL::1 1 (-2.8e+03,2.3e+03) +DEAL::1 2 (-3.3e+03,2.9e+03) +DEAL::1 3 (-3.8e+03,3.5e+03) +DEAL::1 4 (-4.3e+03,4.2e+03) +DEAL::2 0 (-2.8e+03,2.3e+03) +DEAL::2 1 (-3.4e+03,3.2e+03) +DEAL::2 2 (-4.0e+03,4.0e+03) +DEAL::2 3 (-4.6e+03,4.8e+03) +DEAL::2 4 (-5.1e+03,5.6e+03) +DEAL::3 0 (-3.3e+03,3.0e+03) +DEAL::3 1 (-4.0e+03,4.0e+03) +DEAL::3 2 (-4.6e+03,5.0e+03) +DEAL::3 3 (-5.3e+03,6.1e+03) +DEAL::3 4 (-6.0e+03,7.1e+03) +DEAL::4 0 (-3.8e+03,3.6e+03) +DEAL::4 1 (-4.6e+03,4.8e+03) +DEAL::4 2 (-5.3e+03,6.1e+03) +DEAL::4 3 (-6.1e+03,7.4e+03) +DEAL::4 4 (-6.9e+03,8.6e+03) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (-1.8e+03,1.1e+03) +DEAL::0 1 (-2.2e+03,1.5e+03) +DEAL::0 2 (-2.6e+03,1.9e+03) +DEAL::0 3 (-3.0e+03,2.3e+03) +DEAL::0 4 (-3.4e+03,2.7e+03) +DEAL::1 0 (-2.3e+03,1.7e+03) +DEAL::1 1 (-2.8e+03,2.3e+03) +DEAL::1 2 (-3.3e+03,2.9e+03) +DEAL::1 3 (-3.8e+03,3.6e+03) +DEAL::1 4 (-4.3e+03,4.2e+03) +DEAL::2 0 (-2.8e+03,2.3e+03) +DEAL::2 1 (-3.4e+03,3.2e+03) +DEAL::2 2 (-3.9e+03,4.0e+03) +DEAL::2 3 (-4.5e+03,4.8e+03) +DEAL::2 4 (-5.1e+03,5.7e+03) +DEAL::3 0 (-3.3e+03,3.0e+03) +DEAL::3 1 (-3.9e+03,4.0e+03) +DEAL::3 2 (-4.6e+03,5.1e+03) +DEAL::3 3 (-5.3e+03,6.1e+03) +DEAL::3 4 (-6.0e+03,7.2e+03) +DEAL::4 0 (-3.8e+03,3.6e+03) +DEAL::4 1 (-4.5e+03,4.9e+03) +DEAL::4 2 (-5.3e+03,6.1e+03) +DEAL::4 3 (-6.1e+03,7.4e+03) +DEAL::4 4 (-6.8e+03,8.7e+03) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (-1.8e+03,1.0e+03) +DEAL::0 1 (-2.2e+03,1.4e+03) +DEAL::0 2 (-2.6e+03,1.9e+03) +DEAL::0 3 (-3.0e+03,2.3e+03) +DEAL::0 4 (-3.4e+03,2.7e+03) +DEAL::1 0 (-2.3e+03,1.7e+03) +DEAL::1 1 (-2.8e+03,2.3e+03) +DEAL::1 2 (-3.3e+03,2.9e+03) +DEAL::1 3 (-3.8e+03,3.5e+03) +DEAL::1 4 (-4.3e+03,4.2e+03) +DEAL::2 0 (-2.8e+03,2.3e+03) +DEAL::2 1 (-3.4e+03,3.2e+03) +DEAL::2 2 (-4.0e+03,4.0e+03) +DEAL::2 3 (-4.6e+03,4.8e+03) +DEAL::2 4 (-5.1e+03,5.6e+03) +DEAL::3 0 (-3.3e+03,3.0e+03) +DEAL::3 1 (-4.0e+03,4.0e+03) +DEAL::3 2 (-4.6e+03,5.0e+03) +DEAL::3 3 (-5.3e+03,6.1e+03) +DEAL::3 4 (-6.0e+03,7.1e+03) +DEAL::4 0 (-3.8e+03,3.6e+03) +DEAL::4 1 (-4.6e+03,4.8e+03) +DEAL::4 2 (-5.3e+03,6.1e+03) +DEAL::4 3 (-6.1e+03,7.4e+03) +DEAL::4 4 (-6.9e+03,8.6e+03) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_44.cc b/tests/bits/complex_complex_full_matrix_44.cc new file mode 100644 index 0000000000..f0c61c51d2 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_44.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_44/output"; + + +template +void +check () +{ + FullMatrix > m, n, o; + make_complex_matrix (m); + make_complex_matrix (n); + make_complex_matrix (o); + + m.Tmmult (n, o, true); + print_matrix (m); + print_matrix (n); + print_matrix (o); + + m.Tmmult (n, o, false); + print_matrix (m); + print_matrix (n); + print_matrix (o); +} + diff --git a/tests/bits/complex_complex_full_matrix_44/cmp/generic b/tests/bits/complex_complex_full_matrix_44/cmp/generic new file mode 100644 index 0000000000..9079ce4533 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_44/cmp/generic @@ -0,0 +1,301 @@ + +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (-1.9e+03,1.4e+03) +DEAL::0 1 (-2.4e+03,1.9e+03) +DEAL::0 2 (-2.8e+03,2.4e+03) +DEAL::0 3 (-3.2e+03,2.9e+03) +DEAL::0 4 (-3.7e+03,3.4e+03) +DEAL::1 0 (-2.4e+03,1.9e+03) +DEAL::1 1 (-2.9e+03,2.6e+03) +DEAL::1 2 (-3.4e+03,3.3e+03) +DEAL::1 3 (-3.9e+03,3.9e+03) +DEAL::1 4 (-4.4e+03,4.6e+03) +DEAL::2 0 (-2.8e+03,2.4e+03) +DEAL::2 1 (-3.4e+03,3.3e+03) +DEAL::2 2 (-4.0e+03,4.1e+03) +DEAL::2 3 (-4.6e+03,5.0e+03) +DEAL::2 4 (-5.1e+03,5.8e+03) +DEAL::3 0 (-3.2e+03,2.9e+03) +DEAL::3 1 (-3.9e+03,3.9e+03) +DEAL::3 2 (-4.6e+03,5.0e+03) +DEAL::3 3 (-5.2e+03,6.0e+03) +DEAL::3 4 (-5.9e+03,7.0e+03) +DEAL::4 0 (-3.6e+03,3.4e+03) +DEAL::4 1 (-4.4e+03,4.6e+03) +DEAL::4 2 (-5.1e+03,5.8e+03) +DEAL::4 3 (-5.9e+03,7.0e+03) +DEAL::4 4 (-6.6e+03,8.2e+03) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (-1.9e+03,1.4e+03) +DEAL::0 1 (-2.4e+03,1.9e+03) +DEAL::0 2 (-2.8e+03,2.4e+03) +DEAL::0 3 (-3.2e+03,2.9e+03) +DEAL::0 4 (-3.7e+03,3.4e+03) +DEAL::1 0 (-2.4e+03,1.9e+03) +DEAL::1 1 (-2.9e+03,2.6e+03) +DEAL::1 2 (-3.4e+03,3.2e+03) +DEAL::1 3 (-3.9e+03,3.9e+03) +DEAL::1 4 (-4.4e+03,4.6e+03) +DEAL::2 0 (-2.8e+03,2.4e+03) +DEAL::2 1 (-3.4e+03,3.2e+03) +DEAL::2 2 (-4.0e+03,4.1e+03) +DEAL::2 3 (-4.6e+03,4.9e+03) +DEAL::2 4 (-5.2e+03,5.8e+03) +DEAL::3 0 (-3.2e+03,2.9e+03) +DEAL::3 1 (-3.9e+03,3.9e+03) +DEAL::3 2 (-4.6e+03,4.9e+03) +DEAL::3 3 (-5.2e+03,6.0e+03) +DEAL::3 4 (-5.9e+03,7.0e+03) +DEAL::4 0 (-3.7e+03,3.4e+03) +DEAL::4 1 (-4.4e+03,4.6e+03) +DEAL::4 2 (-5.2e+03,5.8e+03) +DEAL::4 3 (-5.9e+03,7.0e+03) +DEAL::4 4 (-6.7e+03,8.2e+03) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (-1.9e+03,1.4e+03) +DEAL::0 1 (-2.4e+03,1.9e+03) +DEAL::0 2 (-2.8e+03,2.4e+03) +DEAL::0 3 (-3.2e+03,2.9e+03) +DEAL::0 4 (-3.7e+03,3.4e+03) +DEAL::1 0 (-2.4e+03,1.9e+03) +DEAL::1 1 (-2.9e+03,2.6e+03) +DEAL::1 2 (-3.4e+03,3.3e+03) +DEAL::1 3 (-3.9e+03,3.9e+03) +DEAL::1 4 (-4.4e+03,4.6e+03) +DEAL::2 0 (-2.8e+03,2.4e+03) +DEAL::2 1 (-3.4e+03,3.3e+03) +DEAL::2 2 (-4.0e+03,4.1e+03) +DEAL::2 3 (-4.6e+03,5.0e+03) +DEAL::2 4 (-5.1e+03,5.8e+03) +DEAL::3 0 (-3.2e+03,2.9e+03) +DEAL::3 1 (-3.9e+03,3.9e+03) +DEAL::3 2 (-4.6e+03,5.0e+03) +DEAL::3 3 (-5.2e+03,6.0e+03) +DEAL::3 4 (-5.9e+03,7.0e+03) +DEAL::4 0 (-3.6e+03,3.4e+03) +DEAL::4 1 (-4.4e+03,4.6e+03) +DEAL::4 2 (-5.1e+03,5.8e+03) +DEAL::4 3 (-5.9e+03,7.0e+03) +DEAL::4 4 (-6.6e+03,8.2e+03) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) +DEAL::0 0 (-1.9e+03,1.4e+03) +DEAL::0 1 (-2.4e+03,1.9e+03) +DEAL::0 2 (-2.8e+03,2.4e+03) +DEAL::0 3 (-3.2e+03,2.9e+03) +DEAL::0 4 (-3.7e+03,3.4e+03) +DEAL::1 0 (-2.4e+03,1.9e+03) +DEAL::1 1 (-2.9e+03,2.6e+03) +DEAL::1 2 (-3.4e+03,3.2e+03) +DEAL::1 3 (-3.9e+03,3.9e+03) +DEAL::1 4 (-4.4e+03,4.6e+03) +DEAL::2 0 (-2.8e+03,2.4e+03) +DEAL::2 1 (-3.4e+03,3.2e+03) +DEAL::2 2 (-4.0e+03,4.1e+03) +DEAL::2 3 (-4.6e+03,4.9e+03) +DEAL::2 4 (-5.2e+03,5.8e+03) +DEAL::3 0 (-3.2e+03,2.9e+03) +DEAL::3 1 (-3.9e+03,3.9e+03) +DEAL::3 2 (-4.6e+03,4.9e+03) +DEAL::3 3 (-5.2e+03,6.0e+03) +DEAL::3 4 (-5.9e+03,7.0e+03) +DEAL::4 0 (-3.7e+03,3.4e+03) +DEAL::4 1 (-4.4e+03,4.6e+03) +DEAL::4 2 (-5.2e+03,5.8e+03) +DEAL::4 3 (-5.9e+03,7.0e+03) +DEAL::4 4 (-6.7e+03,8.2e+03) +DEAL::0 0 (2.0,12.) +DEAL::0 1 (3.0,15.) +DEAL::0 2 (4.0,18.) +DEAL::0 3 (5.0,21.) +DEAL::0 4 (6.0,24.) +DEAL::1 0 (4.0,16.) +DEAL::1 1 (6.0,20.) +DEAL::1 2 (8.0,24.) +DEAL::1 3 (10.,28.) +DEAL::1 4 (12.,32.) +DEAL::2 0 (6.0,20.) +DEAL::2 1 (9.0,25.) +DEAL::2 2 (12.,30.) +DEAL::2 3 (15.,35.) +DEAL::2 4 (18.,40.) +DEAL::3 0 (8.0,24.) +DEAL::3 1 (12.,30.) +DEAL::3 2 (16.,36.) +DEAL::3 3 (20.,42.) +DEAL::3 4 (24.,48.) +DEAL::4 0 (10.,28.) +DEAL::4 1 (15.,35.) +DEAL::4 2 (20.,42.) +DEAL::4 3 (25.,49.) +DEAL::4 4 (30.,56.) diff --git a/tests/bits/complex_complex_full_matrix_45.cc b/tests/bits/complex_complex_full_matrix_45.cc new file mode 100644 index 0000000000..0462394c0a --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_45.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_45/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + Vector > v, w; + make_complex_vector (v); + make_complex_vector (w); + + m.vmult (v, w, true); + print_vector (v); + print_vector (w); + + m.vmult (v, w, false); + print_vector (v); + print_vector (w); +} + diff --git a/tests/bits/complex_complex_full_matrix_45/cmp/generic b/tests/bits/complex_complex_full_matrix_45/cmp/generic new file mode 100644 index 0000000000..b80b7e6cab --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_45/cmp/generic @@ -0,0 +1,41 @@ + +DEAL::0 (-4.1e+02,4.1e+02) +DEAL::1 (-5.0e+02,6.2e+02) +DEAL::2 (-5.9e+02,8.4e+02) +DEAL::3 (-6.8e+02,1.0e+03) +DEAL::4 (-7.6e+02,1.3e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (-4.1e+02,4.1e+02) +DEAL::1 (-5.0e+02,6.2e+02) +DEAL::2 (-5.9e+02,8.3e+02) +DEAL::3 (-6.8e+02,1.0e+03) +DEAL::4 (-7.7e+02,1.2e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (-4.1e+02,4.1e+02) +DEAL::1 (-5.0e+02,6.2e+02) +DEAL::2 (-5.9e+02,8.4e+02) +DEAL::3 (-6.8e+02,1.0e+03) +DEAL::4 (-7.6e+02,1.3e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (-4.1e+02,4.1e+02) +DEAL::1 (-5.0e+02,6.2e+02) +DEAL::2 (-5.9e+02,8.3e+02) +DEAL::3 (-6.8e+02,1.0e+03) +DEAL::4 (-7.7e+02,1.2e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) diff --git a/tests/bits/complex_complex_full_matrix_46.cc b/tests/bits/complex_complex_full_matrix_46.cc new file mode 100644 index 0000000000..4bb1646365 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_46.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_46/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + Vector > v, w; + make_complex_vector (v); + make_complex_vector (w); + + m.vmult_add (v, w); + print_vector (v); + print_vector (w); +} + diff --git a/tests/bits/complex_complex_full_matrix_46/cmp/generic b/tests/bits/complex_complex_full_matrix_46/cmp/generic new file mode 100644 index 0000000000..7dc66d43f4 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_46/cmp/generic @@ -0,0 +1,21 @@ + +DEAL::0 (-4.1e+02,4.1e+02) +DEAL::1 (-5.0e+02,6.2e+02) +DEAL::2 (-5.9e+02,8.4e+02) +DEAL::3 (-6.8e+02,1.0e+03) +DEAL::4 (-7.6e+02,1.3e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (-4.1e+02,4.1e+02) +DEAL::1 (-5.0e+02,6.2e+02) +DEAL::2 (-5.9e+02,8.4e+02) +DEAL::3 (-6.8e+02,1.0e+03) +DEAL::4 (-7.6e+02,1.3e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) diff --git a/tests/bits/complex_complex_full_matrix_47.cc b/tests/bits/complex_complex_full_matrix_47.cc new file mode 100644 index 0000000000..bc6d03cc89 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_47.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_47/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + Vector > v, w; + make_complex_vector (v); + make_complex_vector (w); + + m.Tvmult (v, w, true); + print_vector (v); + print_vector (w); + + m.Tvmult (v, w, false); + print_vector (v); + print_vector (w); +} + diff --git a/tests/bits/complex_complex_full_matrix_47/cmp/generic b/tests/bits/complex_complex_full_matrix_47/cmp/generic new file mode 100644 index 0000000000..1cd3e94d66 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_47/cmp/generic @@ -0,0 +1,41 @@ + +DEAL::0 (-4.3e+02,5.1e+02) +DEAL::1 (-5.1e+02,6.8e+02) +DEAL::2 (-5.9e+02,8.6e+02) +DEAL::3 (-6.7e+02,1.0e+03) +DEAL::4 (-7.4e+02,1.2e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (-4.3e+02,5.1e+02) +DEAL::1 (-5.1e+02,6.8e+02) +DEAL::2 (-5.9e+02,8.5e+02) +DEAL::3 (-6.7e+02,1.0e+03) +DEAL::4 (-7.5e+02,1.2e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (-4.3e+02,5.1e+02) +DEAL::1 (-5.1e+02,6.8e+02) +DEAL::2 (-5.9e+02,8.6e+02) +DEAL::3 (-6.7e+02,1.0e+03) +DEAL::4 (-7.4e+02,1.2e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (-4.3e+02,5.1e+02) +DEAL::1 (-5.1e+02,6.8e+02) +DEAL::2 (-5.9e+02,8.5e+02) +DEAL::3 (-6.7e+02,1.0e+03) +DEAL::4 (-7.5e+02,1.2e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) diff --git a/tests/bits/complex_complex_full_matrix_48.cc b/tests/bits/complex_complex_full_matrix_48.cc new file mode 100644 index 0000000000..0ce11e5fda --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_48.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_48/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + Vector > v, w; + make_complex_vector (v); + make_complex_vector (w); + + m.Tvmult_add (v, w); + print_vector (v); + print_vector (w); +} + diff --git a/tests/bits/complex_complex_full_matrix_48/cmp/generic b/tests/bits/complex_complex_full_matrix_48/cmp/generic new file mode 100644 index 0000000000..0b75c0ef91 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_48/cmp/generic @@ -0,0 +1,21 @@ + +DEAL::0 (-4.3e+02,5.1e+02) +DEAL::1 (-5.1e+02,6.8e+02) +DEAL::2 (-5.9e+02,8.6e+02) +DEAL::3 (-6.7e+02,1.0e+03) +DEAL::4 (-7.4e+02,1.2e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (-4.3e+02,5.1e+02) +DEAL::1 (-5.1e+02,6.8e+02) +DEAL::2 (-5.9e+02,8.6e+02) +DEAL::3 (-6.7e+02,1.0e+03) +DEAL::4 (-7.4e+02,1.2e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) diff --git a/tests/bits/complex_complex_full_matrix_49.cc b/tests/bits/complex_complex_full_matrix_49.cc new file mode 100644 index 0000000000..73f3b889c4 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_49.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_49/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + Vector > v, w; + make_complex_vector (v); + make_complex_vector (w); + + m.precondition_Jacobi (v, w, 3.141); + print_vector (v); + print_vector (w); +} + diff --git a/tests/bits/complex_complex_full_matrix_49/cmp/generic b/tests/bits/complex_complex_full_matrix_49/cmp/generic new file mode 100644 index 0000000000..22a520b353 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_49/cmp/generic @@ -0,0 +1,21 @@ + +DEAL::0 (0.81,-0.13) +DEAL::1 (0.66,-0.12) +DEAL::2 (0.56,-0.090) +DEAL::3 (0.48,-0.070) +DEAL::4 (0.42,-0.054) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (0.81,-0.13) +DEAL::1 (0.66,-0.12) +DEAL::2 (0.56,-0.090) +DEAL::3 (0.48,-0.070) +DEAL::4 (0.42,-0.054) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) diff --git a/tests/bits/complex_complex_full_matrix_50.cc b/tests/bits/complex_complex_full_matrix_50.cc new file mode 100644 index 0000000000..c04f9cc0fc --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_50.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_50/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + Vector > v, w, x; + make_complex_vector (v); + make_complex_vector (w); + make_complex_vector (x); + + m.residual (v, w, x); + print_vector (v); + print_vector (w); + print_vector (x); +} + diff --git a/tests/bits/complex_complex_full_matrix_50/cmp/generic b/tests/bits/complex_complex_full_matrix_50/cmp/generic new file mode 100644 index 0000000000..c219d84978 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_50/cmp/generic @@ -0,0 +1,31 @@ + +DEAL::0 (4.1e+02,-4.1e+02) +DEAL::1 (5.0e+02,-6.2e+02) +DEAL::2 (5.9e+02,-8.2e+02) +DEAL::3 (6.8e+02,-1.0e+03) +DEAL::4 (7.8e+02,-1.2e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (4.1e+02,-4.1e+02) +DEAL::1 (5.0e+02,-6.2e+02) +DEAL::2 (5.9e+02,-8.2e+02) +DEAL::3 (6.8e+02,-1.0e+03) +DEAL::4 (7.8e+02,-1.2e+03) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) diff --git a/tests/bits/complex_complex_full_matrix_51.cc b/tests/bits/complex_complex_full_matrix_51.cc new file mode 100644 index 0000000000..0692a74df7 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_51.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_51/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + Vector > x, y; + make_complex_vector (x); + make_complex_vector (y); + + m.forward (x, y); + print_vector (x); + print_vector (y); +} + diff --git a/tests/bits/complex_complex_full_matrix_51/cmp/generic b/tests/bits/complex_complex_full_matrix_51/cmp/generic new file mode 100644 index 0000000000..1d024aaa59 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_51/cmp/generic @@ -0,0 +1,21 @@ + +DEAL::0 (0.26,-0.041) +DEAL::1 (0.0069,-0.014) +DEAL::2 (0.0045,0.0060) +DEAL::3 (0.0028,0.0029) +DEAL::4 (0.0018,0.0016) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (0.26,-0.041) +DEAL::1 (0.0069,-0.014) +DEAL::2 (0.0045,0.0060) +DEAL::3 (0.0028,0.0029) +DEAL::4 (0.0018,0.0016) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) diff --git a/tests/bits/complex_complex_full_matrix_52.cc b/tests/bits/complex_complex_full_matrix_52.cc new file mode 100644 index 0000000000..9c99b51f91 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_52.cc @@ -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; this time we actually store complex values +// in them + + +#include "../tests.h" +#include "full_matrix_common.h" + + +std::string output_file_name = "complex_complex_full_matrix_52/output"; + + +template +void +check () +{ + FullMatrix > m; + make_complex_matrix (m); + Vector > x, y; + make_complex_vector (x); + make_complex_vector (y); + + m.backward (x, y); + print_vector (x); + print_vector (y); +} + diff --git a/tests/bits/complex_complex_full_matrix_52/cmp/generic b/tests/bits/complex_complex_full_matrix_52/cmp/generic new file mode 100644 index 0000000000..0280c0d8f5 --- /dev/null +++ b/tests/bits/complex_complex_full_matrix_52/cmp/generic @@ -0,0 +1,21 @@ + +DEAL::0 (0.0049,0.0082) +DEAL::1 (0.0029,0.0033) +DEAL::2 (0.0017,0.0015) +DEAL::3 (0.0011,0.00080) +DEAL::4 (0.13,-0.017) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0) +DEAL::0 (0.0049,0.0082) +DEAL::1 (0.0029,0.0033) +DEAL::2 (0.0017,0.0015) +DEAL::3 (0.0011,0.00080) +DEAL::4 (0.13,-0.017) +DEAL::0 (1.0,3.0) +DEAL::1 (2.0,4.0) +DEAL::2 (3.0,5.0) +DEAL::3 (4.0,6.0) +DEAL::4 (5.0,7.0)