Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Updates: Cc: matt...@gmail.com Labels: WrongResult Comment #20 on issue 2607 by asmeurer: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 Regarding the cancelation strategy from as_numer_denom4(): - First, it seems that my fears of non-cancelati

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #21 on issue 2607 by smi...@gmail.com: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 Your method 4 is the same as the one I tested out and didn't post results for because it had more significant failures. Besides the one you noted involving non-com

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #22 on issue 2607 by smi...@gmail.com: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 Would it make sense to pull out the gcd of the denoms? >>> gcd([x*y**2, x*y, 3*x*y]) x*y This slows things down for that step but there may be a speed up since

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #23 on issue 2607 by asmeurer: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 The expression with z is the one that causes he hyperexpand failure (I don't know if you got that). I didn't dig down to see what caused the Matrix error, but I suspect it

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #24 on issue 2607 by asmeurer: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 Would it make sense to pull out the gcd of the denoms? Well, as I said above, the only way to tell is to construct some expressions and benchmark it. But probably it wil

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #25 on issue 2607 by asmeurer: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 By the way, what's the best way to tell if an expression contains non-commutatives? Call is_commutative on the entire expression? -- You received this message because you

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #26 on issue 2607 by matt...@gmail.com: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 terms_gcd() also uses expand(). What you should look into is gcd_terms() from exprtools.py. This is what together() uses and for a few examples I tried, together(

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #27 on issue 2607 by matt...@gmail.com: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 Call is_commutative on the entire expression? Yes (up to bugs). -- You received this message because you are subscribed to the Google Groups "sympy-issues" group

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #28 on issue 2607 by asmeurer: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 for a few examples I tried, together() is only a few percent slower than as_numer_denom4() Not quite: In [4]: numers, denoms = zip(*((Symbol('n%d'%i),Symbol('d%d'%i)) for

Re: Issue 2624 in sympy: Sympy 0.7.1 can't integrate Gaussians

2011-08-05 Thread sympy
Comment #4 on issue 2624 by ness...@googlemail.com: Sympy 0.7.1 can't integrate Gaussians http://code.google.com/p/sympy/issues/detail?id=2624 As I already told Matthew, integrals code does not do any simplification. Perhaps this should be changed (maybe calling cancel if an integral was

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #29 on issue 2607 by matt...@gmail.com: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 I'm not sure what I'm missing here, but: In [1]: numers, denoms = zip(*((Symbol('n%d'%i),Symbol('d%d'%i)) for i in xrange(50))) In [2]: a = Add(*(n/d for n, d in

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #30 on issue 2607 by matt...@gmail.com: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 btw. I created a pull request fixing printing non-commutative objects (https://github.com/sympy/sympy/pull/530). -- You received this message because you are subs

Re: Issue 2618 in sympy: Solve fails on expressions containing finite symbols

2011-08-05 Thread sympy
Comment #12 on issue 2618 by mrock...@gmail.com: Solve fails on expressions containing finite symbols http://code.google.com/p/sympy/issues/detail?id=2618 So a temporary workaround is posted above which has check_assumptions check old assumptions as it checks the new assumptions. Intended

Re: Issue 2625 in sympy: Imaginary unit in R, ordering of complex numbers

2011-08-05 Thread sympy
Comment #3 on issue 2625 by alfr...@gmail.com: Imaginary unit in R, ordering of complex numbers http://code.google.com/p/sympy/issues/detail?id=2625 Actually, this is not even true. We should probably raise an error with that, just as Python does Lexicographic ordering is a natural choic

Re: Issue 2625 in sympy: Imaginary unit in R, ordering of complex numbers

2011-08-05 Thread sympy
Comment #4 on issue 2625 by matt...@gmail.com: Imaginary unit in R, ordering of complex numbers http://code.google.com/p/sympy/issues/detail?id=2625 Just to note this: we use the following ordering when sorting expressions (involving imaginary unit): In [1]: sorted([2 + 3*I, 2 + I, I, -I,

Re: Issue 1746 in sympy: Let bin/test --random accept a seed input

2011-08-05 Thread sympy
Comment #8 on issue 1746 by matt...@gmail.com: Let bin/test --random accept a seed input http://code.google.com/p/sympy/issues/detail?id=1746 btw. bin/test already supports --seed option, e.g.: $ bin/test --random --seed 68466956 == test process starts === executable: /usr/bin/pyth

Re: Issue 1994 in sympy: repr(Matrix) does not give re-executable string

2011-08-05 Thread sympy
Updates: Status: WontFix Owner: --- Comment #2 on issue 1994 by matt...@gmail.com: repr(Matrix) does not give re-executable string http://code.google.com/p/sympy/issues/detail?id=1994 http://groups.google.com/group/sympy/browse_thread/thread/1ec2a3a2c92bb0f5/101ae11736d3f6b2

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #31 on issue 2607 by smi...@gmail.com: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 If a Horner representation is acceptable (like as_numer_denom3) then using a binary approach to build the numer and denom is about 100X faster for the all-differen

Re: Issue 2618 in sympy: Solve fails on expressions containing finite symbols

2011-08-05 Thread sympy
Comment #13 on issue 2618 by nicolas@gmail.com: Solve fails on expressions containing finite symbols http://code.google.com/p/sympy/issues/detail?id=2618 "Intended behavior of Q.is_bounded(x) seems to be that it returns False on a standard Symbol if it has no information. This seems mat

Re: Issue 2618 in sympy: Solve fails on expressions containing finite symbols

2011-08-05 Thread sympy
Comment #14 on issue 2618 by mrock...@gmail.com: Solve fails on expressions containing finite symbols http://code.google.com/p/sympy/issues/detail?id=2618 In that case how about changing the False below to None? sympy/assumptions/handlers/calculus.py:AskBoundedHandler.Symbol(expr, assumptio

Re: Issue 1732 in sympy: integrate(sin(x).series(x)) gives traceback

2011-08-05 Thread sympy
Comment #3 on issue 1732 by matt...@gmail.com: integrate(sin(x).series(x)) gives traceback http://code.google.com/p/sympy/issues/detail?id=1732 Now it doesn't fail, it just return Integral(): In [1]: integrate(sin(x).series(x)) Out[1]: ⌠ ⎮ ⎛ 3 5 ⎞ ⎮ ⎜x x ⎟

Re: Issue 2532 in sympy: Non-commutative inverses should not print as fractions

2011-08-05 Thread sympy
Updates: Status: Started Labels: NeedsReview Comment #1 on issue 2532 by asmeurer: Non-commutative inverses should not print as fractions http://code.google.com/p/sympy/issues/detail?id=2532 https://github.com/sympy/sympy/pull/530 -- You received this message because you are

Re: Issue 2618 in sympy: Solve fails on expressions containing finite symbols

2011-08-05 Thread sympy
Comment #15 on issue 2618 by nicolas@gmail.com: Solve fails on expressions containing finite symbols http://code.google.com/p/sympy/issues/detail?id=2618 "Under what conditions should this return a False? I'm not very familiar with new assumptions." Neither am I. It's at least clear t

Re: Issue 2625 in sympy: Imaginary unit in R, ordering of complex numbers

2011-08-05 Thread sympy
Comment #5 on issue 2625 by asmeurer: Imaginary unit in R, ordering of complex numbers http://code.google.com/p/sympy/issues/detail?id=2625 You can define any number of total orderings on the set of complex numbers, but none of these will be algebraic orders. This is true even if you limit

Re: Issue 1746 in sympy: Let bin/test --random accept a seed input

2011-08-05 Thread sympy
Updates: Cc: ness...@googlemail.com Comment #9 on issue 1746 by asmeurer: Let bin/test --random accept a seed input http://code.google.com/p/sympy/issues/detail?id=1746 Does that seed actually affect the random ordering, or is it just for random tests? -- You received this message

Re: Issue 1994 in sympy: repr(Matrix) does not give re-executable string

2011-08-05 Thread sympy
Comment #3 on issue 1994 by asmeurer: repr(Matrix) does not give re-executable string http://code.google.com/p/sympy/issues/detail?id=1994 Actually, it's interesting that str(Matrix) doesn't give a 1-D representation like it does for all other SymPy objects. -- You received this message b

Re: Issue 2618 in sympy: Solve fails on expressions containing finite symbols

2011-08-05 Thread sympy
Comment #16 on issue 2618 by asmeurer: Solve fails on expressions containing finite symbols http://code.google.com/p/sympy/issues/detail?id=2618 The three way logic of the new and old assumptions are exactly the same. It returns True if it's known to be True in all cases, False if it's kno

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #32 on issue 2607 by asmeurer: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 (those benchmarks are done with cache off, so maybe this is the difference?). Ah, that's it. Here's what I get with the cache off: In [1]: numers, denoms = zip(*((Symbol

Re: Issue 2625 in sympy: Imaginary unit in R, ordering of complex numbers

2011-08-05 Thread sympy
Comment #6 on issue 2625 by alfr...@gmail.com: Imaginary unit in R, ordering of complex numbers http://code.google.com/p/sympy/issues/detail?id=2625 Yeah, you can't make the complex numbers into an ordered field. That's really only a problem if the comparison operator is required to be co

Issue 2626 in sympy: Piecewise should use a different syntax for "otherwise"

2011-08-05 Thread sympy
Status: Accepted Owner: asmeurer CC: andy.ter...@gmail.com Labels: Type-Defect Priority-High Milestone-Release0.7.2 New issue 2626 by asmeurer: Piecewise should use a different syntax for "otherwise" http://code.google.com/p/sympy/issues/detail?id=2626 This is something that's been bothering

Re: Issue 2625 in sympy: Imaginary unit in R, ordering of complex numbers

2011-08-05 Thread sympy
Comment #7 on issue 2625 by asmeurer: Imaginary unit in R, ordering of complex numbers http://code.google.com/p/sympy/issues/detail?id=2625 That's really only a problem if the comparison operator is required to be compatible with the field operations for other reasons. It is a problem. T

Re: Issue 2624 in sympy: Sympy 0.7.1 can't integrate Gaussians

2011-08-05 Thread sympy
Comment #5 on issue 2624 by asmeurer: Sympy 0.7.1 can't integrate Gaussians http://code.google.com/p/sympy/issues/detail?id=2624 I'm not sure actually. Calling simplify() or even cancel() can make integrate() slow. On the other hand, it seems like most of the time there is a very simple ans

Re: Issue 2625 in sympy: Imaginary unit in R, ordering of complex numbers

2011-08-05 Thread sympy
Comment #8 on issue 2625 by alfr...@gmail.com: Imaginary unit in R, ordering of complex numbers http://code.google.com/p/sympy/issues/detail?id=2625 That's the kind of thing I mean =). I am not really advocating either way. I was more pointing out that OP's suggestion is both reasonable and

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #33 on issue 2607 by smi...@gmail.com: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 The unevaluated Mul is a very nice way to go unless you think that `Mul(2,x, evaluate=False)/2` should not be `x`. If you don't like the final touching iwth `+0` o

Re: Issue 1459 in sympy: Summing symbolic elements of a list

2011-08-05 Thread sympy
Comment #2 on issue 1459 by kemellie...@gmail.com: Summing symbolic elements of a list http://code.google.com/p/sympy/issues/detail?id=1459 I tried to run the commands from Ondrej, but they did not work. Did the syntax change since then? Thanks! -- You received this message because you a

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #34 on issue 2607 by candleb...@gmail.com: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 Last thought...what about building up the multiplications from a binary lookup of precomputed muls? e.g. if you had 8 terms x1-x8 then computed x12=x1*x2, x

Re: Issue 1459 in sympy: Summing symbolic elements of a list

2011-08-05 Thread sympy
Comment #3 on issue 1459 by asmeurer: Summing symbolic elements of a list http://code.google.com/p/sympy/issues/detail?id=1459 Yes, sum() was renamed to summation(), as it was overriding the builtin function. -- You received this message because you are subscribed to the Google Groups "symp

Re: Issue 2626 in sympy: Piecewise should use a different syntax for "otherwise"

2011-08-05 Thread sympy
Comment #1 on issue 2626 by matt...@gmail.com: Piecewise should use a different syntax for "otherwise" http://code.google.com/p/sympy/issues/detail?id=2626 Why not to simply allow Piecewise((expr_i, cond_i), ..., otherwise_expr)? The last tuple seems redundant to me. We could also introduce

Re: Issue 1746 in sympy: Let bin/test --random accept a seed input

2011-08-05 Thread sympy
Comment #10 on issue 1746 by matt...@gmail.com: Let bin/test --random accept a seed input http://code.google.com/p/sympy/issues/detail?id=1746 It does affect the order: $ bin/test --random test process starts executable: /usr/bin/python (2.6.6-final-0) architecture: 64-bit cache

Re: Issue 2624 in sympy: Sympy 0.7.1 can't integrate Gaussians

2011-08-05 Thread sympy
Comment #6 on issue 2624 by matt...@gmail.com: Sympy 0.7.1 can't integrate Gaussians http://code.google.com/p/sympy/issues/detail?id=2624 We saw during the conference that the integrator can return very complicated results and conditions for very simple cases. I would try to use simplify(

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #35 on issue 2607 by asmeurer: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 Does your as_numer_denom5() satisfy the requirements of together() that Mateusz mentioned in comment 29? Last thought...what about building up the multiplications from a b

Re: Issue 2626 in sympy: Piecewise should use a different syntax for "otherwise"

2011-08-05 Thread sympy
Comment #2 on issue 2626 by asmeurer: Piecewise should use a different syntax for "otherwise" http://code.google.com/p/sympy/issues/detail?id=2626 Would it ever make sense for Piecewise to hold Tuples? If so, then we have to wrap the last condition around a tuple. Otherwise, there is no w

Re: Issue 1746 in sympy: Let bin/test --random accept a seed input

2011-08-05 Thread sympy
Updates: Status: Fixed Comment #11 on issue 1746 by asmeurer: Let bin/test --random accept a seed input http://code.google.com/p/sympy/issues/detail?id=1746 Oh, I didn't realize that. This can be closed then. -- You received this message because you are subscribed to the Google Grou

Re: Issue 2624 in sympy: Sympy 0.7.1 can't integrate Gaussians

2011-08-05 Thread sympy
Comment #7 on issue 2624 by asmeurer: Sympy 0.7.1 can't integrate Gaussians http://code.google.com/p/sympy/issues/detail?id=2624 1) Is there a chance to merge this branch before August 23rd? If you help review, it can. :) Already, there is a branch for review, which has some integration code

Re: Issue 2624 in sympy: Sympy 0.7.1 can't integrate Gaussians

2011-08-05 Thread sympy
Comment #8 on issue 2624 by matt...@gmail.com: Sympy 0.7.1 can't integrate Gaussians http://code.google.com/p/sympy/issues/detail?id=2624 I can help review it. Yes, it's EuroSciPy. -- You received this message because you are subscribed to the Google Groups "sympy-issues" group. To post to

Re: Issue 2624 in sympy: Sympy 0.7.1 can't integrate Gaussians

2011-08-05 Thread sympy
Comment #9 on issue 2624 by asmeurer: Sympy 0.7.1 can't integrate Gaussians http://code.google.com/p/sympy/issues/detail?id=2624 By the way, that would indeed be very cool if we could do the integral from that paper. The author wants to know how Mathematica does the integral. But with Tom'

Re: Issue 516 in sympy: fix sys.path mess

2011-08-05 Thread sympy
Updates: Status: Fixed Comment #9 on issue 516 by matt...@gmail.com: fix sys.path mess http://code.google.com/p/sympy/issues/detail?id=516 Pyglet isn't a part of SymPy anymore, thus this issue isn't relevant anymore (I mark this a Fixed although more precise status would be Irrelevan

Re: Issue 925 in sympy: sys.path hacks, how to solve it

2011-08-05 Thread sympy
Updates: Status: Fixed Comment #22 on issue 925 by matt...@gmail.com: sys.path hacks, how to solve it http://code.google.com/p/sympy/issues/detail?id=925 Pyglet isn't a part of SymPy anymore, thus this issue isn't relevant anymore (I mark this a Fixed although more precise status wo

Re: Issue 1946 in sympy: Recursion error with SYMPY_GROUND_TYPES=sympy (caching problem)

2011-08-05 Thread sympy
Updates: Status: Fixed Comment #19 on issue 1946 by matt...@gmail.com: Recursion error with SYMPY_GROUND_TYPES=sympy (caching problem) http://code.google.com/p/sympy/issues/detail?id=1946 (No comment was entered for this change.) -- You received this message because you are subscribe

Re: Issue 360 in sympy: Arithmetic with complex infinity

2011-08-05 Thread sympy
Updates: Status: Fixed Comment #10 on issue 360 by matt...@gmail.com: Arithmetic with complex infinity http://code.google.com/p/sympy/issues/detail?id=360 (No comment was entered for this change.) -- You received this message because you are subscribed to the Google Groups "sympy-is

Re: Issue 825 in sympy: SYMPY_USE_CACHE=no py.test sympy/ results

2011-08-05 Thread sympy
Updates: Status: Duplicate Mergedinto: 2495 Comment #5 on issue 825 by matt...@gmail.com: SYMPY_USE_CACHE=no py.test sympy/ results http://code.google.com/p/sympy/issues/detail?id=825 (No comment was entered for this change.) -- You received this message because you are subsc

Re: Issue 2495 in sympy: Test failures with the cache turned off

2011-08-05 Thread sympy
Comment #5 on issue 2495 by matt...@gmail.com: Test failures with the cache turned off http://code.google.com/p/sympy/issues/detail?id=2495 Issue 825 has been merged into this issue. -- You received this message because you are subscribed to the Google Groups "sympy-issues" group. To post t

Re: Issue 2104 in sympy: canonical ordering of terms

2011-08-05 Thread sympy
Updates: Status: Fixed Owner: --- Cc: -matt...@gmail.com Comment #15 on issue 2104 by matt...@gmail.com: canonical ordering of terms http://code.google.com/p/sympy/issues/detail?id=2104 This was fixed in 0.7.0. -- You received this message because you are subscribed to t

Re: Issue 2607 in sympy: as_numer_denom() is too slow

2011-08-05 Thread sympy
Comment #36 on issue 2607 by smi...@gmail.com: as_numer_denom() is too slow http://code.google.com/p/sympy/issues/detail?id=2607 Does your as_numer_denom5() satisfy the requirements of together() that Mateusz mentioned in comment 29? No, it's just a faster version of the existing as_numer_de