forked from tsai/mathnet-numerics
11 changed files with 190 additions and 332 deletions
@ -1,101 +0,0 @@ |
|||||
// <copyright file="DenseMatrix.fs" company="Math.NET"> |
|
||||
// Math.NET Numerics, part of the Math.NET Project |
|
||||
// http://mathnet.opensourcedotnet.info |
|
||||
// http://numerics.mathdotnet.com |
|
||||
// http://github.com/mathnet/mathnet-numerics |
|
||||
// http://mathnetnumerics.codeplex.com |
|
||||
// |
|
||||
// Copyright (c) 2009 Math.NET |
|
||||
// |
|
||||
// Permission is hereby granted, free of charge, to any person |
|
||||
// obtaining a copy of this software and associated documentation |
|
||||
// files (the "Software"), to deal in the Software without |
|
||||
// restriction, including without limitation the rights to use, |
|
||||
// copy, modify, merge, publish, distribute, sublicense, and/or sell |
|
||||
// copies of the Software, and to permit persons to whom the |
|
||||
// Software is furnished to do so, subject to the following |
|
||||
// conditions: |
|
||||
// |
|
||||
// The above copyright notice and this permission notice shall be |
|
||||
// included in all copies or substantial portions of the Software. |
|
||||
// |
|
||||
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, |
|
||||
// EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES |
|
||||
// OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND |
|
||||
// NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT |
|
||||
// HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, |
|
||||
// WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING |
|
||||
// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR |
|
||||
// OTHER DEALINGS IN THE SOFTWARE. |
|
||||
// </copyright> |
|
||||
|
|
||||
namespace MathNet.Numerics.LinearAlgebra.Double |
|
||||
open MathNet.Numerics.LinearAlgebra |
|
||||
open MathNet.Numerics.LinearAlgebra.Generic |
|
||||
|
|
||||
/// A module which implements functional dense vector operations. |
|
||||
module DenseMatrix = |
|
||||
|
|
||||
/// Initialize a matrix by calling a construction function for every element. |
|
||||
let inline init (n: int) (m: int) f = |
|
||||
let A = new DenseMatrix(n,m) |
|
||||
for i=0 to n-1 do |
|
||||
for j=0 to m-1 do |
|
||||
A.[i,j] <- f i j |
|
||||
A |
|
||||
|
|
||||
/// Create a matrix from a list of float lists. Every list in the master list specifies a row. |
|
||||
let inline ofList (fll: float list list) = |
|
||||
let n = List.length fll |
|
||||
let m = List.length (List.head fll) |
|
||||
let A = DenseMatrix(n,m) |
|
||||
fll |> List.iteri (fun i fl -> |
|
||||
if (List.length fl) <> m then failwith "Each subrow must be of the same length." else |
|
||||
List.iteri (fun j f -> A.[i,j] <- f) fl) |
|
||||
A |
|
||||
|
|
||||
/// Create a matrix from a list of sequences. Every sequence in the master sequence specifies a row. |
|
||||
let inline ofSeq (fss: #seq<#seq<float>>) = |
|
||||
let n = Seq.length fss |
|
||||
let m = Seq.length (Seq.head fss) |
|
||||
let A = DenseMatrix(n,m) |
|
||||
fss |> Seq.iteri (fun i fs -> |
|
||||
if (Seq.length fs) <> m then failwith "Each subrow must be of the same length." else |
|
||||
Seq.iteri (fun j f -> A.[i,j] <- f) fs) |
|
||||
A |
|
||||
|
|
||||
/// Create a matrix from a 2D array of floating point numbers. |
|
||||
let inline ofArray2 (arr: float[,]) = new DenseMatrix(arr) |
|
||||
|
|
||||
/// Create a matrix with the given entries. |
|
||||
let inline initDense (n: int) (m: int) (es: #seq<int * int * float>) = |
|
||||
let A = new DenseMatrix(n,m) |
|
||||
Seq.iter (fun (i,j,f) -> A.[i,j] <- f) es |
|
||||
A |
|
||||
|
|
||||
/// Create a square matrix with constant diagonal entries. |
|
||||
let inline constDiag (n: int) (f: float) = |
|
||||
let A = new DenseMatrix(n,n) |
|
||||
for i=0 to n-1 do |
|
||||
A.[i,i] <- f |
|
||||
A |
|
||||
|
|
||||
/// Create a square matrix with the vector elements on the diagonal. |
|
||||
let inline diag (v: #Vector<float>) = |
|
||||
let n = v.Count |
|
||||
let A = new DenseMatrix(n,n) |
|
||||
for i=0 to n-1 do |
|
||||
A.[i,i] <- v.Item(i) |
|
||||
A |
|
||||
|
|
||||
/// Initialize a matrix by calling a construction function for every row. |
|
||||
let inline initRow (n: int) (m: int) (f: int -> #Vector<float>) = |
|
||||
let A = new DenseMatrix(n,m) |
|
||||
for i=0 to n-1 do A.SetRow(i, f i) |
|
||||
A |
|
||||
|
|
||||
/// Initialize a matrix by calling a construction function for every column. |
|
||||
let inline initCol (n: int) (m: int) (f: int -> #Vector<float>) = |
|
||||
let A = new DenseMatrix(n,m) |
|
||||
for i=0 to m-1 do A.SetColumn(i, f i) |
|
||||
A |
|
||||
@ -1,69 +0,0 @@ |
|||||
// <copyright file="DenseVector.fs" company="Math.NET"> |
|
||||
// Math.NET Numerics, part of the Math.NET Project |
|
||||
// http://numerics.mathdotnet.com |
|
||||
// http://github.com/mathnet/mathnet-numerics |
|
||||
// http://mathnetnumerics.codeplex.com |
|
||||
// |
|
||||
// Copyright (c) 2009 Math.NET |
|
||||
// |
|
||||
// Permission is hereby granted, free of charge, to any person |
|
||||
// obtaining a copy of this software and associated documentation |
|
||||
// files (the "Software"), to deal in the Software without |
|
||||
// restriction, including without limitation the rights to use, |
|
||||
// copy, modify, merge, publish, distribute, sublicense, and/or sell |
|
||||
// copies of the Software, and to permit persons to whom the |
|
||||
// Software is furnished to do so, subject to the following |
|
||||
// conditions: |
|
||||
// |
|
||||
// The above copyright notice and this permission notice shall be |
|
||||
// included in all copies or substantial portions of the Software. |
|
||||
// |
|
||||
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, |
|
||||
// EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES |
|
||||
// OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND |
|
||||
// NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT |
|
||||
// HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, |
|
||||
// WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING |
|
||||
// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR |
|
||||
// OTHER DEALINGS IN THE SOFTWARE. |
|
||||
// </copyright> |
|
||||
|
|
||||
namespace MathNet.Numerics.LinearAlgebra.Double |
|
||||
|
|
||||
open MathNet.Numerics.LinearAlgebra |
|
||||
|
|
||||
/// A module which implements functional dense vector operations. |
|
||||
module DenseVector = |
|
||||
|
|
||||
/// Initialize a vector by calling a construction function for every element. |
|
||||
let inline init (n: int) f = |
|
||||
let v = new Double.DenseVector(n) |
|
||||
for i=0 to n-1 do |
|
||||
v.[i] <- f i |
|
||||
v |
|
||||
|
|
||||
/// Create a vector from a float list. |
|
||||
let inline ofList (fl: float list) = |
|
||||
let n = List.length fl |
|
||||
let v = Double.DenseVector(n) |
|
||||
fl |> List.iteri (fun i f -> v.[i] <- f) |
|
||||
v |
|
||||
|
|
||||
/// Create a vector from a sequences. |
|
||||
let inline ofSeq (fs: #seq<float>) = |
|
||||
let n = Seq.length fs |
|
||||
let v = DenseVector(n) |
|
||||
fs |> Seq.iteri (fun i f -> v.[i] <- f) |
|
||||
v |
|
||||
|
|
||||
/// Create a vector with evenly spaced entries: e.g. rangef -1.0 0.5 1.0 = [-1.0 -0.5 0.0 0.5 1.0] |
|
||||
let inline rangef (start: float) (step: float) (stop: float) = |
|
||||
let n = (int ((stop - start) / step)) + 1 |
|
||||
let v = new DenseVector(n) |
|
||||
for i=0 to n-1 do |
|
||||
v.[i] <- (float i) * step + start |
|
||||
v |
|
||||
|
|
||||
/// Create a vector with integer entries in the given range. |
|
||||
let inline range (start: int) (stop: int) = |
|
||||
new DenseVector([| for i in [start .. stop] -> float i |]) |
|
||||
@ -1,10 +1,10 @@ |
|||||
// <copyright file="Extensions.fs" company="Math.NET"> |
// <copyright file="LinearAlgebra.fs" company="Math.NET"> |
||||
// Math.NET Numerics, part of the Math.NET Project |
// Math.NET Numerics, part of the Math.NET Project |
||||
// http://numerics.mathdotnet.com |
// http://numerics.mathdotnet.com |
||||
// http://github.com/mathnet/mathnet-numerics |
// http://github.com/mathnet/mathnet-numerics |
||||
// http://mathnetnumerics.codeplex.com |
// http://mathnetnumerics.codeplex.com |
||||
// |
// |
||||
// Copyright (c) 2009 Math.NET |
// Copyright (c) 2009-2012 Math.NET |
||||
// |
// |
||||
// Permission is hereby granted, free of charge, to any person |
// Permission is hereby granted, free of charge, to any person |
||||
// obtaining a copy of this software and associated documentation |
// obtaining a copy of this software and associated documentation |
||||
@ -1,75 +0,0 @@ |
|||||
// <copyright file="SparseMatrix.fs" company="Math.NET"> |
|
||||
// Math.NET Numerics, part of the Math.NET Project |
|
||||
// http://numerics.mathdotnet.com |
|
||||
// http://github.com/mathnet/mathnet-numerics |
|
||||
// http://mathnetnumerics.codeplex.com |
|
||||
// |
|
||||
// Copyright (c) 2009 Math.NET |
|
||||
// |
|
||||
// Permission is hereby granted, free of charge, to any person |
|
||||
// obtaining a copy of this software and associated documentation |
|
||||
// files (the "Software"), to deal in the Software without |
|
||||
// restriction, including without limitation the rights to use, |
|
||||
// copy, modify, merge, publish, distribute, sublicense, and/or sell |
|
||||
// copies of the Software, and to permit persons to whom the |
|
||||
// Software is furnished to do so, subject to the following |
|
||||
// conditions: |
|
||||
// |
|
||||
// The above copyright notice and this permission notice shall be |
|
||||
// included in all copies or substantial portions of the Software. |
|
||||
// |
|
||||
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, |
|
||||
// EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES |
|
||||
// OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND |
|
||||
// NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT |
|
||||
// HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, |
|
||||
// WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING |
|
||||
// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR |
|
||||
// OTHER DEALINGS IN THE SOFTWARE. |
|
||||
// </copyright> |
|
||||
|
|
||||
namespace MathNet.Numerics.LinearAlgebra.Double |
|
||||
|
|
||||
open MathNet.Numerics.LinearAlgebra |
|
||||
open MathNet.Numerics.LinearAlgebra.Generic |
|
||||
/// A module which implements functional sparse vector operations. |
|
||||
module SparseMatrix = |
|
||||
|
|
||||
/// Create a matrix from a list of float lists. Every list in the master list specifies a row. |
|
||||
let inline ofList (rows: int) (cols: int) (fll: list<int * int * float>) = |
|
||||
let A = new Double.SparseMatrix(rows, cols) |
|
||||
fll |> List.iter (fun (i, j, x) -> A.[i,j] <- x) |
|
||||
A |
|
||||
|
|
||||
/// Create a matrix from a list of sequences. Every sequence in the master sequence specifies a row. |
|
||||
let inline ofSeq (rows: int) (cols: int) (fss: #seq<int * int * float>) = |
|
||||
let A = new Double.SparseMatrix(rows, cols) |
|
||||
fss |> Seq.iter (fun (i, j, x) -> A.[i,j] <- x) |
|
||||
A |
|
||||
|
|
||||
/// Create a square matrix with constant diagonal entries. |
|
||||
let inline constDiag (n: int) (f: float) = |
|
||||
let A = new Double.SparseMatrix(n,n) |
|
||||
for i=0 to n-1 do |
|
||||
A.[i,i] <- f |
|
||||
A |
|
||||
|
|
||||
/// Create a square matrix with the vector elements on the diagonal. |
|
||||
let inline diag (v: #Vector<float>) = |
|
||||
let n = v.Count |
|
||||
let A = new Double.SparseMatrix(n,n) |
|
||||
for i=0 to n-1 do |
|
||||
A.[i,i] <- v.Item(i) |
|
||||
A |
|
||||
|
|
||||
/// Initialize a matrix by calling a construction function for every row. |
|
||||
let inline initRow (n: int) (m: int) (f: int -> #Vector<float>) = |
|
||||
let A = new Double.SparseMatrix(n,m) |
|
||||
for i=0 to n-1 do A.SetRow(i, f i) |
|
||||
A |
|
||||
|
|
||||
/// Initialize a matrix by calling a construction function for every column. |
|
||||
let inline initCol (n: int) (m: int) (f: int -> #Vector<float>) = |
|
||||
let A = new Double.SparseMatrix(n,m) |
|
||||
for i=0 to m-1 do A.SetColumn(i, f i) |
|
||||
A |
|
||||
@ -1,48 +0,0 @@ |
|||||
// <copyright file="SparseVector.fs" company="Math.NET"> |
|
||||
// Math.NET Numerics, part of the Math.NET Project |
|
||||
// http://numerics.mathdotnet.com |
|
||||
// http://github.com/mathnet/mathnet-numerics |
|
||||
// http://mathnetnumerics.codeplex.com |
|
||||
// |
|
||||
// Copyright (c) 2009 Math.NET |
|
||||
// |
|
||||
// Permission is hereby granted, free of charge, to any person |
|
||||
// obtaining a copy of this software and associated documentation |
|
||||
// files (the "Software"), to deal in the Software without |
|
||||
// restriction, including without limitation the rights to use, |
|
||||
// copy, modify, merge, publish, distribute, sublicense, and/or sell |
|
||||
// copies of the Software, and to permit persons to whom the |
|
||||
// Software is furnished to do so, subject to the following |
|
||||
// conditions: |
|
||||
// |
|
||||
// The above copyright notice and this permission notice shall be |
|
||||
// included in all copies or substantial portions of the Software. |
|
||||
// |
|
||||
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, |
|
||||
// EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES |
|
||||
// OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND |
|
||||
// NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT |
|
||||
// HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, |
|
||||
// WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING |
|
||||
// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR |
|
||||
// OTHER DEALINGS IN THE SOFTWARE. |
|
||||
// </copyright> |
|
||||
|
|
||||
namespace MathNet.Numerics.LinearAlgebra.Double |
|
||||
|
|
||||
open MathNet.Numerics.LinearAlgebra |
|
||||
|
|
||||
/// A module which implements functional sparse vector operations. |
|
||||
module SparseVector = |
|
||||
|
|
||||
/// Create a sparse vector with a given dimension from a list of entry, value pairs. |
|
||||
let inline ofList (dim: int) (fl: list<int * float>) = |
|
||||
let v = new Double.SparseVector(dim) |
|
||||
fl |> List.iter (fun (i, f) -> v.[i] <- f) |
|
||||
v |
|
||||
|
|
||||
/// Create a sparse vector with a given dimension from a sequence of entry, value pairs. |
|
||||
let inline ofSeq (dim: int) (fs: #seq<int * float>) = |
|
||||
let v = new Double.SparseVector(dim) |
|
||||
fs |> Seq.iter (fun (i, f) -> v.[i] <- f) |
|
||||
v |
|
||||
Loading…
Reference in new issue