Add OfJaggedArray for dense storage.

Author: wo80

CSparse.Tests/Double/DenseMatrixTest.cs

@@ -304,6 +304,54 @@ public void TestOfMatrix(int rows, int columns)
             Assert.IsTrue(denseA.Equals(denseB));
         }
 
+        [TestCase(2, 2)]
+        [TestCase(2, 3)]
+        public void TestOfArray(int rows, int columns)
+        {
+            var denseData = MatrixHelper.LoadDense(rows, columns);
+
+            var denseA = denseData.A;
+
+            var jArray = new double[rows, columns];
+
+            for (int i = 0; i < rows; i++)
+            {
+                for (int j = 0; j < columns; j++)
+                {
+                    jArray[i, j] = denseA.At(i, j);
+                }
+            }
+
+            var denseB = DenseMatrix.OfArray(jArray);
+
+            Assert.IsTrue(denseA.Equals(denseB));
+        }
+
+        [TestCase(2, 2)]
+        [TestCase(2, 3)]
+        public void TestOfJaggedArray(int rows, int columns)
+        {
+            var denseData = MatrixHelper.LoadDense(rows, columns);
+
+            var denseA = denseData.A;
+
+            var jArray = new double[rows][];
+
+            for (int i = 0; i < rows; i++)
+            {
+                var r = jArray[i] = new double[columns];
+
+                for (int j = 0; j < columns; j++)
+                {
+                    r[j] = denseA.At(i, j);
+                }
+            }
+
+            var denseB = DenseMatrix.OfJaggedArray(jArray);
+
+            Assert.IsTrue(denseA.Equals(denseB));
+        }
+
         [TestCase(2, 2)]
         [TestCase(2, 3)]
         public void TestOfIndexed(int rows, int columns)

CSparse/Storage/DenseColumnMajorStorage.cs

@@ -92,6 +92,27 @@ public static DenseColumnMajorStorage<T> OfArray(T[,] array)
             return A;
         }
 
+        /// <summary>
+        /// Create a new dense matrix as a copy of the given two-dimensional array.
+        /// </summary>
+        public static DenseColumnMajorStorage<T> OfJaggedArray(T[][] array)
+        {
+            int rows = array.Length;
+            int columns = array[0].Length;
+
+            var A = Create(rows, columns);
+
+            for (int i = 0; i < rows; i++)
+            {
+                for (int j = 0; j < columns; j++)
+                {
+                    A.At(i, j, array[i][j]);
+                }
+            }
+
+            return A;
+        }
+
         /// <summary>
         /// Create a new dense matrix as a copy of the given indexed enumerable.
         /// </summary>