オクターブが遅い。提案

Question

Octave 4.0.0とMATLAB 2014の両方で次のコードを実行しました。時差は愚かです。すなわち、2桁以上です。 Windowsラップトップで実行。 Octaveの計算速度を向上させるために何ができるのですか？

startTime = cputime; 
iter = 1; % iter is the current iteration of the loop 
itSum = 0; % itSum is the sum of the iterations 
stopCrit = sqrt(275); % stopCrit is the stopping criteria for the while loop 
while itSum < stopCrit 
    itSum = itSum + 1/iter; 
    iter = iter + 1; 
    if iter > 1e7, break, end 
end 
iter-1 
totTime = cputime - startTime 

オクターブ：totTime〜112

MATLAB：totTime < 0.4

Answer 1

それはあなたのコードで結果を計算するためにループでの反復がかかります。コードをベクトル化することで、スピードアップに役立ちます。私の次のコードはあなたがやったことを正確に行いますが、計算をかなりベクトル化します。それが役立つかどうかを見てください。

startTime = cputime; 
iter = 1; % iter is the current iteration of the loop 
itSum = 0; % itSum is the sum of the iterations 
stopCrit = sqrt(275); % stopCrit is the stopping criteria for the while loop 
step=1000; 
while(itSum < stopCrit && iter <= 1e7) 
    itSum=itSum+sum(1./(iter:iter+step)); 
    iter = iter + step+ 1; 
end 
iter=iter-step-1; 
itSum=sum(1./(1:iter)); 
for i=(iter+1):(iter+step) 
    itSum=itSum+1/i; 
    if(itSum+1/i>stopCrit) 
     iter=i-1; 
     break; 
    end 
end 
totTime = cputime - startTime 

実行時間は、上記のコードを使用して約0.6秒です。ループが停止したとき、あなたは正確に気にしない場合は、次のコードでは、さらに高速である：

startTime = cputime; 
iter = 1; % iter is the current iteration of the loop 
itSum = 0; % itSum is the sum of the iterations 
stopCrit = sqrt(275); % stopCrit is the stopping criteria for the while loop 
step=1000; 
while(itSum < stopCrit && iter <= 1e7) 
    itSum=itSum+sum(1./(iter:iter+step)); 
    iter = iter + step+ 1; 
end 
iter=iter-step-1; 
totTime = cputime - startTime 

私のランタイムでは、後者の場合にのみ、約0.35秒です。

Answer 2

また、試すことができます。

itSum = sum(1./(1:exp(stopCrit))); 
    %start the iteration 
    iter = exp(stopCrit-((stopCrit-itSum)/abs(stopCrit-itSum))*(stopCrit-itSum)); 
    itSum = sum(1./(1:iter)) 

このMethodeのではあなたが唯一の1または2反復を持つことになります。もちろん、配列全体を合計します。