Fix backend::transpose() for block-valued matrices

jamisonlu · Feb 2, 2022 · 67ed91e · 67ed91e
1 parent ab0817a
commit 67ed91e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/amgcl/backend/builtin.hpp b/amgcl/backend/builtin.hpp
@@ -366,7 +366,7 @@ std::shared_ptr< crs<V,C,P> > transpose(const crs<V, C, P> &A)
             P head = T->ptr[A.col[j]]++;
 
             T->col[head] = static_cast<C>(i);
-            T->val[head] = A.val[j];
+            T->val[head] = math::adjoint(A.val[j]);
         }
     }