Fixed bug in SparseDirectUMFPACK #12916

diff --git a/source/lac/sparse_direct.cc b/source/lac/sparse_direct.cc
index 44a3232fbc5352c44a6f49772fe71ebefb19efa2..c34239b21a70831f69111f88e081da5b275ce97c 100644 (file)
--- a/source/lac/sparse_direct.cc
+++ b/source/lac/sparse_direct.cc
@@ -199,7 +199,7 @@ SparseDirectUMFPACK::sort_arrays(const BlockSparseMatrix<number> &matrix)
 
               std::swap(Ax[element], Ax[element + 1]);
               if (numbers::NumberTraits<number>::is_complex == true)
-                std::swap(Az[cursor], Az[cursor + 1]);
+                std::swap(Az[element], Az[element + 1]);
 
               ++element;
             }