Mercurial > repos > public > sbplib

comparison +blockmatrix/multiply.m @ 427:a613960a157b feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
merged with feature/beams
author Ylva Rydin <ylva.rydin@telia.com>
date Thu, 26 Jan 2017 15:59:25 +0100
parents 288138d9586d
children
comparison
equal deleted inserted replaced
426:29944ea7674b 427:a613960a157b
1 function C = multiply(A, B)
2 [n_a, m_a] = size(A);
3 [n_b, m_b] = size(B);
4
5 assert(m_a == n_b, 'Dimensions do not agree.')
6
7 C = cell(n_a, m_b);
8
9 for i = 1:n_a
10 for j = 1:m_b
11 C{i,j} = sparse(size(A{i,1},1), size(B{1,j},2));
12 for k = 1:n_b
13 C{i,j} = C{i,j} + A{i,k}*B{k,j};
14 end
15 end
16 end
17 end