Hello guys again,

I'm not sure if you have received my former email but I'm gonna reply to that so you can read it as well:

While ago I was trying to use polynomial power series to solve a system of partial differential and algebraic equations, when realized there is no implementation of the idea. There is only Mathematica's AsymptoticDSolveValue which is just for ODEs. So I decided to implement it myself. Thanks to the Sympy community we now have some progress. I have one implementation over here by Nicolas CELLIER and some ideas over here by me and some of SymPy developers. I was wondering if we could join forces to come of with a general algorithm, then implementionation on diffrent languages shouldn't be that difficult. I was wondering if you could take a look at this question and help us out. 

FYI I have made a issue on SymPy's github over here which you can follow. 

Best,
Foad


On Thu, Aug 2, 2018 at 9:40 AM Foad Sojoodi Farimani <f.s.farimani@gmail.com> wrote:
Hello everyone,

I want to generate a multivariate polynomial given an array of nonnegative integers D=[d1,...,dm].  I want to replicate what I have been suggested in SymPy, in GiNaC, to see if I can get a better performance. (see here, here and here) I would appreciate if you could help me with this.

Best,
Foad