On the approximation order of triangular Shepard interpolation