LU: fix array length bug I've introduced into MILU(0) a few commits ago

13 years ago · 07725b4cc1
1 changed files with 2 additions and 2 deletions
--- a/src/Numerics/LinearAlgebra/Double/Solvers/MILU0Preconditioner.cs
+++ b/src/Numerics/LinearAlgebra/Double/Solvers/MILU0Preconditioner.cs
@ -100,8 +100,8 @@ namespace MathNet.Numerics.LinearAlgebra.Double.Solvers
            int[] ja = csr.ColumnIndices;
            int[] ia = csr.RowPointers;

-            _alu = new double[ia.Length];
-            _jlu = new int[ia.Length];
+            _alu = new double[ia[n] + 1];
+            _jlu = new int[ia[n] + 1];
            _diag = new int[n];

            int code = Compute(n, a, ja, ia, _alu, _jlu, _diag, UseModified);