## Where communities thrive

• Join over 1.5M+ people
• Join over 100K+ communities
• Free without limits
##### Activity
• 08:54
Noor-Kalibbala closed #23503
• 08:53
Noor-Kalibbala commented #23503
• 08:40
Noor-Kalibbala edited #23503
• 08:20
sympy-bot commented #23504
• 08:19
user202729 opened #23504
• 08:15
Noor-Kalibbala opened #23503
• 03:36
user202729 opened #23502
• 03:22
hEzoooo starred sympy/sympy
• 02:18
jon-deng forked
• 02:05
bertiewooster commented #23462
• 02:04
bertiewooster commented #23462
• 01:41
jon-deng starred sympy/sympy
• 00:25
sympy-bot commented #23462
• 00:24
sympy-bot commented #23462
• 00:24
bertiewooster synchronize #23462
• 00:11
sympy-bot commented #23501
• 00:11
sympy-bot commented #23501
• 00:11
oscarbenjamin synchronize #23501
• 00:06
sympy-bot commented #23501
• 00:06
sympy-bot commented #23501
The lack of solutions is definitely a precision issue - if you input all entries as Rational (e.g. Rational(1,2) instead of 1/2), then linsolve returns the same as Maple (though obviously that's awkward)

A potential solution could be to define

rationalize = lambda x: [nsimplify(e) for e in x]

(this will turn numbers from a list into rationals). And then do

PI = Matrix(map(rationalize, [
[-11/12, 1/24, 0, 0, -1],
[ 17/24, -9/8, 1/24, 0 ,5],
[3/8, 9/8, -9/8, 1/24, -10],
[-5/24, -1/24, 9/8, -9/8, 10],
[1/24, 0, -1/24, 9/8, -5],
[0, 0, 0, -1/24, 1]
]))

and similarly for c. I've just checked and this returns the right solutions

@jksuom I cant remove any variable unfortunately.
Kalevi Suominen
@jksuom
Do not remove any variables. There are too few already. Remove one of the equations.
There are 6 rows and only 5 columns. One row could be omitted.
@jksuom I can't either, they are strict contraints in my solution space.
Kalevi Suominen
@jksuom
If there is a solution, then the same solution is obtained after one equation is omitted.
Kalevi Suominen
@jksuom
In [7]: linsolve(Pc[:5, :])
Out[7]:
{(1.12549831271091, 0.749071392102014, 1.16564957508439, 0.96140585029474, -0.
000620954287620833)}
@valglad Thanks. rationalize function worked (in fact I was using Fractions.from_float() before, this work better). But I could not apply your lambda function to a 2-d matrix, just 1-d.
I have created mrationalize = lambda x: [nsimplify(e) for e in [r for r in x]] that I do not like much the fell.
Doing map(rationalize, list_of_lists) inside the Matrix call (like I did above when creating PI) would probably be easiest

I have a problem with calculating eigenvectors for a simple 3x3 matrix in sympy. Trying to execute

q = symbols("q", positive=True)
m = Matrix([[-2, exp(-q), 1], [exp(q), -2, 1], [1, 1, -2]])
#m.eigenvects(simplify=True)
m.eigenvals()

results in very complicated expressions for the eigenvalues, trying to get the eigenvectors fails with a NotImplementedError. I expect the latter to be due to the former. However, when trying to do the same thing in Mathematica, I get much simpler expressions for everything. Is there some option / flag that I can set to have sympy compute these eigenvalues and eigenvectors? Is there currently some form of limitation within sympy that prevents me from doing this?

Het
@hetp111
what are all the prerequisite i need to have before contributing to opensource (sympy) ? I have a pretty good hold of python now. Github is very confusing for me.
Aaron Meurer
@asmeurer
It helps to understand git/github, but generally the best way to learn them is to use them
@adoa does Mathematica give a simpler expression for the eigenvalues?
Oh apparently they can be simplified if you call simplify
Het
@hetp111
@asmeurer Ty, do i need to learn git or github? is learning github enough? i have never done opensource before, should i start with sympy?
Aaron Meurer
@asmeurer
You need to learn the basics of git, to download the code and push it up to GitHub
@adoa it seems there is a simplify flag, but it doesn't work
Het
@hetp111
@asmeurer can't i do it manually through github website?
Kalevi Suominen
@jksuom
I suspect that there is not enough support for simplifying hyperbolic functions.
Malkhan Singh
@Malkhan52
Aaron Meurer
@asmeurer
Het
@hetp111
@Malkhan52 Thank ypu,
Aaron Meurer
@asmeurer
Het
@hetp111
@asmeurer I am using windows, will i have any problems? or do i need linux?
Aaron Meurer
@asmeurer
You can use windows.
You should use whichever you feel more comfortable with
Het
@hetp111
I am aiming for GSoC 2019, what is the procedure?
(sorry for these silly questions... I am trying my best ,i just need some guidance )
Malkhan Singh
@Malkhan52

I'm getting this error '''>>> expr = cos(x)

expr.subs(x-1,y)
cos(x)
expr = cos(x+1)
expr.subs(x-1,y)
cos(y + 2)
'''

where will it be in code?

Aaron Meurer
@asmeurer
Somewhere in the core
I guess it's because Add does the smart substitution but Symbol does not
The cosine is irrelevant. It's the same with x.subs(x - 1, y) vs. (x + 1).subs(x - 1, y)
Malkhan Singh
@Malkhan52
yeah that's why i'm searching subs()
Malkhan Singh
@Malkhan52
got that in basic.py
Aaron Meurer
@asmeurer
Subs is implemented on the different classes in _eval_subs
So Symbol._eval_subs should do the same thing as Add._eval_subs.
@asmeurer: I am somewhat confused regarding the simplification. Apparently there is a difference between simplify(a_list) and [simplify(element) for element in a_list]. With the latter, I get the same expressions as in Mathematica, the former does not simplify at all o_O. In any case: thanks for opening the issue :-)
Aaron Meurer
@asmeurer
I don't know why simplify returns a list unchanged. In general, though, you shouldn't rely on SymPy functions automatically vectorizing over lists
Ashwin Shenoy
@Ultra-Instinct-05
In the ellipse class, do hradius & vradius represent the semi-major axis length & semi-minor axis length respectively ?
It's a bit unclear in the docstring. I take it they are being represented as radius vectors from the origin along the major & minor axis ?
Akshit Srivastava
@akshitt
Hi I am Akshit, fresher in CS at Indian Institute of Technology, Bombay. Could someone guide me to a basic level bug to start with?
I have some basic level knowledge in Python, FLINT, NumPy
Het
@hetp111
@akshit1511 Hi, i am also new, i'll tell you what i'll be doing to get started, first of all i'll learn to use git and github and then https://github.com/sympy/sympy/wiki has a lot of information on how to start
@akshit1511 Tell me if i am missing something.
Akshit Srivastava
@akshitt