Computing bivariate splines in scattered data fitting and the finite-element method