Find all fixed points of exp[b]
(10/07/2014, 01:50 AM)mike3 Wrote: Addendum: You can replace Wapprox with this:

Wapprox(k, z) = {
   local(brl = brlog(k, z));
   local(lbrl = log(brl));

   return(brl - lbrl + sum(kk=0,5,sum(m=1,5,c(kk, m)*lbrl^m * brl^(-kk-m))));

to get even more speed (saves some repeated log calculations).
Hi Mike -

very well, thank you! I've also made the coefficients c() constants in a small matrix, because your routine uses them only to small arguments but repeatedly, also made Pi2I a constant. I've got 120 frixpoints to base 4 with 200 digits precision in 656 ms on a 2 Ghz machine. Great!

Gottfried Helms, Kassel

Possibly Related Threads…
Thread Author Replies Views Last Post
  Kneser-iteration on n-periodic-points (base say \sqrt(2)) Gottfried 11 9,487 05/05/2021, 04:53 AM
Last Post: Gottfried
  On n-periodic points of the exp() - A discussion with pictures and methods Gottfried 1 4,152 06/10/2020, 09:34 AM
Last Post: Gottfried
  fixed point formula sheldonison 6 21,848 05/23/2015, 04:32 AM
Last Post: mike3
  Attempt to find a limit point but each step needs doubling the precision... Gottfried 15 39,134 11/09/2014, 10:25 PM
Last Post: tommy1729
  sexp(strip) is winding around the fixed points Kouznetsov 8 24,209 06/29/2009, 10:05 AM
Last Post: bo198214
  An error estimate for fixed point computation of b^x bo198214 0 4,835 05/31/2008, 04:11 PM
Last Post: bo198214

Users browsing this thread: 1 Guest(s)