From cee5c7eb97f991c9ec94cc6b9f00baa345da6665 Mon Sep 17 00:00:00 2001 From: bangerth Date: Thu, 7 Jan 2010 14:40:13 +0000 Subject: [PATCH] Use correct (?) sign. git-svn-id: https://svn.dealii.org/trunk@20330 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-42/step-42.cc | 16 ++++------------ 1 file changed, 4 insertions(+), 12 deletions(-) diff --git a/deal.II/examples/step-42/step-42.cc b/deal.II/examples/step-42/step-42.cc index 55891b9d7a..d5bcacb38b 100644 --- a/deal.II/examples/step-42/step-42.cc +++ b/deal.II/examples/step-42/step-42.cc @@ -401,7 +401,9 @@ void SchurComplement::vmult (Vector &dst, A_inverse->name = "in schur"; A_inverse->vmult (tmp2, tmp1); system_matrix->block(1,0).vmult (dst, tmp2); + dst *= -1; system_matrix->block(1,1).vmult_add (dst, src); + dst *= -1; } @@ -798,8 +800,7 @@ vmult (BlockVector &dst, rhs = src; direct_solver.vmult(solution, rhs); dst = solution; - - /* +/* const InverseMatrix,InnerPreconditioner> A_inverse (system_matrix->block(0,0), *A_preconditioner); Vector tmp (dst.block(0).size()); @@ -829,15 +830,6 @@ vmult (BlockVector &dst, << schur_complement.m() << " unknowns" << std::endl; #endif -// FullMatrix full_schur(schur_complement.n(),schur_complement.m()); -// copy(schur_complement, full_schur); -// std::ostringstream filename; -// filename << "schur_matrix"; -// std::ofstream output (filename.str().c_str()); -// full_schur.print_formatted(output, 1, true,0,"0",1,0); -// std::cout << full_schur.relative_symmetry_norm2 () << std::endl; -// std::cout << full_schur.frobenius_norm () << std::endl; -// abort(); try { cg.solve (schur_complement, dst.block(1), schur_rhs, @@ -873,7 +865,7 @@ vmult (BlockVector &dst, #ifdef STEP_42_TEST std::cout << "Exiting smoother with " << dst.size() << " unknowns" << std::endl; #endif - */ +*/ } -- 2.39.5