<div dir="ltr">Hello guys again,<div><br></div><div>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:</div><div><div style="margin:0px;padding:0px;border:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial"><br></div><div style="margin:0px;padding:0px;border:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial">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<span> </span><a href="http://nbviewer.jupyter.org/gist/celliern/b38158d04d9dc3d8079dc44e3b747ac8" target="_blank" rel="nofollow" style="margin:0px;padding:0px;border:0px;text-decoration:none;color:rgb(102,17,204)">here</a> by <a href="https://github.com/celliern" target="_blank" rel="nofollow" style="margin:0px;padding:0px;border:0px;text-decoration:none;color:rgb(102,17,204)">Nicolas CELLIER</a> and some ideas over<span> </span><a href="https://cs.stackexchange.com/questions/95886/algorithm-for-using-power-series-to-numerically-solve-a-partial-differential-equ" target="_blank" rel="nofollow" style="margin:0px;padding:0px;border:0px;text-decoration:none;color:rgb(102,17,204)">here</a><span> </span>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<span> </span><a href="https://cs.stackexchange.com/questions/95886/algorithm-for-using-power-series-to-numerically-solve-a-partial-differential-equ" target="_blank" rel="nofollow" style="margin:0px;padding:0px;border:0px;text-decoration:none;color:rgb(102,17,204)">this question</a><span> </span>and help us out. </div><div style="margin:0px;padding:0px;border:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial"><br></div><div style="margin:0px;padding:0px;border:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial">FYI I have made a issue on SymPy's github over<span> </span><a href="https://github.com/sympy/sympy/issues/15015#" target="_blank" rel="nofollow" style="margin:0px;padding:0px;border:0px;text-decoration:none;color:rgb(102,17,204)">here</a><span> </span>which you can follow. </div><div style="margin:0px;padding:0px;border:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial"><br></div><div style="margin:0px;padding:0px;border:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial">Best,</div><div style="margin:0px;padding:0px;border:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial">Foad</div><br></div><br><div class="gmail_quote"><div dir="ltr">On Thu, Aug 2, 2018 at 9:40 AM Foad Sojoodi Farimani <<a href="mailto:f.s.farimani@gmail.com">f.s.farimani@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello everyone,<div><br></div><div>I want to generate a multivariate polynomial given an array of nonnegative integers D=[d1,...,dm]. 

<span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">I want to replicate what I have been suggested in SymPy, in GiNaC, to see if I can get a better performance. (see<span> </span></span><a href="https://groups.google.com/forum/#!topic/sympy/V4KRzA-3doI" style="color:rgb(17,85,204);font-size:small;background-color:rgb(255,255,255)" target="_blank">here</a><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">,<span> </span></span><a href="https://gist.github.com/Foadsf/07d4649f94f4d71876612a0ec6d9939b" style="color:rgb(17,85,204);font-size:small;background-color:rgb(255,255,255)" target="_blank">here<span> </span></a><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">and<span> </span></span><a href="https://stackoverflow.com/questions/51628270/how-to-generate-a-symbolic-multivariate-polynomial-of-a-given-dimension-in-sympy" style="color:rgb(17,85,204);font-size:small;background-color:rgb(255,255,255)" target="_blank">here</a><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">) I would appreciate if you could help me with this.</span><br></div><div><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></div><div><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Best,</span></div><div><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Foad</span></div></div>
</blockquote></div></div>