Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

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


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

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]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.