Mercurial > repos > public > sbplib
view +blockmatrix/isDivision.m @ 1320:cd0bfcaef0ba feature/poroelastic
Add some interface forcing penalties in ViscoElastic2d
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Tue, 28 Jul 2020 21:59:41 -0700 |
parents | a5f1b0267dba |
children |
line wrap: on
line source
function b = isDivision(div) % Make sure it is a cellarray if ~iscell(div) b = false; return end % Make sure it has the right shape if numel(div) ~= 2 b = false; return end if ~isDivisionVector(div{1}) || ~isDivisionVector(div{2}) b = false; return end b = true; end function b = isDivisionVector(v) if isempty(v) b = true; return end if any(v < 0) b = false; return end b = true; end