oldk1331 wrote:
>
> count(zero?, [real acosh((x*0.001)::COMPLEX Float) for x in 0..1000])
> returns 314, after sqrt$Float patch, that number becomes 622.
> After "acosh x == log(x + sqrt(x^2 - 1))" change, that number
> becomes 761.
>
> I think "acosh x == log(x + sqrt(x^2 - 1))" is better than
> "acosh x == 2*log(sqrt((x+1)/2::%) + sqrt((x-1)/2::%))",
> it's more efficient and more accurate.
There is question of branches: AFAIK the later formula
gives correct branch. The one you propose is likely
to lead to wrong branches. The factor 2 in the later
formula means that they produce different sets of values...
--
Waldek Hebisch
--
You received this message because you are subscribed to the Google Groups
"FriCAS - computer algebra system" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at https://groups.google.com/group/fricas-devel.
For more options, visit https://groups.google.com/d/optout.
