From: guido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
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<MATRIX>::Tvmult_add (BlockVector<VECTOR>& 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);
 	}
     }
 }