From cb58405b853388f0a1923dd10225fcfe054af3a1 Mon Sep 17 00:00:00 2001
From: Christoph Ruegg <git@cdrnet.ch>
Date: Tue, 17 Sep 2013 21:55:55 +0200
Subject: [PATCH] LA: add direct (factorization-based) Solve methods to matrix
 type

---
 src/Numerics/LinearAlgebra/Matrix.Solve.cs | 69 +++++++++++++++++++++-
 1 file changed, 68 insertions(+), 1 deletion(-)

diff --git a/src/Numerics/LinearAlgebra/Matrix.Solve.cs b/src/Numerics/LinearAlgebra/Matrix.Solve.cs
index c41994ab..70ad8a2c 100644
--- a/src/Numerics/LinearAlgebra/Matrix.Solve.cs
+++ b/src/Numerics/LinearAlgebra/Matrix.Solve.cs
@@ -81,7 +81,73 @@ namespace MathNet.Numerics.LinearAlgebra
         public abstract Evd<T> Evd();
 
 
-        // Iterative Solvers: Try
+
+        // Direct Solvers: Full
+
+        /// <summary>
+        /// Solves a system of linear equations, <b>Ax = b</b>, with A QR factorized.
+        /// </summary>
+        /// <param name="input">The right hand side vector, <b>b</b>.</param>
+        /// <param name="result">The left hand side <see cref="Matrix{T}"/>, <b>x</b>.</param>
+        public void Solve(Vector<T> input, Vector<T> result)
+        {
+            if (ColumnCount == RowCount)
+            {
+                LU().Solve(input, result);
+                return;
+            }
+
+            QR().Solve(input, result);
+        }
+
+        /// <summary>
+        /// Solves a system of linear equations, <b>AX = B</b>, with A QR factorized.
+        /// </summary>
+        /// <param name="input">The right hand side <see cref="Matrix{T}"/>, <b>B</b>.</param>
+        /// <param name="result">The left hand side <see cref="Matrix{T}"/>, <b>X</b>.</param>
+        public void Solve(Matrix<T> input, Matrix<T> result)
+        {
+            if (ColumnCount == RowCount)
+            {
+                LU().Solve(input, result);
+                return;
+            }
+
+            QR().Solve(input, result);
+        }
+
+
+
+        // Direct Solvers: Simple
+
+        /// <summary>
+        /// Solves a system of linear equations, <b>AX = B</b>, with A QR factorized.
+        /// </summary>
+        /// <param name="input">The right hand side <see cref="Matrix{T}"/>, <b>B</b>.</param>
+        /// <returns>The left hand side <see cref="Matrix{T}"/>, <b>X</b>.</returns>
+        public Matrix<T> Solve(Matrix<T> input)
+        {
+            var matrixX = CreateMatrix(ColumnCount, input.ColumnCount);
+            Solve(input, matrixX);
+            return matrixX;
+        }
+
+
+        /// <summary>
+        /// Solves a system of linear equations, <b>Ax = b</b>, with A QR factorized.
+        /// </summary>
+        /// <param name="input">The right hand side vector, <b>b</b>.</param>
+        /// <returns>The left hand side <see cref="Vector{T}"/>, <b>x</b>.</returns>
+        public Vector<T> Solve(Vector<T> input)
+        {
+            var x = CreateVector(ColumnCount);
+            Solve(input, x);
+            return x;
+        }
+
+
+
+        // Iterative Solvers: Full
 
         /// <summary>
         /// Solves the matrix equation Ax = b, where A is the coefficient matrix (this matrix), b is the solution vector and x is the unknown vector.
@@ -197,6 +263,7 @@ namespace MathNet.Numerics.LinearAlgebra
         }
 
 
+
         // Iterative Solvers: Simple
 
         public Vector<T> SolveIterative(Vector<T> input, IIterativeSolver<T> solver, IPreconditioner<T> preconditioner, Iterator<T> iterator)