Избегание for-loops в матричных операциях с элементами ячеек

Question

Размеры этой проблемы: model.nlf = 4. Каждый {r} Куу или Куу Гамма - это 500 × 500 матриц данных.

Как подавить цикл? моя интуиция идет на использование cellfun с функцией для logdet.

logDetKuu = 0; 
for r=1:model.nlf, 
    if isfield(model, 'gamma') && ~isempty(model.gamma) 
     [model.Kuuinv{r}, model.sqrtKuu{r}] = pdinv(model.KuuGamma{r}); 
     model.logDetKuu{r} = logdet(model.KuuGamma{r}, model.sqrtKuu{r}); 
    else 
     [model.Kuuinv{r}, model.sqrtKuu{r}] = pdinv(model.Kuu{r}); 
     model.logDetKuu{r} = logdet(model.Kuu{r}, model.sqrtKuu{r}); 
    end 
    logDetKuu = logDetKuu + model.logDetKuu{r}; 
end 

Благодарен за указатели. Спасибо

Следующий вопрос: Можно ли векторизовать следующие типы петель на элементах ячейки? nlf = 4; nout = 16; каждый KuuinvKuy {1,1} является 150x650

for k =1: model.nout, 
    for r =1: model.nlf, 
     model.KuuinvKuy{r,k} = model.Kuuinv{r}*model.Kyu{k,r}'; 
    end 
end 

Answer 1

Если все матрицы являются настолько большими, и вы выполняете свои только 4 для цикла раз, то нет никаких причин, чтобы удалить для-петли, так как это будет не приводят к ускорению. Мое единственное наблюдение заключается в том, что состояние if, по-видимому, не зависит от цикла, поэтому чистнее перемещать это значение if перед циклом. Что-то вроде этого:

if isfield(model, 'gamma') && ~isempty(model.gamma) 
    myKuu = model.KuuGamma; 
else 
    myKuu = model.Kuu; 
end 

logDetKuu = 0;  
for r=1:model.nlf, 
    [model.Kuuinv{r}, model.sqrtKuu{r}] = pdinv(myKuu{r}); 
    model.logDetKuu{r} = logdet(myKuu{r}, model.sqrtKuu{r}); 
    logDetKuu = logDetKuu + model.logDetKuu{r}; 
end