[sympy] Re: Feynman diagrams

2013-06-30 Thread mario
In https://github.com/sympy/sympy/pull/1699 in group_factors the theoretical group factors for SU(N) are computed; for instance in test_group_factors it is computed the group factor for the four-loop gluon Feynman graph in eq.(1.1) of P. Cvitanovic "Group Theory" On Saturday, June 29, 2013 7:18:30

Re: [sympy] Circle and it's equation:subs() and solve()

2013-06-30 Thread Aaron Meurer
On Sat, Jun 29, 2013 at 7:34 PM, Amit Saha wrote: > Hello, > > This is more of a note of my findings and also just implicitly > verifying whether I am doing the right thing. > > Problem: I want to find whether a point lies on a circle. > > I was hoping that like Linear entity objects, the circle/e

Re: [sympy] Diophantine Module: Solution Representation

2013-06-30 Thread Aaron Meurer
If it's a recurrence, just return the recurrence relation, and any initial conditions to that relation (i.e., the same form as the input to rsolve). It might be useful to have a helper function that could take such a relation and generate the nth term of each function, by the way. Aaron Meurer O

Re: [sympy] exp function on Matrix

2013-06-30 Thread Aaron Meurer
It seems I misremembered _eval_power. It handles when the object is the base, not the exponent. So I'm actually not sure if there is a way to override this. This ties in to the dispatch thing that is often discussed here. Aaron Meurer On Wed, Jun 26, 2013 at 5:48 PM, F. B. wrote: > I don't under

Re: [sympy] Solving a quadratic symbolically (Sympy vs. Mathematica)

2013-06-30 Thread Aaron Meurer
It's also possible that Mathematica automatically returns numeric solutions when explicit ones can't be found. Aaron Meurer On Sun, Jun 30, 2013 at 5:48 PM, Stefan Krastanov wrote: > For the moment mathematica is more advanced than sympy both in finding > symbolic solutions and in finding _multi

Re: [sympy] Feynman diagrams

2013-06-30 Thread Brian Granger
I don't know the physics well enough to know from your pseudo code. sympy.physics.quantum knows about quantum mechanics though: * Operators, states (bras, kets), commutators, anticommutators, outer projects, tensor projects, inner projects * The representation of those entities in different bases

Re: [sympy] Solving a quadratic symbolically (Sympy vs. Mathematica)

2013-06-30 Thread Stefan Krastanov
For the moment mathematica is more advanced than sympy both in finding symbolic solutions and in finding _multiple_ numeric solutions. This explains most of the differences in results. -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from

Re: [sympy] Solving a quadratic symbolically (Sympy vs. Mathematica)

2013-06-30 Thread Justin Carden
So, I switched to the 'sympy.nsolve' method and got the approximate solutions similar to Mathematica. Yes, the system is nonlinear. While the method in the Mathematica code was 'NSolve', the solution set returned was the same for the 'Solve' method. Regardless, this helped a bunch. Thanks! On S

Re: [sympy] Feynman diagrams

2013-06-30 Thread Stefan Krastanov
@Brian, what would be the correct way to write something like gluon(lorentz_index, su3_index) * group_generator(su3_index, antispinor_index, spinor_index) * top_quark(spinor_index) in the physics.quantum module? On 30 June 2013 21:31, Brian Granger wrote: > Extra bonus points if this stuff can

Re: [sympy] Solving a quadratic symbolically (Sympy vs. Mathematica)

2013-06-30 Thread Stefan Krastanov
The equations are nonlinear and it seems that sympy is unable to solve them. Maybe there is a way to instruct sympy how to deal with them but I do not know how. However, there is a very big difference between your mathematica code and your sympy code. In mathematica you are calling `NSolve`, which

Re: [sympy] Solving a quadratic symbolically (Sympy vs. Mathematica)

2013-06-30 Thread Justin Carden
Hi Stefan, Both of these solutions return an empty array. Where I should get an array of possible solutions similar to the output of the Mathematica version. Such as: {{x -> -1.4282, y -> -1.08259, z -> -1401.51}, {x -> -6.67923, y -> -598.932, z -> -0.500032}, {x -> -2.17279 - 53.7467 I, y

Re: [sympy] Solving a quadratic symbolically (Sympy vs. Mathematica)

2013-06-30 Thread Stefan Krastanov
f_1 = x + s1 + s2 + s5 - t1 this is all that you need to do if I understand your question correctly Or if you wish, you can create 'Eq(right_hand, left_hand)' instances. On 30 June 2013 21:09, Justin Carden wrote: > Hi all, > > I'm trying to make the switch from Mathematica to Python in the lab

Re: [sympy] Feynman diagrams

2013-06-30 Thread Brian Granger
Extra bonus points if this stuff can be integrated with sympy.physics or even better sympy.physics.quantum...otherwise we just keep creating new corners of sympy that don't work together. I know it is easier, but in the long run it is a horrible situation. The stuff in sympy.physics.quantum is co

[sympy] Solving a quadratic symbolically (Sympy vs. Mathematica)

2013-06-30 Thread Justin Carden
Hi all, I'm trying to make the switch from Mathematica to Python in the lab, but I'm running into a small problem binding the result space to a boundary. Specifically, I'm trying to assign equality to each polynomial equation in a system to a total t1,t2 and t3 to solve the system symbolically.

Re: [sympy] Re: Release Notes Request

2013-06-30 Thread Ondřej Čertík
I've added all changes from this thread into: https://github.com/sympy/sympy/wiki/Release-Notes-for-0.7.3 Feel free to polish it there. For myself I think I only contributed the Gauss-Legendre and Gauss-Laguerre points and weights, so I put it there as well. On Sat, Jun 29, 2013 at 8:04 AM, Aar

Re: [sympy] Diophantine Module: Solution Representation

2013-06-30 Thread Ondřej Čertík
On Sun, Jun 30, 2013 at 3:00 AM, Thilina Rathnayake wrote: > Looks like we will run into more and more trouble representing the > solutions. > > The latest is that when solving quadratic Diophantine equation, > A*x**2 + B*x*y + C*y**2 + D*x + E*y + F = 0, for the case B**2 - 4AC > 0, > when we kno

Re: [sympy] Feynman diagrams

2013-06-30 Thread Ondřej Čertík
Hi Stefan, This is cool. I would be interested in having this in sympy. Long time ago I wrote some diagram generating code for scalar diagrams: https://github.com/certik/sympy/blob/wick/t.py but I didn't have time to polish it up and send a PR. As you correctly mentioned, there are lots of steps

Re: [sympy] planet.sympy

2013-06-30 Thread Ondřej Čertík
On Sat, Jun 29, 2013 at 12:14 PM, Matthew Rocklin wrote: > Ondrej got planet.sympy.org up and running again. > > There are now a lot of great blogposts by our GSoC students. Take a look! > > Kudos to Ondrej. No problem. For reference, there were 3 issues why it didn't work: 1) The repository at

Re: [sympy] Unbundle mpmath

2013-06-30 Thread Ondřej Čertík
On Sat, Jun 29, 2013 at 1:52 AM, Sergey B Kirpichev wrote: > On Fri, Jun 28, 2013 at 09:38:35PM -0500, Ondřej Čertík wrote: >> >> That being said, if projects like binstar (https://binstar.org/), >> >> which was announced *yesterday*, take off and allow easy >> >> installation on all platforms (in

Re: [sympy] Diophantine Module: Solution Representation

2013-06-30 Thread Thilina Rathnayake
Thank you Sergey for your reply. On Sun, Jun 30, 2013 at 4:27 PM, Sergey Kirpichev wrote: > > On Sunday, June 30, 2013 2:49:20 PM UTC+4, Thilina Rathnayake wrote: >> >> We can solve this difference equation and output the general solution >> perhaps. >> Is there a way to solve recurrences in sy

Re: [sympy] Diophantine Module: Solution Representation

2013-06-30 Thread Sergey Kirpichev
On Sunday, June 30, 2013 2:49:20 PM UTC+4, Thilina Rathnayake wrote: > > We can solve this difference equation and output the general solution > perhaps. > Is there a way to solve recurrences in sympy? I searched and found there > is a > function called rsolve() but couldn't find it's documentat

Re: [sympy] Diophantine Module: Solution Representation

2013-06-30 Thread Thilina Rathnayake
We can solve this difference equation and output the general solution perhaps. Is there a way to solve recurrences in sympy? I searched and found there is a function called rsolve() but couldn't find it's documentation. On Sun, Jun 30, 2013 at 2:30 PM, Thilina Rathnayake wrote: > Looks like we w

[sympy] Diophantine Module: Solution Representation

2013-06-30 Thread Thilina Rathnayake
Looks like we will run into more and more trouble representing the solutions. The latest is that when solving quadratic Diophantine equation, A*x**2 + B*x*y + C*y**2 + D*x + E*y + F = 0, for the case B**2 - 4AC > 0, when we know a basic solution, all the other solutions can be represented as a re