Efficient computation of partition of unity interpolants through a block-based searching technique