[sympy] Re: Quartic formula

Hi, Tuom I'm trying to attach a file answering your questions. I hope it arrives. Kalevi > > -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receivi

[sympy] Are all sympy Floats rational?

Hi, The attribute 'is_rational' seems to give surprising results: >>> from sympy import nfloat, sqrt >>> x = nfloat(sqrt(2)) >>> x 1.41421356237310 >>> x.is_rational True This leads to problems if used as a basis for decisions. There might be an explanation I don't see. (Other than 'all truncate

Re: [sympy] Re: Are all sympy Floats rational?

It looks like nobody really planned the semantics of Float.is_rational. I believe it does not make sense and should be removed from Float's dictionary. If this breaks something it should be fixed anyway. Kalevi Suominen -- You received this message because you are subscribed to the G

Re: [sympy] Re: Are all sympy Floats rational?

On Friday, March 21, 2014 10:46:19 PM UTC+2, Aaron Meurer wrote: > > Actually, that's a good point. Try removing it and see what tests fail. > > Aaron Meurer > > On Fri, Mar 21, 2014 at 2:55 PM, Kalevi Suominen > > > wrote: > > It looks like

Re: [sympy] Re: Are all sympy Floats rational?

oat). A missing attribute would reveal every attempted use. I think this is how it should be in production code, to warn off any accidental use. Kalevi Suominen -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this grou

Re: [sympy] Re: Removing the UniversalSet

xists) has contradictory properties, all later axiomatic set theories have restricted the correspondence between sets and statements in one way or another. Kalevi Suominen On Monday, May 12, 2014 10:31:53 PM UTC+3, Aaron Meurer wrote: > > Another argument for keeping it is that it works

[sympy] Re: elliptic_e._eval_conjugate for one arg?

On Sunday, September 28, 2014 4:15:06 AM UTC+3, Chris Smith wrote: > > elliptic_e can take 1 or two args but the _eval_conjugate assumes 2 (and > fails in a branch I am working on). The following mod works without > returning anything for the 1 arg case, but can someone tell me if I should > b

[sympy] Re: elliptic_e._eval_conjugate for one arg?

On Sunday, September 28, 2014 5:19:09 PM UTC+3, Sergey Kirpichev wrote: > > Oh, no. It's just the Legendre complete elliptic integral of the second > kind. Please read the docs. > > Yes, you are right. But denoting the argument by z is very confusing as it really is the modulus. There ar

Re: [sympy] Sympy assumptions

its derivative P_y, if necessary.) > and how do you make cuts in a consistent manner? > Branch cuts should not be needed as the functions are single-valued on the Riemann surface. Kalevi Suominen -- You received this message because you are subscribed to the Google Groups "sympy"

[sympy] Re: Tensors and differential geometry

On Friday, December 19, 2014 8:20:58 AM UTC+2, Alok Gahlot wrote: > > my question is what is the difference in contravariant and covariant >> tensor >> > The simple (and stupid) answer is that they are objects of different type. In more detail, their components behave differently under general

Re: [sympy] Probable bug in meijerint.py

I think it is run in every meijerint test. But the situation it is intended to check should probably never occur in the G-functions produced by the code (in _rewrite_single()). Otherwise there would be an error in its lookup table. -- You received this message because you are subscribed to the

[sympy] Re: Slow inverse_laplace_transformation

On Sunday, February 8, 2015 at 12:06:46 AM UTC+2, Sten Sogaard wrote > > > H = 1/(0.07071*s**2 + 258.0*s + 110.) > inverse_laplace_transform(1/s*H,s,t) > > > Any ideas? > The current implementation leads to floating point numbers with very large positive and negative exponents. After the s

[sympy] Re: IndexeFunction

Technically an index is just another argument, though it often takes only a finite number of values. Would it be possible to use functions like f(i, x)? -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving

[sympy] Re: Algorithm for integration of mixed functions

Sympy currently uses the parallel Risch algorithm for this function. The full Risch algorithm for functions with nonrational algebraic terms is not implemented. -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop

[sympy] Status of RootOf

The RootOf class is explicitly defined as representing numerical complex roots of equations. The only possibility of dealing with roots of polynomials in other domains, e.g., in function fields, is by means of radicals. (Though this may not be a too big impediment in practise.) However, there a

[sympy] Re: Status of RootOf

It should not be necessary to lose any of the current functionality of RootOf. A call will now return a complex value. It could be replaced with a call to rootof() which would do exactly the same. (Some editing would be needed in solvers, of course.) The field of complex numbers is special becau

Re: [sympy] Re: Status of RootOf

> Couldn't we still just number the roots? Or would that lead to > inconsistencies? > It would be fairly complicated to do in full generality. It is necessary to construct one root at a time. (This is essentially so in the complex case, too.) For each root one has to remove those that have al

Re: [sympy] Re: Slow inverse_laplace_transformation

he denominator before rounding. > On Sun, Feb 8, 2015 at 2:57 PM, Kalevi Suominen > wrote: > > > > > > On Sunday, February 8, 2015 at 12:06:46 AM UTC+2, Sten Sogaard wrote > >> > >> > >> H = 1/(0.07071*s**2 + 258.0*s + 110

Re: [sympy] Re: CAD Implementation - GSoC 2015

On Saturday, March 21, 2015 at 3:40:07 AM UTC+2, Luv Agarwal wrote: > > I have added the implementation details of few functions in the draft > > . > > How can I implement refine_root (a

Re: [sympy] Re: CAD Implementation - GSoC 2015

On Sunday, March 22, 2015 at 11:33:43 AM UTC+2, Luv Agarwal wrote: > Thanks, I looked at them. > How can I implement dup_isolate_real_roots_sqf for algebraic domains? > It seems that the algorithm in rootisolation should work for real algebraic numbers as well. The critical points are the ineq

Re: [sympy] Re: CAD Implementation - GSoC 2015

On Tuesday, March 24, 2015 at 2:44:03 AM UTC+2, Luv Agarwal wrote > I have updated my proposal > > . > Please have a look especially at the Timeline part. Is there anything that > I sho

Re: [sympy] Re: CAD Implementation - GSoC 2015

AlgebraicField require revision. Kalevi > > On Tue, Mar 24, 2015 at 7:26 AM, Kalevi Suominen > wrote: > > > > > > On Tuesday, March 24, 2015 at 2:44:03 AM UTC+2, Luv Agarwal wrote > >> > >> I have updated my proposal. > >> Please h

Re: [sympy] Re: CAD Implementation - GSoC 2015

On Wednesday, March 25, 2015 at 1:54:24 AM UTC+2, Luv Agarwal wrote: > > > Is primitive_element doing something similar to simple > ?. > > If yes, will there be difference in ef

Re: [sympy] Re: CAD Implementation - GSoC 2015

On Thursday, March 26, 2015 at 12:27:13 AM UTC+2, Luv Agarwal wrote: > What reconsideration are you talking about? > I was just pondering if you really have time to complete also isolate_roots in week 1 already. Even proj might take some time... Kalevi -- You received this message

[sympy] Re: expanding powers of cos()

-2)*x) + ...) / 2**(n - 1) . (If n is even, the last term in the parentheses is C(n, n/2)/2.) Kalevi Suominen -- 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 sym

[sympy] Re: Sine product integral

n), (0, True)) > > In [17]: Imn2 > Out[17]: Integral(sin(pi*m*(L/2 + z)/L)*sin(pi*n*(L/2 + z)/L), (z, -L/2, > L/2)) > > Best regards, > > Julien > It seems that SymPy can compute the integral if you write sin(pi*m/2 + pi*m*z/L) instead of sin(pi*m*(L/2 + z)/L). Apparently

[sympy] Re: Help on Sympy HyperGeometric functions

> But the values generated by sympy's hyp2f1() function are different than the ones given in the example. Can you explain in more detail in which way the values are different? -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this g

[sympy] Re: Tons of LaTeX warnings on Travis

On Friday, June 5, 2015 at 9:02:44 PM UTC+3, Jason Moore wrote: > > I think we found out timeout issue but the warnings are still odd, > especially since we don't get them locally. That's be nice to fix too. > > > Jason > moorepants.info > +01 530-601-9791 > > On Fri, Jun 5, 2015 at 10:51 AM, Ja

Re: [sympy] Re: Tons of LaTeX warnings on Travis

On Friday, June 5, 2015 at 10:09:19 PM UTC+3, Jason Moore wrote: > > But if I `make clean && make latex` locally I do not get these warnings. > I'm not sure where latex aux files would be other than in the directories > generated under those deleted by clean. > > Perhaps there is some differenc

[sympy] Re: Sympy zoo

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 and log(z) as z tends to 0. The default domain in SymPy is the complex field, so the limits are computed i

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

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 in SymPy. The inverse transform involves integration along a vertical line in the comp

[sympy] Re: Introducing package maintainers

On Monday, July 13, 2015 at 4:11:30 AM UTC+3, Jason Moore wrote: > > After discussions at SciPy, we decided to try out a new idea to help > delegate maintenance responsibilities to a broader group of people. In > particular, we'd like to have a volunteer maintainer for each substantial > packa

[sympy] Re: analytical integrate

On Thursday, July 23, 2015 at 8:31:41 PM UTC+3, Denis Akhiyarov wrote: > > Sympy cannot do this? > > integrate(sqrt(x+sqrt(x))) > > BTW, this is computed by SAGE: > > -1/12*(3*(sqrt(x) + 1)^(5/2)/x^(5/4) - 8*(sqrt(x) + 1)^(3/2)/x^(3/4) - > 3*sqrt(sqrt(x) + 1)/x^(1/4))/((sqrt(x) + 1)^3/x^(3/2) -

Re: [sympy] Re: analytical integrate

) - 3*sqrt(sqrt(x) + 1)/x^(1/4))/((sqrt(x) + 1)^3/x^(3/2) - >>> 3*(sqrt(x) + 1)^2/x + 3*(sqrt(x) + 1)/sqrt(x) - 1) + 1/8*log(sqrt(sqrt(x) + >>> 1)/x^(1/4) + 1) - 1/8*log(sqrt(sqrt(x) + 1)/x^(1/4) - 1)") >>> >>> In [7]: print(simplify(a.diff(x))) >>>

[sympy] Re: bug #9869

On Saturday, August 29, 2015 at 9:05:34 PM UTC+3, Preet Singh Khalsa wrote: > > Hey! > I am new to open source. I would like to fix bug #9869 (label: easy to > fix). Can someone please get me started? > Thanking in anticipation, > Preet. > It seems to be already taken care of by https://github

[sympy] Re: Distribute derivative of sum with constant term after substitution

On Thursday, September 3, 2015 at 8:54:54 AM UTC+3, Carl Sandrock wrote: > > If I do > > import sympy > k, V, Vprime = sympy.symbols('k, V, Vprime') > print sympy.diff(k + V(t), t) > > I get Derivative(V(t), t) as I expect - the derivative distributes and > the constant term has zero derivative.

[sympy] Re: Equivalent to Sage's assume: inverse Laplace of second order transfer function

On Thursday, September 3, 2015 at 7:19:38 PM UTC+3, Carl Sandrock wrote: > > I am trying to build a workbook to illustrate the effect of various > parameters of second order transfer functions. The full workbook is on > GitHub >

[sympy] Re: Equivalent to Sage's assume: inverse Laplace of second order transfer function

It appears possible to compute the transform for ``zeta > 1`` by replacing ``zeta`` with ``1 + zeta``. In [22]: G = K/(tau**2*s**2 + 2*tau*(1 + zeta)*s + 1) In [24]: print(inverse_laplace_transform(G, s, t)) K*(exp(t*sqrt(zeta)*sqrt(zeta + 2)/tau) - 1)*(exp(t*sqrt(zeta)*sqrt(zeta + 2)/tau) + 1)

[sympy] Re: Result depends on the sign of I

On Monday, September 7, 2015 at 9:37:58 AM UTC+3, Paul Royik wrote: > > I see this error when try to calculate limit of (-1)^x /sqrt(x) as x > approaches infinity. > Can somebody explain me what is the problem here? > ``(-1)^x`` is equal to ``exp(pi*I*x)``. As the function has no serial expans

[sympy] Re: Problems with integrating ln(x-y) wrt. y

On Thursday, October 22, 2015 at 4:56:35 PM UTC+3, Filip wrote: > > I'm trying to integrate ln(x-y) dy, but the result I get from SymPy is > very different from Maple's answer. > Here is what I get when I integrate in SymPy: > In [2]: > > integrate(ln(x-y), y) > > Out[2]: > > -x*log(-x + y) + y*

[sympy] Re: Singular Value

On Monday, November 2, 2015 at 1:47:52 AM UTC+2, Trevor Hutchinson wrote: > > Hello, > > I have what I think is a simple problem but I am struggling. I seek the > non-trivial solution to a homogenous linear system of equations given by Ax > = b, where A is matrix of symbolic variables, x is coe

[sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

On Thursday, November 5, 2015 at 3:56:09 PM UTC+2, Imran Ali wrote: > > If I simply hard code each metric, every thing works properly. For > example, take egg carton surface, > > from sympy.diffgeom import Manifold, Patch, CoordSystem, TensorProduct > from sympy import sin,cos > dim

[sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

l(BaseScalarField(CoordSystem(Symbol('egg_carton'), >> Patch(Symbol('P'), Manifold(Symbol('M'), Integer(2))), Tuple(Symbol('u'), >> Symbol('v'))), Integer(0))) >> >> I see now what the error is! My code used u as a sympy expression,

[sympy] Re: How do I convert products and powers of sin, cos to sin, cos of sums and multiple angles?

On Wednesday, December 2, 2015 at 5:22:24 AM UTC+2, Andrey Grozin wrote: > > вторник, 1 декабря 2015 г., 12:58:34 UTC+6 пользователь Andrey Grozin > написал: >> >> I need the action inverse to expand_trig (I think it is needed more often >> than expand_trig; suppose I want to multiply two trunc

[sympy] Re: Help: simple univariate polynomial, sympy solve gives no results

On Wednesday, January 20, 2016 at 4:37:18 AM UTC+2, Junwei Huang wrote: > > Hello, I am new to sympy and try to solve the following equation > > import sympy as sy > A,B,C,D,x=sy.var('A,B,C,D,x',positive=True) > sy.solve(A*x**5+B*x**4+C*x-D,x) > > but got no result. There are no roots, or I used

[sympy] Re: Rewriting expressions as a Meijer G-function for integration

On Thursday, January 28, 2016 at 5:14:41 AM UTC+2, Ondřej Čertík wrote: > > On Wed, Jan 27, 2016 at 4:30 PM, Ondřej Čertík > wrote: > > On Wed, Jan 27, 2016 at 3:34 PM, Ondřej Čertík > wrote: > >> Hi, > >> > >> Our Meijer G-function integrator can find integrals (definite and > >> indefini

[sympy] Re: GSOC16 : Risch algorithm for symbolic integration

On Tuesday, March 1, 2016 at 5:01:21 PM UTC+2, Ashwani Gautam wrote: > > Hi, the idea "Risch algorithm for symbolic integration" seems interesting > to me and thus want to contribute to it, anyone who can give me the > direction to proceed further, the latest information about the project, > w

Re: [sympy] Bug in factor() and simplify()

On Tuesday, March 1, 2016 at 5:53:11 PM UTC+2, Oscar wrote: > > On 26 February 2016 at 17:15, Aaron Meurer > wrote: > >> I traced this to here: > >> > >> In [1]: expr = x*sin(y)**2 + x*cos(y)**2 + 0.1*sin(y)**2 + > 0.1*cos(y)**2 > >> > >> In [2]: p = Poly(expr, x, cos(y), sin(y), domain='R

[sympy] Re: Remainder sequences

On Sunday, March 27, 2016 at 6:20:57 PM UTC+3, ecir...@gmail.com wrote: > > Hello! > > I would like to compute the remainder sequence of two polynomials as used > by GCD. If I understood the Wikipedia article about [Pseudo-remainder > sequence]( > https://en.wikipedia.org/wiki/Polynomial_greate

[sympy] Re: polynomials over GF(2)

On Tuesday, May 24, 2016 at 3:37:42 AM UTC+3, Phillip M. Feldman wrote: > > I would like to perform operations on polynomials over GF(2), i.e., > polynomials with binary coefficients. Is there a way to do this with SymPy? > > There is a sympy.polys.galoistools module, but I haven't found any us

[sympy] Re: Italics in math equations

On Tuesday, May 24, 2016 at 12:09:09 PM UTC+3, Björn Dahlgren wrote: > > I was recently made aware that only variables should be in italics in > equations (see e.g. [1] > , > > [2] >

Re: [sympy] Re: polynomials over GF(2)

on: SympifyError: x**5 + 1 mod 2 > > Any advice will be appreciated. > div(p, q) is intended for expressions. p and q are polynomials, not expressions. Division of polynomials is defined as a method: In [4]: p.div(q) Out[4]: (x**4 + x**3 + x**2 + x + 1 mod 2, 0 mod 2) > On Mon

Re: [sympy] Re: polynomials over GF(2)

On Wednesday, May 25, 2016 at 7:55:46 AM UTC+3, Kalevi Suominen wrote: > > > > On Wednesday, May 25, 2016 at 3:31:23 AM UTC+3, Phillip M. Feldman wrote: >> >> I tried the following: >> >> In [1]: from sympy import * >> In [2]: R, x= ring('x', G

[sympy] Re: On Parametric Polynomial

On Thursday, June 2, 2016 at 8:11:26 AM UTC+3, Yiming Yang wrote: > > Hello All, > > My question may be naive, as I'm new to Sympy. I tried to find how > Sympy defines parametric polynomials, i.e. polynomials with indeterminate > coefficients (K[U][X] with coefficients in K[U]), but could

Re: [sympy] Sphinx Documentation

On Friday, July 1, 2016 at 12:10:33 PM UTC+3, SAMPAD SAHA wrote: > > I was trying to add the docstrings of Singularity Function module in the > documentation in PR 11237 (the > last commit). But there was an error with message "* WARNING: Inline > li

[sympy] Re: Function and Indexed objects behavior with .has( )

On Tuesday, July 19, 2016 at 11:52:43 AM UTC+3, Francesco Bonazzi wrote: > > > Concerning functions, a similar problem could arise: > > In [15]: f(x).diff(f(x)) > Out[15]: 1 > > In [17]: f(x).diff(f(y)) > Out[17]: 0 > > > I think that output 17 should be a Dirac delta function: DiracDelta(x - >

Re: [sympy] Re: Function and Indexed objects behavior with .has( )

On Friday, July 22, 2016 at 12:02:39 PM UTC+3, Francesco Bonazzi wrote: > > > > Kalevi, I think that generalizing the Dirac delta to complex numbers is a > bit out of scope. Besides, do SymPy users really need it? > My thoughts, exactly. I only wanted to point out that Dirac delta should not b

[sympy] Re: Galois Fields

On Thursday, August 18, 2016 at 9:30:23 PM UTC+3, clemens novak wrote: > > Hello, > > I try to perform calculations in GF(2^3); e.g. calculate (x^2+x)(x+1) = > x^3 + 2x^2 + x = x^3 + x as polynomial coefficients are modulo-2. Using the > irreducible polynomial x^3 + x + 1 I arrive at (x^2+x)(x+

[sympy] Re: Results for integrate(Rational(1,4)*(1/(4-t)))

There are several integrators in SymPy. The latter form is produced by ratint. >>> from sympy.integrals.rationaltools import ratint >>> from sympy.abc import t >>> ratint(1/(4*(-t + 4)), t) -log(t - 4)/4 Kalevi Suominen On Monday, September 19, 2016 at 5:06:22 AM

Re: [sympy] QEPCAD

On Saturday, October 8, 2016 at 3:40:47 PM UTC+3, Ferran Pujol Camins wrote: > > In sympy 1.0 docs I see a general n-dimensional point: > http://docs.sympy.org/latest/modules/geometry/points.html > > Anyway, let's say we create our custom class. What's the appropiate class > to store each coo

Re: [sympy] QEPCAD

4 UTC+2, brombo va escriure: >> >> You might find the sympy modules at github.com/brombo/galgebra useful >> since they include an n-dimensional vector space with operations. Look at >> galgebra.pdf in the doc directory for a complete description. The >> following lin

Re: [sympy] QEPCAD

needed for CAD. On Tuesday, November 1, 2016 at 7:08:09 PM UTC+2, Ferran Pujol Camins wrote: > > Thank you Kalevi, is something like this what solve does? > > El diumenge, 30 octubre de 2016 19:30:44 UTC+1, Kalevi Suominen va > escriure: >> >> >> The implementation

[sympy] Re: Compute minimal polynomial

On Wednesday, January 4, 2017 at 6:25:02 PM UTC+2, ecir...@gmail.com wrote: > > I have two algebraic numbers defined by minimal polynomials: > > x, z = symbols('x, z') > p = Poly((x-1)*(x-2)*(x-3)*(x-4)) > q = Poly((x-5)*(x-6)*(x-7)*(x-8)) > > and I would like to compute the sum of th

Re: [sympy] Compute minimal polynomial

On Thursday, January 5, 2017 at 1:08:22 PM UTC+2, ecir...@gmail.com wrote: > > No, not a typo but an oversight from me. > > Is a polynomial "irreducible" if `factor_list(Poly(...))` is exactly of > length 1? > It is irreducible if there is only one factor *and* its exponent is 1. The square of

[sympy] Re: trunc () result

On Wednesday, January 11, 2017 at 6:29:40 PM UTC+2, swapnil sharma wrote: > > trunc(2*x**3 + 3*x**2 + 5*x + 7, 3) gives result : -x**3 - x + 1. is there > some way to get coefficents positive > instead of negative (i.e. here soln be like: 2*x**3+2*x+1)? > > One possibility is to write trunc(2

[sympy] Re: trunc () result

will force it to be the smallest field containing the coefficients, the field of rational numbers, in this case. > > On Wednesday, January 11, 2017 at 11:51:37 PM UTC+5:30, Kalevi Suominen > wrote: >> >> >> >> On Wednesday, January 11, 2017 at 6:29:40 PM UTC+2, swapnil

[sympy] Re: Problem with PR#12130

On Sunday, February 5, 2017 at 6:04:54 PM UTC+2, Yathartha Joshi wrote: > > Can anyone help me with [PR#12130]( > https://github.com/sympy/sympy/pull/12130). I don't know where I went > wrong. > All builds have been failing for the last couple of days because of a bug in Travis. -- You recei

[sympy] Re: SO(3) and SU(2) lie algebras

mplex simple Lie algebras have been dealt with. It would be useful to have an implementation of real simple Lie algebras as well. You are welcome to contribute. Kalevi Suominen -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubsc

[sympy] Re: Mentor for Possible GSoC 2017 project : Cylindrical Algebraic Decomposition

algebraic numbers has not been developed. The implementation should be capable of finding the real roots of a polynomial with real algebraic coefficients and ordering them properly. I would start by investigating that subject. Kalevi Suominen On Thursday, March 2, 2017 at 11:29:24 AM UTC+2, Arif

Re: [sympy] GSOC Project: Expand the number theory class

ambitious. Have you considered how much you would be able to implement in three months? Kalevi Suominen > > I believe it is a great addition to Sympy as it would definitely serve as > a complement to the current crypto module, specifically the RSA method. > > > On Tuesda

[sympy] Re: GSoC 2017: Lie groups in SymPy

There is a module on complex Lie algebras in SymPy that also deals with Dynkin diagrams. The possibility of adding real Lie algebras was discussed in this thread <https://groups.google.com/forum/#!topic/sympy/lbTwNgehbWk> some time ago. Kalevi Suominen On Saturday, March 25, 2017 at 4

[sympy] Re: integrate Heaviside give not correct answer

rg(((3, 1), ()), ((), (1, 0)), x), > True)) > > Manual integration will give the correct answer: > >>> integrate(Heaviside(x-1)*(x-1) + 1, (x, 0, 2), manual=True) > 5/2 > So the bug is in the evaluation of the indefinite integral. Kalevi Suominen -- You received th

[sympy] Re: Mellin transform modification

On Monday, May 8, 2017 at 8:13:42 PM UTC+3, Chris Smith wrote: > > Could someone please check to see if the modification of the Mellin > transformation at line 117 of *sympy/integrals/tests/test_transforms.py* >

Re: [sympy] roots of unity

trivial sympification routine. >>> class MyMatrix(Matrix): ... _sympify = staticmethod(lambda x: x) ... Kalevi Suominen On Friday, May 26, 2017 at 9:39:46 PM UTC+3, Calvin McPhail-Snyder wrote: > > Have there been any relevant updates since this post? I sometimes have to > do m

[sympy] Re: Huge subscripts in latex rendering of vector expressions

ne a subclass of CoordSys3D that redefines vector_names and variable_names in its __init__ function. (Some experimentation may be necessary to see exactly which attributes need to be changed.) Kalevi Suominen -- You received this message because you are subscribed to the Google Groups "sympy"

[sympy] Re: performance issue of sympy 1.1.1

t; or ray intersect with polygon),encloses_point.* > > * Json.dou* > > Can you find out how `primerange` is called in your code? Kalevi Suominen -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this gro

[sympy] Re: Rational with Symbol

ould probably use n/sympy.Integer(261), or briefly n/S(261) where S = sympy.sympify. Kalevi Suominen -- 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 sym

[sympy] Re: Rational function factor cancellation in simplify

t; from sympy import * >>> x = Symbol('x') >>> x != -1 True It is possible to substitute a number for x, but then then the expression will also be different. Kalevi Suominen -- You received this message because you are subscribed to the Google Groups "sympy

[sympy] Re: Perimeter of Ellipse

This looks like a computation of the arithmetic-geometric mean (https://en.wikipedia.org/wiki/Arithmetic%E2%80%93geometric_mean). There is an implementation called agm in mpmath. You could also use that for a complete elliptic integral of the first kind. Kalevi Suominen On Saturday, November

[sympy] Re: Introducing myself to community

in the project . > Thank you :-) > Hi, This link contains useful information: https://github.com/sympy/sympy/wiki/Introduction-to-contributing. Kalevi Suominen -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this

[sympy] Re: Can someone help me with this command?

ng my night! > > I get 683772⋅t^783, but I have not checked that for correctness. Kalevi Suominen -- 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 em

[sympy] Re: failures in master at Travis

Sphinx updated 1.6.7 -> 1.7.0 On Monday, February 12, 2018 at 3:58:19 PM UTC+2, Chris Smith wrote: > > Runs are failing the TEST_SPHINX tests with the following error: > > Warning, treated as error: > /home/travis/build/sympy/sympy/sympy/assumptions/handlers/ntheory.py:docstring > of sympy.assum

[sympy] Re: failures in master at Travis

Maybe this would help. --- a/sympy/assumptions/handlers/ntheory.py +++ b/sympy/assumptions/handlers/ntheory.py @@ -136,7 +136,7 @@ def Mul(expr, assumptions): Odd * Odd -> Odd Even * Even -> Even Integer * Integer -> Even if Integer + Integer = Odd -

[sympy] Re: failures in master at Travis

It seems that another weak point was found: Warning, treated as error: /home/travis/build/sympy/sympy/sympy/codegen/cfunctions.py:docstring of sympy.codegen.cfunctions.exp2.eval:15:Unexpected indentation. make: *** [html] Error 2 On Monday, February 12, 2018 at 6:52:48 PM UTC+2, Francesco Bonaz

[sympy] Re: GSoc 2018 - Ordinary Differential Equations

or the members are requested to fix this. > > Thankyou. > Hi, Google currently gives the following link: http://rubendebeerst.de/master/. It probably leads to the same master thesis as the broken link originally did. Kalevi Suominen -- You received this message because you are subs

[sympy] Re: GSoC 2018 - Univariate polynomials over algebraic domains

ensions of Q can be created by QQ.frac_field(self, *symbols, **kwargs) and their elements are of type `FracElement`. Transcendental extensions of F_p can be created similarly. Algebraic extensions of Q are created by QQ.algebraic_field(self, *extension) and their elements are of type ANP. There

[sympy] Re: Inverse Fourier transform of a partial fraction decomposition looks to give a wrong answer

variable assumed by SymPy. The transforms are also differently normalized. That may explain the differing results. Kalevi Suominen On Tuesday, April 10, 2018 at 9:16:57 AM UTC+3, Janko Slavič wrote: > > I am trying to replicate this: > > https://math.stackexchange.com/questions/228

[sympy] Re: I found some wrong integration

olframe alpha, symbolab answer: > 6 * x**(Rational(1,6)) - 6 *atan(x**(Rational(1,6))) > Both answers are correct. They only differ by a constant as atan(1/x**(1/6)) is equal to π/2 - atan(x**(1/6)). Kalevi Suominen -- You received this message because you are subscribed to the Google G

[sympy] Re: using Sympy for solving PDAEs as taylor series?

r equations with analytic coefficients. For the algorithm, you should look into the proof of the theorem. Kalevi Suominen On Friday, July 27, 2018 at 7:13:02 PM UTC+3, foadsf wrote: > > I posted this question here on Reddit > <https://www.reddit.com/r

[sympy] Re: How to generate a symbolic multivariate polynomial of a given dimension in SymPy?

variables this could be done: indices = [i for i in itertools.product(range(4), repeat=3) if sum(i) < 4] a = IndexBase('a') coeffs = {i: a[i] for i in indices} vars = symbols('x:3') Poly(coeffs, *vars) Kalevi Suominen On Wednesday, August 1, 2018 at 2:01:13 PM UTC+3, foadsf w

[sympy] Re: How to generate a symbolic multivariate polynomial of a given dimension in SymPy?

Correction: IndexedBase instead of IndexBase On Wednesday, August 1, 2018 at 4:10:35 PM UTC+3, Kalevi Suominen wrote: > > I would first generate the list of monomial indices by using e.g. > itertool.product, then create a dictionary containing the indexed > coefficients, and fina

[sympy] Re: Sympy returning weird polynomial roots

mpy/sympy/issues/14690> for more information.) Kalevi Suominen On Wednesday, October 3, 2018 at 9:03:22 AM UTC+3, Cesar Gomes wrote: > > Hello everyone! > > While trying to find the zeros of the following polynomial: > > > q = (-2/15)*x**3 + (23/10)*x**2 - (47/30)*x - 21 &

[sympy] Re: Introduction and Question

, for example). Kalevi Suominen On Thursday, November 8, 2018 at 10:15:52 PM UTC+2, Sara J Zelenberg wrote: > > Hello, > > My name is Sara and I'm a post-doc in the Applied Mathematics department > at IIT (in Chicago). > > One of my thesis projects, I tinkered qui

[sympy] Re: polynomial division, irrational coefficients

back to an expression if desired. >>> p = Poly(16*sqrt(3)*h**3 - 36*h**2 + 9, extension=True) >>> q = Poly(6*(4*sqrt(3)*h**2 - 12*h + 3*sqrt(3)), extension=True) >>> p.quo(q) Poly(2/3*h + sqrt(3)/6, h, domain='QQ') >>> _.as_expr() 2*h/3 + sqrt(3)/6 Ka

[sympy] Re: Pyglet plotting

"license" for more information. >>> from sympy.plotting.pygletplot import PygletPlot as Plot >>> from sympy import symbols >>> from sympy.abc import x, y, z >>> Plot(x*y**3-y*x**3) [0]: -x**3*y + x*y**3, 'mode=cartesian' >>> Kalevi Suomin

[sympy] Re: Unrelated PRs failing Travis tests due to test_sympify

.6 to 3.6.7. The change breaks sympification of `lambda` because the code expects only two tokens instead of three. Kalevi Suominen On Wednesday, November 21, 2018 at 1:16:03 PM UTC+2, Oscar wrote: > > Hi all, > > Travis tests are failing in unrelated PRs due to the same issue.

[sympy] Re: sympy inverse_laplace_transform, how to get the result quickly

- 15625*exp(-t/5 + 11/10)*Heaviside(t - 1/2)/48 + 1944*exp(-t/6 + 13/12)*Heaviside(t - 1/2)/5 - 117649*exp(-t/7 + 15/14)*Heaviside(t - 1/2)/720 + E*Heaviside(t - 1/2) Kalevi Suominen On Thursday, December 13, 2018 at 6:06:06 PM UTC+2, rupu.y...@gmail.com wrote: > > Dear all, > Whe

Re: [sympy] Re: Azure pipelines

I have found Travis logs quite useful as the browser can show them. Azure logs come in zip files which are rather inconvenient to deal with. I don't think that dropping Azure would be a loss for me. Kalevi Suominen On Monday, January 28, 2019 at 10:01:51 PM UTC+2, Aaron Meurer wrote: &

[sympy] Re: GSoC'19 Project in Multivariate polynomials and factorization and proposed plan of action.

I think that most of those topics have already been implemented. You can find the implementations in factortools.py and galoistools.py to see what remains to be done. Kalevi Suominen On Saturday, February 16, 2019 at 7:15:15 AM UTC+2, shrivastavaavi...@gmail.com wrote: > > Hello, >

[sympy] Re: Status on Risch Algorithm for Symbolic Integration

Hi, The latest update on the status is this project report: https://github.com/gxyd/GSoC-2017-Report. Kalevi Suominen On Sunday, March 24, 2019 at 9:35:13 AM UTC+2, Anirudh Chandramouli wrote: > > Hello all, > > I am a second year undergrad Computer Science student at I

  1   2   >