good thought -- I must be missing an expected function
julia> #Pkg.clone("git://github.com/JuliaArbTypes/ArbFloats.jl.git")
julia> using GenericSVD, ArbFloats
julia> setprecision(ArbFloat,53);
julia> n,m = 4,2;
julia> srand(1);
julia> X = randn(n,m)
4×2 Array{Float64,2}:
0.297288 -0.839027
0.382396 0.311111
-0.597634 2.29509
-0.0104452 -2.26709
julia> A = reshape(ArbFloat(vec(X)),(n,m))
4×2 Array{ArbFloats.ArbFloat{53},2}:
0.2972879845354616 -0.839026854388764
0.3823959677906078 0.3111113384983338
-0.5976344767282311 2.29508782383731
-0.01044524463737564 -2.267086348800531
julia> svdvals(X)
2-element Array{Float64,1}:
2.87552
0.791637
julia> svdvals(A)
ERROR: MethodError: no method matching
svd!(::Bidiagonal{ArbFloats.ArbFloat{53}}, ::Void, ::Void, ::Float64)
On Wednesday, August 24, 2016 at 1:08:43 PM UTC-4, Tony Kelman wrote:
>
> Try https://github.com/simonbyrne/GenericSVD.jl ?
