Three-dimensional indexing in GPU for numerical approximation of solutions of the Laplace equation