There seems to be a bug in Complex with domains Float or DoubleFloat:
Working:
(68) -> (-1 :: Complex(MachineFloat))^(3/2)
(68) - %i
Type:
Expression(Complex(MachineFloat))
But for Float it's not working as expected:
(65) -> (-1 :: Complex(Float))^(3/2)
(65) %i
Type:
Complex(Float)
(this is fine:)
(64) -> (-1 :: Complex(Float))^(1/2)
(64) %i
Type:
Complex(Float)
Neither working for DoubleFloat:
(69) -> (-1 :: Complex(DoubleFloat))^(3/2)
(69) %i
Type:
Complex(DoubleFloat)
I looked at the Complex implementation and there is some special case
"if R is Float or R is DoubleFloat then" but I don't see how this would
lead to wrong results.
A temporary workaround for me is to convert the exponent to Float first.
Thanks and best wishes,
Tobias
--
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 view this discussion on the web visit
https://groups.google.com/d/msgid/fricas-devel/62c39e1b-cb07-4b41-8f48-6086fa653f43n%40googlegroups.com.
