]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Bug in transposed removed
authorguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 31 Oct 2001 13:47:26 +0000 (13:47 +0000)
committerguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 31 Oct 2001 13:47:26 +0000 (13:47 +0000)
git-svn-id: https://svn.dealii.org/trunk@5178 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/lac/include/lac/block_matrix_array.h

index 347ef511b299c8dc0bdd613ed3f02954d2f633af..1ce794b928ef1ae4c01af95aca844c103464a5e2 100644 (file)
@@ -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);
        }
     }
 }

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.