Interpret numpy arrays as quaternionic arrays with numba acceleration