How can i speed up my function?

2 次查看（过去 30 天）
Gerrit
Gerrit 2020-1-8
评论： Gerrit 2020-1-12
Hello everybody,
I'm trying to increase the performance of my Matlab function. Input is a table with the same columns as airfoil_update.
The function splits the X-Coordinates of one AeroDataId at the minimum of X/C and does this for every AeroDataId in the table. Airfoil, the input table, will have a size of 131123421x4. But the fourth column "Upper" is empty and should be filled up with a 1, if the coordinate belongs to the upper airfoil and 0, if it belongs to the lower airfoil.
Thank you very much in advance.
function [airfoil_update] = splitAirfoil(airfoil)
%% create table for upper and lower airfoil
airfoil_update=array2table(zeros(0,4),'VariableNames',{'AeroDataId','X_C','CP','Upper'});
%% vector with unique AeroDataIds
DataId_unique=unique(airfoil.AeroDataId);
for i=1:1:size(DataId_unique)
%% temp vector filled with every point of specific AeroDataId
disp(DataId_unique(i,:));
AeroDataId=airfoil(airfoil.AeroDataId==DataId_unique(i,:),:);
%% find the minimum of X/C
[~, min_idx]=min(AeroDataId.X_C);
%% every point with same ID before minimum -> lower_airfoil -> Upper =0
airfoil_new=AeroDataId(1:min_idx,:);
airfoil_new.Upper(:,:)=0;
airfoil_update=[airfoil_update; airfoil_new];
%% every point with same ID after minimum = upper_airfoil -> Upper =1
airfoil_new=AeroDataId(min_idx+1:end,:);
airfoil_new.Upper(:,:)=1;
%% add new rows
airfoil_update=[airfoil_update; airfoil_new];
end
end
  4 个评论
显示 2更早的评论
Daniel Ribeiro Gonçalves
Daniel Ribeiro Gonçalves 2020-1-9
They wont be deleted.. Its only for the matter of performance.. It will lower your program computacional time
Gerrit
Gerrit 2020-1-9
编辑：Gerrit 2020-1-9
Oh okay I understand. But where would I put it?
And I just saw sparse does not work for tables.

请先登录，再进行评论。

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

采纳的回答

Philippe Lebel
Philippe Lebel 2020-1-8
编辑：Philippe Lebel 2020-1-8
In order to help you in a significant way, it would be useful to have a sample data set on which we could run the function.
On the top of my head, i'd remove the disp() function call. Printing info to the console takes a significant amount of time for no real benefit.
I would also pre-alocate the size of
airfoil_update
since you know it's maximum size, you can pre alocate it to that size and keep track of the number of elements you put in it then chop the unused space only once at the end of the script.
  9 个评论
显示 7更早的评论
Gerrit
Gerrit 2020-1-9
Thank you veryyy much Philippe!!
here is the link: https://de.mathworks.com/matlabcentral/answers/499639-how-can-i-speed-up-my-code
Gerrit
Gerrit 2020-1-12
Philippe, do you have any idea for my other problem?

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Airfoil tools 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by