diff --git a/src/Numerics/LinearAlgebra/Complex/DiagonalMatrix.cs b/src/Numerics/LinearAlgebra/Complex/DiagonalMatrix.cs
index 9f58d067..5a86975d 100644
--- a/src/Numerics/LinearAlgebra/Complex/DiagonalMatrix.cs
+++ b/src/Numerics/LinearAlgebra/Complex/DiagonalMatrix.cs
@@ -301,11 +301,11 @@ namespace MathNet.Numerics.LinearAlgebra.Complex
             Matrix<Complex> result;
             if (other is DiagonalMatrix)
             {
-                result = new DenseMatrix(RowCount, ColumnCount);
+                result = new DiagonalMatrix(RowCount, ColumnCount);
             }
             else
             {
-                result = new DiagonalMatrix(RowCount, ColumnCount);
+                result = new DenseMatrix(RowCount, ColumnCount);
             }
 
             Add(other, result);
diff --git a/src/Numerics/LinearAlgebra/Complex32/DiagonalMatrix.cs b/src/Numerics/LinearAlgebra/Complex32/DiagonalMatrix.cs
index d765e2c8..b1a3e7d6 100644
--- a/src/Numerics/LinearAlgebra/Complex32/DiagonalMatrix.cs
+++ b/src/Numerics/LinearAlgebra/Complex32/DiagonalMatrix.cs
@@ -306,11 +306,11 @@ namespace MathNet.Numerics.LinearAlgebra.Complex32
             Matrix<Complex32> result;
             if (other is DiagonalMatrix)
             {
-                result = new DenseMatrix(RowCount, ColumnCount);
+                result = new DiagonalMatrix(RowCount, ColumnCount);
             }
             else
             {
-                result = new DiagonalMatrix(RowCount, ColumnCount);
+                result = new DenseMatrix(RowCount, ColumnCount);
             }
 
             Add(other, result);
@@ -381,11 +381,11 @@ namespace MathNet.Numerics.LinearAlgebra.Complex32
             Matrix<Complex32> result;
             if (other is DiagonalMatrix)
             {
-                result = new DenseMatrix(RowCount, ColumnCount);
+                result = new DiagonalMatrix(RowCount, ColumnCount);
             }
             else
             {
-                result = new DiagonalMatrix(RowCount, ColumnCount);
+                result = new DenseMatrix(RowCount, ColumnCount);
             }
 
             Subtract(other, result);
diff --git a/src/Numerics/LinearAlgebra/Double/DiagonalMatrix.cs b/src/Numerics/LinearAlgebra/Double/DiagonalMatrix.cs
index ffcdb900..6922b101 100644
--- a/src/Numerics/LinearAlgebra/Double/DiagonalMatrix.cs
+++ b/src/Numerics/LinearAlgebra/Double/DiagonalMatrix.cs
@@ -300,11 +300,11 @@ namespace MathNet.Numerics.LinearAlgebra.Double
             Matrix<double> result;
             if (other is DiagonalMatrix)
             {
-                result = new DenseMatrix(RowCount, ColumnCount);
+                result = new DiagonalMatrix(RowCount, ColumnCount);
             }
             else
             {
-                result = new DiagonalMatrix(RowCount, ColumnCount);
+                result = new DenseMatrix(RowCount, ColumnCount);
             }
 
             Add(other, result);
@@ -375,11 +375,11 @@ namespace MathNet.Numerics.LinearAlgebra.Double
             Matrix<double> result;
             if (other is DiagonalMatrix)
             {
-                result = new DenseMatrix(RowCount, ColumnCount);
+                result = new DiagonalMatrix(RowCount, ColumnCount);
             }
             else
             {
-                result = new DiagonalMatrix(RowCount, ColumnCount);
+                result = new DenseMatrix(RowCount, ColumnCount);
             }
 
             Subtract(other, result);
diff --git a/src/Numerics/LinearAlgebra/Single/DiagonalMatrix.cs b/src/Numerics/LinearAlgebra/Single/DiagonalMatrix.cs
index 5bac5ed6..55c5e45c 100644
--- a/src/Numerics/LinearAlgebra/Single/DiagonalMatrix.cs
+++ b/src/Numerics/LinearAlgebra/Single/DiagonalMatrix.cs
@@ -300,11 +300,11 @@ namespace MathNet.Numerics.LinearAlgebra.Single
             Matrix<float> result;
             if (other is DiagonalMatrix)
             {
-                result = new DenseMatrix(RowCount, ColumnCount);
+                result = new DiagonalMatrix(RowCount, ColumnCount);
             }
             else
             {
-                result = new DiagonalMatrix(RowCount, ColumnCount);
+                result = new DenseMatrix(RowCount, ColumnCount);
             }
 
             Add(other, result);
@@ -375,11 +375,11 @@ namespace MathNet.Numerics.LinearAlgebra.Single
             Matrix<float> result;
             if (other is DiagonalMatrix)
             {
-                result = new DenseMatrix(RowCount, ColumnCount);
+                result = new DiagonalMatrix(RowCount, ColumnCount);
             }
             else
             {
-                result = new DiagonalMatrix(RowCount, ColumnCount);
+                result = new DenseMatrix(RowCount, ColumnCount);
             }
 
             Subtract(other, result);