Re: Issue 1720 in sympy: integrate(asin(x/2),x) crashes

2009-12-03 Thread sympy
Updates: Labels: NeedsReview Comment #10 on issue 1720 by plaes800: integrate(asin(x/2),x) crashes http://code.google.com/p/sympy/issues/detail?id=1720 Ahh.. thanks for pointing out the solution. I was also worried about all these simplify() calls. :) I went through all the code that

Re: Issue 1694 in sympy: solve has many issues with fractions

2009-12-03 Thread sympy
Comment #34 on issue 1694 by nicolas.pourcelot: solve has many issues with fractions http://code.google.com/p/sympy/issues/detail?id=1694 Chris, I agree with you, but the two issues are linked, since improving fraction handling will remove some singularities from the solutions, as you

Issue 1739 in sympy: Constants 'C_i' in differential equations confuse solver

2009-12-03 Thread sympy
Status: Accepted Owner: jorn.baayen CC: asmeurer Labels: Type-Defect Priority-Medium WrongResult New issue 1739 by jorn.baayen: Constants 'C_i' in differential equations confuse solver http://code.google.com/p/sympy/issues/detail?id=1739 The last call to dsolve() in the following snippet does

Re: Issue 1739 in sympy: Constants 'C_i' in differential equations confuse solver

2009-12-03 Thread sympy
Comment #2 on issue 1739 by jorn.baayen: Constants 'C_i' in differential equations confuse solver http://code.google.com/p/sympy/issues/detail?id=1739 My example was a very much simplified version of an actual problem I'm working with, which requires solving two different ODEs, where the

Re: Issue 1722 in sympy: atomic mass unit (amu, u) is defined wrong

2009-12-03 Thread sympy
Updates: Labels: -NeedsReview NeedsBetterPatch Comment #4 on issue 1722 by Vinzent.Steinberg: atomic mass unit (amu, u) is defined wrong http://code.google.com/p/sympy/issues/detail?id=1722 I think 'u' should indeed use the unit kg. BTW, Chris, to credit a reporter you can use

Issue 1740 in sympy: prime(symbol('x')) never returns

2009-12-03 Thread sympy
Status: New Owner: Labels: Type-Defect Priority-Medium New issue 1740 by Nerdbeard: prime(symbol('x')) never returns http://code.google.com/p/sympy/issues/detail?id=1740 sympy/ntheory/generate.py, line 118, in prime() is: assert n 0 if 'n' is a symbol, 'n 0' is a symbolic

Re: Issue 1733 in sympy: Remove useless reimplementation of 2 unused methods in functions.

2009-12-03 Thread sympy
Updates: Status: Fixed Labels: -Type-Defect -NeedsReview Type-Enhancement Comment #1 on issue 1733 by Vinzent.Steinberg: Remove useless reimplementation of 2 unused methods in functions. http://code.google.com/p/sympy/issues/detail?id=1733 Thank you, this is in. -- You

Re: Issue 1734 in sympy: Simplify Basic.rewrite

2009-12-03 Thread sympy
Updates: Status: Fixed Labels: -NeedsReview Comment #1 on issue 1734 by Vinzent.Steinberg: Simplify Basic.rewrite http://code.google.com/p/sympy/issues/detail?id=1734 Thanks, I pushed it in. -- You received this message because you are listed in the owner or CC fields of this

Re: Issue 1735 in sympy: Use expr.func to get the head of expressions.

2009-12-03 Thread sympy
Comment #2 on issue 1735 by Vinzent.Steinberg: Use expr.func to get the head of expressions. http://code.google.com/p/sympy/issues/detail?id=1735 Don't you think you could even use 'is' instead of '=='? For things like if arg.func == asinh: ... it makes sense in my opinion. This could really

Re: Issue 1736 in sympy: Refactor pattern matching and substitution

2009-12-03 Thread sympy
Updates: Status: Started Labels: -NeedsReview Comment #6 on issue 1736 by Vinzent.Steinberg: Refactor pattern matching and substitution http://code.google.com/p/sympy/issues/detail?id=1736 Some comments to your wip branch: Why do you use Rational(1) instead of S.one? What is

Issue 1741 in sympy: py.bench broken

2009-12-03 Thread sympy
Status: Accepted Owner: Labels: Type-Defect Priority-Medium New issue 1741 by Vinzent.Steinberg: py.bench broken http://code.google.com/p/sympy/issues/detail?id=1741 $ bin/py.bench Traceback (most recent call last): File bin/py.bench, line 15, in module from sympy.utilities import

Re: Issue 1731 in sympy: fraction can be more flexible

2009-12-03 Thread sympy
Updates: Labels: -NeedsReview NeedsBetterPatch Comment #7 on issue 1731 by Vinzent.Steinberg: fraction can be more flexible http://code.google.com/p/sympy/issues/detail?id=1731 (No comment was entered for this change.) -- You received this message because you are listed in the owner or

Re: Issue 1719 in sympy: implement from sympy.abc import lower, greek

2009-12-03 Thread sympy
Updates: Labels: -NeedsReview NeedsBetterPatch Comment #7 on issue 1719 by Vinzent.Steinberg: implement from sympy.abc import lower, greek http://code.google.com/p/sympy/issues/detail?id=1719 (No comment was entered for this change.) -- You received this message because you are listed

Re: Issue 1720 in sympy: integrate(asin(x/2),x) crashes

2009-12-03 Thread sympy
Updates: Labels: -NeedsBetterPatch Comment #12 on issue 1720 by Vinzent.Steinberg: integrate(asin(x/2),x) crashes http://code.google.com/p/sympy/issues/detail?id=1720 (No comment was entered for this change.) -- You received this message because you are listed in the owner or CC

Issue 1742 in sympy: how to implement functions for sympy

2009-12-03 Thread sympy
Status: New Owner: Vinzent.Steinberg CC: ondrej.certik Labels: Type-Defect Priority-Medium New issue 1742 by Vinzent.Steinberg: how to implement functions for sympy http://code.google.com/p/sympy/issues/detail?id=1742 From doc/src/guide.txt: [...] if isinstance(arg, Basic.NaN): return

Re: Issue 1741 in sympy: py.bench broken

2009-12-03 Thread sympy
Comment #1 on issue 1741 by asmeurer: py.bench broken http://code.google.com/p/sympy/issues/detail?id=1741 It works for me, but I have 0.9.2. So something in py must have broken it. -- You received this message because you are listed in the owner or CC fields of this issue, or because you

Re: Issue 1215 in sympy: fix html documentation

2009-12-03 Thread sympy
Updates: Summary: fix html documentation Comment #9 on issue 1215 by Vinzent.Steinberg: fix html documentation http://code.google.com/p/sympy/issues/detail?id=1215 While looking at the docs I found several issues: The ODE doc formatting is broken. ('== Usage ==' nor rendered properly

Re: Issue 1215 in sympy: fix html documentation

2009-12-03 Thread sympy
Comment #10 on issue 1215 by Vinzent.Steinberg: fix html documentation http://code.google.com/p/sympy/issues/detail?id=1215 We need to break the modules into categories, there are currently too many. -- You received this message because you are listed in the owner or CC fields of this issue, or

Re: Issue 1215 in sympy: fix html documentation

2009-12-03 Thread sympy
Comment #13 on issue 1215 by asmeurer: fix html documentation http://code.google.com/p/sympy/issues/detail?id=1215 I did not intend for the == Titles == to be rendered with Sphinx. I would like to do Title = but Sphinx will not allow it. I just wrote them so they would look nice in

Re: Issue 1215 in sympy: fix html documentation

2009-12-03 Thread sympy
Comment #14 on issue 1215 by asmeurer: fix html documentation http://code.google.com/p/sympy/issues/detail?id=1215 Also, one thing that I like about the way that I did the math in the ODE docs is that I made them doctests, so that way if the printer is improved, they can also be improved, and

Re: Issue 1736 in sympy: Refactor pattern matching and substitution

2009-12-03 Thread sympy
Comment #7 on issue 1736 by ronan.l...@gmail.com: Refactor pattern matching and substitution http://code.google.com/p/sympy/issues/detail?id=1736 Why do you use Rational(1) instead of S.one? For no good reason. It was just convenient to replace Rational(sign) with Rational(1). This part of

Re: Issue 1741 in sympy: py.bench broken

2009-12-03 Thread sympy
Comment #2 on issue 1741 by ronan.l...@gmail.com: py.bench broken http://code.google.com/p/sympy/issues/detail?id=1741 It's kinda written on the tin that this would break with any new release. Surely, a module named '__' is meant to be private! NB: I get the same traceback as Vinzent with

Issue 1743 in sympy: Vectorize subs and make a dummy function for it

2009-12-03 Thread sympy
Status: Accepted Owner: asmeurer Labels: Type-Defect Priority-Medium EasyToFix New issue 1743 by asmeurer: Vectorize subs and make a dummy function for it http://code.google.com/p/sympy/issues/detail?id=1743 The following should work: In [8]: s = x*pi In [9]: s.subs(x, range(10))

Re: Issue 1671 in sympy: runtests modifications and docstring corrections

2009-12-03 Thread sympy
Comment #61 on issue 1671 by asmeurer: runtests modifications and docstring corrections http://code.google.com/p/sympy/issues/detail?id=1671 It seems you were right. It works now. The tests runner seems to work fine for me, though I haven't looked at the code. Is there anything else that

Re: factoring upgrade

2009-12-03 Thread Mateusz Paprocki
Hi, On Wed, Dec 02, 2009 at 11:49:24PM -0800, smichr wrote: Is gcdfactor() any different from the following? In [1]: f = x**2*exp(x)+exp(x+y)*x/y In [2]: f Out[2]:            x + y  2  x   x⋅ℯ     x ⋅ℯ  +            y     In [3]: factor(f) Out[3]:   ⎛     y⎞  

Re: factoring upgrade

2009-12-03 Thread Vinzent Steinberg
There are several doctests failing on your polys2 branch. Vinzent -- You received this message because you are subscribed to the Google Groups sympy-patches group. To post to this group, send email to sympy-patc...@googlegroups.com. To unsubscribe from this group, send email to

Re: factoring upgrade

2009-12-03 Thread Aaron S. Meurer
The failures are all in the Sphinx file polynomials.txt. Aaron Meurer On Dec 3, 2009, at 11:21 PM, Ondrej Certik wrote: On Thu, Dec 3, 2009 at 4:38 PM, Vinzent Steinberg vinzent.steinb...@googlemail.com wrote: There are several doctests failing on your polys2 branch. Which ones? All

Re: factoring upgrade

2009-12-03 Thread Aaron S. Meurer
There are also failures that will not show up until Chris's runtest branch is merged in if you do: bin/doctest sympy/polys/galoistools.py bin/doctest sympy/polys/factortools.py bin/doctest sympy/polys/monomialtools.py bin/doctest sympy/polys/polyroots.py They are all based on the name of the

[sympy] Re: symbolic one and 2d tensors

2009-12-03 Thread Scott
My integral has several pieces like this: int^x1_x0 int ^t0_t1 (dot(del_a(t)).T* R *b ) dt dx del_a is 3*1, b is 3*1 and R is 3*3 rotation tensor( or a non orthogonal velocity transformation tensor H) R,H and b are are constant over the interval of integration. R=H inv(H.T) My goal is to

Re: [sympy] Re: symbolic one and 2d tensors

2009-12-03 Thread Alan Bromborsky
Scott wrote: My integral has several pieces like this: int^x1_x0 int ^t0_t1 (dot(del_a(t)).T* R *b ) dt dx del_a is 3*1, b is 3*1 and R is 3*3 rotation tensor( or a non orthogonal velocity transformation tensor H) R,H and b are are constant over the interval of integration. R=H inv(H.T)

Re: [sympy] Re: symbolic one and 2d tensors

2009-12-03 Thread Alan Bromborsky
Alan Bromborsky wrote: Scott wrote: My integral has several pieces like this: int^x1_x0 int ^t0_t1 (dot(del_a(t)).T* R *b ) dt dx del_a is 3*1, b is 3*1 and R is 3*3 rotation tensor( or a non orthogonal velocity transformation tensor H) R,H and b are are constant over the interval

[sympy] Re: printing exponential format vs float numbers

2009-12-03 Thread Vinzent Steinberg
On 2 déc, 22:37, archeryguru2000 archeryguru2...@gmail.com wrote: Is there any way to convert a sympy Matrix (or even numpy array, etc) to an mpmath matrix?  I have tried patching sympy.mpmath, but now I get errors trying to convert sympy to mpmath! Actually, is there a way to tell N(x,

[sympy] Re: printing exponential format vs float numbers

2009-12-03 Thread Vinzent Steinberg
Sorry for double-posting, sent previous message by mistake. On 2 déc, 22:37, archeryguru2000 archeryguru2...@gmail.com wrote: Is there any way to convert a sympy Matrix (or even numpy array, etc) to an mpmath matrix?  I have tried patching sympy.mpmath, but now I get errors trying to convert

[sympy] Re: symbolic one and 2d tensors

2009-12-03 Thread Scott
.T turns the column vector into a row vector. The integrands that involve the rotation tensors are functions of t or x , not both. int^x1_x0 int ^t0_t1 (dot(del_a(t)).T* R *b ) dt dx = (x1-x0)( del_a(t1).T-del_a(t0).T )* R* b = scalar I will test drive the galgebra module tonight. Is your

Re: [sympy] Re: symbolic one and 2d tensors

2009-12-03 Thread Alan Bromborsky
Scott wrote: .T turns the column vector into a row vector. The integrands that involve the rotation tensors are functions of t or x , not both. int^x1_x0 int ^t0_t1 (dot(del_a(t)).T* R *b ) dt dx = (x1-x0)( del_a(t1).T-del_a(t0).T )* R* b = scalar I will test drive the galgebra module

[sympy] intersection of Line and Line with Real points (floating point values)

2009-12-03 Thread sebastian
Hi, sympy.geometry.intersection doesn't work for me if the two lines are based on float values. For example, l1 Line(Point(Real('491.476998', prec=15), Real ('591.027601', prec=15)), Point(Real('656.726', prec=15), Real ('-2.46040501', prec=15))) l2

Re: [sympy] intersection of Line and Line with Real points (floating point values)

2009-12-03 Thread Ronan Lamy
Le jeudi 03 décembre 2009 à 03:52 -0800, sebastian a écrit : What am I doing wrong? Nothing, this is a known bug. See http://code.google.com/p/sympy/issues/detail?id=1649 for more information and a workaround. If you wish to contribute a fix, that would be most welcome. Ronan -- You received

[sympy] Re: printing exponential format vs float numbers

2009-12-03 Thread archeryguru2000
Ok, I think I have it. I have (hopefully) compiled a set of patches that implements what I've done. I have modified the following files: .../sympy/mpmath/mptypes.py .../sympy/mpmath/matrices.py .../sympy.mpmath/libmpf.py Basically I've just modified the patch that is found on the

Re: [sympy] Re: printing exponential format vs float numbers

2009-12-03 Thread Aaron S. Meurer
We use the Google Code issues page at http://code.google.com/p/sympy/issues/list. Create an issue and upload the patch to it. Aaron Meurer On Dec 3, 2009, at 1:06 PM, archeryguru2000 wrote: Ok, I think I have it. I have (hopefully) compiled a set of patches that implements what I've done.

[sympy] Using symbolic expressions for numeric calculations

2009-12-03 Thread twk
Hi, I am relatively new to sympy but like it and try to use it for my calculations. Now, I would like to use symbolic terms for further numeric calculations. I tried to use something like .subs(x,range(20)) but it did not work. Is there a way to treat symbols in a symbolic expression like a list

Re: [sympy] Using symbolic expressions for numeric calculations

2009-12-03 Thread Aaron S. Meurer
I also made this issue for this: http://code.google.com/p/sympy/issues/detail?id=1743 What do you think about a subs function? I just realized that would have to be a lambda for map anyway because it takes multiple arguments. So I could go either way on that. Aaron Meurer On Dec 3, 2009, at