scipy.special.orthogonal.p_roots with cache