Decomposition and reconstruction algorithms for bivariate spline wavelets on the unit square