Re: [sage-devel] SageMath ArraryFire and afnumpy

2016-09-20 Thread Vincent Delecroix
Hello, Just for information. A lot of efforts are currently done to have Sage more system friendly. In particular - on gentoo distribution https://github.com/cschwan/sage-on-gentoo - on debian distribution (experimental) https://wiki.debian.org/DebianScience/Sage I think this is the only viabl

Re: [sage-devel] Re: AlgebraicReal.minpoly slow

2016-09-20 Thread Clemens Heuberger
Am 2016-09-20 um 20:22 schrieb Marc Mezzarobba: > Clemens Heuberger wrote: >> x = polygen(QQ) >> equation = -96000*x^7 + 41600*x^6 - 6640*x^5 + 560*x^4 >> - 28*x^3 + 8400*x^2 - 140*x + 1 roots = equation.roots(QQbar) >> a_root = roots[-1][0] >> abs_root = abs(a_root) > [...] >>

[sage-devel] Sage's references: new policy?

2016-09-20 Thread John H Palmieri
As discussed in another thread [1]_ on sage-devel recently, I propose changing our policy toward references: - all references should be put into a master bibliography file, and - all references should be, insofar as possible, in a standard form: for a work by a single author "Author" published i

Re: [sage-devel] Re: How to check that something is a real number?

2016-09-20 Thread Bill Page
Admittedly, the author of this thread did not consider the case in the list of targets but when comparing 'conjugate' to 'imag_part' in Sage I was surprized to see the following results. sage: Q. = QuaternionAlgebra(SR, -1, -1) sage: a,b,c=var('a,b,c') sage: q = x + a*i + b*j + c*k sage: conjugate

[sage-devel] SageMath ArraryFire and afnumpy

2016-09-20 Thread Bill Page
SageMath is a computer algebra system http://www.sagemath.org/ that provides an interface to a large number of open source packages for both symbolic and numeric computations, including especially numpy. But unlike a lot of open source SageMath bundles most of these software packages in a sing

Re: [sage-devel] Re: How to check that something is a real number?

2016-09-20 Thread Emmanuel Charpentier
Le mardi 20 septembre 2016 04:02:16 UTC+2, Bill Page a écrit : > > In keeping with Richard's suggestion, in Sage I think a good > _algebraic_ definition of 'real' is > > bool(x/2+conjugate(x)/2 == x) > why not bool(SR(x).imag_part()==0) ? > > So > > sage: def RN(x): > : try:

[sage-devel] Re: Blocker tickets needs review

2016-09-20 Thread 'Martin R' via sage-devel
I think that the emacs ticket should not be a blocker anymore. stakemori did a great job! Martin Am Dienstag, 20. September 2016 15:19:20 UTC+2 schrieb Jeroen Demeyer: > > Hello, this is a reminder that currently 3 blocker tickets need review. > All of them have some discussion, but they all s

Re: [sage-devel] Sage contributor in AAUP journal

2016-09-20 Thread William Stein
On Mon, Sep 19, 2016 at 6:37 PM, kcrisman wrote: > Some of you may be interested in this article by Paul-Olivier Dehaye: > https://www.aaup.org/article/mooc-platforms-surveillance-and-control > It is a rant against the main trends in modern computing. Paul perhaps underestimates the value of "b

Re: [sage-devel] Re: AlgebraicReal.minpoly slow

2016-09-20 Thread Jonas Jermann
Hi Factoring first and then calculating the roots of the factors (you don't need the multiplicity of the factor) worked fast for me... Side note: In AA I noticed that most operations become much faster after a num.simplify(). But in this case the simplify operation is slow itself... Best

Re: [sage-devel] doctesting and sys.path

2016-09-20 Thread Michael Orlitzky
On 09/20/2016 01:40 AM, Clemens Heuberger wrote: > I have two .py files (not in the library) in the same directory where one > imports the other. I'd like to run doctests, but the local directory does not > seem to be in the search path. Is there any way of doing this apart from > fiddling around w

[sage-devel] Blocker tickets needs review

2016-09-20 Thread Jeroen Demeyer
Hello, this is a reminder that currently 3 blocker tickets need review. All of them have some discussion, but they all seem stalled: * Set JUPYTER_CONFIG_DIR https://trac.sagemath.org/ticket/21430 * Old installed version of Cython is used https://trac.sagemath.org/ticket/21441 * CoinBackend: _

Re: [sage-devel] Re: AlgebraicReal.minpoly slow

2016-09-20 Thread Vincent Delecroix
Hi, This thread would better be on ask.sagemath.org or sage-support. Vincent -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-devel+unsubscr...@googlegroups.

[sage-devel] Re: AlgebraicReal.minpoly slow

2016-09-20 Thread Marc Mezzarobba
Clemens Heuberger wrote: > x = polygen(QQ) > equation = -96000*x^7 + 41600*x^6 - 6640*x^5 + 560*x^4 > - 28*x^3 + 8400*x^2 - 140*x + 1 roots = equation.roots(QQbar) > a_root = roots[-1][0] > abs_root = abs(a_root) [...] > Is this expected behaviour? Well, QQbar has a number of w

Re: [sage-devel] Re: How to check that something is a real number?

2016-09-20 Thread Jeroen Demeyer
On 2016-09-19 19:51, Thierry wrote: you should first remove SR(2.3) and 2*pi from your first list That's a problem with the Symbolic Ring. Obviously, it doesn't coerce to a real field (since it contains things like "sin(x)") but these constants *should* coerce. Anyway, I can always special-

Re: [sage-devel] AlgebraicReal.minpoly slow

2016-09-20 Thread Vincent Delecroix
Yes: have a look at composed_op from #18356 -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-devel+unsubscr...@googlegroups.com. To post to this group, send em

[sage-devel] Re: Re: How to check that something is a real number?

2016-09-20 Thread Marc Mezzarobba
Thierry wrote: > My favourite is "in RealField(2)" just in case you fall into something > with smaller precision. But if you prefer to work with coercion, then > you should first remove SR(2.3) and 2*pi from your first list, and you > could do something like: > > is_real = lambda self: > get_coerc

[sage-devel] AlgebraicReal.minpoly slow

2016-09-20 Thread Clemens Heuberger
I am running x = polygen(QQ) equation = -96000*x^7 + 41600*x^6 - 6640*x^5 + 560*x^4 - 28*x^3 + 8400*x^2 - 140*x + 1 roots = equation.roots(QQbar) a_root = roots[-1][0] abs_root = abs(a_root) %time abs_root.minpoly() and am surprised that it takes approx. three minutes. Conse

[sage-devel] Re: broken %attach and %debug ?

2016-09-20 Thread Frédéric Chapoton
I have opened a ticket https://trac.sagemath.org/ticket/21545 to upgrade to ipython 5.1, in the hope that this will solve this *damn* %attach issue Frederic Le vendredi 12 août 2016 14:22:51 UTC+2, Frédéric Chapoton a écrit : > > is this only me, or the recent update to ipython 5.0 has broken s

[sage-devel] doctesting and sys.path

2016-09-20 Thread Clemens Heuberger
I have two .py files (not in the library) in the same directory where one imports the other. I'd like to run doctests, but the local directory does not seem to be in the search path. Is there any way of doing this apart from fiddling around with sys.path in the .py files themselves? Minimal exampl

[sage-devel] SageManifolds 0.9.1 released

2016-09-20 Thread Eric Gourgoulhon
Hi, A new stable version of SageManifolds has just been released. Among the novelties are the computation of the Schouten and Cotton tensors of a pseudo-Riemannian metric and the parallelization (via multiprocessing) of the vector field plots. More details in th