Simpler example:
(7) -> op3c(f : POLY FRAC INT) : POLY FRAC INT == 1/2*z1^2*D(f, z1, 2)
Function declaration op3c : Polynomial(Fraction(Integer)) ->
Polynomial(Fraction(Integer)) has been added to workspace.
Type: Void
(8) -> op3(f : POLY INT) : POLY INT == 1/2*z1^2*D(f, z1, 2)
Function declaration op3 : Polynomial(Integer) -> Polynomial(Integer
) has been added to workspace.
Type: Void
(9) -> pol16g := 17121201828914589238*z1^15
15
(9) 17121201828914589238z1
Type: Polynomial(Integer)
(10) -> op3c(pol16g) - op3(pol16g)
Compiling function op3c with type Polynomial(Fraction(Integer)) ->
Polynomial(Fraction(Integer))
Compiling function op3 with type Polynomial(Integer) -> Polynomial(
Integer)
15
(10) - 2z1
Type: Polynomial(Fraction(Integer))
Tracing suggest that op3 is using floating point arithmetic, that
would explain error.
--
Waldek Hebisch
[email protected]
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"FriCAS - computer algebra system" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/fricas-devel?hl=en
-~----------~----~----~----~------~----~------~--~---