How can i speed up this part of my code?

1 次查看（过去 30 天）

amanita 2013-4-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/72662-how-can-i-speed-up-this-part-of-my-code

I have this part of code that i apply for different indns (which are different vectors each time):

    phi1=phi(indn1);
    p1=p(indn1,indn1);
    rtn1=(phi1'*(p1*phi1));
    Ctn1=((rtn1*ep2)/(st*(1+rtn1)));
    Ctn1=max(Ctn1,0);
    Sn1=gammainc(Ctn1*0.5,0.5) ;
    Sn1=1-Sn1;
    phi2=phi(indn2);
    ....
    ....

I tried cellfun and structfun but it was slower. Is there something more i can do to improve the speed?

5 个评论
显示 3更早的评论隐藏 3更早的评论

amanita 2013-4-19

thanx Cedric

amanita 2013-6-16

The problem was solved by creating mex files and calling them through the GA btw.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Roger Stafford 2013-6-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/72662-how-can-i-speed-up-this-part-of-my-code#answer_88972

在 MATLAB Online 中打开

In your code the multiplication part of doing

phi1'*p1*phi1

is being repeated over and over again as you select different index combinations. Perhaps you could speed up things, even in your mex method, by doing things in the following way. (I am assuming phi is a column vector: one column, with the same number of rows as there are rows and columns in p.)

First, compute, one time only, the matrix

A = (phi*phi').*p;

Then for each new indn vector as it comes along do only this

rtn = sum(sum(A(indn,indn)));

rtn = sum(reshape(A(indn,indn),[],1));

whichever is faster. This step requires only addition operations rather than the multiplications plus additions required for matrix multiplication, and gives the same result.

Also in gammainc use the 'upper' option instead of doing Sn1=1-Sn1.