Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

generalized_series_solutions() fails with some recurrences #28

Open
mezzarobba opened this issue Aug 27, 2021 · 0 comments
Open

generalized_series_solutions() fails with some recurrences #28

mezzarobba opened this issue Aug 27, 2021 · 0 comments
Labels
bug Something isn't working

Comments

@mezzarobba
Copy link
Collaborator

(Might try fixing it myself someday, but not now.)

sage: Pol.<z> = QQ[]
sage: Rec.<Sz> = OreAlgebra(Pol)
sage: rec = (-25165779*z^3 - 1031796939*z^2 - 14092836240*z - 64122404892)*Sz^13 + (1702884379*z^3 + 652213105
....: 75*z^2 + 835034101592*z + 3573825830162)*Sz^12 + (-22196217078*z^3 - 871105051492*z^2 - 11362215001012*z
....:  - 49294576910526)*Sz^11 + (108363844374*z^3 + 4401360529612*z^2 + 59150888473364*z + 263562562419066)*S
....: z^10 + (-197744302789*z^3 - 8602728613511*z^2 - 121160273897408*z - 556950616254958)*Sz^9 + (20375892397
....: *z^3 + 3501356775235*z^2 + 74072014386184*z + 418316114610808)*Sz^8 + (350408306796*z^3 + 9333484115520*
....: z^2 + 72659777977448*z + 140897660149620)*Sz^7 + (-350408306796*z^3 - 11691014292240*z^2 - 1198103815118
....: 48*z - 382170707559340)*Sz^6 + (-20375892397*z^3 + 2278803231415*z^2 + 41531185746816*z + 17441139780512
....: 8)*Sz^5 + (197744302789*z^3 + 3261929553829*z^2 + 14344292703768*z + 7117838600802)*Sz^4 + (-10836384437
....: 4*z^3 - 2100470132828*z^2 - 13133080537684*z - 25821750195414)*Sz^3 + (22196217078*z^3 + 460667973188*z^
....: 2 + 3153473434932*z + 7077439136914)*Sz^2 + (-1702884379*z^3 - 36951752165*z^2 - 269642933392*z - 661407
....: 003678)*Sz + 25165779*z^3 + 478149801*z^2 + 3019893480*z + 6341776308
sage: rec.generalized_series_solutions()
---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
<ipython-input-5-1e6f1da48059> in <module>
----> 1 rec.generalized_series_solutions()

~/docs/code/sage/ore_algebra/ore_algebra/src/ore_algebra/ore_operator_1_1.py in generalized_series_solutions(self, n, dominant_only, real_only, infolevel)
   4486                             while killer.is_zero():
   4487                                 dl += 1
-> 4488                                 adjustment = operator_applied_to_term(k, l + dl)
   4489                                 killer = adjustment[l + dl][ram*prec - degdrop - k]
   4490                             # update solution

~/docs/code/sage/ore_algebra/ore_algebra/src/ore_algebra/ore_operator_1_1.py in operator_applied_to_term(k, l)
   4458                             if i != l and j == 0: # [log(n)^i] log(n)^l
   4459                                 continue
-> 4460                             contrib += ((coeffs[j]*log_tails[j][l - i]).shift(-ram*prec)* \
   4461                                         (poly_tails[j][k]*spoly_tails[j]).shift(-ram*prec)).shift(-ram*prec - k)
   4462                         out.append(_binomial(l, i)*contrib)

IndexError: list index out of range
@mezzarobba mezzarobba added the bug Something isn't working label Oct 17, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

1 participant