Opened 10 years ago

Last modified 10 years ago

#6397 needs_work enhancement

[with patch, needs work] implement general Newton's method root finding for power series rings

Reported by: ncalexan Owned by: malb
Priority: major Milestone: sage-feature
Component: commutative algebra Keywords: power series root newton method
Cc: robertwb, was Merged in:
Authors: Nick Alexander Reviewers:
Report Upstream: Work issues:
Branch: Commit:
Dependencies: Stopgaps:

Description

Extracting a square root of a power series is implemented in power_series_ring_element.pyx. Could we have the more general "improving a root of a polynomial" Newton's method?

My use case is calculating Puiseaux expansions around points of algebraic curves.

Attachments (1)

trac_6397-power-series-newton-ncalexan-1.pyx​ (4.3 KB) - added by ncalexan 10 years ago.

Download all attachments as: .zip

Change History (2)

comment:1 Changed 10 years ago by ncalexan

  • Authors set to Nick Alexander
  • Summary changed from implement general Newton's method root finding for power series rings to [with patch, needs work] implement general Newton's method root finding for power series rings

Here's a stand-alone implementation that needs to be plugged into the hell that is polynomial.roots().

Changed 10 years ago by ncalexan

Note: See TracTickets for help on using tickets.