[sympy] Sympy Integral

2014-11-30 Thread Paul Royik
Is there any way to determine whether integral can be calculated in terms of elementary functions? I.e. without erf, airy, etc. -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an

[sympy] Custom ordering in Sympy

2014-11-30 Thread Paul Royik
I want to rewrite Integral of 10x^7-xy+y/x+2x^5 as sum of integrals. Something like Int(10x^7-xy+y/x+2x^5) = Int(10x^7)-Int(xy)+Int(y/x)+Int(2x^5) But sympy makes its own ordering and outputs Int(10x^7-xy+y/x+2x^5) = Int(2x^5)+Int(10x^7)+Int(y/x)-Int(xy) Where can I override this behaviour

[sympy] Sympy assumptions

2014-12-04 Thread Paul Royik
Suppose that I have an expression e = x**2+sqrt(a**2*b**2) What is the best way to assign to expression assumption that all variables are Positive, so that equivalent of e is x**2+ab? -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from

Re: [sympy] Sympy assumptions

2014-12-04 Thread Paul Royik
I know this way, but is there any way to change this after symbol creation? On Thursday, December 4, 2014 6:29:27 PM UTC+2, Ondřej Čertík wrote: On Thu, Dec 4, 2014 at 7:38 AM, Paul Royik distan...@gmail.com javascript: wrote: Suppose that I have an expression e = x**2+sqrt(a**2*b

[sympy] Factor options

2014-12-07 Thread Paul Royik
How should I use factor to factor expression over irrational numbers? For example, x^2-4 produces (x-2)(x+2) x^2-2 produces (x-sqrt(2))(x+sqrt(2)) x^4+1 produces (x^2-sqrt(2) x+1) (x^2+sqrt(2) x+1) x^2+1 produces x^2+1 (only complex roots) x^4-9 x^2-22 produces (x^2+2)(x-sqrt(11))(x+sqrt(11)) --

[sympy] Re: Factor options

2014-12-09 Thread Paul Royik
Exactly. Thanks. -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to this group, send email to sympy@googlegroups.com.

[sympy] Subs behavior

2014-12-09 Thread Paul Royik
from sympy import * a = symbols('a',positive=True) (1/(sec(a)**2-1)**(3/2)).subs(sec(a)**2-1,tan(a)**2) last statement should output something close to 1/tan^3(a), or abs(tan(a))/tan^4(a) But instead it outputs something unrelated: abs(tan(a))/tan^2(a) Is this a bug or I'm missing something?

[sympy] Re: Subs behavior

2014-12-09 Thread Paul Royik
or Rational(3, 2) On Tuesday, December 9, 2014 6:22:00 AM UTC-6, Paul Royik wrote: from sympy import * a = symbols('a',positive=True) (1/(sec(a)**2-1)**(3/2)).subs(sec(a)**2-1,tan(a)**2) last statement should output something close to 1/tan^3(a), or abs(tan(a))/tan^4(a) But instead

[sympy] How find works?

2014-12-09 Thread Paul Royik
Can somebody explain how find works. For example, a = sympy.Wild('a', exclude=[x]) b = sympy.Wild('b', exclude=[x]) c = sympy.Wild('c', exclude=[x]) pattern = sympy.sqrt(a*x**2+b*x+c) 1/((x^2+a^2)+sqrt(a^2-x^2)).find(pattern) # - set([sqrt(a^2-x^2)]) 1/((x^2+a^2)*sqrt(a^2-x^2)).find(pattern) # -

[sympy] Sympy rewriting

2014-12-09 Thread Paul Royik
What function should I use to rewrite 2/(a*(u**2 + 1) - b*(u**2 - 1)) into 2/((a-b)*u**2+(a+b))? Thanks. -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to

Re: [sympy] Re: Subs behavior

2014-12-13 Thread Paul Royik
live.sympy.org gives Abs(tan(a))/tan(a)**2 -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to this group, send email to

[sympy] How to get the last version of symy on windows

2014-12-24 Thread Paul Royik
I need the freshest version of sympy. Developers told me that they solved issue with inverse laplace trhansform (1/(s-2)+3/s^7 raised an exception in 0.7.6) How can I install it from git? -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe

Re: [sympy] How to get the last version of symy on windows

2014-12-24 Thread Paul Royik
it). Aaron Meurer On Wed, Dec 24, 2014 at 5:55 AM, Paul Royik distan...@gmail.com javascript: wrote: I need the freshest version of sympy. Developers told me that they solved issue with inverse laplace trhansform (1/(s-2)+3/s^7 raised an exception in 0.7.6) How can I install

[sympy] Sympy multiprocessing

2015-02-18 Thread Paul Royik
I'm using sympy in my web-based project. Specifically, I need to calculate integral. However, I need to limit time and interrupt function. Code, that I use. import multiprocessing class RunableProcessing(multiprocessing.Process): def __init__(self, func, *args, **kwargs): self.queue

[sympy] Bug in parse expression

2015-02-18 Thread Paul Royik
Just noticed, that official example works: from sympy.parsing.sympy_parser import (parse_expr, standard_transformations, implicit_multiplication_application) parse_expr(10sin**2 x**2 + 3xyz + tan theta,transformations=(standard_transformations +(implicit_multiplication_application,))) But

[sympy] Stop long-running function

2015-02-20 Thread Paul Royik
Is there any way to stop long running function? Simplify sometimes runs too long. I want to terminate it after 30 seconds. Thank you. -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it,

Re: [sympy] Stop long-running function

2015-02-20 Thread Paul Royik
I mean in program, not in command window. On Saturday, February 21, 2015 at 1:08:10 AM UTC+2, Jason Moore wrote: ctrl-C will terminate a running command. Jason moorepants.info +01 530-601-9791 On Fri, Feb 20, 2015 at 2:53 PM, Paul Royik distan...@gmail.com javascript: wrote

Re: [sympy] Concurrency in Sympy

2015-03-12 Thread Paul Royik
Can you explain more deeply concerning integration? How many people can access integration simultaneously? On Thursday, March 12, 2015 at 2:25:53 PM UTC+2, Joachim Durchholz wrote: Am 12.03.2015 um 13:17 schrieb Paul Royik: How many concurrent instances can Sympy handle? I.e, how many

Re: [sympy] Concurrency in Sympy

2015-03-12 Thread Paul Royik
schrieb Paul Royik: Can you explain more deeply concerning integration? How many people can access integration simultaneously? Can you explain what scenario you mean? I.e. how many people accessing SymPy via what means on what machine(s)? I suspect a misunderstanding somewhere, but I

Re: [sympy] Concurrency in Sympy

2015-03-12 Thread Paul Royik
Thank you will look. On Thursday, March 12, 2015 at 10:50:28 PM UTC+2, Joachim Durchholz wrote: Am 12.03.2015 um 21:48 schrieb Paul Royik: That's what I asked. How to configure sympy? SymPy does not manage multiple processes at all, it's a single-threaded application. You'd have

Re: [sympy] Concurrency in Sympy

2015-03-12 Thread Paul Royik
That's what I asked. How to configure sympy? On Thursday, March 12, 2015 at 10:30:56 PM UTC+2, Joachim Durchholz wrote: Am 12.03.2015 um 21:14 schrieb Paul Royik: But if integral is easy (take not much memory) is number of people is restricted only by available memory? Can Sympy serve

Re: [sympy] Concurrency in Sympy

2015-03-12 Thread Paul Royik
But if integral is easy (take not much memory) is number of people is restricted only by available memory? Can Sympy serve as much as memory allows? On Thursday, March 12, 2015 at 8:44:44 PM UTC+2, Joachim Durchholz wrote: Am 12.03.2015 um 17:39 schrieb Paul Royik: I'm running my site

[sympy] Concurrency in Sympy

2015-03-12 Thread Paul Royik
How many concurrent instances can Sympy handle? I.e, how many people can simultenously access diff function? -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to

[sympy] sympy.Eq, equality versus is

2015-03-24 Thread Paul Royik
Maybe this is property of Python, but I'm not aware about it. expr = 1/(2*sin(x/2 + pi/4)*cos(x/2 + pi/4)) equation = Eq(expr,expr) equation == True # outputs True equation is True # outputs False Why True equation is not True? Moreover, in IDE in which I'm working, when I write equation !=

[sympy] Difference between expr.is_number and expr.is_Number

2015-03-25 Thread Paul Royik
What is the difference between above two properties? -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to this group, send

Re: [sympy] Difference between expr.is_number and expr.is_Number

2015-03-25 Thread Paul Royik
, 2015 at 7:32 AM, Paul Royik distan...@gmail.com javascript: wrote: What is the difference between above two properties? -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from

[sympy] Power simplification

2015-03-31 Thread Paul Royik
Is there any way I can rewrite x^2y^4 to (xy^2)^2 or x^3(x+y)^6 to (x(x+y)^2)^3 , i.e. combine base? Exponents are guaranteed to be integers. Thank you. -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop

[sympy] Re: Keeping symbolic expressions raw

2015-03-31 Thread Paul Royik
No, there is no way. You can use evaluate=False, when creating expressions, but it will help partly. For example, Add(4*x*y,-6*x*y, evaluate=False) will output -6xy+4xy (notice order changed), but Add(4*x*y,-6*x*y) will output -2xy On Tuesday, March 31, 2015 at 5:40:31 PM UTC+3, Milos Bogataj

[sympy] Difference between fraction and as_numer_denom

2015-03-26 Thread Paul Royik
I used as_numer_denom earlier and now found fraction. Is there any difference between them? -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to

Re: [sympy] Difference between fraction and as_numer_denom

2015-03-27 Thread Paul Royik
of fractions (e.g., fraction(1/x + 1/y) gives (1/x + 1/y, 1), whereas (1/x + 1/y).as_numer_denom() gives (x + y, x*y)). Aaron Meurer On Thu, Mar 26, 2015 at 12:09 PM, Paul Royik distan...@gmail.com javascript: wrote: I used as_numer_denom earlier and now found fraction. Is there any

Re: [sympy] Re: Combine sqrt?

2015-02-27 Thread Paul Royik
, 2015 at 6:20 PM, Paul Royik distan...@gmail.com javascript: wrote: Thank you. Got it! On Friday, February 27, 2015 at 1:52:32 AM UTC+2, Ondřej Čertík wrote: On Thu, Feb 26, 2015 at 4:05 PM, Paul Royik distan...@gmail.com wrote: It doesn't work if x is positive

Re: [sympy] Re: Combine sqrt?

2015-02-26 Thread Paul Royik
Thank you. Got it! On Friday, February 27, 2015 at 1:52:32 AM UTC+2, Ondřej Čertík wrote: On Thu, Feb 26, 2015 at 4:05 PM, Paul Royik distan...@gmail.com javascript: wrote: It doesn't work if x is positive. Indeed, looks like a bug. As a workaround, you can always substitute general x

[sympy] Issue with positive variable assumption

2015-02-25 Thread Paul Royik
Hello. I recently needed to make all variables in expression positive. I do this like following: f = log(x) a=Wild('a') symbols = f.free_symbols for symbol in symbols: f = f.subs(symbol, dummy) f = f.subs(symbol, sympy.Symbol(symbol.name, positive=True)) Now, when I do f.replace(log(a),

[sympy] Re: Issue with positive variable assumption

2015-02-25 Thread Paul Royik
It appeared, that problem is with cache. Clearing cache solves the problem. Now, I want to ask, how cache is working in multithreaded environment? Do all people share same cache? On Wednesday, February 25, 2015 at 4:11:03 PM UTC+2, Paul Royik wrote: Hello. I recently needed to make all

Re: [sympy] Re: Issue with positive variable assumption

2015-02-25 Thread Paul Royik
symbols = f.free_symbols for symbol in symbols: f = f.subs(symbol, sympy.Symbol(symbol.name)) print f.replace(log(a),log(Abs(a))) # OOPS! prints log(x) instead of log(|x|) On Wednesday, February 25, 2015 at 6:38:12 PM UTC+2, Joachim Durchholz wrote: Am 25.02.2015 um 17:04 schrieb Paul Royik

Re: [sympy] Re: Issue with positive variable assumption

2015-02-25 Thread Paul Royik
a new cache)? And if so, are you using fastcache? Aaron Meurer On Wed, Feb 25, 2015 at 10:04 AM, Paul Royik distan...@gmail.com javascript: wrote: It appeared, that problem is with cache. Clearing cache solves the problem. Now, I want to ask, how cache is working in multithreaded

Re: [sympy] Re: Issue with positive variable assumption

2015-02-25 Thread Paul Royik
25.02.2015 um 18:10 schrieb Paul Royik: I mean in multithreaded environment, do people share same cache (and same symbols) If person1 created x=Symbol('x') and this is cached, and person2 created x=Symbol('x', positive=True), do they have separate x's or they mixed, so person1 can get

[sympy] Combine sqrt?

2015-02-26 Thread Paul Royik
What is the best way to convert sqrt(x)*sqrt(x-4) to sqrt(x^2-4x) or (x^2+5x+4)/sqrt(x)/sqrt(x-4) to (x^2+5x+4)/sqrt(x^2-4x) I tried replace, but it doesn't work in second case. -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this

Re: [sympy] Combine sqrt?

2015-02-26 Thread Paul Royik
I have assumption that x is positive. a, b = Wild... I use like f.replace(sqrt(a)*sqrt(b), sqrt((a*b).expand())) On Thursday, February 26, 2015 at 7:51:17 PM UTC+2, Ondřej Čertík wrote: Hi Paul, On Thu, Feb 26, 2015 at 8:32 AM, Paul Royik distan...@gmail.com javascript: wrote: What

[sympy] Re: Combine sqrt?

2015-02-26 Thread Paul Royik
I need something more general On Thursday, February 26, 2015 at 7:50:15 PM UTC+2, John Peterson wrote: On Thursday, February 26, 2015 at 8:32:11 AM UTC-7, Paul Royik wrote: What is the best way to convert sqrt(x)*sqrt(x-4) to sqrt(x^2-4x) or (x^2+5x+4)/sqrt(x)/sqrt(x-4) to (x^2+5x+4)/sqrt

Re: [sympy] Stop long-running function

2015-02-21 Thread Paul Royik
...@gmail.com javascript: wrote: Maybe this: http://stackoverflow.com/questions/366682/how-to-limit-execution-time-of-a-function-call-in-python Jason moorepants.info +01 530-601-9791 On Fri, Feb 20, 2015 at 3:12 PM, Paul Royik distan...@gmail.com javascript

Re: [sympy] Stop long-running function

2015-02-21 Thread Paul Royik
I see. On Saturday, February 21, 2015 at 4:45:49 PM UTC+2, Joachim Durchholz wrote: Am 21.02.2015 um 11:32 schrieb Paul Royik: So, this feature is closed? We'd really like to have that, we're just not sure how to implement it. -- You received this message because you are subscribed

Re: [sympy] Re: Combine sqrt?

2015-02-26 Thread Paul Royik
=True) will do what you want. Aaron Meurer On Thu, Feb 26, 2015 at 12:18 PM, Paul Royik distan...@gmail.com javascript: wrote: I need something more general On Thursday, February 26, 2015 at 7:50:15 PM UTC+2, John Peterson wrote: On Thursday, February 26, 2015 at 8:32:11 AM

[sympy] Re: expanding powers of cos()

2015-03-27 Thread Paul Royik
Hello, Jeremy. Try this: from sympy.simplify.fu import TR7 TR7(cos(x)**2) On Friday, March 27, 2015 at 8:36:43 PM UTC+2, Jeremy Holleman wrote: Hi, I'm relatively new to SymPy, and enjoying it so far. I'm trying to analyze the harmonics that result when a sinusoid goes through a

Re: [sympy] Restore Poly from dict

2015-03-27 Thread Paul Royik
*y**2*x, y, x, domain='ZZ') The order of the generators is the same as the order in the Poly. I don't remember if there is a better method on Poly for reordering generators, but in the above example Poly(p, x, y) will reorder them. Aaron Meurer On Fri, Mar 27, 2015 at 2:50 PM, Paul

[sympy] Re: expanding powers of cos()

2015-03-27 Thread Paul Royik
What do you mean how I would do that in a programmatic way? On Friday, March 27, 2015 at 11:33:27 PM UTC+2, Jeremy Holleman wrote: Update - Some combination of TR7, TR8, and expand() seems to work for powers that I've tried. For example: In [5]: TR8(expand(TR8(cos(x)**7))) Out[5]:

Re: [sympy] Difference between fraction and as_numer_denom

2015-03-27 Thread Paul Royik
Meurer On Fri, Mar 27, 2015 at 1:47 PM, Paul Royik distan...@gmail.com javascript: wrote: This dosn't seem reasonable. Expanding goes on the line Mul(*numer). Mul(*numer, evaluate=False) fixes things. On Friday, March 27, 2015 at 4:31:58 PM UTC+2, Chris Smith wrote

Re: [sympy] Re: expanding powers of cos()

2015-03-28 Thread Paul Royik
It seems like TR7 reduces only cos(x)**2 On Saturday, March 28, 2015 at 3:16:53 AM UTC+2, Jeremy Holleman wrote: I'll do that. Thanks to both of you for the help! On Friday, March 27, 2015 at 8:30:03 PM UTC-4, Aaron Meurer wrote: I would take a look at the fu paper for more references on

[sympy] Re: Smart subs

2015-04-02 Thread Paul Royik
-sqrt(x+1), x)[0]) Out[2]: 2 ╱ 2 y + ╲╱ y + 2 If *y* is declared to be positive, you get the right expression. On Wednesday, April 1, 2015 at 7:01:29 PM UTC+2, Paul Royik wrote: Is there a function that will convert sqrt(x+1)+x+3 into y^2+y+2 with y=sqrt

[sympy] Re: Smart subs

2015-04-03 Thread Paul Royik
OK. Thank you. On Friday, April 3, 2015 at 4:06:29 AM UTC+3, Francesco Bonazzi wrote: Maybe you should read the documentation. Have a look at unevaluated expressions. -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group

Re: [sympy] What is the fastest way to calculate CDF of T-distribution?

2015-04-23 Thread Paul Royik
Thanks. On Thursday, April 23, 2015 at 8:04:46 PM UTC+3, Matthew Brett wrote: Hi, On Wed, Apr 22, 2015 at 11:51 PM, Paul Royik distan...@gmail.com javascript: wrote: I'm looking for numerical answer. No need to find analytical solution. On Wednesday, April 22, 2015 at 6:38:29 PM

Re: [sympy] What is the fastest way to calculate CDF of T-distribution?

2015-04-23 Thread Paul Royik
. Aaron Meurer On Wed, Apr 22, 2015 at 6:17 AM, Paul Royik distan...@gmail.com javascript: wrote: I need to calculate P(X1.2) where X follows T distribution with 5 degrees of freedom. cdf(StudentT('x',5))(1.2) takes too much time. Is there another way? -- You

[sympy] What is the fastest way to calculate CDF of T-distribution?

2015-04-22 Thread Paul Royik
I need to calculate P(X1.2) where X follows T distribution with 5 degrees of freedom. cdf(StudentT('x',5))(1.2) takes too much time. Is there another way? -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop

Re: [sympy] Sympy solve on interval

2015-04-28 Thread Paul Royik
How? On Tuesday, April 28, 2015 at 2:04:07 PM UTC+3, Joachim Durchholz wrote: Am 28.04.2015 um 12:42 schrieb Paul Royik: How can I solve equation on the interval? For example sin(x)=0, 2pi=x=4pi Apply an assumption to x. -- You received this message because you are subscribed

Re: [sympy] Sympy solve on interval

2015-04-28 Thread Paul Royik
site:docs.sympy.org Am 28.04.2015 um 13:08 schrieb Paul Royik: How? On Tuesday, April 28, 2015 at 2:04:07 PM UTC+3, Joachim Durchholz wrote: Am 28.04.2015 um 12:42 schrieb Paul Royik: How can I solve equation on the interval? For example sin(x)=0, 2pi=x=4pi Apply an assumption to x

[sympy] Sympy solve on interval

2015-04-28 Thread Paul Royik
How can I solve equation on the interval? For example sin(x)=0, 2pi=x=4pi -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to

[sympy] How evalf chop works?

2015-05-05 Thread Paul Royik
expr = 2682469501914725726236807146020965585039500052171867725591732162347*exp(-60)/48648670709168952944186714683693506574127 N(expr) outputs 0.482830727370627 but N(expr, chop=True) outputs 0 In some other expressions chop doesn't make 0. So, how chop works? -- You received this message

[sympy] Sympy limit - result depends on the sign

2015-04-09 Thread Paul Royik
Hello. I'm trying to calulate limit as x approaches infinity of e^((t-1)x) Naturally, I got Notimplemented error Result depends on the sign of sign(t - 1). My question how can I in general make assumption that this limit is finite? I mean in general, because expression for limit can be

Re: [sympy] factor with sqrt

2015-04-06 Thread Paul Royik
Thank you for your response. I made a mistake, wanted (sqrt(x)+2)(sqrt(x)-1). But, in general, I've got the idea. On Sunday, April 5, 2015 at 7:45:51 AM UTC+3, Ondřej Čertík wrote: On Sat, Apr 4, 2015 at 7:15 AM, Paul Royik distan...@gmail.com javascript: wrote: Hwow can I make sympy

[sympy] factor with sqrt

2015-04-04 Thread Paul Royik
Hwow can I make sympy factor x+sqrt(x)-2 into (sqrt(x)+1)(sqrt(x)-2) ? -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to

[sympy] Smart subs

2015-04-01 Thread Paul Royik
Is there a function that will convert sqrt(x+1)+x+3 into y^2+y+2 with y=sqrt(x+1)? Subs is smart enough to make (x**4+x**2+2).subs(x**2, y) into y**2+y+2. But it doesn't work with initial example. Even sqrt(x+1)+x+1 is not converted into y^2+y -- You received this message because you are

[sympy] Re: Smart subs

2015-04-02 Thread Paul Royik
Your aaproach fails when we want to substitute sqrt(xy+1) for u. What to solve here? On Thursday, April 2, 2015 at 4:17:54 PM UTC+3, Paul Royik wrote: Another smart solution. Thanks a ton. Maybe you can tell me why the following code doesn't work: https://github.com/sympy/sympy/issues/9233

[sympy] Re: Smart subs

2015-04-02 Thread Paul Royik
[2]: 2 ╱ 2 y + ╲╱ y + 2 If *y* is declared to be positive, you get the right expression. On Wednesday, April 1, 2015 at 7:01:29 PM UTC+2, Paul Royik wrote: Is there a function that will convert sqrt(x+1)+x+3 into y^2+y+2 with y=sqrt(x+1)? Subs is smart

[sympy] Re: Power simplification

2015-04-02 Thread Paul Royik
]: Pow(mt[w], 2, evaluate=False) Out[6]: 2 ⎛ 2⎞ ⎝x⋅y ⎠ On Tuesday, March 31, 2015 at 4:50:59 PM UTC+2, Paul Royik wrote: Is there any way I can rewrite x^2y^4 to (xy^2)^2 or x^3(x+y)^6 to (x(x+y)^2)^3 , i.e. combine base? Exponents are guaranteed to be integers

[sympy] Sympy zoo

2015-06-07 Thread Paul Royik
Why 1/0 is complex infinity and log(0) is complex infinity? I also found a bug with oo**zoo. It is recursion error. -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to

[sympy] Sympy exponential equation

2015-06-05 Thread Paul Royik
Hello. Why solve(-(21/20)**(12*x) + 120, x) gives empty set? -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to this group,

Re: [sympy] Sympy exponential equation

2015-06-05 Thread Paul Royik
/sympy/sympy/issues/2999. Aaron Meurer On Fri, Jun 5, 2015 at 10:56 AM, Paul Royik distan...@gmail.com javascript: wrote: Hello. Why solve(-(21/20)**(12*x) + 120, x) gives empty set? -- You received this message because you are subscribed to the Google Groups sympy group

Re: [sympy] Sympy exponential equation

2015-06-06 Thread Paul Royik
. I think the solution should be log(120)/log(S(21)/20)/12. This is related to this issue https://github.com/sympy/sympy/issues/2999. Aaron Meurer On Fri, Jun 5, 2015 at 10:56 AM, Paul Royik distan...@gmail.com javascript: wrote: Hello. Why solve(-(21/20)**(12*x) + 120, x

[sympy] Re: Sympy zoo

2015-06-08 Thread Paul Royik
Thank you. How to make it work in real field? On Monday, June 8, 2015 at 12:05:24 PM UTC+3, Kalevi Suominen wrote: On Sunday, June 7, 2015 at 2:52:52 PM UTC+3, Paul Royik wrote: Why 1/0 is complex infinity and log(0) is complex infinity? They are shorthand notations for the limits of 1/z

Re: [sympy] Sympy CPU usage

2015-06-02 Thread Paul Royik
at 11:53 AM, Paul Royik distan...@gmail.com javascript: wrote: Hello. I noticed that sympy uses much CPU. Is there any way I can control its usage or control the code where it might occur. For now my code uses match extensively, but I still don't know what functions use much CPU

Re: [sympy] Sympy CPU usage

2015-06-02 Thread Paul Royik
to them. Aaron Meurer On Tue, Jun 2, 2015 at 11:53 AM, Paul Royik distan...@gmail.com javascript: wrote: Hello. I noticed that sympy uses much CPU. Is there any way I can control its usage or control the code where it might occur. For now my code uses match extensively

[sympy] 1/0

2015-06-09 Thread Paul Royik
Is it possible to make sympy to return oo instead of zoo for expressions like 1/0 and tan(pi/2). What I need is real field. -1/0 is -oo, not zoo. -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails

[sympy] System Error

2015-06-09 Thread Paul Royik
SystemError: Parent module 'sympy.core' not loaded, cannot perform relative import When I try to run the following: expr=sqrt(pi)*(sqrt(2)*(1 - I)*erf(sqrt(-I)) + sqrt(2)*(-1 + I)*erf(sqrt(2)*(-1/2 + I/2)) + sqrt(2)*(1 + I)*erf((-1)**(1/4)) + sqrt(2)*(1 + I)*erf(sqrt(2)*(1/2 + I/2)))/16

Re: [sympy] 1/0

2015-06-09 Thread Paul Royik
at infinities than using the infinites directly, as at best an expression simplifying itself at infinities naively will lead to nan when there is finite limit, and at worst it will lead to a wrong result. Aaron Meurer On Tue, Jun 9, 2015 at 9:23 AM, Paul Royik distan...@gmail.com javascript

[sympy] Bug with limit

2015-06-22 Thread Paul Royik
limit(sin(x)**15,x,0,'+') works limit(sin(x)**15,x,0,'-') hangs -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to this

[sympy] Powsimp

2015-06-22 Thread Paul Royik
Can't rewrite sin(x)**17/x**17 as (sin(x)/x)**17. Is there any way to do that? -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To

[sympy] Does sympy has two-sided limit?

2015-06-11 Thread Paul Royik
Why limit(1/x,x,0) gives me infinity? I digged in the code and saw that default direction is right. But how to make sympy to compute two-sided limit? limit(1/x,x,0) should return something like nana, because one-sided limits are not equal. -- You received this message because you are

[sympy] Extracting highest power and dividing

2015-06-17 Thread Paul Royik
I'm curently try to implement a little tutorial concerning limits. On this stage I'm trying to explain inifnite limits, solved by factoring. For example, lim x-oo (sqrt(x^2+x))/(2x+5)=(xsqrt(1+1/x))/(x(2+5/x))=sqrt(1+0)/(2+5)=1/2 So, basically, I need to extract highest power from expression

Re: [sympy] Re: How evalf chop works?

2015-06-02 Thread Paul Royik
if you specify chop smaller than that you will not get 0 for the result: expr.n(chop=1e-27) 0.482830727370613 expr.n(chop=1e-26) 0 On Thursday, May 21, 2015 at 1:49:34 PM UTC-5, Paul Royik wrote: OK. Thank you. Could you move it to issues? On Thursday, May 21, 2015 at 6:45:41 PM

[sympy] Sympy CPU usage

2015-06-02 Thread Paul Royik
Hello. I noticed that sympy uses much CPU. Is there any way I can control its usage or control the code where it might occur. For now my code uses match extensively, but I still don't know what functions use much CPU: simplify, integrate, match etc. Also, is clear_cache impacts CPU usage? Thank

[sympy] Re: Is it possible to calculate inverse Laplace Transform of exponent?

2015-07-03 Thread Paul Royik
OK. Thanks. On Friday, July 3, 2015 at 2:47:19 PM UTC+3, Kalevi Suominen wrote: On Friday, July 3, 2015 at 1:58:06 PM UTC+3, Paul Royik wrote: I want to calculate inverse laplace of e^(-2s) which is dirac(t-2), but sympy gives unevaluated. It is currently not possible to compute

[sympy] Is it possible to calculate inverse Laplace Transform of exponent?

2015-07-03 Thread Paul Royik
I want to calculate inverse laplace of e^(-2s) which is dirac(t-2), but sympy gives unevaluated. -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving emails from it, send an email to

Re: [sympy] Re: How evalf chop works?

2015-05-21 Thread Paul Royik
at 8:53 AM, Paul Royik distan...@gmail.com javascript: wrote: Nobody could tell how chop works? On Tuesday, May 5, 2015 at 11:41:24 AM UTC+3, Paul Royik wrote: expr = 2682469501914725726236807146020965585039500052171867725591732162347*exp(-60

[sympy] Re: How evalf chop works?

2015-05-21 Thread Paul Royik
Nobody could tell how chop works? On Tuesday, May 5, 2015 at 11:41:24 AM UTC+3, Paul Royik wrote: expr = 2682469501914725726236807146020965585039500052171867725591732162347*exp(-60)/48648670709168952944186714683693506574127 N(expr) outputs 0.482830727370627 but N(expr, chop=True

Re: [sympy] mod_wsgi with sympy

2015-08-03 Thread Paul Royik
Cache is cleared every request. On Saturday, August 1, 2015 at 2:11:12 PM UTC+3, Denis Akhiyarov wrote: Have you checked caching in sympy? I have to force clean the sympy cache for long-term simulations. -- You received this message because you are subscribed to the Google Groups sympy

[sympy] mod_wsgi with sympy

2015-07-24 Thread Paul Royik
Again I'm asking about mod_wsgi. I'm running some sort of webservice that uses sympy. Site is running under apache mod_wsgi. Users access it an everything is fine. But sometimes when I access the page, it loads too long (even simple example). It looks like memory used and another query just

Re: [sympy] mod_wsgi with sympy

2015-07-24 Thread Paul Royik
It never happens on the development server. On Saturday, July 25, 2015 at 6:36:42 AM UTC+3, Sudhanshu Mishra wrote: Hi, Does it happen on development server too? On Jul 25, 2015 1:39 AM, Paul Royik distan...@gmail.com javascript: wrote: Again I'm asking about mod_wsgi. I'm running some

[sympy] Re: Float precision

2015-07-14 Thread Paul Royik
str(a).rstrip('0') trick won't won't work with complex numbers, like 0.3+0.500i On Monday, July 13, 2015 at 7:08:27 PM UTC+3, Vinzent Steinberg wrote: On Sunday, July 12, 2015 at 6:18:26 PM UTC+2, Paul Royik wrote: Is it possible to set for float precision of 50

Re: [sympy] Re: Float precision

2015-07-15 Thread Paul Royik
Thank you. On Tuesday, July 14, 2015 at 11:10:38 PM UTC+3, Aaron Meurer wrote: The printer has an option to do this. Use sstr() and set the full_prec flag. Aaron Meurer On Tue, Jul 14, 2015 at 2:53 PM, Paul Royik distan...@gmail.com javascript: wrote: str(a).rstrip('0') trick won't

Re: [sympy] don't understand piece of code

2015-11-09 Thread Paul Royik
OK. Thanks. On Monday, November 9, 2015 at 12:42:51 PM UTC+2, Mateusz Paprocki wrote: > > Hi, > > On 9 November 2015 at 09:49, Paul Royik <distan...@gmail.com > > wrote: > > In Mul class I see this piece of code: > > > > def _eval_is_real(self): &g

[sympy] don't understand piece of code

2015-11-09 Thread Paul Royik
In Mul class I see this piece of code: def _eval_is_real(self): real = True zero = one_neither = False for t in self.args: if not t.is_complex: return t.is_complex Don't understand last two lines. How is it possible that if t is not

Re: [sympy] Powers and complex numbers

2015-11-04 Thread Paul Royik
tom Pow subclass that evaluates like > you want, and replace instances of Pow with your class before doing a > substitution. > > Aaron Meurer > > On Wed, Nov 4, 2015 at 1:19 PM, Paul Royik <distan...@gmail.com > > wrote: > > So, there is no way to do it us

Re: [sympy] Powers and complex numbers

2015-11-04 Thread Paul Royik
> Out[3]: -1 > > In [4]: real_root(-1, 3)**2 > Out[4]: 1 > > SymPy, like most math libraries, uses complex roots (i.e., principal > roots) because they have nicer mathematical properties. > > Aaron Meurer > > On Wed, Nov 4, 2015 at 3:30 AM, Paul Royik <distan...

[sympy] Powers and complex numbers

2015-11-04 Thread Paul Royik
I have the following expresssion: f=x**(Rational(2,3)) How can I get 1, when substituting (-1) instead of complex number? For now, I got complex number when run f.subs(x,-1).evalf() -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from

Re: [sympy] Powers and complex numbers

2015-11-06 Thread Paul Royik
(replace -1 with whatever value you want to replace). A "cleaner" > solution would be to create a custom Pow subclass that evaluates like > you want, and replace instances of Pow with your class before doing a > substitution. > > Aaron Meurer > > On Wed, Nov

Re: [sympy] Powers and complex numbers

2015-11-06 Thread Paul Royik
> > Aaron Meurer > > On Fri, Nov 6, 2015 at 1:42 PM, Paul Royik <distan...@gmail.com > > wrote: > > I decided to use "cleaner" solution: create a custom Pow subclass that > > evaluates. What method should I override? > > > > On Wed

[sympy] Re: Logcombine

2015-10-19 Thread Paul Royik
simplify(2log(5)) returns log(25) How to revert this? On Monday, October 19, 2015 at 10:35:45 AM UTC+3, Francesco Bonazzi wrote: > > > > On Monday, 19 October 2015 07:12:46 UTC+2, Paul Royik wrote: >> >> Hello. >> I noticed that simplify moves coefficient of logarith

[sympy] Re: Logcombine

2015-10-19 Thread Paul Royik
s: > > In [6]: l = log(25) > > In [7]: l.replace(lambda expr: isinstance(expr, log), lambda expr: factor( > expr)) > Out[7]: 2*log(5) > > > > On Monday, 19 October 2015 11:32:34 UTC+2, Paul Royik wrote: >> >> simplify(2log(5)) returns log(25) >> How to re

Re: [sympy] TypeError: cannot determine truth value of

2015-07-11 Thread Paul Royik
But this can't be reproduced. It is a bug from time to time. How can this be possible? On Friday, July 10, 2015 at 1:49:56 AM UTC+3, Aaron Meurer wrote: This looks like a bug. Numerical values should be evaluated in inequalities. Aaron Meurer On Thu, Jul 9, 2015 at 5:23 PM, Paul Royik

[sympy] TypeError: cannot determine truth value of

2015-07-08 Thread Paul Royik
I'm using sympy on my site together with django and mod_wsgi. From time to time (i.e. this error is not reproducible) I've got the following error: TypeError: cannot determine truth value of 3*pi/8 = 7*pi/8 In other words, error is raised when I try to compare numerical values. Can someone

  1   2   >