Mercurial > repos > public > sbplib

diff +blockmatrix/fromMatrix.m @ 338:3974dccff55b feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Made block matrices created from a matrix always have matrix blocks, never empty.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 24 Oct 2016 17:01:05 +0200
parents 40dda96c8c9c
children
line wrap: on
line diff
--- a/+blockmatrix/fromMatrix.m	Fri Oct 21 12:59:07 2016 +0200
+++ b/+blockmatrix/fromMatrix.m	Mon Oct 24 17:01:05 2016 +0200
@@ -13,7 +13,7 @@
         for j = 1:length(d2)
             Asub = A(I:(I + d1(i)-1), J:(J + d2(j)-1));
             if nnz(Asub) == 0
-                bm{i,j} = [];
+                bm{i,j} = sparse(d1(i), d2(j));
             else
                 bm{i,j} = Asub;
             end