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
-~----------~----~----~----~------~----~------~--~---

Reply via email to