MATH-1389: Performance improvement for Array2DRowRealMatrix.getSubMatrix()
diff --git a/src/changes/changes.xml b/src/changes/changes.xml
index 1e43272..ac6d23a 100644
--- a/src/changes/changes.xml
+++ b/src/changes/changes.xml
@@ -54,6 +54,9 @@
     </release>
 
     <release version="4.0" date="XXXX-XX-XX" description="">
+      <action dev="ebourg" type="update" issue="MATH-1389" due-to="Christoph Dibak">
+        Performance improvement for Array2DRowRealMatrix.getSubMatrix()
+      </action>
       <action dev="erans" type="add" issue="MATH-1383">
         Dependency toward the "Commons Rng" component.
       </action>
diff --git a/src/main/java/org/apache/commons/math4/linear/Array2DRowRealMatrix.java b/src/main/java/org/apache/commons/math4/linear/Array2DRowRealMatrix.java
index 3778db5..ed1ca9f 100644
--- a/src/main/java/org/apache/commons/math4/linear/Array2DRowRealMatrix.java
+++ b/src/main/java/org/apache/commons/math4/linear/Array2DRowRealMatrix.java
@@ -390,6 +390,24 @@
 
     /** {@inheritDoc} */
     @Override
+    public RealMatrix getSubMatrix(final int startRow, final int endRow, 
+                                   final int startColumn, final int endColumn)
+            throws OutOfRangeException, NumberIsTooSmallException {
+        MatrixUtils.checkSubMatrixIndex(this, startRow, endRow, startColumn, endColumn);
+        final int rowCount = endRow - startRow + 1;
+        final int columnCount = endColumn - startColumn + 1;
+        final double[][] outData = new double[rowCount][columnCount];
+        for (int i = 0; i < rowCount; ++i) {
+            System.arraycopy(data[startRow + i], startColumn, outData[i], 0, columnCount);
+        }
+        
+        Array2DRowRealMatrix subMatrix = new Array2DRowRealMatrix();
+        subMatrix.data = outData;
+        return subMatrix;
+    }
+
+    /** {@inheritDoc} */
+    @Override
     public double walkInRowOrder(final RealMatrixChangingVisitor visitor) {
         final int rows    = getRowDimension();
         final int columns = getColumnDimension();