Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

svdl: results differ from those of LinearAlgebra.svd #278

Open
severinson opened this issue Sep 6, 2020 · 0 comments
Open

svdl: results differ from those of LinearAlgebra.svd #278

severinson opened this issue Sep 6, 2020 · 0 comments

Comments

@severinson
Copy link

severinson commented Sep 6, 2020

I'm getting some surprising values from svdl (IterativeSolvers.jl 0.8.4). Anyone that knows what's going on here?

Random.seed!(123);
X = randn(10, 6);
Σ, L = svdl(X, vecs=:both);
F = svd(X);

Singular values computed by svdl

Σ.S
 6-element Array{Float64,1}:
 320.9211937151166
 314.34224845686225
 213.92773975361044
 209.19233874206984
 142.11384510704383
 139.89924384072287

Singular values computed by LinearAlgebra.svd

F.S
 6-element Array{Float64,1}:
 5.753436432294746
 4.414616158400172
 3.369985724208918
 1.707195384548486
 1.6125442621275017
 1.0507184615303307

Platform info

versioninfo()
Julia Version 1.5.1
Commit 697e782ab8 (2020-08-25 20:08 UTC)
Platform Info:
  OS: Linux (x86_64-pc-linux-gnu)
  CPU: AMD EPYC 7451 24-Core Processor
  WORD_SIZE: 64
  LIBM: libopenlibm
  LLVM: libLLVM-9.0.1 (ORCJIT, znver1)
Environment:
  JULIA_NUM_THREADS = 48
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant