From aba11fa4a7ae7b87ad931d6b972dca0a48bd95d5 Mon Sep 17 00:00:00 2001 From: guido Date: Wed, 31 Oct 2001 13:47:26 +0000 Subject: [PATCH] Bug in transposed removed git-svn-id: https://svn.dealii.org/trunk@5178 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/lac/include/lac/block_matrix_array.h | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/deal.II/lac/include/lac/block_matrix_array.h b/deal.II/lac/include/lac/block_matrix_array.h index 347ef511b2..1ce794b928 100644 --- a/deal.II/lac/include/lac/block_matrix_array.h +++ b/deal.II/lac/include/lac/block_matrix_array.h @@ -314,20 +314,20 @@ BlockMatrixArray::Tvmult_add (BlockVector& dst, if (m->prefix==1.) { if (m->transpose) - m->matrix->vmult_add(dst.block(m->row), - src.block(m->col)); + m->matrix->vmult_add(dst.block(m->col), + src.block(m->row)); else - m->matrix->Tvmult_add(dst.block(m->row), - src.block(m->col)); + m->matrix->Tvmult_add(dst.block(m->col), + src.block(m->row)); } else { - aux.reinit(dst.block(m->row)); + aux.reinit(dst.block(m->col)); if (m->transpose) m->matrix->vmult(aux, - src.block(m->col)); + src.block(m->row)); else m->matrix->Tvmult(aux, - src.block(m->col)); - dst.block(m->row).add (m->prefix, aux); + src.block(m->row)); + dst.block(m->col).add (m->prefix, aux); } } } -- 2.39.5