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 ]
Francesco Perrone

Posts: 39
Registered: 5/2/12
Re: Lapack spptrf function
Posted: Apr 11, 2013 3:07 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

Dear Mr. Tursa,

I have a certain n*n matrix. Then calling

my_matrix = nn_matrix(itril(size(nn_matrix)));

I am able to wrap it into a lower triangular factor.

Afterwards, I manipulate my_matrix by means of a coherence function; the outcome, for only one step of the above mentioned FOR loop, is

Coh_u = exp(-a.*my_matrix.*sqrt((f(ii)/Uhub).^2 + (0.12/Lc).^2)).*(df.*psd(ii,1));

Therefore, Coh_u is alread packed as spptrf input.

Then I call the lapack wrapper as follows:

[C1u,C2u,HH_u,C3u] = lapack('spptrf','L',size(nn_matrix,1),Coh_u,0);

HH_u is the expected result.

I would more than appreciate if you could write a mex file of the spptrf function for me. Unfortunately I started only now learning some C++ and it would require a long time for me to catch up with mex file.

Look really forward to hearing from you.

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.