The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Software » comp.soft-sys.matlab

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Lapack spptrf function
Replies: 10   Last Post: Apr 12, 2013 11:08 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Steven Lord

Posts: 18,038
Registered: 12/7/04
Re: Lapack spptrf function
Posted: Apr 10, 2013 3:47 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

"Francesco Perrone" <> wrote in message
> Dear Mr. Lord,
> yes I am aware of the function CHOL and I am already using it.
> But CHOL gets tremendously slow when factorizing, e.g., a 1945*1945
> matrix.

Please define "tremendously slow." When I use a sample SPD matrix on my
machine, it takes about a tenth of a second. YTMV.

A = gallery('minij', 1945);
[R, p] = chol(A);

> Within my code I call CHOL, within a 2048 iterations FOR loop, for three
> different matrices with the above mentioned size.
> The elapsed time is about an hour.

So 3*2048 = 6144 calls to CHOL, averaging about half a second per call,
assuming that's the only thing in the loop? Depending on the hardware
configuration of your machine, that may not be too bad.

> For the same case, after downloading the LAPACK wrapper
> (, the
> elapsed time was around 17 min less than before: quite a bargain I'd urge.
> Btw, after a bunch of tests, MATLAB started crashing at any time I tried
> calling the function: blatantly, I don't even have a clue why that
> happened.

By "the function" do you mean CHOL or the LAPACK wrapper function from the
File Exchange? If the former, please send the crash log, a MAT-file
containing the matrix, and the exact command you used to Technical Support
so we can investigate the cause of the crash. If you mean the LAPACK wrapper
function, that suggests that it's not operating correctly. I'd be cautious
about using it until you and/or the developer can determine the cause of the
crash and correct it.

Steve Lord
To contact Technical Support use the Contact Us link on
YTMV = 'Your Timing May Vary"

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.