Speeding up the Scipy Bessel function when applied over large arrays

Question

I am working on a numerical problem that requires me to calculate the value of the Bessel function of the first kind of order 1/4 over a very large 4-dimensional array. To g