From: "Jim Love" <Jim.Love@asml.com>
To: <jt@lanl.gov>
Cc: <gsl-discuss@sources.redhat.com>
Subject: Re: Problem with Singular Value Decomposition Algorithm
Date: Wed, 19 Dec 2001 13:20:00 -0000 [thread overview]
Message-ID: <sb9f4602.003@wiltonhub.svg.com> (raw)
In every text that I have read dealing with the process of SVD, they have always placed the eigenvalues in a decreasing order. The stated reason is to provide a standard solution and methodology. Why would anybody build an API that does something different?? If you are using the SVD to find the least squares fit of a plane or a cylinder for example, order matters (if you want the right answer hehehe). Yes I can add to my code to re-order the output, but so will most everybody else that uses this function. So it should be done in the API. Just a suggestion.
James A. Love
X4477
Pager 1-800-286-1188 Pin# 400659
>>> James Theiler <jt@lanl.gov> 09/12/01 10:51AM >>>
On Wed, 12 Sep 2001, Jim Love wrote:
] I have downloaded the latest beta release and the SVD algorithm
] produces the wrong answers. It appears that columns are swapped
] in the output and possibly a sign problem.
]
] Here was my test array:
]
] 1 1 0.975
] 1 -1 0.975
] -1 -1 -0.925
] -1 1 -1.025
]
] The correct answer for the S vector is: 2.7940 2.0000 0.0358
] The gls output was: 2.0000 2.7940 0.0358
]
] The Correct Q matrix is:
]
] -0.7155 0.0256 -0.6981
] 0.0183 0.9997 0.0179
] -0.6983 -0.0000 0.7158
]
] The gls output was:
]
] -0.025633 -0.715538 -0.698103
] -0.999671 0.018347 0.017900
] -0.000000 -0.698332 0.715774
]
] A similar problem was seen in the U matrix.
]
] Any ideas? Is this caused by my implementation or is it a real bug?
I don't believe this is a bug at all.
Both answers are correct. For SVD, the algorithm is
asked to find matrices U,S,Q such that U.S.Q' equals
the original matrix. If you swap the columns of Q
and swap the equivalent rows of U, and also swap the
corresponding elements of S, you have another solution.
Also if you multiply one of the columns of Q by -1,
and multiply the corresponding row of U by -1, you
will get another equivalent solution.
Sometimes you want the solution with the eigenvalues (the
diagonal elements of the S matrix -- represented as a vector)
sorted numerically, as the one you cite as "correct" but
that is a convenience that can be performed after the fact.
And in fact, I notice that the documentation says that
eigenvalues are "generally chosen" to form a non-decreasing
sequence. Maybe worth a sentence to say that these algorithms
do not always do that, or else add a utility that does this
postprocessing step.
jt
---------------------------------------------
James Theiler jt@lanl.gov
MS-D436, NIS-2, LANL tel: 505/665-5682
Los Alamos, NM 87545 fax: 505/665-4414
----- Space and Remote Sensing Sciences -----
next reply other threads:[~2001-12-19 13:20 UTC|newest]
Thread overview: 15+ messages / expand[flat|nested] mbox.gz Atom feed top
2001-12-19 13:20 Jim Love [this message]
2001-12-19 13:20 ` Brian Gough
2001-12-19 13:20 ` Brian Gough
2001-12-19 13:20 Jim Love
2001-12-19 13:20 ` James Theiler
2001-12-19 13:20 Jim Love
2001-12-19 13:20 Jim Love
2001-12-19 13:20 ` James Theiler
2001-12-19 13:20 ` Brian Gough
2001-12-19 13:20 ` Peter Schmitteckert
2001-12-19 13:20 ` Brian Gough
2001-12-19 13:20 Jim Love
[not found] <sba09ba1.005@wiltonhub.svg.com>
2001-12-19 13:20 ` Brian Gough
2001-12-19 13:20 Jim Love
2001-12-19 13:20 ` Brian Gough
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=sb9f4602.003@wiltonhub.svg.com \
--to=jim.love@asml.com \
--cc=gsl-discuss@sources.redhat.com \
--cc=jt@lanl.gov \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).