Abstract: We show that stable local bases for the spaces of polynomial splines on a triangulation of a bivariate polygonal domain can be efficiently computed by using either singular value decomposition or pivoted QR-decompositon of certain small matrices of nodal smoothness conditions.
Preprint version available: pdf