Changeset 5428:e4e8a90e85d2 for sage/rings/power_series_ring_element.pyx

Timestamp:

06/07/07 22:53:49 (6 years ago)

Author:

Robert Bradshaw <robertwb@…>

Branch:

default

Message:

optimize sqrt some more

File:

• sage/rings/power_series_ring_element.pyx (modified) (2 diffs)

sage/rings/power_series_ring_element.pyx

@@ -358 +358 @@
         """
         return self._prec
-
+        
     def _repr_(self):
         """
@@ -1096 +1096 @@
             a = a.change_ring(R)
         half = ~R(2)
-            
-        for i in range (ceil(log(prec, 2))):
+        
+        s = a.parent()([s])
+        for cur_prec in sage.misc.misc.newton_method_sizes(prec)[1:]:
+            (<PowerSeries>s)._prec = cur_prec
             s = half * (s + a/s)
-            
-        ans = P.gen(0)**(val/2) * s
+        
+        ans = s
+        if val != 0:
+            ans *= P.gen(0)**(val/2)
 
         if all: