Mercurial > repos > public > sbplib
view +blockmatrix/getDivision.m @ 1213:43f1cd11e8e8 feature/poroelastic
Add physical normals to AnisotropicCurvilinear
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Mon, 14 Oct 2019 13:54:50 -0700 |
parents | a5f1b0267dba |
children |
line wrap: on
line source
function div = getDivision(bm) if ~blockmatrix.isBlockmatrix(bm) error('blockmatrix:getDivision:NotABlockmatrix', 'Input is not a blockmatrix'); end if isempty(bm) div = {[],[]}; return end div = {row_height(bm),col_width(bm)}; end function m = col_width(C) m = zeros(1,size(C,2)); for j = 1:size(C,2) for i = 1:size(C,1) if isNullMatrix(C{i,j}) continue end m(j) = size(C{i,j},2); end end end function n = row_height(C) n = zeros(1,size(C,1)); for i = 1:size(C,1) for j = 1:size(C,2) if isNullMatrix(C{i,j}) continue end n(i) = size(C{i,j},1); end end end function b = isNullMatrix(A) [n, m] = size(A); b = n == 0 && m == 0; end