On the fast solution of Toeplitz-block linear systems arising in multivariate approximation theory