Is there a way to reduce compilation time for this segment?

Question

Souheil Araoud 2019-9-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479415-is-there-a-way-to-reduce-compilation-time-for-this-segment

评论： Souheil Araoud 2019-9-8

tic
    iscorr = 1;
    myunc = zeros(1,length(Kurzschluss.Plateaus_I));
    myfunction = '(X2/X1)*(sqrt(1 - X3^2/(X2^2 * X1^2)))';  
    for i=1:length(mystruct.Plateaus_I)
         x1 = mystruct.Plateaus_I{1,i};
         x2 = mystruct.Plateaus_U{1,i};
         x3 = mystruct.Plateaus_P{1,i};
         myunc(i) =  calculateEndUncertainty_3(x1,x2,x3,iscorr,myfunction);
    end 
    
    toc

4 个评论
显示 2更早的评论隐藏 2更早的评论

Souheil Araoud 2019-9-8

在 MATLAB Online 中打开

function [unc] = CalculateEndUncertainty_3(x1,x2,x3,iscorr,myfunction) 
N = 3; 
if iscorr
Enable_correlation  = 1;
else
Enable_correlation = 0;
end
StrModelFunction = myfunction;
StrInput = '';
StrMean = '';
for c1 = 1:N
    StrInput = [StrInput 'X' num2str(c1) ','];    
    StrMean = [StrMean 'mean(x' num2str(c1) '),']; 
end
StrInput(end) = [];
StrMean(end) = [];
eval(['syms f(' StrInput ')'])
eval(['f(' StrInput ') ='  StrModelFunction ';'])
u2 = 0;
for c1 = 1:N
    
    temp = diff(f,eval(['X' num2str(c1)]));
    df_dx = double(eval(['temp(' StrMean ')']));
    
    u2 = u2 + df_dx^2*std(eval(['x' num2str(c1)]))^2/length(x1) ;
    
end
 if Enable_correlation
    for c1 = 1:N-1
        for c2 = c1+1:N    
            temp_c1 = diff(f,eval(['X' num2str(c1)]));
            df_dx_c1 = double(eval(['temp_c1(' StrMean ')']));
            temp_c2 = diff(f,eval(['X' num2str(c2)]));
            df_dx_c2 = double(eval(['temp_c2(' StrMean ')']));
            temp = cov(eval(['x' num2str(c1)]),eval(['x' num2str(c2)]));
            u2 = u2 + 2*df_dx_c1*df_dx_c2*temp(2)/length(x1); 
        end
    end
end
unc = sqrt(abs(u2));
res = double(eval(['f(' StrMean ')']));
end

darova 2019-9-8

You have a lot of eval() that is why script is slow

The function below depends on f(), temp(), temp_c1(), temp_c2()

Maybe you can attach whole your project? Any description?

Souheil Araoud 2019-9-8

The whole project is way too long. This function is supposed to take 3 inputs of a function as vectors (in my case containing about 5 lightly diffefrent values) and calculate the uncertainty of the result based on the mean value of each input vector etc.. The function has been applied on 3 unidimensional cells named Plateaus_I, Plateaus_U and Plateaus_P containing a vector in each cell.

请先登录，再进行评论。

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

Answer 1

darova 2019-9-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479415-is-there-a-way-to-reduce-compilation-time-for-this-segment#answer_390898

CalculateEndUncertainty_3.m

See if it works

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Souheil Araoud 2019-9-8

Great! Thanks a lot

请先登录，再进行评论。

Answer 2

Jackson Burns 2019-9-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479415-is-there-a-way-to-reduce-compilation-time-for-this-segment#answer_390882

在 MATLAB Online 中打开

Hi Souheil!

If you have access to the parallel computing toolbox, you can improve execution time with a parfor loop.

tic
iscorr = 1;
myunc = zeros(1,length(Kurzschluss.Plateaus_I));
myfunction = '(X2/X1)*(sqrt(1 - X3^2/(X2^2 * X1^2)))';  
parfor i=1:length(mystruct.Plateaus_I)
     x1 = mystruct.Plateaus_I{1,i};
     x2 = mystruct.Plateaus_U{1,i};
     x3 = mystruct.Plateaus_P{1,i};
     myunc(i) =  calculateEndUncertainty_3(x1,x2,x3,iscorr,myfunction);
end 
toc

Other optimizations will likely come from improving your functions, such as calculateEndUncertainty_3. This page from MathWorks has some great advice on how to do this. I would reccomend looking at vectorization in particular, as MATLAB is extremely powerful when properly vectorized.

Hope this helps, good luck!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Is there a way to reduce compilation time for this segment?

4 个评论
显示 2更早的评论隐藏 2更早的评论

采纳的回答

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

更多回答（1 个）

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Is there a way to reduce compilation time for this segment?

4 个评论 显示 2更早的评论隐藏 2更早的评论

采纳的回答

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

更多回答（1 个）

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

4 个评论
显示 2更早的评论隐藏 2更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

0 个评论
显示 -2更早的评论隐藏 -2更早的评论