matlab оптимизация использования factorial

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 

В обоих случаях линия "стойкие фары" соответственно. «global facs» теперь занимает большую часть времени в профилировщике.

Есть ли другой способ сделать это? Могу ли я определить что-то вроде переменной, доступной во всем мире, без ее загрузки?

Answer 1

Избегайте вызывать функции и использовать if..else как можно больше. Векторизация - лучший подход. вы можете сделать факторный матрицу, как это:

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

Предположим, что вы хотите назвать Факториал 5 раз со значением n с являются

ns = [3 0 56 23 456]; 

так что вы можете получить факториал, как показано ниже:

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

, то вы продолжаете вычисления со значениями f_ns.

Если вы не можете предвычисление f_ns вы можете создать встроенную функцию facn внутри функции, которая нуждается в Mutiple вычисления факториала и вызвать facn несколько раз:

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