SymPy计算特定无穷级数返回错误结果，咨询其计算正确性边界及异常原因

我最近在用SymPy做数值计算，尤其是无穷求和的场景，碰到了一个很费解的问题。运行下面的代码后，返回结果居然是0：

import sympy as sp
from sympy.abc import k

sp.summation(1 / ((k + 1) ** 4 + 1), (k, -sp.oo, sp.oo))

但这个级数里的每一项都是正数，结果怎么可能是0呢？根据SymPy的文档说明，如果求和无法计算，应该会直接返回原式，但这次它却给出了一个明显错误的结果。这让我不禁疑惑：SymPy在哪些情况下无法保证计算的正确性？

补充说明：这个级数实际上是可计算的（比如通过留数定理就能求解，不需要用到特殊函数），收敛结果是一个正数。我本来接受SymPy可能算不出这个结果，但至少它应该承认自己无法计算，而不是返回一个错误的答案。

备注：内容来源于stack exchange，提问作者Ilya Kuleshov