Scratch this patch. I see what Waldek mean in the other thread:
> There is also another issue: if some coefficient mathematically
> should be 0 then it is good to make it zero without calculation.
> This helps if we want to reuse routine for say floating point
> computations. AFAICS your 'pseudoDivide' does not have this
> property. Version from NSUP uses 'fmecg' which should ensure
> this.
y :DoubleFloat := (0.1)^200
(y*y=0)@Boolean -- returns true
So we have to check c2 in loop.
But doesn't that causes trouble in other places:
places that don't check if coeff is zero, for example:
(16) -> z := y*t
(16) 1.e-200 t
Type: Polynomial(DoubleFloat)
(17) -> z*y
(17)
Type: Polynomial(DoubleFloat)
(18) -> (z*y) pretend SEX
(18) (1 t (1 0 . 0.0))
Type: SExpression
(19) -> z*y + 1
(19) + 1.0
Type: Polynomial(DoubleFloat)
--
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.
