Re: Issue 3311 in sympy: SymPy is lacking a expressiontree facility / creator

2012-07-06 Thread sympy
Comment #6 on issue 3311 by jari-pek...@kolumbus.fi: SymPy is lacking a expressiontree facility / creator http://code.google.com/p/sympy/issues/detail?id=3311 SymPy interface to Axiom CAS would also be possible to use in integrations and symbolic manipulation if SymPy has one. -- You

Re: [sympy] sympy.plotting.experimental_lambdify vs. sympy.utilities.lambdify

2012-07-06 Thread Vinzent Steinberg
Am Dienstag, 3. Juli 2012 11:09:27 UTC+2 schrieb Stefan Krastanov: The one exception is code inside sympy which was agreed at some point, should not require more than a python interpreter. However, if you are using nsolve with lambdify you are again doing it wrong, because the only

Re: [sympy] sympy.plotting.experimental_lambdify vs. sympy.utilities.lambdify

2012-07-06 Thread Vinzent Steinberg
Am Dienstag, 3. Juli 2012 22:00:41 UTC+2 schrieb Ondřej Čertík: There is a difference between calling matplotlib, numpy or scipy, and having a working fortran or C compiler working at your machine. In particular, in scipy ODE, you can provide a simple Python function to get the ODE

[sympy] Re: autowrap a matrix expression

2012-07-06 Thread Gustavo
Is it possible to accomplish this using an out parameter? This would still be useful to me. On Monday, January 23, 2012 2:20:50 PM UTC-5, Gustavo wrote: Right now it returns a float. I'd like an NumPy array, so maybe I can create the NumPy array with the correct size and row/column

Re: [sympy] Re: autowrap a matrix expression

2012-07-06 Thread Matthew Rocklin
It would be nice if SymPy could do this. I'm unfortunately not very familiar with the code generation module. If it's not possible now then we should make an issue for it. Alternatively you could take a look at the theano projecthttp://deeplearning.net/software/theano/. They optimize and generate

Re: [sympy] Re: autowrap a matrix expression

2012-07-06 Thread Matthew Rocklin
On a semi-related note I think it'd be awesome if SymPy created a Theano code generator. Our two projects are related and I think that there is something to be gained by cooperation. We optimize and simplify expressions in semi-orthogonal ways. On Fri, Jul 6, 2012 at 10:45 AM, Matthew Rocklin

[sympy] can not do integral of sqrt(sin(x))

2012-07-06 Thread pallab
It seems sympy can not integrate sqrt(sin(x)). I did the following: import sympy as sm from sympy.abc import x,y,z tointegrate=sm.sqrt(sm.sin(y)) sm.integrate(tointegrate) output is : Integral(sqrt(sin(y)), y) After a simple change of variable the integral is doable: def ytoz(z):

Re: [sympy] Test failures in master: Hash randomization

2012-07-06 Thread Vinzent Steinberg
Am Dienstag, 3. Juli 2012 19:16:37 UTC+2 schrieb Joachim Durchholz: Am 01.07.2012 22:48, schrieb Sergiu Ivanov: [0] http://pypi.python.org/pypi/ordereddict We want key order. OrderedDict remembers insertion order, so it's close but doesn't match. There is a pure Python implementation

[theano-dev] Re: [sympy] Re: autowrap a matrix expression

2012-07-06 Thread Frédéric Bastien
Hi, I agree that more cooperation would be great. Up to now the problem on our side is time to look at SymPy. I'm going to SciPy in ~1 week in Texas. I was under the impression there was a tutorial/talk about SimPy There, but I don't find it anymore. Is it a sprint? The sprints list is not on

Re: [theano-dev] Re: [sympy] Re: autowrap a matrix expression

2012-07-06 Thread Aaron Meurer
On Jul 6, 2012, at 12:41 PM, Frédéric Bastien no...@nouiz.org wrote: Hi, I agree that more cooperation would be great. Up to now the problem on our side is time to look at SymPy. I'm going to SciPy in ~1 week in Texas. I was under the impression there was a tutorial/talk about SimPy There,

Re: [sympy] unable to solve system with more equations than unknowns?

2012-07-06 Thread Aaron Meurer
On Thu, Jul 5, 2012 at 7:48 PM, Chris Smith smi...@gmail.com wrote: On Thu, Jul 5, 2012 at 5:56 PM, Andrew zxwand...@gmail.com wrote: Hi, I am very new to sympy and python, so pardon my simple question I have a matrix right here like this (4 equations, 3 unknowns) mat1= [ 0,15.0, 10.0,

Re: [sympy] Test failures in master: Hash randomization

2012-07-06 Thread Joachim Durchholz
Am 03.07.2012 20:34, schrieb Sergiu Ivanov: On Tue, Jul 3, 2012 at 8:16 PM, Joachim Durchholz j...@durchholz.org wrote: Am 01.07.2012 22:48, schrieb Sergiu Ivanov: [0] http://pypi.python.org/pypi/ordereddict We want key order. OrderedDict remembers insertion order, so it's close but

Re: [sympy] can not do integral of sqrt(sin(x))

2012-07-06 Thread Aaron Meurer
WolframAlpha gives an answer in terms of an elliptic integral: http://www.wolframalpha.com/input/?i=integrate(sqrt(sin(x)),%20x) So I imagine that if we added those that we would be able to compute this. Aaron Meurer On Fri, Jul 6, 2012 at 10:00 AM, pallab pallabb...@gmail.com wrote: It

Re: [sympy] can not do integral of sqrt(sin(x))

2012-07-06 Thread rl
The problem I see is *how* to add these to the integration routines. There is no nice Meijer-G representation. And I suppose that Risch can not handle these. At least not w/o major extensions. -- You received this message because you are subscribed to the Google Groups sympy group. To view this

Re: [sympy] can not do integral of sqrt(sin(x))

2012-07-06 Thread Aaron Meurer
But if Pallab's answer is correct, then it can be, right, at least for this specific one? Aaron Meurer On Fri, Jul 6, 2012 at 3:03 PM, rl someb...@bluewin.ch wrote: The problem I see is *how* to add these to the integration routines. There is no nice Meijer-G representation. And I suppose that

Re: [sympy] can not do integral of sqrt(sin(x))

2012-07-06 Thread Aaron Meurer
According to http://en.wikipedia.org/wiki/Elliptic_integral, every elliptic integral can be reduced to an expression containing only the three Legendre canonical forms. So I wonder if that is algorithmic. Aaron Meurer On Fri, Jul 6, 2012 at 4:08 PM, Aaron Meurer asmeu...@gmail.com wrote: But

[sympy] Re: can not do integral of sqrt(sin(x))

2012-07-06 Thread pallab
The answer seems to be wrong. I do the following to numerically integral evaluation: import numpy as np import scipy as sp import scipy.integrate f=lambda x:np.sqrt(np.sin(x)) sp.integrate.quad(f,0.6,0.7) get the following value: 0.07776347731181982 which matches with mathematica. whereas

[sympy] Re: can not do integral of sqrt(sin(x))

2012-07-06 Thread pallab
By answer wrong I mean the answer given by sympy, On Friday, July 6, 2012 10:06:57 PM UTC-4, pallab wrote: The answer seems to be wrong. I do the following to numerically integral evaluation: import numpy as np import scipy as sp import scipy.integrate f=lambda x:np.sqrt(np.sin(x))