Comment #14 on issue 4133 by someb...@bluewin.ch: Suggestion: implement
classical orthogonal polynomials
http://code.google.com/p/sympy/issues/detail?id=4133
Ok. I have started work on this. The current eval in the default case
is getting a pass instead of the poly eval. I will add a new method
for conversion to monomial expansions.
You can use this for example: "def _eval_rewrite_as_polynomial(self, ...):"
I'm not sure that jacobi and gegenbauer polynomials should
automatically convert to other types when the parameters correspond
to legendre and chebyshev etc.
No, they should not. In the new model, nothing should change to
simpler/other representation if not explicitly requested by
either simplify/expand/expand_func or rewrite.
I will also write a dictionary for numerical evaluation. NumPy
doesn't implement chebyshev of the second kind, gegenbauer or jacobi,
but I will write a routine for these.
Look at scipy. Maybe there is an implementation there.
Furthermore, I would suggest new classes for polynomial series, which
could be converted from one type to another.
Should I create a new fork for this?
Probably, yes.
--
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 sympy-issues+unsubscr...@googlegroups.com.
To post to this group, send email to sympy-issues@googlegroups.com.
Visit this group at http://groups.google.com/group/sympy-issues.
For more options, visit https://groups.google.com/groups/opt_out.
