Opened 11 years ago
Last modified 5 weeks ago
#4618 new enhancement
Puiseux series
Reported by:  ljpk  Owned by:  AlexGhitza 

Priority:  minor  Milestone:  sage8.2 
Component:  algebra  Keywords:  Puiseux 
Cc:  mmezzarobba, vdelecroix, dkrenn  Merged in:  
Authors:  Reviewers:  
Report Upstream:  N/A  Work issues:  
Branch:  public/4618 (Commits)  Commit:  a922e2acef3c871a23600de9d299d210673fcee4 
Dependencies:  #24420, #24431  Stopgaps: 
Description (last modified by )
We provide an implementation of Puiseux series, that is power series in x^(1/n)
where n
is an arbitrary integer.
When the base ring is an algebraically closed field, this is an algebraically closed field. In other words, any polynomial in QQ[X,Y]
has a solution in Y
as a Puiseux series in X
over QQbar
.
Change History (33)
comment:1 Changed 11 years ago by
comment:2 Changed 10 years ago by
 Type changed from defect to enhancement
comment:3 Changed 9 years ago by
 Component changed from linear algebra to algebra
 Owner changed from was to AlexGhitza
 Report Upstream set to N/A
comment:4 Changed 8 years ago by
See also #9555.
comment:5 Changed 6 years ago by
 Milestone changed from sage5.11 to sage5.12
comment:6 Changed 5 years ago by
 Milestone changed from sage6.1 to sage6.2
comment:7 Changed 5 years ago by
 Keywords Puiseux added
comment:8 Changed 5 years ago by
 Description modified (diff)
comment:9 Changed 5 years ago by
 Description modified (diff)
comment:10 Changed 5 years ago by
 Milestone changed from sage6.2 to sage6.3
comment:11 Changed 5 years ago by
 Milestone changed from sage6.3 to sage6.4
comment:12 Changed 4 years ago by
See also #9289.
comment:13 Changed 4 years ago by
 Milestone changed from sage6.4 to sagewishlist
comment:14 Changed 3 years ago by
some code is available here:
https://github.com/abelfunctions/abelfunctions/tree/master/abelfunctions
comment:15 Changed 3 years ago by
As mentioned in an email correspondence with chapoton I believe my implementation of PuiseuxSeriesRing
and PuiseuxSeriesRingElement
to be pretty complete in regards to fitting into Sage's coercion model. Hopefully it will, at the very least, be a good starting point.
Thanks to chapoton for offering to plug this code into Sage!
On a side note, I also have some code for computing Puiseux series representations of places on a plane algebraic curve. The entry function is abelfunctions.puiseux_series.puiseux_series
. I'm sure someone with better Sage skills than myself can improve the algorithm and plug that in as well.
comment:16 Changed 3 years ago by
 Branch set to public/4618
 Commit set to 7c6e9fbd65e4b19eed62d8cc8d8694659b49e996
just made a git branch, not tested at all
New commits:
7c6e9fb  trac #4618 creating a branch from Chris Swierczewski code

comment:17 Changed 3 years ago by
 Commit changed from 7c6e9fbd65e4b19eed62d8cc8d8694659b49e996 to a828db5df2268a20bea62b83033c2325f7cfeaca
Branch pushed to git repo; I updated commit sha1. New commits:
a828db5  trac #4618 details, not working yet

comment:18 Changed 3 years ago by
 Commit changed from a828db5df2268a20bea62b83033c2325f7cfeaca to 0e38712d192a35f9cf35d1c92e053b412ad8e940
Branch pushed to git repo; I updated commit sha1. New commits:
0e38712  trac #4618 now working

comment:19 Changed 3 years ago by
 Commit changed from 0e38712d192a35f9cf35d1c92e053b412ad8e940 to 8fe475bada98e8b7a2912693525f56d036dcaa10
Branch pushed to git repo; I updated commit sha1. New commits:
8fe475b  trac #4618, now working a little bit more, and added some doc

comment:20 Changed 3 years ago by
 Commit changed from 8fe475bada98e8b7a2912693525f56d036dcaa10 to 2e05bc8e741b8f943f4fedff664e24a1f8ba9bc7
comment:21 Changed 3 years ago by
 Cc mmezzarobba added
comment:22 Changed 3 years ago by
 Commit changed from 2e05bc8e741b8f943f4fedff664e24a1f8ba9bc7 to 9af42db0e6db06161e484f318e456e664cbf9ebc
Branch pushed to git repo; I updated commit sha1. New commits:
9af42db  trac #4618 some more doc and tests

comment:23 Changed 3 years ago by
 Commit changed from 9af42db0e6db06161e484f318e456e664cbf9ebc to 98afe229b00ae434a5a2bd4849bf94081ff39bca
comment:24 Changed 2 years ago by
This is a nice addition. Is it ready for review?
A few comments:
 Is it planned to extend it to the multivariate case?
 It would be nice to add a more detailed mathematical explanation about what a Puiseux series is in the docstring.
comment:25 Changed 2 years ago by
It is probably not ready, there seems to be a bug lurking around.
No plan to implement the multivariate case.
and the coverage is not 100%
comment:26 Changed 17 months ago by
 Cc vdelecroix added
comment:27 Changed 17 months ago by
This is needed to allow extend=True
for the nth_root
method in #10720 (see also this sagedevel thread).
comment:28 Changed 17 months ago by
An innocent operation like the following will multiply the memory footprint by 24 with no change in information
sage: p = prod(1  q^n + O(q^100) for n in range(1,100)) # fine: 100 * coeff size sage: q^(1/24) * p # bad: 100 * 24 * coeff size
The above example is the qseries expansion of the Dedekind eta function. There might be a more clever data structure to use as this kind of series is frequently encountered when dealing with modular forms.
comment:29 Changed 17 months ago by
 Dependencies set to #24420, #24431
comment:30 Changed 17 months ago by
 Description modified (diff)
 Milestone changed from sagewishlist to sage8.2
 Summary changed from Request for Puiseux series in SAGE to Puiseux series
comment:31 Changed 17 months ago by
 Cc dkrenn added
comment:32 Changed 5 weeks ago by
 Commit changed from 98afe229b00ae434a5a2bd4849bf94081ff39bca to a922e2acef3c871a23600de9d299d210673fcee4
comment:33 Changed 5 weeks ago by
I do the following changes:
 Adaptation to current version 8.8.beta3 in order to have the code compile again.
 I replace the method
_cmp_
ofPuiseuxSeries
by_richcmp_
and let it completely rely on the corresponding method of classLaurentSeries
. The reason for this is not only modernisation, but also that the former method doesn't work correctly (for example comparison with zero returns wrong results).  I add a specification of the representative of the Puiseux series inside the constructor of
PuiseuxSeries
in such a way that the ramification index is minimized. This improves the methodslaurent_series
andpower_series
(the examples given there would not work elsewise).  I add more doctests.
There is still work to do, for example:
 There are several workarounds on other bugs in Sage. In general I think, these workarounds should be removed and the corresponding bugs should be treated in separate tickets. Examples:
 In the methods
_repr_
,exponents
andcoefficients
there is the following comment: NOTE:self.__l.coefficients()
is bugged when the coefficients are in QQbar but coerced into SR .... I have no idea how far it makes sense to consider Puiseux series (Laurent series, polynomial rings, ...) over theSymbolicRing
. But since these constructions are possible, they should work (or should be blocked). There are simple examples where this is not the case:sage: PS = PolynomialRing(SR,'x') sage: P = PolynomialRing(QQ,'x') sage: q = P((1,1,5)); q 5*x^2 + x + 1 sage: p = PS(q) sage: p.coefficients() [5*x^2 + x + 1] sage: p in SR True
Is this a known bug? Concerning the methods in question, I think that they should rely more directly on the according methods ofLaurentSeries
.  In the method
add_bigoh
the following error is caught:sage: L.<x> = LaurentSeriesRing(QQ) sage: q = x^2 + x^3 sage: q.add_bigoh(1) Traceback (most recent call last): ... ValueError: prec (= 3) must be nonnegative
This should be fixed inLaurentSeries
.
 In the methods
 I think we should clarify the following behavior of the method
add_bigoh
:sage: R.<x> = PuiseuxSeriesRing(ZZ) sage: p = x**(1/3) + 2*x**(1/5) sage: p.add_bigoh(1/2) x^(1/3) + 2*x^(1/5) + O(x^(7/15))
is this acceptable?
 The method
_repr_
needs work:sage: R.<x> = PuiseuxSeriesRing(Zp(5)) sage: x**(1/2) + 5 * x^(1/3) 5 + O(5^21)*x^(1/3) + (1 + O(5^20))*x^(1/2)
 Further doctests are needed.
 Integration into documentation.
I will continue to work on this ticket according to feedback!
Hi ljpk,
requests like this shouls always go to [sagedevel] too since people generally do not look for things to do in trac. Sending the same request to [sagedevel] will make people aware of the problem and might get some design discussion going. Obviously if you are working on code this is different, but if you expect someone else to do the dirty work a little advertisement cannot hurt :)
Cheers,
Michael