python中bigfloat变量的trigamma和digamma函数

我需要计算python中bigfloat类型的变量的

scipy.special

trigamma和digamma函数，但我得到以下错误消息：

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'

为了精确起见，我需要将变量保持为bigfloat。 有人知道我如何计算这些变量的digamma和trigamma函数吗？

---

谢谢。

和中提供了digamma的多精度实现。我不知道trigamma的任何实现

---

免责声明：我坚持认为

gmpy2

trigamma是polygamma的特例