general/basic/lincomb_sequence2.m

function res = lincomb_sequence2(seq,sigma1,sigma2)
%function res = lincomb_sequence2(seq,sigma1,sigma2)
%
% function performing a linear combination of the elements in the
% 2d cell array seq with coefficients in sigma1 and sigma2 result is a
% vector/matrix the same size as the entries of seq.
% size of seq is length(sigma1) x length(sigma2)
% if some sigma = 0, the component is not touched, i.e. the component
% may also be an empty matrix. 

% Bernard Haasdonk 15.5.2007

  Q1 = length(sigma1);
  Q2 = length(sigma2);
  res = zeros(size(seq{1,1}));
  for q1=1:Q1
    if sigma1(q1)~=0      
      for q2=1:Q2
        if sigma2(q2)~=0      
          res = res + sigma1(q1)*sigma2(q2)*seq{q1,q2};
        end;
      end;
    end;
  end;