From 8fa7e1211fb7d40396a6355e84b4f7ffeb76feb0 Mon Sep 17 00:00:00 2001
From: kalmarek <kalmar@amu.edu.pl>
Date: Wed, 5 Sep 2018 10:40:08 +0200
Subject: [PATCH] move orthSVD close to OrbitData

---
 src/Orbit-wise.jl         | 6 ++++++
 src/OrbitDecomposition.jl | 7 -------
 2 files changed, 6 insertions(+), 7 deletions(-)

diff --git a/src/Orbit-wise.jl b/src/Orbit-wise.jl
index 167d612..f36cf95 100644
--- a/src/Orbit-wise.jl
+++ b/src/Orbit-wise.jl
@@ -107,6 +107,12 @@ function save_preps(fname::String, preps)
     save(fname, "perms_d", [preps[elt].d for elt in elements(autS)])
 end
 
+function orthSVD(M::AbstractMatrix{T}) where {T<:AbstractFloat}
+    M = full(M)
+    fact = svdfact(M)
+    M_rank = sum(fact[:S] .> maximum(size(M))*eps(T))
+    return fact[:U][:,1:M_rank]
+end
 
 ###############################################################################
 #
diff --git a/src/OrbitDecomposition.jl b/src/OrbitDecomposition.jl
index 562f1f4..0589180 100644
--- a/src/OrbitDecomposition.jl
+++ b/src/OrbitDecomposition.jl
@@ -72,13 +72,6 @@ function matrix_repr(x::GroupRingElem, mreps::Dict)
     return sum(x[i].*mreps[parent(x).basis[i]] for i in nzeros)
 end
 
-function orthSVD(M::AbstractMatrix{T}) where {T<:AbstractFloat}
-    M = full(M)
-    fact = svdfact(M)
-    M_rank = sum(fact[:S] .> maximum(size(M))*eps(T))
-    return fact[:U][:,1:M_rank]
-end
-
 function matrix_reps(preps::Dict{T,perm{I}}) where {T<:GroupElem, I<:Integer}
     kk = collect(keys(preps))
     mreps = Vector{SparseMatrixCSC{Float64, Int}}(length(kk))