From: Wolfgang Bangerth Date: Mon, 9 Feb 2004 23:02:51 +0000 (+0000) Subject: New tests. X-Git-Tag: v8.0.0~15810 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4056307af673ff133142323f267b1201ddbb65b5;p=dealii.git New tests. git-svn-id: https://svn.dealii.org/trunk@8445 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/bits/petsc_sparse_matrix_vector_01.cc b/tests/bits/petsc_sparse_matrix_vector_01.cc new file mode 100644 index 0000000000..2adbe046ba --- /dev/null +++ b/tests/bits/petsc_sparse_matrix_vector_01.cc @@ -0,0 +1,99 @@ +//---------------------------- petsc_sparse_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_sparse_matrix_vector_01.cc --------------------------- + + +// check SparseMatrix::vmult + +#include "../tests.h" +#include +#include +#include +#include +#include + + +void test (petsc_wrappers::Vector &v, + petsc_wrappers::Vector &w) +{ + petsc_wrappers::SparseMatrix m(v.size(),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::SparseMatrix m(v.size(),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::SparseMatrix m(v.size(),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::SparseMatrix m(v.size(),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::SparseMatrix m(v.size(),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::SparseMatrix m(v.size(),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::SparseMatrix m(v.size(),v.size(),v.size()); + for (unsigned int i=0; i