Scipy 1.14.0 changed the upcasting rule of several special functions when passing integers

[ricardoV94]

Description

Functions like scipy.special.psi now always upcast to float64, even if the input is a low precision integer like int8. We need to handle these types, but:

1. float64 may be a bit too wasteful
2. We can't have the dtype depend on the scipy version, and lower pinning to 1.14.0 may be too restrictive?

[ricardoV94]

I posed a vague question on the scipy discourse: https://discuss.scientific-python.org/t/several-special-functions-started-upcasting-low-precision-integer-inputs-to-float64-in-v1-14-0/1266