Issue 1438 in sympy: bug in intersection: segment - segment

2009-05-26 Thread codesite-noreply
Status: New Owner: Labels: Type-Defect Priority-Medium New issue 1438 by flor...@mickler.org: bug in intersection: segment - segment http://code.google.com/p/sympy/issues/detail?id=1438 entity1 = Segment(Point(-10,10),Point(10,10)) entity2 = Segment(Point(-5,-5),Point(-5,5)) print

Issue 1390 in sympy: setup.py --record: manpage of isympy.1 is not listed correctly

2009-05-26 Thread codesite-noreply
Comment #2 on issue 1390 by werner.ho: setup.py --record: manpage of isympy.1 is not listed correctly http://code.google.com/p/sympy/issues/detail?id=1390 I'm sorry. The error is not inside sympy. It's a rpm packaging error. Explanation: The rpm build script of openSUSE compresses all man

Issue 1438 in sympy: bug in intersection: segment - segment

2009-05-26 Thread codesite-noreply
Comment #1 on issue 1438 by flor...@mickler.org: bug in intersection: segment - segment http://code.google.com/p/sympy/issues/detail?id=1438 attaching now test for vertical segment, which fails with current implementation and a patch which solves the containment tests for vertical lines a

Issue 1439 in sympy: polynomial solution missing real roots

2009-05-26 Thread codesite-noreply
Status: New Owner: Labels: Type-Defect Priority-Medium New issue 1439 by smichr: polynomial solution missing real roots http://code.google.com/p/sympy/issues/detail?id=1439 The equation -1936 - 5056*x - 7592*x**2 + 2704*x**3 - 49*x**4 has two real roots at about 3.6 and 52 but only

Issue 1441 in sympy: integrate takes forever

2009-05-26 Thread codesite-noreply
Status: New Owner: asmeurer Labels: Type-Defect Priority-Medium Integration New issue 1441 by asmeurer: integrate takes forever http://code.google.com/p/sympy/issues/detail?id=1441 I ran the below for several hours, and when after that it had still not returned anything, I gave up. It should

Issue 1127 in sympy: integrate(2*a/((((2*a)**2+x**2))**(3/2)),x) gives traceback

2009-05-26 Thread codesite-noreply
Updates: Labels: Integration Comment #4 on issue 1127 by asmeurer: integrate(2*a/2*a)**2+x**2))**(3/2)),x) gives traceback http://code.google.com/p/sympy/issues/detail?id=1127 I have found another integral that does this: integrate(-(1 + (-x + x**2)**(1/2))/(-x + (1 + (-x +

Issue 1442 in sympy: mathml patches

2009-05-26 Thread codesite-noreply
Status: Accepted Owner: Labels: Type-Defect Priority-Medium NeedsReview Milestone-Release0.6.5 New issue 1442 by ondrej.certik: mathml patches http://code.google.com/p/sympy/issues/detail?id=1442 http://groups.google.com/group/sympy-patches/browse_thread/thread/af189fc4af1e9748 -- You

Issue 1390 in sympy: setup.py --record: manpage of isympy.1 is not listed correctly

2009-05-26 Thread codesite-noreply
Updates: Status: Invalid Comment #3 on issue 1390 by ondrej.certik: setup.py --record: manpage of isympy.1 is not listed correctly http://code.google.com/p/sympy/issues/detail?id=1390 Ok, excellent. If you find more bugs, just report them as well, until we figure out whose bug it is

Issue 1438 in sympy: bug in intersection: segment - segment

2009-05-26 Thread codesite-noreply
Updates: Labels: NeedsReview Milestone-Release0.6.5 Comment #2 on issue 1438 by ondrej.certik: bug in intersection: segment - segment http://code.google.com/p/sympy/issues/detail?id=1438 (No comment was entered for this change.) -- You received this message because you are listed in

Re: Resending logcombine patch.

2009-05-26 Thread Alan Bromborsky
asmeurer wrote: I see that the gods do not want me to get this in. Apparently the files did not make it to the last message. I am just going to do what I know works: below is the text copied and pasted from the 4 patches. If it is easier, you can still pull from my github repo with the

Re: Resending logcombine patch.

2009-05-26 Thread Aaron S. Meurer
What are you using for an email program? My email program is not the problem. I use Apple Mail, but those send just fine. The problem was in sending patches through git using the send-email command. The files were corrupted because they were in the wrong character encoding. For

Re: [PATCH] Fix for ISSUE 626

2009-05-26 Thread Ondrej Certik
The second patch is +1 and I think I reviewed the first patch already and I had some comments to it. Ondrej On Fri, May 22, 2009 at 11:42 AM, Luke hazelnu...@gmail.com wrote: Issue 626 is fixed by these two patches.  Also, solve can now Symbol, Function, and Derivative types. I implemented

Re: Updates and Fixeds to Geometric Algebra - module, tests, examples, and documentation

2009-05-26 Thread Ondrej Certik
It fails the following tests: $ bin/test sympy/utilities/ = test process starts == executable: /usr/bin/python (2.6.2-final-0) sympy/utilities/tests/test_code_quality.py[2] F. [FAIL]

Re: 2 MathML Patches

2009-05-26 Thread Ondrej Certik
Let's wait for Fabian's review. http://code.google.com/p/sympy/issues/detail?id=1442 Thanks for the patches! Ondrej On Mon, May 25, 2009 at 4:37 PM, Fabian Pedregosa fab...@fseoane.net wrote: Thanks for the patches, thanks for making the appropriate changes. I am quite busy right now with

Re: Fixes and tests for Issue 1434

2009-05-26 Thread Ondrej Certik
Ok. Fabian discovered though that the patch makes some other (untested yet) things fail. So we need to fix it the right way --- Fabian, are you still working on that? Ondrej On Sun, May 24, 2009 at 2:01 PM, Luke hazelnu...@gmail.com wrote: I put it in there, but commented out, so that it

Re: [PATCH] Fix for ISSUE 626, cleaner implementation, all as one patch

2009-05-26 Thread Ondrej Certik
The patch looks good, it's in. Thanks! On Mon, May 25, 2009 at 4:44 PM, Luke hazelnu...@gmail.com wrote: Attached is the patch that Fixes issue 626 and allows for solve() to handle Function and Derivative instances, implements tests for both cases, and has extensive commenting. All tests

Re: Fix shebang line in test scripts to /usr/bin/env python

2009-05-26 Thread Ondrej Certik
Hi Aaron -- I wanted to push this in, but git is telling me: $ git am ~/Desktop/0001-Fixed-shebang-lines-in-test-scripts-to-usr-bin-env\(3\).patch Patch does not have a valid e-mail address. I am too busy now to figure out what went wrong. I don't see why the email address should be wrong it

Re: Resending logcombine patch.

2009-05-26 Thread Ondrej Certik
On Tue, May 26, 2009 at 10:08 AM, Aaron S. Meurer asmeu...@gmail.com wrote: What are you using for an email program? My email program is not the problem.  I use Apple Mail, but those send just fine.  The problem was in sending patches through git using the send-email command.  The files were

Re: Resending logcombine patch.

2009-05-26 Thread Aaron S. Meurer
OK. Those were actually on my github account in the odes branch. It is possible to pull the commits from the hashes, no? I provided them in the main thread. Did you think the logcombine function was ok? Aaron Meurer On May 26, 2009, at 11:02 PM, Ondrej Certik wrote: On Tue, May 26,

[sympy] Re: trigsimp

2009-05-26 Thread Luke
Thanks. On May 25, 9:02 pm, Aaron S. Meurer asmeu...@gmail.com wrote: You need to use the exclude option I was telling you about.  Do,   a = Wild('a', exclude=[pi]) and   b = Wild('a', exclude=[pi]). and you get   e1.match(a*pi + b) {a_: 1, b_: x}   b = Wild('b', exclude=[pi])  

[sympy] Best way to determine equivalence of 1 and numbers.One

2009-05-26 Thread Aaron S. Meurer
I wrote a function for dsolve that determine the homogeneous order of an expression. A function f(x,y) is or order n if f(xt,yt) == f(x,y)*t**n. The way that I did this was to pull the additive parts and determine the order of each part. It is only homogeneous if each additive part has

[sympy] Re: Best way to determine equivalence of 1 and numbers.One

2009-05-26 Thread Abderrahim Kitouni
On Tue, 26 May 2009 10:38:17 -0600 Aaron S. Meurer asmeu...@gmail.com wrote: [...] The problem is that sometimes, I will get something like set(0, numbers.Zero), which then fails. What would be the best way to determine the equivalence of things like 1, -1, 0, and Rational(1,2)

[sympy] Re: Best way to determine equivalence of 1 and numbers.One

2009-05-26 Thread Aaron S. Meurer
Thanks. I think that should do the trick. Aaron Meurer On May 26, 2009, at 11:40 AM, Abderrahim Kitouni wrote: On Tue, 26 May 2009 10:38:17 -0600 Aaron S. Meurer asmeu...@gmail.com wrote: [...] The problem is that sometimes, I will get something like set(0, numbers.Zero), which then

[sympy] Lambdify does not work for Ylm (spherical harmonics)

2009-05-26 Thread Sanjay
I'm trying to lambdify an expression that involves Ylm and I get the following error: x=Symbol('x') y=Symbol('y') g=lambdify((x,y),Ylm(3,2,x,y)) g(1.0,1.0) Traceback (most recent call last): File input, line 1, in module File string, line 1, in lambda TypeError: a float is required Am I

[sympy] Strange simplification

2009-05-26 Thread Neal Becker
Out[80]: (-2*s**2*w**2 + w**4)/(s**4 + w**4) In [81]: _.subs(s,2*pi*I*f) Out[81]: (8*pi**2*f**2*w**2 + w**4)/(w**4 + 16*pi**4*f**4) In [82]: simplify(_) Out[82]: (8*pi**2*f**2*w**2 + w**4)/(w**4 + 16*pi**4*f**4) Notice that this integrand should be strictly positive. It also should give a

[sympy] Re: Strange simplification

2009-05-26 Thread Luke
I get the same behavior on my machine. It seems like a bug to me. Not sure if the bug is in integrate or in simplify, but could you file this in the issues? Thanks, ~Luke On May 26, 12:16 pm, Neal Becker ndbeck...@gmail.com wrote: Out[80]: (-2*s**2*w**2 + w**4)/(s**4 + w**4) In [81]:

[sympy] Re: Strange simplification

2009-05-26 Thread Robert Kern
On Tue, May 26, 2009 at 15:08, Luke hazelnu...@gmail.com wrote: I get the same behavior on my machine. Hmm, I don't. In [1]: from sympy import * In [2]: %sym f w Adding variables: f w In [3]: s = 2*pi*I*f In [4]: ia = (-2*s**2*w**2 + w**4)/(s**4 + w**4) In [5]: simplify(integrate(ia,

[sympy] Re: trigsimp

2009-05-26 Thread Luke
Aaron, It seems like this only works for Symbol instances, not FunctionClass instances, i.e: from sympy import * x = Function('x') a = Wild('a', exclude=[pi]) b = Wild('b') e = pi/S(2) + x e.match(a*pi + b) Traceback (most recent call last): File stdin, line 1, in module File

[sympy] Re: Strange simplification

2009-05-26 Thread Neal Becker
On Tuesday 26 May 2009, Robert Kern wrote: On Tue, May 26, 2009 at 15:08, Luke hazelnu...@gmail.com wrote: I get the same behavior on my machine. Hmm, I don't. I'm using 0.6.4, are you using something newer perhaps? --~--~-~--~~~---~--~~ You received this

[sympy] Re: Strange simplification

2009-05-26 Thread Robert Kern
On Tue, May 26, 2009 at 17:04, Neal Becker ndbeck...@gmail.com wrote: On Tuesday 26 May 2009, Robert Kern wrote: On Tue, May 26, 2009 at 15:08, Luke hazelnu...@gmail.com wrote: I get the same behavior on my machine. Hmm, I don't. I'm using 0.6.4, are you using something newer perhaps?

[sympy] Re: Strange simplification

2009-05-26 Thread Neal Becker
On Tuesday 26 May 2009, Robert Kern wrote: On Tue, May 26, 2009 at 17:04, Neal Becker ndbeck...@gmail.com wrote: On Tuesday 26 May 2009, Robert Kern wrote: On Tue, May 26, 2009 at 15:08, Luke hazelnu...@gmail.com wrote: I get the same behavior on my machine. Hmm, I don't. I'm

[sympy] Re: Lambdify does not work for Ylm (spherical harmonics)

2009-05-26 Thread Sanjay
Strange...it's working now. Sorry for wasting your time! --~--~-~--~~~---~--~~ You received this message because you are subscribed to the Google Groups sympy group. To post to this group, send email to sympy@googlegroups.com To unsubscribe from this group, send

[sympy] Re: Lambdify does not work for Ylm (spherical harmonics)

2009-05-26 Thread Sanjay
OK...now...sorry once again. There seems to be a bug in the lambdify. When I tried lambdifying Ylm(4,0,t,p) it worked but it does not work with Ylm(4,1,t,p). I see that Ylm(4,1,t,p) has a exp(I*p) term. So, when I then try lambdifying exp(I*p) I get an error Traceback (most recent call last):

[sympy] Re: Strange simplification

2009-05-26 Thread Luke
I'm using the latest pull from git://git.sympy.org/sympy.git, and this is the response I get: In [1]: from sympy import * In [2]: f, w = symbols('fw') In [3]: s = 2*pi*I*f In [4]: ia = (-2*s**2*w**2 + w**4)/(s**4 + w**4) In [5]: simplify(integrate(ia, (f, 0, infty)))

[sympy] Re: Strange simplification

2009-05-26 Thread Robert Kern
On Tue, May 26, 2009 at 19:22, Luke hazelnu...@gmail.com wrote: I'm using the latest pull from git://git.sympy.org/sympy.git, and this is the response I get: In [1]: from sympy import * In [2]: f, w = symbols('fw') In [3]: s = 2*pi*I*f In [4]: ia = (-2*s**2*w**2 + w**4)/(s**4 + w**4)

[sympy] Re: Strange simplification

2009-05-26 Thread Luke
Ok, so there is some sort of bug. Here is what Matlab (maple) 2008a gives: int((8*pi^2*f^2*w^2+w^4)/(16*pi^4*f^4+w^4), f, 0, inf) ans = PIECEWISE([NaN, And(0 (w^4)^(1/4)*2^(1/2)*pi+pi*(-2*csgn(w^2)*w^2)^ (1/2),(w^4)^(1/4)*2^(1/2)*pi-pi*(-2*csgn(w^2)*w^2)^(1/2) 0,(w^4)^

[sympy] Re: Strange simplification

2009-05-26 Thread Neal Becker
On Tuesday 26 May 2009, Luke wrote: Ok, so there is some sort of bug. Here is what Matlab (maple) 2008a gives: int((8*pi^2*f^2*w^2+w^4)/(16*pi^4*f^4+w^4), f, 0, inf) ans = PIECEWISE([NaN, And(0 (w^4)^(1/4)*2^(1/2)*pi+pi*(-2*csgn(w^2)*w^2)^

[sympy] Function vs. Symbol

2009-05-26 Thread Luke
I'm trying to better understand how Sympy is structured with regard to Function and Symbol. In most problems I have encountered with ODE's (and PDE's, but I'll limit my discussion to ODE's) of the form: dx/dt = f(x, t) x \in R^n, f: R^n x R --- R^n there simply is no closed form solution

[sympy] Re: Function vs. Symbol

2009-05-26 Thread Aaron S. Meurer
What I wrote in the other thread is wrong. You can make x an implicit function of t. Do t = Symbol('t') x = Function('x')(t) x_ = diff(x, t) x x(t) x_ d ──(x(t)) dt See the docstring for dsolve. I agree that typing f(x) and f(x).diff(x) is very time consuming if you have a lot