> I recently had a need for Jacobi polynomials with the property that > > int(P(i, a, b, x)*P(j, a, b, x)*(1-x)^a*(1+b)^b, (x, -1, 1) = > delta_ij > > with the key element being the normalization factor. Would it be > possible to get this upstream (it is a pain to code up!), for example > as a norm=True kwarg?
I think that this should go upstream and I'm preparing a PR for this. But there are a few open points: - What about the *symbolic* representation? (Unevaluated P(j,a,b,x) with j,a,b unspecified.) Are they normalized or not? Does the flag take action here too? How to get from one to the other form? (Something like "rewrite") And how to differ them in pretty prints, latex etc? - What if we have more different normalizations? F.e. for Hermite polynomials we have the probabilists' and the physicists' Hermite polynomials. [1] For spherical harmonics there are even more conventions [2]. How should we call the flag? And what values should it take? Maybe "normalization" with enumerated values "none", "geodesic", "magnetic" etc. But these values could be different for each polynomial type! I'd like to see a general solution we can apply to all orthogonal polynomials and related functions. [1]: http://en.wikipedia.org/wiki/Hermite_polynomial [2]: http://en.wikipedia.org/wiki/Spherical_harmonics#Orthogonality_and_normalization -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to this group, send email to sympy@googlegroups.com. Visit this group at http://groups.google.com/group/sympy?hl=en-US. For more options, visit https://groups.google.com/groups/opt_out.
