Mercurial > repos > public > sbplib
comparison +blockmatrix/toMatrix.m @ 886:8894e9c49e40 feature/timesteppers
Merge with default for latest changes
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Thu, 15 Nov 2018 16:36:21 -0800 |
| parents | 2501067f2fc7 |
| children |
comparison
equal
deleted
inserted
replaced
| 816:b5e5b195da1e | 886:8894e9c49e40 |
|---|---|
| 1 function A = toMatrix(bm) | |
| 2 if ~blockmatrix.isBlockmatrix(bm) | |
| 3 error('blockmatrix:toMatrix:NotABlockmatrix', 'Input is not a blockmatrix'); | |
| 4 end | |
| 5 | |
| 6 div = blockmatrix.getDivision(bm); | |
| 7 n = div{1}; | |
| 8 m = div{2}; | |
| 9 | |
| 10 N = sum(n); | |
| 11 M = sum(m); | |
| 12 | |
| 13 A = sparse(N,M); | |
| 14 | |
| 15 for i = 1:size(bm,1) | |
| 16 for j = 1:size(bm,2) | |
| 17 if isempty(bm{i,j}) | |
| 18 bm{i,j} = sparse(n(i),m(j)); | |
| 19 end | |
| 20 end | |
| 21 end | |
| 22 A = cell2mat(bm); | |
| 23 end |