Comment #2 on issue 4096 by [email protected]: Poly.coeff works only
after .simply().expand()
http://code.google.com/p/sympy/issues/detail?id=4096
Thanks for the feedback. I am using it as an Expr.
But if I do the following, using a Poly, I also appear to get bad behavior:
(Pdb) poly.as_poly().coeff(f0)
*** NotImplementedError: Either convert to Expr with `as_expr` method to
use Expr's coeff method or else use the `coeff_monomial` method of Polys.
(Pdb) poly.as_poly().coeff_monomial(f0)
0
Regards,
Ben
In case you're interested in trying it yourself with my exact use case,
here is the exact polynomial for one of the simplest cases:
(Pdb) poly.as_poly()
Poly(2*x0**2*x*f0*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3
- x1*x3**2 - x2*x3**2 + x3**3) + 2*x0**2*x*f0*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) -
x0**2*x1*f0*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) - x0**2*x1*f0*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) -
x0**2*x2*f0*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) - x0**2*x3*f0*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) -
3*x0*x**2*f0*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) - 3*x0*x**2*f0*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) +
2*x0*x*x1*f1*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) + 2*x0*x*x1*f1*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) +
2*x0*x*x2*f0*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) - 2*x0*x*x2*f1*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) -
2*x0*x*x2*f1*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) + 2*x0*x*x2*f2*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) +
2*x0*x*x3*f0*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) - 2*x0*x*x3*f2*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) -
2*x0*x*f0*1/(x1**2 - x1*x2 - x1*x3 + x2*x3) + x0*x1**2*f0*1/(-x0*x1*x2 +
x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) +
x0*x1**2*f0*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) - x0*x1**2*f1*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) -
x0*x1**2*f1*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) + x0*x1*x2*f1*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) -
x0*x1*x2*f2*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) - x0*x1*x3*f1*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) +
x0*x1*x3*f2*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) + x0*x2**2*f1*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) -
x0*x2**2*f2*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) + x0*x2*x3*f1*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) +
x0*x2*x3*f2*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) + x0*x2*f0*1/(x1**2 - x1*x2 - x1*x3 + x2*x3) +
x0*x3*f0*1/(x1**2 - x1*x2 - x1*x3 + x2*x3) + 3*x**2*x1*f0*1/(-x0*x1*x2 +
x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) +
3*x**2*x1*f0*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) - 3*x**2*x1*f1*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) -
3*x**2*x1*f1*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) + 3*x**2*x2*f1*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) +
3*x**2*x2*f1*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) - 3*x**2*x2*f2*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) +
3*x**2*x3*f2*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) + 3*x**2*f0*1/(x1**2 - x1*x2 - x1*x3 + x2*x3)
- 2*x*x1**2*f0*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) - 2*x*x1**2*f0*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) +
2*x*x1**2*f1*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) + 2*x*x1**2*f1*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) -
2*x*x1*x2*f0*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) - 2*x*x1*x2*f1*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) +
2*x*x1*x2*f2*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) - 2*x*x1*x3*f0*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) +
2*x*x1*x3*f1*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) - 2*x*x1*x3*f2*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) -
2*x*x2**2*f1*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) + 2*x*x2**2*f2*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) -
2*x*x2*x3*f1*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) - 2*x*x2*x3*f2*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) -
2*x*x2*f0*1/(x1**2 - x1*x2 - x1*x3 + x2*x3) - 2*x*x3*f0*1/(x1**2 - x1*x2 -
x1*x3 + x2*x3) + x1**2*x2*f0*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2
+ x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) - x1**2*x2*f1*1/(-x0*x1*x2 +
x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) +
x1**2*x3*f0*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) - x1**2*x3*f1*1/(x0*x1*x2 - x0*x1*x3 -
x0*x2**2 + x0*x2*x3 - x1*x2**2 + x1*x2*x3 + x2**3 - x2**2*x3) +
x1*x2**2*f1*1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 -
x1*x3**2 - x2*x3**2 + x3**3) - x1*x2**2*f2*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) +
x1*x2*x3*f1*1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2 +
x1*x2*x3 + x2**3 - x2**2*x3) + x1*x2*x3*f2*1/(-x0*x1*x2 + x0*x1*x3 +
x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 - x2*x3**2 + x3**3) +
x2*x3*f0*1/(x1**2 - x1*x2 - x1*x3 + x2*x3), x0, x, x1, x2, x3, f0, f1, f2,
1/(-x0*x1*x2 + x0*x1*x3 + x0*x2*x3 - x0*x3**2 + x1*x2*x3 - x1*x3**2 -
x2*x3**2 + x3**3), 1/(x0*x1*x2 - x0*x1*x3 - x0*x2**2 + x0*x2*x3 - x1*x2**2
+ x1*x2*x3 + x2**3 - x2**2*x3), 1/(x1**2 - x1*x2 - x1*x3 + x2*x3),
domain='ZZ')
--
You received this message because this project is configured to send all
issue notifications to this address.
You may adjust your notification preferences at:
https://code.google.com/hosting/settings
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" 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 http://groups.google.com/group/sympy-issues.
For more options, visit https://groups.google.com/groups/opt_out.
