Mercurial > repos > public > sbplib
view +blockmatrix/multiply.m @ 1249:1f82ad27616d feature/dirac_discr
Close branch feature/dirac_discr
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Wed, 20 Nov 2019 22:24:06 +0000 |
parents | 288138d9586d |
children |
line wrap: on
line source
function C = multiply(A, B) [n_a, m_a] = size(A); [n_b, m_b] = size(B); assert(m_a == n_b, 'Dimensions do not agree.') C = cell(n_a, m_b); for i = 1:n_a for j = 1:m_b C{i,j} = sparse(size(A{i,1},1), size(B{1,j},2)); for k = 1:n_b C{i,j} = C{i,j} + A{i,k}*B{k,j}; end end end end