I have numpy arrays of coordinates X, Y, Z, at which I want to evaluate the obtained solution u.
I know that one thing I is probing u(x, y, z) point by point. Sadly it takes more time than calculation of the solution (5 minutes and 3 seconds respectively). Is there any more efficient approach?