Mercurial > repos > public > sbplib
view +blockmatrix/toMatrix.m @ 1012:1e437c9e5132 feature/advectionRV
Create residual viscosity package +rv and generalize the ResidualViscosity class
- Generalize residual viscosity, by passing user-defined flux and calculating the time derivative outside of the update.
- Create separate RungekuttaRV specifically using interior RV updates
- Separate the artifical dissipation operator from the scheme AdvectionRV1D so that the same scheme can be reused for creating the diff op used by the ResidualViscosity class
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Wed, 05 Dec 2018 13:44:10 +0100 |
parents | 2501067f2fc7 |
children |
line wrap: on
line source
function A = toMatrix(bm) if ~blockmatrix.isBlockmatrix(bm) error('blockmatrix:toMatrix:NotABlockmatrix', 'Input is not a blockmatrix'); end div = blockmatrix.getDivision(bm); n = div{1}; m = div{2}; N = sum(n); M = sum(m); A = sparse(N,M); for i = 1:size(bm,1) for j = 1:size(bm,2) if isempty(bm{i,j}) bm{i,j} = sparse(n(i),m(j)); end end end A = cell2mat(bm); end