I need to compute the scipy.special
trigamma and digamma functions of a variable with type bigfloat in python but I get the following error message:
TypeError: ufunc 'psi' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule 'safe'
I need to keep my variables as bigfloat for precision purposes. Does any body know how I can compute digamma and trigamma functions for such variables? Thanks.