From: Peter Schmitteckert <support@schmitteckert.com>
To: 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: <01091407482400.01118@baxter> (raw)
In-Reply-To: <15265.6178.106122.979002@debian>
Dear All,
please excuse me, I have used gsl only once and havn't read the documents,
but I got curious by the follwing comment:
On Thursday 13 September 2001 22:33, Brian Gough wrote:
> Yes, the right singular vectors are given by the columns of V. The
> third column is the appropriate one.
>
Isn't that counter-intuitive fo a C programmer ?
In the original post I found no hints on a desired
Fortran semantic for the memory layout. But I would
expect that the default semantics of a C matrix library is
C style, so eigen vectors / singular vectors should be in
the rows, otherwise they are not in a contiguous memory
block.
Best wishes
Peter
P.S:
I just looked up the documentation:
> Matrices are stored in row-major order, meaning that each row of elements
> forms a contiguous block in memory. This is the standard "C-language
> ordering" of two-dimensional arrays. Note that FORTRAN stores arrays in
> column-major order. The number of rows is size1. The range of valid
therefore, returnning vectors in the coulmn should be considered
as a bug. Have I missed something ?
Note, that the documentation says: "The matrix Q contains the elements of Q
in untransposed form".
next prev parent 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
2001-12-19 13:20 ` James Theiler
2001-12-19 13:20 ` Brian Gough
2001-12-19 13:20 ` Peter Schmitteckert [this message]
2001-12-19 13:20 ` Brian Gough
-- strict thread matches above, loose matches on Subject: below --
2001-12-19 13:20 Jim Love
2001-12-19 13:20 ` Brian Gough
[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 Jim Love
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
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=01091407482400.01118@baxter \
--to=support@schmitteckert.com \
--cc=gsl-discuss@sources.redhat.com \
/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).