From: wolf Date: Mon, 23 Feb 2004 17:14:49 +0000 (+0000) Subject: Implement PETSc full matrices. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7e42d40a2e5b18c537ec65ac80e9b87e007c3e7b;p=dealii-svn.git Implement PETSc full matrices. git-svn-id: https://svn.dealii.org/trunk@8491 0785d39b-7218-0410-832d-ea1e28bc413d --- 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