Non-linear Multi-block Partial Least Squares via Univariate and Bivariate B-splines