clear all;

ni = 2;
no = 3;
nr = 2;
np = 2;

x = ones(np, ni);

w = reshape(1:nr*np, nr, 1, np)
c = bsxfun(@times, permute(1:nr, [1 3 2]), ones(no, ni))

f = zeros(no, ni, np);
for ii=1:np
 for jj=1:nr
   f(:,:,ii) = f(:,:,ii) + w(jj, 1, ii)*c(:,:,jj);
  end
end
f