From 8ef531ff8e5e025e2ff249427ca6998782b7e7c5 Mon Sep 17 00:00:00 2001 From: ESeNonFossiIo Date: Tue, 16 Jun 2015 22:33:47 +0200 Subject: [PATCH] removing test that requires lower triangular matrix --- tests/lac/block_linear_operator_03.cc | 101 ------------------ ...ck_linear_operator_03.with_cxx11=on.output | 40 ------- 2 files changed, 141 deletions(-) delete mode 100644 tests/lac/block_linear_operator_03.cc delete mode 100644 tests/lac/block_linear_operator_03.with_cxx11=on.output diff --git a/tests/lac/block_linear_operator_03.cc b/tests/lac/block_linear_operator_03.cc deleted file mode 100644 index dad87eb770..0000000000 --- a/tests/lac/block_linear_operator_03.cc +++ /dev/null @@ -1,101 +0,0 @@ -// --------------------------------------------------------------------- -// -// Copyright (C) 2015 by the deal.II authors -// -// This file is part of the deal.II library. -// -// The deal.II library is free software; you can use it, redistribute -// it, and/or modify it under the terms of the GNU Lesser General -// Public License as published by the Free Software Foundation; either -// version 2.1 of the License, or (at your option) any later version. -// The full text of the license can be found in the file LICENSE at -// the top level of the deal.II distribution. -// -// --------------------------------------------------------------------- - -// Test lower_triangular_operator: - -#include "../tests.h" - -#include -#include -#include -#include -#include -#include - -#define PRINTME(title, name, var) \ - deallog << title << std::endl; \ - deallog << "Block vector: " name << ":" << std::endl; \ - for (unsigned int i = 0; i < var.n_blocks(); ++i) \ - deallog << "[block " << i << " ] " << var.block(i); - - -using namespace dealii; - -int main() -{ - initlog(); - deallog << std::setprecision(2); - - // BlockSparseMatrix: - { - BlockDynamicSparsityPattern dsp(3, 3); - for (unsigned int i = 0; i < 3; ++i) - for (unsigned int j = 0; j < 3; ++j) - dsp.block(i, j).reinit (1, 1); - dsp.collect_sizes (); - - BlockSparsityPattern sparsity_pattern; - sparsity_pattern.copy_from(dsp); - sparsity_pattern.compress(); - - // | 2 1 0 | | 0 0 0 | - // | 3 2 1 | -> | 3 0 0 | - // | 4 3 2 | | 4 3 0 | - BlockSparseMatrix a (sparsity_pattern); - for (unsigned int i = 0; i < 3; ++i) - for (unsigned int j = 0; j < 3; ++j) - a.block(i,j).set(0, 0, 2 + i - j); - - auto op_a = linear_operator>(a); - auto triangular_block_op = lower_triangular_operator<3, BlockVector, BlockVector, BlockSparseMatrix >(a); - - BlockVector u; - op_a.reinit_domain_vector(u, false); - BlockVector v; - op_a.reinit_range_vector(v, false); - - - // vmult: - for(unsigned int i = 0; i<3; ++i) - { - u.block(i)[0] = i+1; - v.block(i)[0] = 1; - } - - triangular_block_op.vmult(v, u); - PRINTME(" -- vmult --", "v", v); - - // vmult_add: - for(unsigned int i = 0; i<3; ++i) - v.block(i)[0] = 1; - - triangular_block_op.vmult_add(v, u); - PRINTME(" -- vmult_add --", "v", v); - - // Tvmult - for(unsigned int i = 0; i<3; ++i) - v.block(i)[0] = i+1; - - triangular_block_op.Tvmult(u, v); - PRINTME(" -- Tvmult --", "u", u); - - // Tvmult_add - for(unsigned int i = 0; i<3; ++i) - u.block(i)[0] = 1; - - triangular_block_op.Tvmult_add(u, v); - PRINTME(" -- Tvmult_add --", "u", u); - } -} diff --git a/tests/lac/block_linear_operator_03.with_cxx11=on.output b/tests/lac/block_linear_operator_03.with_cxx11=on.output deleted file mode 100644 index f4633ca2e2..0000000000 --- a/tests/lac/block_linear_operator_03.with_cxx11=on.output +++ /dev/null @@ -1,40 +0,0 @@ - -DEAL:: -- Matrix -- -DEAL::Block vector: v: -DEAL::[block 0 ] 1.00000000000 -DEAL::[block 1 ] 10.0000000000 -DEAL::[block 2 ] 10.0000000000 -DEAL::Block vector: v: -DEAL::[block 0 ] 0.00000000000 -DEAL::[block 1 ] 3.00000000000 -DEAL::[block 2 ] 10.0000000000 -DEAL::Block vector: v: -DEAL::[block 0 ] 0.00000000000 -DEAL::[block 1 ] 0.00000000000 -DEAL::[block 2 ] 5.00000000000 -DEAL:: -- Inverse -- -DEAL::Block vector: v: -DEAL::[block 0 ] 1.00000000000 -DEAL::[block 1 ] -3.33333333333 -DEAL::[block 2 ] 4.66666666667 -DEAL::Block vector: v: -DEAL::[block 0 ] 0.00000000000 -DEAL::[block 1 ] 0.333333333333 -DEAL::[block 2 ] -0.666666666667 -DEAL::Block vector: v: -DEAL::[block 0 ] 0.00000000000 -DEAL::[block 1 ] 0.00000000000 -DEAL::[block 2 ] 0.200000000000 -DEAL:: -- Identity -- -DEAL::Block vector: v: -DEAL::[block 0 ] 1.00000000000 -DEAL::[block 1 ] 0.00000000000 -DEAL::[block 2 ] 0.00000000000 -DEAL::Block vector: v: -DEAL::[block 0 ] 0.00000000000 -DEAL::[block 1 ] 1.00000000000 -DEAL::[block 2 ] 0.00000000000 -DEAL::Block vector: v: -DEAL::[block 0 ] 0.00000000000 -DEAL::[block 1 ] 0.00000000000 -DEAL::[block 2 ] 1.00000000000 -- 2.39.5