Improve performance of this index creation and accessing code

Hi @eldee I am trying to 3D plot Black Hole simulation data. I am converting Kerr-Schild coordinates to Cartesian coordinates.
Bp is magnetic field ratio(poloidal/toroidal). I need log.(r) in interpolate because r(geometric series) is not uniform so i have to make it uniform.

I could not notice improvement in plot on making changes in ϕ_eval.
a is spin of Black hole. Yes, We need extrapolate as we are interpolating on cx_all which are location where color is assigned and calculating cxe_all. Although extrapolation can be ignored if i tolerate little bit error/approximation in color position and interpolate directly on cx_all.

That i too don’t understand but it gives correct plot as i have matched it with scatter plot and every feature looks same.