Re: [sympy] Why does subs() not work as expected - a related question ?

2021-04-07 Thread Oscar Benjamin
On Wed, 7 Apr 2021 at 15:27, 'Bruce Allen' via sympy wrote: > > On the downside, this makes it hard to work on a calculation which > requires large numbers of CPU hours and needs to be checkpointed and saved. > > Do you know of any workarounds? Or is there an alternative to > pickle.save() and

Re: [sympy] Why does subs() not work as expected - a related question ?

2021-04-07 Thread Oscar Benjamin
On Wed, 7 Apr 2021 at 14:43, 'Bruce Allen' via sympy wrote: > > David, Oscar, > > Thank you for your help. > > Oscar, the list 'u' was created in the course of a calculation, and > saved as a .pkl file. I then reloaded it and want to manipulate the > saved equations, of which u[7] is an example.

Re: [sympy] Why does subs() not work as expected - a related question ?

2021-04-07 Thread Oscar Benjamin
On Wed, 7 Apr 2021 at 14:55, Oscar Benjamin wrote: > > On Wed, 7 Apr 2021 at 14:43, 'Bruce Allen' via sympy > wrote: > > > > David, Oscar, > > > > Thank you for your help. > > > > Oscar, the list 'u' was created in the course of a calculation,

Re: [sympy] Profiling functions

2021-04-09 Thread Oscar Benjamin
On Fri, 9 Apr 2021 at 13:35, Shri Keshavinee wrote: > > I'm writing a proposal on the following project "Classical Mechanics: > Efficient Equation of Motion Generation with Python". As a first step, I > would like to contribute to it by finding where the functions are getting > slower through

[sympy] SymPy 1.8 released

2021-04-09 Thread Oscar Benjamin
to this release! - 彭于斌* - Thomas Aarholt* - Abhay_Dhiman* - Sachin Agarwal - Shubham Agrawal* - Sakirul Alam* - AlexCQY* - Abhinav Anand - Mohit Balwani - Soumi Bardhan* - Elias Basler* - Nijso Beishuizen* - Mark Bell* - Oscar Benjamin - Vaibhav Bhat* - Akshansh Bhatt* - Mohammed Bilal* - Ayush

Re: [sympy] Re: How to change/influence the dependent/independent variables from solve

2021-04-09 Thread Oscar Benjamin
On Sat, 10 Apr 2021 at 01:17, brandon...@gmail.com wrote: > > This almost worked for me as is. My desired dependent variables are [w, y]. > If I just pass those, then I don't get the solution for z=1. So I did two > solves: > > 1. Solve with my desired dependent variables (e.g. [w, y]) > 2.

Re: [sympy] Names that aren't unique

2021-04-12 Thread Oscar Benjamin
On Mon, 12 Apr 2021 at 12:03, 'Bruce Allen' via sympy wrote: > > Hi David, > > > However, my feeling is that some proportion of SymPy users will work > > interactively - in one scope - without defining any Python functions. So > > they might calculate a polynomial without regard to any

Re: [sympy] Names that aren't unique

2021-04-12 Thread Oscar Benjamin
On Mon, 12 Apr 2021 at 13:23, 'Bruce Allen' via sympy wrote: > > Hi Oscar, > > > It wouldn't be hard to make any new definition of a Symbol with the > > same name as a previously created symbol raise an error but it would > > break the assumption that it is okay to define a symbol that is only >

Re: [sympy] Sympy apart() with complex variables

2021-04-12 Thread Oscar Benjamin
On Mon, 12 Apr 2021 at 14:14, s5s wrote: > > I am trying to do partial fractions involving a complex variable. I was > wondering if one can do this with sympy apart function. Below is the code to > attempt this but it does not return a decomposed expression. > > > from sympy import * > x, y, z

[sympy] Automated releases and the authors file

2021-04-10 Thread Oscar Benjamin
Hi all, Over the past couple of releases I've done a bunch of work to try and automate the release processes. This is based on the release script that Aaron originally wrote although I have rewritten it so that it is no longer based on rever and is easier to debug locally. Instead of being based

Re: [sympy] Automated releases and the authors file

2021-04-11 Thread Oscar Benjamin
On Sun, 11 Apr 2021 at 07:51, Aaron Meurer wrote: > > On Sat, Apr 10, 2021 at 3:27 AM Oscar Benjamin > wrote: > > > > That script builds all of the release artifacts and then installs the > > SymPy wheel in a bunch of Python versions and runs sympy.test() in >

Re: [sympy] Why does subs() not work as expected - a related question ?

2021-04-07 Thread Oscar Benjamin
On Wed, 7 Apr 2021 at 12:24, David Bailey wrote: > > On 07/04/2021 11:59, Oscar Benjamin wrote: > > On Wed, 7 Apr 2021 at 10:28, 'Bruce Allen' via sympy > > wrote: > >> > >> Why does subs() work in one case, and not in the other? > > I'm not sure. I can'

Re: [sympy] Names that aren't unique

2021-04-11 Thread Oscar Benjamin
On Sun, 11 Apr 2021 at 16:26, David Bailey wrote: > > Dear group, > > Recently Bruce Allen discussed a problem that he had after pickling a > list. However, I think he revealed a deeper problem that is nothing to > do with pickling. > > w=sin(Symbol("x", positive=True))*cos(Symbol("x")) > >

Re: [sympy] GSoC Idea Discussion

2021-02-12 Thread Oscar Benjamin
On Fri, 12 Feb 2021 at 13:13, cheshta babbar wrote: > > While going through ideas , I found this idea , Multivariate polynomials and > factorization, using Gao's partial differential equations approach, very > appealing . I wanted to know if this idea is outdated or is being considered > for

Re: [sympy] Solvers improvement roadmap

2021-02-15 Thread Oscar Benjamin
On Mon, 15 Feb 2021 at 21:01, Sayandip Halder wrote: > > Hello! I hope to discuss the possible action plan for improving solvers and > solveset here. The GitHub wiki is old and needs to be updated. Please do update it! > The following are what I can think of right now. Also, I want to know

Re: [sympy] Solvers improvement roadmap

2021-02-16 Thread Oscar Benjamin
On Tue, 16 Feb 2021 at 07:34, Sayandip Halder wrote: > > On Tuesday, February 16, 2021 at 5:21:21 AM UTC+5:30 Oscar wrote: >> >> It might make more sense to fix these things in solveset and then make >> use of solveset in solve rather than trying to fix things in solve >> itself. > > There are

Re: [sympy] Linear Inequality Simplifier

2021-02-18 Thread Oscar Benjamin
On Thu, 18 Feb 2021 at 11:32, Oscar Gustafsson wrote: > > After currently using Mathematica for similar things, I would just like to > encourage you to provide some nice method to simplify constraints of > piecewise functions using your simplifier, including additional constraints > on the

Re: [sympy] Linear Inequality Simplifier

2021-02-17 Thread Oscar Benjamin
‪On Tue, 16 Feb 2021 at 17:01, ‫אוריאל מליחי‬‎ wrote:‬ > > For your first question, I intend to create a new function that when given a > set of linear inequalities and a target inequality it would output True if > the target is implied by this set, False if it is not, and Unknown otherwise. >

Re: [sympy] GSoC Organization Applications Open

2021-02-17 Thread Oscar Benjamin
On Wed, 17 Feb 2021 at 22:28, Aaron Meurer wrote: > > I have updated our application for this year > https://github.com/sympy/sympy/wiki/GSoC-2021-Organization-Application. > If anyone has any comments on it, let me know. If you edit that page, > please make sure to ping me, as I will need to

Re: [sympy] GSoC Organization Applications Open

2021-02-17 Thread Oscar Benjamin
On Wed, 17 Feb 2021 at 23:24, Aaron Meurer wrote: > > On Wed, Feb 17, 2021 at 3:39 PM Oscar Benjamin > wrote: > > > > On Wed, 17 Feb 2021 at 22:28, Aaron Meurer wrote: > > > > > > I have updated our application for this year > > > https://gi

Re: [sympy] Unable to see my message

2021-02-04 Thread Oscar Benjamin
On Thu, 4 Feb 2021 at 18:48, mayank gaur wrote: > > Now I have posted several times, so now I will not be moderated right? When a message is in moderation the person who allows the message through needs to tick a box to say "let this message through *and* all future messages from the same

Re: [sympy] We need some community input on the desired behavior of an equation class...

2021-02-06 Thread Oscar Benjamin
Hi Jonathan, I just looked at the survey and I think the question as posed is something of a distraction from the real issue. If I set that as an exercise to students then it would be a trick question. The correct answer is "it depends on information not given" :) I've set this to students many

Re: [sympy] [SymPy] Is this a bug or a feature?

2021-02-07 Thread Oscar Benjamin
On Sun, 7 Feb 2021 at 00:07, David Bailey wrote: > > Dear Group, > > While thinking about Jonathon's question, I came across this oddity: > > x=symbols('x') > > f=symbols('f',cls=Function) > > diff(f,x) > > 1 > > Why 1? I think I would have expected it to generate a TypeError, just like > f+x

Re: [sympy] AttributeError: 'LessThan' object has no attribute 'as_poly'

2021-02-12 Thread Oscar Benjamin
On Fri, 12 Feb 2021 at 14:33, mpierro3 wrote: > > I tried using nsolve like you mentioned and got it to work. I then tried > again, except adding a bound, and it produced the following error: > > from sympy import * > > M_3a = symbols('M_3a') > eqn = nsolve([Eq(0.56*(0.33*M_3a**2 +

Re: [sympy] Is piece of code in StrPrinter ever reachable?

2021-02-13 Thread Oscar Benjamin
On Sat, 13 Feb 2021 at 12:34, Paul Royik wrote: > > In StrPrinter there is > def emptyPrinter(self, expr): > if isinstance(expr, str): > return expr > elif isinstance(expr, Basic): > return repr(expr) > else: > return str(expr) > > Is `return repr(expr)` reachable if `StrPrinter` defines

Re: [sympy] GSoc Ideas

2021-02-13 Thread Oscar Benjamin
Yes, that project is definitely still relevant and needed. I would supervise it if the project goes ahead. The description links to a Github issue which describes the idea also has links to other pull requests where some of this work has been done including some unifinished PRs. Oscar On Sat,

Re: [sympy] refine: Doc error in refine.py

2021-02-13 Thread Oscar Benjamin
You seem to presume that someone reading your message understands what it is that you are doing. It isn't clear what question you are asking or what it is you are doing that does not work. Please take the time to explain yourself more clearly. On Sat, 13 Feb 2021 at 03:02, Ayush Bisht wrote: >

Re: [sympy] Gsoc Idea Discussion

2021-02-20 Thread Oscar Benjamin
On Sat, 20 Feb 2021 at 07:40, Kartik Sethi wrote: > > Hey, my name is Kartik Sethi. I am interested in taking part in Gsoc this > year. I have been contributing to sympy for the past couple of months. Hi Kartik, > I have come up with a few ideas for Gsoc. I am focusing primarily on the >

Re: [sympy] Gsoc Idea Discussion

2021-02-20 Thread Oscar Benjamin
On Sat, 20 Feb 2021 at 11:33, Oscar Benjamin wrote: > This is not listed anywhere on the ideas page or really documented > anywhere yet but there is a new mostly internal implementation of > matrices in the sympy.polys.matrices module. This new implementation > is called DomainMatrix

Re: [sympy] Gsoc Idea Discussion

2021-02-20 Thread Oscar Benjamin
On Sat, 20 Feb 2021 at 14:28, Kartik Sethi wrote: > > Thanks Oscar Benjamin for the wonderful suggestions. I will start looking > into algorithms and methods that have not been implemented for the domain > Matrix class yet. > I would be happy to help in building Domain Matrix cl

Re: [sympy] AttributeError: 'LessThan' object has no attribute 'as_poly'

2021-02-11 Thread Oscar Benjamin
On Thu, 11 Feb 2021 at 22:54, mpierro3 wrote: > > Yes, thank you. The rest is as follows: > > from sympy import * > import numpy as np > > M_e = 1 > gamma_1 = 1.667 > gamma_4 = 1.667 > MW_1 = 39.948 > MW_4 = 4.0026 > D_4 = 5 / 39.37 > D_1 = 3 / 39.37 > A_4 = (np.pi / 4) * Pow(D_4, 2) > A_1 =

Re: [sympy] GSoc Ideas

2021-02-13 Thread Oscar Benjamin
On Sat, 13 Feb 2021 at 16:59, Suryam Kalra wrote: > I would be glad to take the project ahead from here under your supervision > . I have read the description and the other related PR's. Should I start by > finishing those unmerged PR's or shall I create a new PR for each of the > solver methods

Re: [sympy] AttributeError: 'LessThan' object has no attribute 'as_poly'

2021-02-11 Thread Oscar Benjamin
On Thu, 11 Feb 2021 at 21:42, mpierro3 wrote: > > New to Sympy. I am trying to solve an equation for a variable with given > boundaries and I am reviving the following error: > > M_3a = symbols('M_3a', positive=true, nonzero=true) > eqn = nonlinsolve([Eq((M_e / M_3a) * Pow((2 + (gamma_4 - 1) *

Re: [sympy] Is it possible to get polynomial multiplyiers?

2021-02-22 Thread Oscar Benjamin
On Mon, 22 Feb 2021 at 15:39, Michał Pawłowski wrote: > > Thank you so much. But is it also possible to get formula from such array? > > I.e. > getFormula("[2, 0, 2]") > 2*x**2+2 Why is the "array" given as a string (in quotes)? If you have a list of the coefficients then you can construct the

[sympy] Supporting older Python versions

2021-08-17 Thread Oscar Benjamin
Hi all, This issue discusses SymPy's policy for supporting older Python versions: https://github.com/sympy/sympy/issues/21884 Currently at the point SymPy is released it is expected to support all Python versions that are not yet EOL:

[sympy] Cylindrical polar unit vectors

2021-08-13 Thread Oscar Benjamin
Hi all, Does anyone know how to use cylindrical polars with SymPy's vector module? I just saw this SO question here: https://stackoverflow.com/questions/68767945/how-do-i-work-with-vectors-in-cylindrical-and-spherical-systems-with-sympy We can define a coordinate system like this: In [25]:

Re: [sympy] Cylindrical polar unit vectors

2021-08-13 Thread Oscar Benjamin
On Fri, 13 Aug 2021 at 14:06, Davide Sandona' wrote: >> >> Although the scalars are r, theta and z the unit vectors are i, j and >> k. How would I get vectors like e_r and e_theta to work with those? > > As of now, that's not possible. It is hard-coded and it only represents cartesian unit

Re: [sympy] Re: Keeping the AUTHORS file up to date

2021-08-24 Thread Oscar Benjamin
I've realised another problem with this which is that if a PR was opened before this change then it is not enough to merge with the latest master and rerun the bin/authors_update.py script. The script updates the AUTHORS file in the order of the commits but that order is different when merging

Re: [sympy] Re: Keeping the AUTHORS file up to date

2021-08-24 Thread Oscar Benjamin
On Wed, 25 Aug 2021 at 01:07, Aaron Meurer wrote: > > Are you simulating having the PR run in a "merged with master" state > when running the authors script on CI? Maybe we should update the > script itself so that it can do this. It's not a simulation. When a PR is pushed GitHub actions will

Re: [sympy] Extended sympy interface in SageMath 9.4

2021-08-24 Thread Oscar Benjamin
Hi Matthias, That's excellent. Is there a way to run those same examples from Python or do they need to be run from the SageMath application? Also since SymPy switched from Travis to GitHub Actions for CI I haven't figured out how to make a setup that tests the latest master SymPy as part of

Re: [sympy] Extended sympy interface in SageMath 9.4

2021-08-24 Thread Oscar Benjamin
On Wed, 25 Aug 2021 at 01:54, Matthias Köppe wrote: > On Tuesday, August 24, 2021 at 5:21:02 PM UTC-7 Oscar wrote: > >> >> Also since SymPy switched from Travis to GitHub Actions for CI I >> haven't figured out how to make a setup that tests the latest master >> SymPy as part of Sage. Currently

Re: [sympy] Extended sympy interface in SageMath 9.4

2021-08-24 Thread Oscar Benjamin
I just downloaded SageMath 9.4 on OSX from here: https://github.com/3-manifolds/Sage_macOS/releases I ran most of the examples fine but there was a problem with the Weyl group example: sage: from sage.all import WeylGroup sage: W = WeylGroup(["A",1,1]) Assertion failed: ((STATE_SLOTS_SIZE -

Re: [sympy] Re: Keeping the AUTHORS file up to date

2021-08-25 Thread Oscar Benjamin
On Wed, 25 Aug 2021 at 09:22, Aaron Meurer wrote: > On Tue, Aug 24, 2021 at 6:30 PM Oscar Benjamin > wrote: > > > > On Wed, 25 Aug 2021 at 01:07, Aaron Meurer wrote: > > > > > > Are you simulating having the PR run in a "merged with master" sta

Re: [sympy] Re: Keeping the AUTHORS file up to date

2021-08-25 Thread Oscar Benjamin
On Wed, 25 Aug 2021 at 11:16, prathamesh bhole wrote: > @Oscar the test is still failing with the same error > Look carefully at the diff and check every changed line: https://github.com/sympy/sympy/pull/21797/files There are changes in the authors file from the merge conflict not being fully

Re: [sympy] Re: Keeping the AUTHORS file up to date

2021-08-25 Thread Oscar Benjamin
On Wed, 25 Aug 2021 at 09:43, prathamesh bhole wrote: > Today I made a commit to my pull request which was created way before this > change. Now the code quality test is failing because of exit code 1 by > authors_update.py. Can you please help, so that I can clear the tests Hi Prathamesh,

Re: [sympy] Re: Keeping the AUTHORS file up to date

2021-08-30 Thread Oscar Benjamin
On Mon, 30 Aug 2021 at 22:44, Aaron Meurer wrote: > On Mon, Aug 30, 2021 at 1:46 PM Oscar wrote: > > > > Okay so here's a proposal: > > > > List all authors in .mailmap which can be sorted alphabetically by > primary name/email combination. New contributors have to add their name and > email

[sympy] Keeping the AUTHORS file up to date

2021-08-21 Thread Oscar Benjamin
Hi all, I just merged a PR to update the AUTHORS file: https://github.com/sympy/sympy/pull/21881 This is something that I've needed to do before each release and is quite time consuming so in that PR I also added the scripts that check the AUTHORS file to the CI tests. Now any PR that does not

Re: [sympy] when do imagesets and intersections resolve to finite sets?

2021-09-01 Thread Oscar Benjamin
On Wed, 1 Sept 2021 at 21:16, Aaron Meurer wrote: > I think there's a few things that could be done here: > > - Make top-level intersection work by inverting the function in the > image set. This needs to be done carefully, however. > - Make the inner intersection resolve to {n*pi | n in Z+} > -

Re: [sympy] when do imagesets and intersections resolve to finite sets?

2021-09-01 Thread Oscar Benjamin
On Wed, 1 Sept 2021 at 20:06, Charlie Hartman < charlesmichaelhart...@gmail.com> wrote: > When do image sets resolve to finite sets when intersected with an > interval? For example: > > I'm seeing that this looks great: > > solveset(sin(x),domain = Reals).intersect(Interval(2,50)) >

Re: [sympy] increasing weight of n-size subsets algorithm

2021-09-03 Thread Oscar Benjamin
On Fri, 3 Sept 2021 at 21:48, Aaron Meurer wrote: > On Fri, Sep 3, 2021 at 12:10 PM Chris Smith wrote: > > > > I also realize that there is an aspect to solving systems of non-linear > > equations that I have not fully appreciated. I would normally think of > > working manually through a

[sympy] SymPy 1.9 release

2021-09-09 Thread Oscar Benjamin
Hi all, I want to put out the SymPy 1.9 release ASAP. There have been a large number of changes since 1.8: https://github.com/sympy/sympy/wiki/Release-Notes-for-1.9 I intend to spend a little time editing those notes before the release but I don't have time to go through all of them.

Re: [sympy] Re: Difference between f.is_number and bool(f.free_symbols)

2021-09-13 Thread Oscar Benjamin
Think about things that are literally not numbers: In [9]: Interval(1, 2).is_number Out[9]: False In [10]: ImmutableMatrix([[1, 2], [3, 4]]).is_number Out[10]: False On Mon, 13 Sept 2021 at 13:00, Chris Smith wrote: > To confirm, if you mean that it is free from any Symbol (free or bound) >

Re: [sympy] Solution of nonlinear ODE

2021-09-10 Thread Oscar Benjamin
On Fri, 10 Sept 2021 at 16:40, Nicolas Guarin wrote: > > > I have the following ODE > > z'' = 1 - z² > > That has as solutions z=tanh(C + t) and z=coth(C + t), depending on the initial condition being greater or less than 1. When I use dsolve I get the latter > > > from sympy import * >

Re: [sympy] Lambdify - performance using sympy as the module

2021-09-15 Thread Oscar Benjamin
On Tue, 14 Sept 2021 at 23:12, sandona...@gmail.com < sandona.dav...@gmail.com> wrote: > Hello, > let's say I'd like to numerically evaluate a single sympy function over an > array using sympy as the module. Curiously, passing in regular Python's > float numbers makes the evaluation much faster

Re: [sympy] Lambdify - performance using sympy as the module

2021-09-15 Thread Oscar Benjamin
On Wed, 15 Sept 2021 at 13:41, Oscar Benjamin wrote: > On Tue, 14 Sept 2021 at 23:12, sandona...@gmail.com < > sandona.dav...@gmail.com> wrote: > >> Hello, >> let's say I'd like to numerically evaluate a single sympy function over >> an array using sympy

Re: [sympy] Re: Difference between f.is_number and bool(f.free_symbols)

2021-09-13 Thread Oscar Benjamin
There is also is_Number vs is_number. There is also is_comparable which means can be evalf'd to a real number. There should be a glossary of all the different attributes that are used throughout sympy: In [*1*]: *for* name *in* dir(x): ...: *if* name.startswith('is'): ...:

Re: [sympy] SymPy 1.9 release

2021-09-09 Thread Oscar Benjamin
" On Thu, 9 Sept 2021 at 23:12, Aaron Meurer wrote: > > On Thu, Sep 9, 2021 at 2:25 PM Oscar Benjamin > wrote: > > > > Anyone is welcome to have an opinion on the printing of Heaviside functions: > > https://github.com/sympy/sympy/issues/21945 > > In the a

Re: [sympy] SymPy 1.9 release

2021-09-09 Thread Oscar Benjamin
On Fri, 10 Sept 2021 at 02:02, Aaron Meurer wrote: > > On Thu, Sep 9, 2021 at 6:07 PM Oscar Benjamin > wrote: > > > > " > > On Thu, 9 Sept 2021 at 23:12, Aaron Meurer wrote: > > > > > > On Thu, Sep 9, 2021 at 2:25 PM Oscar Benjamin >

Re: [sympy] SymPy 1.9 release

2021-09-09 Thread Oscar Benjamin
On Fri, 10 Sept 2021 at 02:02, Aaron Meurer wrote: > > On Thu, Sep 9, 2021 at 6:07 PM Oscar Benjamin > wrote: > > > > You can see another example in the release notes where I wrote: > > > > - The (internal) PolyMatrix class has been changed and now requires

Re: [sympy] Request for discussion on kind system

2021-08-02 Thread Oscar Benjamin
On Mon, 2 Aug 2021 at 22:24, Jonathan Gutow wrote: > > Is there a discussion page/tab attached to github wikis? If so, how do I find > it? I do have some questions/thoughts on this. Best place to start the discussion is here. Then the wiki page can be updated following the discussion (it can

Re: [sympy] Q: codegen

2021-08-03 Thread Oscar Benjamin
The issue is here: https://github.com/sympy/sympy/issues/21815 On Tue, 3 Aug 2021 at 04:03, Audrius-St wrote: > > Aaron Meurer . Thank you for your reply and suggestions > > I've now opened an issue in the SymPy issue tracker regarding this matter. > > On Monday, August 2, 2021 at 3:42:55 PM

Re: [sympy] Basic plotting problem

2021-08-07 Thread Oscar Benjamin
Maybe that error message could be improved... On Sat, 7 Aug 2021 at 11:44, Davide Sandona' wrote: > > Hello Davide, > > In order to use Qt5 you need to install the following dependency: > > pip install PyQt5 > > After that it should work fine! > > Davide. > > > Il giorno sab 7 ago 2021 alle ore

Re: [sympy] Using TensorFlow with SymPy

2021-07-29 Thread Oscar Benjamin
de > that make use of them are mostly independent. > > Aaron Meurer > > On Wed, Jul 28, 2021 at 6:36 PM Oscar Benjamin > wrote: > > > > I should say that since migrating to GitHub Actions we no longer test > > tensorflow in CI. Someone could add that back in bu

Re: [sympy] Using TensorFlow with SymPy

2021-07-28 Thread Oscar Benjamin
I should say that since migrating to GitHub Actions we no longer test tensorflow in CI. Someone could add that back in but I removed it because it was messing up the other dependencies in the optional dependencies job. Possibly it should have its own separate job. On Thu, 29 Jul 2021 at 01:10,

Re: [sympy] evalf and Integer

2021-08-01 Thread Oscar Benjamin
On Sun, 1 Aug 2021 at 11:02, Paul Royik wrote: > > For a given number, I want to approximate it to 5 decimal points. > If a number of digits in the integer part is greater than 7, then round to 7 > significant figures. What does "round" mean here? Do you want to round to nearest, or round

Re: [sympy] Request for discussion on kind system

2021-08-03 Thread Oscar Benjamin
On Tue, 3 Aug 2021 at 23:17, Aaron Meurer wrote: > > > Finally, it must be stressed that Kind instance must be singleton, and we > > should not implement complicated set theory logic here. This is to > > facilitate the determination of the kind of associative operation with very > > large

Re: [sympy] Version 1.9 improvement, and indexing a variable

2021-10-20 Thread Oscar Benjamin
On Tue, 19 Oct 2021 at 20:55, David Bailey wrote: > > Dear group, > > First I would like to say how good it was to discover that the online > SymPy documentation now supports copy/paste operations without > constantly switching to SymPyLive. I think this will make the > documentation considerably

[sympy] SymPy 1.9 is now released!

2021-10-08 Thread Oscar Benjamin
Balasubramanyam* - Mohit Balwani - Oscar Benjamin - Akshansh Bhatt - Remco de Boer* - Matt Bogosian* - Francesco Bonazzi - Zach Carmichael* - Kaustubh Chaudhari - Michael Chu* - Alexander Cockburn* - Lorenzo Contento - Miguel Torres Costa* - Björn Dahlgren - Jonathan Daniel* - Gaurav Dhingra - dimasvq

Re: [sympy] emulate a lambda

2021-09-28 Thread Oscar Benjamin
On Tue, 28 Sept 2021 at 10:22, Aaron Meurer wrote: > Why not just Lambda(x, Eq(x, 1))? > > I don't think redefining __eq__ to return a symbolic result is a good > idea. You're liable to get a lot of "TypeError: cannot determine truth > value of Relational" errors from using that object. There

Re: [sympy] emulate a lambda

2021-09-28 Thread Oscar Benjamin
On Tue, 28 Sept 2021 at 04:13, Chris Smith wrote: > I would like to emulate something like this with Basic objects. I am > drawing a blank on how that might be done. Does anyone have any ideas? > >>> f=lambda x: x==1 > >>> f(1) > True > >>> Lambda(x, f(x))(1) # doesn't work > False > Something

Re: [sympy] emulate a lambda

2021-09-30 Thread Oscar Benjamin
I don't know what the OracleGate class is for but I'm pretty sure it's not very useful and it's certainly poorly designed. I'd rather just delete it than try to come up with hacks to make it work. If someone wants to maintain the quantum module then that's great. Until then we shouldn't allow

Re: [sympy] SymPEP: make MatchPy a SymPy dependency

2021-10-06 Thread Oscar Benjamin
On Wed, 6 Oct 2021 at 08:52, Jason Moore wrote: > I think we should relax our stringent no dependency stance for pure python > dependencies. Pure python dependency management is essentially solved now > for the python ecosystem with the various package managers. This dependency > would also be

Re: [sympy] SymPEP: make MatchPy a SymPy dependency

2021-10-06 Thread Oscar Benjamin
On Wed, 6 Oct 2021 at 23:36, Francesco Bonazzi wrote: > > It would be curious if we could rewrite SymPy's core using MatchPy. I > mean, there have been some development with a lot of code developed to > properly handle the subclassing (or .kind property, that has been > introduced now) the *Add*

[sympy] SymPY 1.9rc1 released

2021-09-20 Thread Oscar Benjamin
Hi all, I've just put up SymPy 1.9rc1 which is the release candidate for SymPy 1.9. This is a pre-release intended for testing rather than production use. You can find the release files on GitHub here: https://github.com/sympy/sympy/releases/tag/sympy-1.9rc1 The release notes for 1.9 are here

Re: [sympy] Re: Idiomatic use of derivative replacements

2021-10-05 Thread Oscar Benjamin
I think that the basic problem is just this (which is a bug): In [86]: expr = f(y(t)).diff(y(t)) In [87]: print(expr) Derivative(f(y(t)), y(t)) In [88]: print(expr.subs(t, 0)) Derivative(f(y(0)), y(0)) There is code to detect this in some cases but it doesn't work for your example. If y(t) was

Re: [sympy] SymPEP: make MatchPy a SymPy dependency

2021-10-05 Thread Oscar Benjamin
On Tue, 5 Oct 2021 at 07:48, Francesco Bonazzi wrote: > Hi everyone, > > I have written a draft for a SymPEP (SymPy enhancement proposal) to > include MatchPy as a dependency of SymPy. > > https://github.com/sympy/SymPEPs/pull/3 > For those less familiar to github you need to click the "files

[sympy] 1.9 release is delayed

2021-09-25 Thread Oscar Benjamin
I was intending to put out the release this weekend i.e. nowish but it is held up by these issues: https://github.com/sympy/sympy/issues/14933 https://github.com/sympy/sympy/issues/21494 Christopher Smith, Oscar Gustafsson and I have been investigating fixes but nothing is confirmed yet so the

[sympy] Re: 1.9 release is delayed

2021-09-27 Thread Oscar Benjamin
On Sat, 25 Sept 2021 at 23:55, Oscar Benjamin wrote: > I was intending to put out the release this weekend i.e. nowish but it is > held up by these issues: > > https://github.com/sympy/sympy/issues/14933 > https://github.com/sympy/sympy/issues/21494 > > Christopher Smith, O

Re: [sympy] How to create a vector function in sympy with actual vectors, that can be transformed between coordinate systems?

2021-10-25 Thread Oscar Benjamin
I don't see the exception that you showed: In [9]: O = CoordSys3D('O') ...: r = O.x*O.i + O.y*O.j + O.z*O.k ...: dot(r,r) # this works ...: O.x**2 + O.y**2 + O.z**2 ...: from sympy import symbols ...: a, b, c = symbols("a, b, c") ...: p = a*O.i +b*O.j + c*O.k ...: p ...:

Re: [sympy] Re: Another problem with Sympy's solver.

2021-12-21 Thread Oscar Benjamin
It solves a lot quicker using check=False: In [1]: from sympy import symbols, solve ...: f, q, r, h, i, g = map(symbols, 'fqrhig') ...: e1 = -4*f**2*q - 5*f**2*r - 5*f*r + 4*q**2 + 3*q*r ...: e2 = -4*f**2*h*q - 4*f**2*i*q - 8*f*g*q - 10*f*g*r - 8*f*h*q - 20*f*h*r - 10*f*i*r + 8*f*q**2/5

Re: [sympy] Re: Another problem with Sympy's solver.

2021-12-23 Thread Oscar Benjamin
On Thu, 23 Dec 2021 at 08:33, emanuel.c...@gmail.com wrote: > > Le mardi 21 décembre 2021 à 16:41:06 UTC+1, smi...@gmail.com a écrit : >> >> e1 is quadratic only if f and e2 is linear in g and quadratic in f...it is a >> pretty simple system but there are lots of extraneous symbols. Even paring

Re: [sympy] Poly expands the expression

2021-11-20 Thread Oscar Benjamin
On Sat, 20 Nov 2021 at 07:59, Paul Royik wrote: > > Poly((x-1)**2, x) returns Poly(x**2-2x+1, x). > Is it possbile to retain the expression as it is, i.e. (x-1)**2 ? It is possible to retain the unexpanded expression: don't convert to Poly! The internal representation of Poly cannot represent

Re: [sympy] Self-introduction--Interest: Making SymPy easier for new users via documentation

2021-11-28 Thread Oscar Benjamin
On Sun, 28 Nov 2021 at 11:09, Jeremy Monat wrote: > > Hello, I'm Jeremy Monat. My primary interest in helping with SymPy > development is making it easier for new users to use SymPy by improving the > documentation. It is plenty powerful for the application I needed; learning > how to use

Re: [sympy] Why Poly((x-6)**100*(2-x)**9000,x) is very slow?

2021-11-22 Thread Oscar Benjamin
On Mon, 22 Nov 2021 at 09:18, Paul Royik wrote: > > `(x-2)**9000` takes much time, but `(x-6)**100*(2-x)**9000` takes forever. It's slow because it involves explicit coefficient calculations with very large polynomials. Note that if you don't use Poly and you don't expand the expressions then

Re: [sympy] Open issues for solving Diophantine equations

2021-10-30 Thread Oscar Benjamin
Issues and pull requests in the sympy repo are labelled by topic which is very useful in this situation. There is a label for anything related to diophantine which shows 20 open issues: https://github.com/sympy/sympy/labels/solvers.diophantine Oscar On Sat, 30 Oct 2021 at 21:42, Raphael C

Re: [sympy] How to factorize a noncommuting symbol?

2021-11-02 Thread Oscar Benjamin
On Tue, 2 Nov 2021 at 15:38, Gerardo Suarez wrote: > > So basically I have a lot of objects like H: > > a= sympy.Symbol('a', commutative=False) > beta= sympy.Symbol('beta',commutative=True) > gamma= sympy.Symbol('\gamma',commutative=True) > ad = sympy.Symbol('a^{\dagger}', commutative=False) >

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

2021-11-02 Thread Oscar Benjamin
On Mon, 1 Nov 2021 at 17:17, Andreas Schuldei wrote: > . > Oscar schrieb am Sonntag, 31. Oktober 2021 um 23:50:57 UTC+1: >> >> On Sun, 31 Oct 2021 at 20:10, Andreas Schuldei wrote: >> > >> > Thank you very much, Oscar. That is excellent feedback. >> > >> > You are right. The solver didn't find a

Re: [sympy] Melvin: An algorithm for automatically designing new quantum experiments using SymPy

2021-10-27 Thread Oscar Benjamin
Hi Xuemei, I'm sure we can list this on the webpage. I think you can just open a PR to add this. Oscar On Fri, 15 Oct 2021 at 18:42, Xuemei Gu wrote: > > Hello all, > > I worked on quantum physics and studied computer science. Using AI to > discover new quantum experiments is one of my

Re: [sympy] Doit method in concrete

2021-11-06 Thread Oscar Benjamin
On Sat, 6 Nov 2021 at 11:58, Anderson Bhat wrote: > > Hello guys , I am working on couple of Pr's extending the functionality of > the doit method in the concrete module , I noticed that one inconsistency > leads to couple of errors . Product(1, (n, 1, oo)).doit() returns 1 and 1**oo > returns

Re: [sympy] zero in sympy: Issue #22425

2021-11-04 Thread Oscar Benjamin
On Thu, 4 Nov 2021 at 23:13, Zoufiné Lauer-Baré wrote: > > Dear all, Hi and thanks for reporting this. > I just created an issue in the SymPy github project (Zero in SymPy > #22425 opened 2 minutes ago by zolabar). Here goes the content, may be > someone knows this topic and there are already

Re: [sympy] Doit method in concrete

2021-11-06 Thread Oscar Benjamin
On Sat, 6 Nov 2021 at 16:55, Peter Stahlecker wrote: > > I am just a hobby mathematician, but it seems to me like this: > > 1^oo := lim(1^n) = lim(1) = 1. > > The other 'limits' seem to me to be an inadmissible 'exchange' of limits: > 1 != (1 + 1/n) for any finite n When implementing the

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

2021-10-31 Thread Oscar Benjamin
On Sat, 30 Oct 2021 at 15:41, Andreas Schuldei wrote: > > > Thank you, that seems to have been the issue. I split up the vectors into > their coefficients and entered those vector components as separate equations. > Now the python process is humming along at 1.8Gbyte memory and low to medium >

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

2021-10-31 Thread Oscar Benjamin
On Sun, 31 Oct 2021 at 20:10, Andreas Schuldei wrote: > > Thank you very much, Oscar. That is excellent feedback. > > You are right. The solver didn't find a solution after about 18h of running. > For the solver to work, I had broken up my vector equation into vector > component equations, and

Re: [sympy] How to create a vector function in sympy with actual vectors, that can be transformed between coordinate systems?

2021-10-23 Thread Oscar Benjamin
A SymPy Vector is constructed algebraically from the unit vectors i, j and k of the coordinate system. For a vector field you also use the coordinate system base scalars x, y and z. In [11]: from sympy.vector import CoordSys3D, dot In [12]: O = CoordSys3D('O') In [13]: r = O.x*O.i + O.y*O.j +

Re: [sympy] Self-introduction--Interest: Making SymPy easier for new users via documentation

2021-11-30 Thread Oscar Benjamin
On Tue, 30 Nov 2021 at 02:37, Jeremy Monat wrote: > > Improve the tutorial solvers page's order by first explaining how to solve > equations and inequalities, particularly explaining the strengths of > different solvers, when to use each, and giving examples of each (as > suggested by Oscar)

Re: [sympy] How to factor out a constant?

2021-12-14 Thread Oscar Benjamin
On Tue, 14 Dec 2021 at 13:53, Paul Royik wrote: > > In any expression, how to factor a constant (that doesn't depend on certain > symbols)? > > For example (-2xy+2).factor_function(x,y) should give 2, (1-xy). > > as_coeff_mul doesn't work. Maybe factor_terms is what you want: In [50]:

Re: [sympy] Common shape with lambdify on Matrix expression

2021-07-20 Thread Oscar Benjamin
On Tue, 20 Jul 2021 at 12:28, Hanno Klemm wrote: > > Well, this is what you currently get. Personally, I am not very concerned by > the ordering. However, since the shape of the vector in sympy is (2,1), > changing the shape on the “internal” axis to (1,2) seems like an odd choice. That was a

Re: [sympy] questions about orientation of ParametricRegions

2021-07-19 Thread Oscar Benjamin
Hi Joe, The vector integration functionality is relatively new. It was added last year I think in a GSOC project. It's possible that this is a bug or a corner case that wasn't fully considered. Have you looked into the code at all? Oscar On Sun, 18 Jul 2021 at 00:19, Joe Heafner wrote: > >

Re: [sympy] Common shape with lambdify on Matrix expression

2021-07-19 Thread Oscar Benjamin
This is a common problem that comes up. There are two separate issues. The first is that lamdifying a constant function does not respect the shape of array inputs which is a longstanding still open issue: https://github.com/sympy/sympy/issues/5642 The second issue (which also makes it hard to

Re: [sympy] Common shape with lambdify on Matrix expression

2021-07-19 Thread Oscar Benjamin
t I did overlook (which does not appear > to be the case). > > Currently, I use a workaround and define a spurious symbol ‘null’, instead of > using 0, which I then evaluate with an array of zeros of the appropriate > shape. It isn’t particularly beautiful but works in my case. &g

<    1   2   3   4   5   6   7   8   >