Re: [sage-devel] bug in matrix solve over inexact coefficients?

2020-05-22 Thread Dave Morris
>From the documentation at https://www.mathworks.com/help/matlab/ref/mldivide.html, it appears to me that MATLAB gives a warning: "Warning: Matrix is close to singular or badly scaled. Results may be inaccurate." That seems to me to be a better default behavior than what sage is doing now, but

Re: [sage-devel] bug in matrix solve over inexact coefficients?

2020-05-22 Thread Michael Orlitzky
On 5/22/20 6:40 PM, AlexGhitza wrote: > I would also argue that, despite the validity of the arguments regarding > inexact rings, this is a change in behavior that would have benefited > from a deprecation warning for a short while. We were pretty careful not to break anything in the sage library.

Re: [sage-devel] Re: Could someone point me to the code that supports math on symbolic equations?

2020-05-22 Thread Samuel Lelievre
Le samedi 23 mai 2020 02:14:58 UTC+2, Dima: > > Conda does have Sagemath available. > Not 100% sure how it works on Windows, though. One can install SageMath from Conda on Linux and macOS. Not on Windows. -- You received this message because you are subscribed to the Google Groups "sage-devel"

Re: [sage-devel] Re: Could someone point me to the code that supports math on symbolic equations?

2020-05-22 Thread dimpase
On Fri, May 22, 2020 at 04:07:47PM -0700, Jonathan wrote: > Emmanuel, > > Thanks, that is one of the places I was starting. It turns out that doesn't > quite pick up the necessary stuff from the `Expr` type. I have had better > luck extending the base type `Expr`. It was not hard to get the ari

[sage-devel] Re: Could someone point me to the code that supports math on symbolic equations?

2020-05-22 Thread Jonathan
Emmanuel, Thanks, that is one of the places I was starting. It turns out that doesn't quite pick up the necessary stuff from the `Expr` type. I have had better luck extending the base type `Expr`. It was not hard to get the arithmetic parts (+, -, /,*, pow) working. I'm still looking for/worki

[sage-devel] Weierstrass Preparation Theorem and Hensel Lemma

2020-05-22 Thread mahsa kazemi
Hi, I have two questions regarding the implementation of multivariate power series is Sage. I was wondering if Weierstrass Preparation Theorem and Hensel Lemma ( which use multivariate power series) have been implemented in Sage. About the implementation of multivariate power series (which use

Re: [sage-devel] bug in matrix solve over inexact coefficients?

2020-05-22 Thread AlexGhitza
Hi, On Saturday, May 23, 2020 at 3:55:06 AM UTC+10, Michael Orlitzky wrote: > > > This was changed to "do what MATLAB does" because of the numerical > issues inherent to inexact rings. While > > m = matrix(SR, [0]) > > is singular and the system `m*x == [1]` has no solutions, with > > m =

Re: [sage-devel] bug in matrix solve over inexact coefficients?

2020-05-22 Thread Michael Orlitzky
On 5/22/20 3:39 PM, Nils Bruin wrote: > I think this might need some work: > > S=RealBallField(100) > M=Matrix(S,2,1,[1,1]) > M.solve_right(vector([1,2])) > > There's enough information here to conclude there is no solution; or in > a rather deranged way, perhaps it should give a rather large bal

[sage-devel] Re: Pre-Announcement: Global Virtual SageDays 109 - May 28, 2020 (all timezones)

2020-05-22 Thread Matthias Koeppe
On Thursday, May 21, 2020 at 10:04:26 AM UTC-7, Matthias Koeppe wrote: > > The Global Virtual SageDays 109 will be held on May 28, 2020 (all > timezones). > > I have updated https://wiki.sagemath.org/days109 and created the stream > #sd109 on https://zulip.sagemath.org/ > > > People interested in

Re: [sage-devel] bug in matrix solve over inexact coefficients?

2020-05-22 Thread Nils Bruin
I think this might need some work: S=RealBallField(100) M=Matrix(S,2,1,[1,1]) M.solve_right(vector([1,2])) There's enough information here to conclude there is no solution; or in a rather deranged way, perhaps it should give a rather large ball back so that the multiplication results in a ve

Re: [sage-devel] bug in matrix solve over inexact coefficients?

2020-05-22 Thread Michael Orlitzky
On 5/21/20 8:44 PM, AlexGhitza wrote: > Hi, > > I'm observing the following with version 9.1 (but not with 9.0 where the > behavior is correct): > > sage: m = matrix(SR, [0]) > sage: b = vector([1]) > sage: m.solve_right(b) > (0) > > This should of course raise > > ValueError: matrix equation h

[sage-devel] Re: Could someone point me to the code that supports math on symbolic equations?

2020-05-22 Thread Emmanuel Charpentier
Well, you might consider working on the expressions -. A quick test with Sympy: Python 3.8.3 (default, May 14 2020, 11:03:12) [GCC 9.3.0] on linux Type "help", "copyright", "credits" or "license" for more information. >>> python.el: native completion setup loaded >>> from sympy import * >>> p,V,