% Usage: A = matmul_bycolumn(B, C)
%
% Computes the matrix product A = B * C, but instead of
% just calling B * C, it multiplies B by each column of
% C in sequence.
function A = matmul_bycolumn(B, C)
A = zeros(size(B,1), size(C,2));
for j = 1:size(C,2)
A(:,j) = B * C(:,j);
end