Mercurial > repos > public > sbplib
view +blockmatrix/fromMatrix.m @ 753:44c46bd6913a feature/d1_staggered
Fix stage rate bug in RK4 discrete data
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Fri, 15 Jun 2018 13:02:29 -0700 |
parents | 3974dccff55b |
children |
line wrap: on
line source
function bm = fromMatrix(A, div) d1 = div{1}; d2 = div{2}; [n, m] = size(A); if n ~= sum(d1) || m ~= sum(d2) error('blockmatrix:fromMatrix:NonMatchingDim','The dimensions in div does not sum to the dimensions in A.'); end bm = cell(length(d1), length(d2)); I = 1; for i = 1:length(d1) J = 1; for j = 1:length(d2) Asub = A(I:(I + d1(i)-1), J:(J + d2(j)-1)); if nnz(Asub) == 0 bm{i,j} = sparse(d1(i), d2(j)); else bm{i,j} = Asub; end J = J + d2(j); end I = I + d1(i); end end