10/07/2014, 11:00 AM

(10/07/2014, 01:50 AM)mike3 Wrote: Addendum: You can replace Wapprox with this:Hi Mike -

Code:`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).

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

Gottfried Helms, Kassel