Skip to content

Conversation

@martinholters
Copy link
Member

@martinholters martinholters added the linear algebra Linear algebra label Jul 3, 2018
@martinholters martinholters merged commit 8b4c7bb into master Jul 4, 2018
@martinholters martinholters deleted the mh/svd_adjoint branch July 4, 2018 13:29
@test_throws ErrorException usv.Z
b = rand(eltya,n)
@test usv\b transpose(a)\b
end
Copy link
Member

@StefanKarpinski StefanKarpinski Jul 5, 2018

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I would think that the most obvious tests here would be these:

svd(a').U == svd(a).V
svd(a').S == svd(a).S
svd(a').V == svd(a).U

Is there some reason not to test this directly? I realize that the precise U and V matrices are not well defined in many cases but it seems unlikely that this will ever break in this case, no?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

linear algebra Linear algebra

Projects

None yet

Development

Successfully merging this pull request may close these issues.

svdfact fails on transposed matrices (julia 0.7.0-DEV.5110)

4 participants