From 7e42d40a2e5b18c537ec65ac80e9b87e007c3e7b Mon Sep 17 00:00:00 2001 From: wolf Date: Mon, 23 Feb 2004 17:14:49 +0000 Subject: [PATCH] Implement PETSc full matrices. git-svn-id: https://svn.dealii.org/trunk@8491 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/bits/petsc_full_matrix_vector_01.cc | 99 ++++++++++++++++ tests/bits/petsc_full_matrix_vector_02.cc | 99 ++++++++++++++++ tests/bits/petsc_full_matrix_vector_03.cc | 102 +++++++++++++++++ tests/bits/petsc_full_matrix_vector_04.cc | 102 +++++++++++++++++ tests/bits/petsc_full_matrix_vector_05.cc | 105 +++++++++++++++++ tests/bits/petsc_full_matrix_vector_06.cc | 96 ++++++++++++++++ tests/bits/petsc_full_matrix_vector_07.cc | 108 ++++++++++++++++++ .../bits/petsc_full_matrix_vector_01.output | 2 + .../bits/petsc_full_matrix_vector_02.output | 2 + .../bits/petsc_full_matrix_vector_03.output | 2 + .../bits/petsc_full_matrix_vector_04.output | 2 + .../bits/petsc_full_matrix_vector_05.output | 2 + .../bits/petsc_full_matrix_vector_06.output | 2 + .../bits/petsc_full_matrix_vector_07.output | 2 + 14 files changed, 725 insertions(+) create mode 100644 tests/bits/petsc_full_matrix_vector_01.cc create mode 100644 tests/bits/petsc_full_matrix_vector_02.cc create mode 100644 tests/bits/petsc_full_matrix_vector_03.cc create mode 100644 tests/bits/petsc_full_matrix_vector_04.cc create mode 100644 tests/bits/petsc_full_matrix_vector_05.cc create mode 100644 tests/bits/petsc_full_matrix_vector_06.cc create mode 100644 tests/bits/petsc_full_matrix_vector_07.cc create mode 100644 tests/results/i686-pc-linux-gnu+gcc3.2/bits/petsc_full_matrix_vector_01.output create mode 100644 tests/results/i686-pc-linux-gnu+gcc3.2/bits/petsc_full_matrix_vector_02.output create mode 100644 tests/results/i686-pc-linux-gnu+gcc3.2/bits/petsc_full_matrix_vector_03.output create mode 100644 tests/results/i686-pc-linux-gnu+gcc3.2/bits/petsc_full_matrix_vector_04.output create mode 100644 tests/results/i686-pc-linux-gnu+gcc3.2/bits/petsc_full_matrix_vector_05.output create mode 100644 tests/results/i686-pc-linux-gnu+gcc3.2/bits/petsc_full_matrix_vector_06.output create mode 100644 tests/results/i686-pc-linux-gnu+gcc3.2/bits/petsc_full_matrix_vector_07.output diff --git a/tests/bits/petsc_full_matrix_vector_01.cc b/tests/bits/petsc_full_matrix_vector_01.cc new file mode 100644 index 0000000000..8878f9b7e4 --- /dev/null +++ b/tests/bits/petsc_full_matrix_vector_01.cc @@ -0,0 +1,99 @@ +//---------------------------- petsc_full_matrix_vector_01.cc --------------------------- +// vector_11.cc,v 1.4 2003/07/03 10:31:46 guido Exp +// Version: +// +// Copyright (C) 2004 by the deal.II authors and Anna Schneebeli +// +// 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. +// +//---------------------------- petsc_full_matrix_vector_01.cc --------------------------- + + +// check FullMatrix::vmult + +#include "../tests.h" +#include +#include +#include +#include +#include + + +void test (petsc_wrappers::Vector &v, + petsc_wrappers::Vector &w) +{ + petsc_wrappers::FullMatrix m(v.size(),v.size()); + for (unsigned int i=0; i +#include +#include +#include +#include + + +void test (petsc_wrappers::Vector &v, + petsc_wrappers::Vector &w) +{ + petsc_wrappers::FullMatrix m(v.size(),v.size()); + for (unsigned int i=0; i +#include +#include +#include +#include + + +void test (petsc_wrappers::Vector &v, + petsc_wrappers::Vector &w) +{ + petsc_wrappers::FullMatrix m(v.size(),v.size()); + for (unsigned int i=0; i +#include +#include +#include +#include + + +void test (petsc_wrappers::Vector &v, + petsc_wrappers::Vector &w) +{ + petsc_wrappers::FullMatrix m(v.size(),v.size()); + for (unsigned int i=0; i +#include +#include +#include +#include + + +void test (petsc_wrappers::Vector &v, + petsc_wrappers::Vector &w) +{ + petsc_wrappers::FullMatrix m(v.size(),v.size()); + for (unsigned int i=0; i + const double s = m.matrix_scalar_product (w,v); + + // make sure we get the expected result + for (unsigned int i=0; i +#include +#include +#include +#include + + +void test (petsc_wrappers::Vector &v) +{ + petsc_wrappers::FullMatrix m(v.size(),v.size()); + for (unsigned int i=0; i + const double s = m.matrix_norm_square (v); + + // make sure we get the expected result + for (unsigned int i=0; i +#include +#include +#include +#include + + +void test (petsc_wrappers::Vector &v, + petsc_wrappers::Vector &w, + petsc_wrappers::Vector &x) +{ + petsc_wrappers::FullMatrix m(v.size(),v.size()); + for (unsigned int i=0; i