Re: [sympy] type of cse(matrix)[1]

2011-12-20 Thread Chris Smith
I've thought this about functions, too. It's easy enough to see what was given and to just return the same thing: iterable in, iterable out; expr in, expr out. What do others think? -- You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group,

Re: [sympy] Integrating over multivariate dirac deltafunctions

2011-12-20 Thread Matthew Rocklin
@Tom - I agree that it makes sense to separate out the common easy case from the very rare very difficult case. @Stephan - In the multivariate case the zeros are likely to consist of continuous lines/surfaces/manifolds of co-dimension 1. We'll need to do an integral over complex surfaces. This coul

Re: [sympy] Integrating over multivariate dirac deltafunctions

2011-12-20 Thread krastanov.ste...@gmail.com
@Tom, the operation is well defined for all sufficiently nice functions. @Matthew, The notation that is used there is a bit unclear (for me). The domain is all zeros of g and the sigma is just some measure theory notation meaning that you should actually take the sum. To solve this in sympy you ca

Re: [sympy] Integrating over multivariate dirac deltafunctions

2011-12-20 Thread Tom Bachmann
I think there are two basic cases to consider: 1) delta functions with *linear* arguments. This is (I believe) by a large margin the most common case. For this deltaintegrate should simply be fixed; one only needs to correctly implement integral(f(x) delta(a*x+b), (x, c, d)) and this shouldn't

[sympy] Integrating over multivariate dirac deltafunctions

2011-12-20 Thread Matthew Rocklin
Hi Everyone, I'd like to compute multivariate integrals that contain Dirac Deltafunctions. I.e. expressions like integrate(exp(-(x**2+y**2))/pi * delta(2*x+3*y), (x,-oo, oo), (y,-oo, oo)) The deltaintegrate function inside sympy fails to compute these correctly, see issue 2630.

[sympy] type of cse(matrix)[1]

2011-12-20 Thread Bastian Weber
Hello, Currently I am dealing with common subexpressions of a matrix. In [1]: from sympy import * In [2]: from sympy.abc import * In [3]: M = Matrix([a**2 + b, a**2 + b+ c]) In [4]: v1, v2 = cse(M) In [5]: type(v2) Out[5]: My question is: wouldn't it be more intuitive if v2 is a matrix

Re: [sympy] Re: need suggestions to the new concept of using SymPy

2011-12-20 Thread Aaron Meurer
On Tue, Dec 6, 2011 at 2:22 PM, Tim Lahey wrote: > On Sat, Nov 26, 2011 at 11:37 PM, Aaron Meurer wrote: >> Let us know how this goes.  I don't have an Android device, though I >> do have iOS devices, so I would be interested if you port to that >> platform. > > Aaron, > > There's a universal app

Re: [sympy] assumptions

2011-12-20 Thread Aaron Meurer
To be sure, I think we should consider keeping the old assumptions syntax (.is_positive, Symbol(positive=True), etc.), but make them call the new assumptions. This way, we will break very little compatibility, and we will still have these useful shortcuts for the simple cases. Aaron Meurer 2011/