Structure-Preserving Matrix Methods for Computations on Univariate and Bivariate Bernstein Polynomials