matlabで階乗の使用を最適化

Question

私は同じ階乗を何度も何度も計算するコードを最適化しています。呼び出しはいくつかの異なる関数からのものですので、関数の引数としてあらかじめ計算された階乗を配列に渡すことを避けたかったのですが、これが最も速い方法だと思います。

グローバル変数を使用する方が効果的ですが、グローバル変数の読み込みがファクタルアルの計算とほぼ同じくらい長くかかるので、それほど高速ではありません。

function [ facn ] = cfactorial(n) 

global facs 

if n > 170 
    facn = Inf; 
elseif n == 0 
    facn = 1; 
else 
    facn = facs(n); 
end 

私はまた永続的に設定しようとしましたが、あまり速くはありません。

function [ facn ] = cfactorial(n) 

persistent facs 

if isempty(facs) 
    load('facs.mat') 
    facs = faccs; 
end 

if n > 170 
    facn = Inf; 
elseif n == 0 
    facn = 1; 
else 
    facn = facs(n); 
end 

どちらの場合も、 "persistent facs"という行が表示されます。 「グローバルなfacs」はプロファイラでほとんどの時間を費やします。

これを行う別の方法はありますか？変数をロードせずにグローバルに利用できるようなものを定義できますか？

Answer 1

可能な限り関数を呼び出したり、if..elseを使用したりしないでください。ベクトル化はより良いアプローチです。あなたはこのよう階乗行列を作ることができます ：

f_ns = factorials(min(ns, 171) + 1) ; 

を：

factorials = [1 cumprod(1:170) Inf]; 

あなたはn sの値と階乗5回呼び出す必要があるとしますので、怒鳴るよう階乗を得ることができます

ns = [3 0 56 23 456]; 

です

次に、f_nsの値で計算を続けます。

あなたがf_nsを事前に計算することができない場合は、階乗のみとめ計算を必要とする関数内でインライン関数facnを作成することとfacn複数回呼び出す：あなたはすべての値の配列を事前に計算いけないのはなぜ

function myfunction 
    factorials = [1 cumprod(1:170) Inf]; 
    facn =@(n) factorials(min(n, 171) + 1) ; 
    facn(3); 
    facn(0); 
    facn(456); 
end