Padua2D: Lagrange Interpolation at Padua Points on Bivariate Domains