[sympy] Re: Area of Surface of Revolution Integral too Hard to be Computed By JULIA’ SymPy and Python’ SymPy

2023-01-23 Thread Kalevi Suominen
EDIT: In general, the square root does *not* simplify. On Monday, January 23, 2023 at 7:22:01 PM UTC+2 Kalevi Suominen wrote: > The derivative should actually be squared in the square root expression: > sqrt(1 + f'(x)^2) (see e.g. > https://en.wikipedia.org/wiki/Surface_of_revolution

[sympy] Re: Area of Surface of Revolution Integral too Hard to be Computed By JULIA’ SymPy and Python’ SymPy

2023-01-23 Thread Kalevi Suominen
should be able to handle it. In general, the square root does simplify. In that case the result will be a hyperelliptic integral, which is non-elementary and cannot be represented by means of common special functions. There is no support in SymPy for such integrals. Kalevi Suominen On Saturday

[sympy] Re: Doubt regarding as_leading_term for multivairate case

2021-11-09 Thread Kalevi Suominen
I think that 1/x should be returned in both cases. Kalevi Suominen On Tuesday, November 9, 2021 at 1:30:36 PM UTC+2 anderso...@gmail.com wrote: > Hey guys ! > > I am working on a Pr based on the as_leading_term() method of core/expr.py > (specifically on multivariable cases) . I w

[sympy] Re: How to narrow down "Invalid argument types for subtraction" in a system of equations?

2021-10-30 Thread Kalevi Suominen
I do not know solve well but I suspect that it cannot handle vector equations. There should probably be separate equations for the components. Kalevi Suominen On Saturday, October 30, 2021 at 9:28:06 AM UTC+3 Andreas Schuldei wrote: > I have a system of longish equations and when calling .so

[sympy] Re: Monitoring resource usage

2021-07-31 Thread Kalevi Suominen
You should send the request to SimPy's mailing list https://groups.google.com/g/python-simpy. This is SymPy's list. Kalevi Suominen On Saturday, July 31, 2021 at 5:04:06 PM UTC+3 Parthasarati Dileepan wrote: > Hi, I am trying to understand how to collect data for computing resource &g

[sympy] Re: What is the difference between is_symbol and is_Symbol?

2021-05-11 Thread Kalevi Suominen
>>> from sympy import Indexed, IndexedBase, Idx, symbols >>> i, j = symbols('i j', cls=Idx) >>> Aij = Indexed('A', i, j) >>> Aij.is_symbol True >>> Aij.is_Symbol False Aij is not a Symbol but it can be used like one in many cases. Kalevi Suominen

[sympy] Re: How to convert erf to decimal and compute normal distribution

2021-04-27 Thread Kalevi Suominen
evalf is a function. You have to give the precision as an argument: >>> Prob.evalf(4) 0.7333 >>> Prob.evalf(10) 0.7333134908 Kalevi Suominen On Monday, April 26, 2021 at 7:56:31 PM UTC+3 areebs...@gmail.com wrote: > Hi everyone, > I am trying to compute the normal di

[sympy] Re: Need help with Fourier transform

2021-04-21 Thread Kalevi Suominen
. Kalevi Suominen On Wednesday, April 21, 2021 at 7:29:50 PM UTC+3 emanuel.c...@gmail.com wrote: > Le mardi 20 avril 2021 à 13:01:13 UTC+2, jks...@gmail.com a écrit : > >> Fourier transform is currently implemented in SymPy only for integrable >> functions. None of those funct

[sympy] Re: Need help with Fourier transform

2021-04-20 Thread Kalevi Suominen
Fourier transform is currently implemented in SymPy only for integrable functions. None of those functions is integrable so SymPy cannot be used find the transform. Kalevi Suominen On Tuesday, April 20, 2021 at 11:38:08 AM UTC+3 aTPer wrote: > I am trying to compute the integral four

[sympy] Re: persistent detached head state

2021-04-06 Thread Kalevi Suominen
. I would consider removing those (after securing sufficient backup). Kalevi Suominen On Tuesday, April 6, 2021 at 6:05:38 PM UTC+3 smi...@gmail.com wrote: > I tried to rebase (and then decided to merge) a branch with master > > git co origin/pr/18489 > git rebase master > git r

[sympy] Re: Problem on dup_zz_zassenhaus?

2020-05-11 Thread Kalevi Suominen
algorithm. There are better bounds for estimating the size of the roots but they are not suitable for Zassenhaus' algorithm. Kalevi Suominen On Monday, May 11, 2020 at 10:38:39 AM UTC+3, Lagamura Chess wrote: > > Hello to everybody, > > Assisting my professor we managed

Re: [sympy] Error with inverse_laplace_transform (easy to solve)

2020-05-07 Thread Kalevi Suominen
returned, but those should be fixed anyway. Kalevi Suominen On Thursday, May 7, 2020 at 12:28:29 AM UTC+3, czgdp1807 wrote: > > May be using `as_int` with the second argument in `RisingFactorial` should > solve the problem. Just an idea. > > With Regards, > Gagandeep Sing

[sympy] Re: LambertW with lambdify

2020-04-23 Thread Kalevi Suominen
You can probable solve this by running `from scipy.special import lambertw` and then adding the function to the local lambdify dict: `f_S_mm = sym.lambdify(t, sol_mm_init.rhs, [{'LambertW': lambertw}, 'numpy'] )` Kalevi Suominen On Thursday, April 23, 2020 at 6:04:19 PM UTC+3, Luis Itzá

[sympy] Re: Redefine derivative of a IndexedBase

2020-03-03 Thread Kalevi Suominen
. As to 2), I believe that your way is correct. Kalevi Suominen On Tuesday, March 3, 2020 at 9:42:14 AM UTC+2, Lorenzo Monacelli wrote: > > Dear all, > I am manipulating an IndexedBase object, and I would like to modify its > behavior when differentiated. > > I thought to do something

[sympy] Re: Problem with subs and Simpify

2020-02-27 Thread Kalevi Suominen
Hi, It looks like solve returns a list containing the solution: 129 after solving for b1 [-(F_1 + G_1*b_{-1}**2 + G_1*b_{-2})/(2*F_1*b_{-1} - 1)] SymPy's Eq expects something that can be converted to an expression. In this case, this should work: Eq(b1, ex129b[0]) Kalevi Suominen

[sympy] Re: Gsoc 2020 : Implementing better integer factorization algorithms

2020-02-25 Thread Kalevi Suominen
Hi, I think that it is enough work for a summer of code (maybe even more than enough). It would be a good addition to the prime testing and factorization modules of SymPy. Kalevi Suominen On Tuesday, February 25, 2020 at 10:02:07 PM UTC+2, ABHINAV ANAND wrote: > > Hey, i hope the wor

[sympy] Re: Gsoc 2020 : Implementing better integer factorization algorithms

2020-02-16 Thread Kalevi Suominen
Hi, Yes, those algorithms are relevant to SymPy. For the first one, you might be interested in completing this PR: https://github.com/sympy/sympy/pull/2449. It should be possible to make it work also with finite fields and rings. Kalevi Suominen On Sunday, February 16, 2020 at 7:52:26 PM

[sympy] Re: Behaviour of symbols with and without assumptions in sympify

2020-02-14 Thread Kalevi Suominen
'a') == a False >>> sympify('a', {'a': a}) == a True Kalevi Suominen On Friday, February 14, 2020 at 5:03:02 PM UTC+2, ludi wrote: > > Hi > > If I do > > a=Symbol('a', real=True) > b=Symbol('b') > > and do > > solve(a**2-1,a) > I get > [−1, 1] > solve(b**2-1,b)

[sympy] Re: can't convert -oo to int

2019-10-12 Thread Kalevi Suominen
Hi, Limit values like -oo should not be directly substituted into complicated expressions. Instead of `g.subs(x, -oo)`, one should compute `limit(g, x, -oo)` (which is -1/2). Kalevi Suominen On Saturday, October 12, 2019 at 5:25:03 PM UTC+3, Rainer Dorsch wrote: > > Hi, > > I hav

[sympy] Re: Quaternions with non-commutative arguments

2019-09-03 Thread Kalevi Suominen
such as quaternions and (finite dimensional) Lie algebras. Kalevi Suominen On Tuesday, September 3, 2019 at 3:33:14 PM UTC+3, Maxence Mayrand wrote: > > The current implementation of quaternions (sympy/algebras/quaternion) > requires that the arguments are commutative expressions. But I can't see

[sympy] Re: Trandformed Field

2019-09-01 Thread Kalevi Suominen
There is currently some work in progress aiming to implement such transformations: https://github.com/sympy/sympy/pull/17507 Kalevi Suominen On Sunday, September 1, 2019 at 8:32:56 PM UTC+3, Sebastián Pérez wrote: > > Hi everybody, I hope you are doing well. > > I am trying to use s

[sympy] Re: Making experimental changes to the sympy sources

2019-08-27 Thread Kalevi Suominen
of the library in a directory where you have all access rights. The sources can be downloaded from https://www.sympy.org/en/download.html. Kalevi Suominen On Tuesday, August 27, 2019 at 1:07:36 PM UTC+3, David Bailey wrote: > > I have sympy running under Windows 10-64. > > Recently

[sympy] Re: Should symbolic integration fail gracefully?

2019-08-21 Thread Kalevi Suominen
We should have, for example, `if (b < 0) == True or ...` here File "C:\SymPyWorkbook\lib\site-packages\sympy\integrals\meijerint.py", line 1684, in _meijerint_indefinite_1 if b < 0 or f.subs(x, 0).has(nan, zoo): Kalevi Suominen On Tuesday, August 20, 2019 at 11:08:59 PM UTC

[sympy] Re: Can't get sympy.functions.atan2 to work with _eval_evalf

2019-08-20 Thread Kalevi Suominen
It seems that evalf of atan2 currently expects that its arguments are real. For symbols like x and y, this means that they have to be initialized with `real=True`. It does not suffice that evalf would return a real value as assumption code does not try to run evalf. Kalevi Suominen On Tuesday

[sympy] Re: Problem with evalf

2019-08-01 Thread Kalevi Suominen
Hi, Decimal(2.7) is not a SymPy type. Instead, you can use sympy.Rational(27, 10). Kalevi Suominen On Thursday, August 1, 2019 at 10:11:20 AM UTC+3, jm allard wrote: > > Hi > > See following code: > > >>> import sympy > > >>> sympy.E.evalf(25) > 2.

[sympy] Re: Creating a new matrix with removing columns in anoter doesn't work

2019-05-03 Thread Kalevi Suominen
Perhaps you should make a copy of q: w = q.copy(). Kalevi Suominen On Friday, May 3, 2019 at 10:50:29 AM UTC+3, murat...@erzurum.edu.tr wrote: > Hello, > > I am new to sympy and I want removed one column of a matrix and I want > to keep the original one also but when I

Re: [sympy] diffgeom: Lie Bracket violates Jacobi identities

2019-04-29 Thread Kalevi Suominen
It looks like there is a bug in Expr.coeff. It cannot find the coefficient of e_r in Commutator(v1, v2), apparently because of non-commutativity. Can you open an issue? Kalevi Suominen On Monday, April 29, 2019 at 7:56:43 PM UTC+3, ric wrote: > > A manual calculation suggests that LieD(L

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

2019-03-24 Thread Kalevi Suominen
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

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

2019-02-16 Thread Kalevi Suominen
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, >

Re: [sympy] Re: Azure pipelines

2019-01-28 Thread Kalevi Suominen
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: sympy inverse_laplace_transform, how to get the result quickly

2018-12-13 Thread Kalevi Suominen
p(-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, > When I try to

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

2018-11-21 Thread Kalevi Suominen
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. Yo

[sympy] Re: Pyglet plotting

2018-11-17 Thread Kalevi Suominen
"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 Suominen On Fri

[sympy] Re: polynomial division, irrational coefficients

2018-11-16 Thread Kalevi Suominen
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 Kalevi Su

[sympy] Re: Introduction and Question

2018-11-08 Thread Kalevi Suominen
, 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 quite a

[sympy] Re: Sympy returning weird polynomial roots

2018-10-03 Thread Kalevi Suominen
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: How to generate a symbolic multivariate polynomial of a given dimension in SymPy?

2018-08-01 Thread Kalevi Suominen
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 fi

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

2018-08-01 Thread Kalevi Suominen
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 wrote: > > I wa

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

2018-07-27 Thread Kalevi Suominen
ations 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/Pyth

[sympy] Re: I found some wrong integration

2018-04-13 Thread Kalevi Suominen
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: Inverse Fourier transform of a partial fraction decomposition looks to give a wrong answer

2018-04-10 Thread Kalevi Suominen
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: GSoc 2018 - Ordinary Differential Equations

2018-03-03 Thread Kalevi Suominen
sted 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 subscribed to the Google Grou

[sympy] Re: failures in master at Travis

2018-02-12 Thread Kalevi Suominen
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

[sympy] Re: failures in master at Travis

2018-02-12 Thread Kalevi Suominen
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

2018-02-12 Thread Kalevi Suominen
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] Re: Can someone help me with this command?

2018-02-12 Thread Kalevi Suominen
t; > 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 email to sympy+unsubscr...

[sympy] Re: Introducing myself to community

2017-12-09 Thread Kalevi Suominen
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: Perimeter of Ellipse

2017-11-04 Thread Kalevi Suominen
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: Rational function factor cancellation in simplify

2017-10-18 Thread Kalevi Suominen
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" group. To unsubscr

[sympy] Re: Rational with Symbol

2017-10-04 Thread Kalevi Suominen
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 sympy+un

[sympy] Re: performance issue of sympy 1.1.1

2017-09-29 Thread Kalevi Suominen
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: Huge subscripts in latex rendering of vector expressions

2017-09-03 Thread Kalevi Suominen
e 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" grou

Re: [sympy] roots of unity

2017-05-27 Thread Kalevi Suominen
. >>> 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 matrix computations whose entri

[sympy] Re: Mellin transform modification

2017-05-08 Thread Kalevi Suominen
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* >

[sympy] Re: integrate Heaviside give not correct answer

2017-04-28 Thread Kalevi Suominen
, ()), ((), (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 this message

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

2017-03-25 Thread Kalevi Suominen
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

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

2017-03-02 Thread Kalevi Suominen
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 Tuesday, February 28, 2017 at 6:11:25

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

2017-03-02 Thread Kalevi Suominen
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

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

2017-02-10 Thread Kalevi Suominen
ave 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 unsubscribe from this group and stop

[sympy] Re: Problem with PR#12130

2017-02-05 Thread Kalevi Suominen
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

[sympy] Re: trunc () result

2017-01-11 Thread Kalevi Suominen
g field=True 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: trunc () result

2017-01-11 Thread Kalevi Suominen
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

Re: [sympy] Compute minimal polynomial

2017-01-05 Thread Kalevi Suominen
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

[sympy] Re: Compute minimal polynomial

2017-01-04 Thread Kalevi Suominen
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

Re: [sympy] QEPCAD

2016-11-01 Thread Kalevi Suominen
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

Re: [sympy] QEPCAD

2016-10-30 Thread Kalevi Suominen
+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 link mig

Re: [sympy] QEPCAD

2016-10-08 Thread Kalevi Suominen
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

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

2016-09-19 Thread Kalevi Suominen
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

[sympy] Re: Galois Fields

2016-08-18 Thread Kalevi Suominen
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

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

2016-07-22 Thread Kalevi Suominen
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

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

2016-07-19 Thread Kalevi Suominen
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] Sphinx Documentation

2016-07-01 Thread Kalevi Suominen
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 >

[sympy] Re: On Parametric Polynomial

2016-06-02 Thread Kalevi Suominen
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] Re: polynomials over GF(2)

2016-05-24 Thread Kalevi Suominen
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', GF(2)

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

2016-05-24 Thread Kalevi Suominen
yError: 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, May

[sympy] Re: Italics in math equations

2016-05-24 Thread Kalevi Suominen
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] >

[sympy] Re: polynomials over GF(2)

2016-05-23 Thread Kalevi Suominen
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

[sympy] Re: Remainder sequences

2016-03-28 Thread Kalevi Suominen
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]( >

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

2016-03-01 Thread Kalevi Suominen
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,

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

2016-03-01 Thread Kalevi Suominen
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, >

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

2016-02-01 Thread Kalevi Suominen
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

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

2016-01-20 Thread Kalevi Suominen
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: How do I convert products and powers of sin, cos to sin, cos of sums and multiple angles?

2015-12-02 Thread Kalevi Suominen
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

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

2015-11-08 Thread Kalevi Suominen
; 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, >> instead of it being one of the coord_functions(). >> >> Thanks! >> >> On Fr

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

2015-11-06 Thread Kalevi Suominen
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: Singular Value

2015-11-02 Thread Kalevi Suominen
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

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

2015-10-22 Thread Kalevi Suominen
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) +

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

2015-09-04 Thread Kalevi Suominen
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

2015-09-04 Thread Kalevi Suominen
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) +

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

2015-09-03 Thread Kalevi Suominen
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

[sympy] Re: bug #9869

2015-08-29 Thread Kalevi Suominen
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

Re: [sympy] Re: analytical integrate

2015-07-23 Thread Kalevi Suominen
/4) - 1)) In [7]: print(simplify(a.diff(x))) x**(1/4)*sqrt(sqrt(x) + 1) (if you distribute sqrt(sqrt(x)) inside the sqrt(sqrt(x) + 1) you get sqrt(x + sqrt(x)) Aaron Meurer On Thu, Jul 23, 2015 at 2:37 PM, Kalevi Suominen jks...@gmail.com javascript: wrote: On Thursday, July 23, 2015

[sympy] Re: Introducing package maintainers

2015-07-13 Thread Kalevi Suominen
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 package

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

2015-07-03 Thread Kalevi Suominen
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 complex

[sympy] Re: Sympy zoo

2015-06-08 Thread Kalevi Suominen
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 in a

[sympy] Re: Tons of LaTeX warnings on Travis

2015-06-05 Thread Kalevi Suominen
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, Jason Moore

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

2015-06-05 Thread Kalevi Suominen
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 difference in

[sympy] Re: Help on Sympy HyperGeometric functions

2015-05-16 Thread Kalevi Suominen
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

[sympy] Re: Sine product integral

2015-03-31 Thread Kalevi Suominen
/2 + pi*m*z/L) instead of sin(pi*m*(L/2 + z)/L). Apparently its pattern matching does not cover the latter form. (fu will not be needed.) Kalevi Suominen -- You received this message because you are subscribed to the Google Groups sympy group. To unsubscribe from this group and stop receiving

[sympy] Re: expanding powers of cos()

2015-03-28 Thread Kalevi Suominen
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 sympy+unsubscr...@googlegroups.com. To post to this group, send email to sympy

  1   2   >