From: guido Date: Wed, 31 Oct 2001 13:47:26 +0000 (+0000) Subject: Bug in transposed removed X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=aba11fa4a7ae7b87ad931d6b972dca0a48bd95d5;p=dealii-svn.git Bug in transposed removed git-svn-id: https://svn.dealii.org/trunk@5178 0785d39b-7218-0410-832d-ea1e28bc413d --- 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); } } }