wolffd@0: function p = isscalarBNT(v)
wolffd@0: % ISSCALAR Returns 1 if all dimensions have size 1.
wolffd@0: % p = isscalar(v)
wolffd@0: 
wolffd@0: p = (prod(size(v))==1);