RuudH's profile - activity

Could FLINT/ARB-functions be run multi-threaded under Sage? A number of FLINT/ARB functions, for instance the integral,

This is not an answer, but a possible clue that it is a problem of ARB and not the SageMath wrapper. When I run this co

Super! That does the trick :-)

Maybe a very naive question, but when I code something like:

RDD=RealBallField(256)
RDD(3.1)*RDD(2.1)

I obtain:

[6.510000000000000461852778244065128704839826135899304117285652827862296732064351 +/- 9.03e-80]

as the output. Where does all the noise after the 15th digit come from and how could I get rid of it?

Or am I not allowed to increase the accuracy above 53 when using ARB in SageMath?

This is not an answer, but a possible clue that it is a problem of ARB and not the SageMath wrapper. When I run this co

Ghost numbers when using ARB Maybe a very naive question, but when I code something like: RDD=RealBallField(256) RDD(3.

Ghost numbers when using ARB Maybe a very naive question, but when I code something like: RDD=RealBallField(256) RDD(3.

Ghost numbers when using ARB Maybe a very naive question, but when I code something like: RDD=RealBallField(256) RDD(3.

Encountered a strange problem when using latex in the titles, labels or legends of plots:

a=point2d(xy, title="$\lceil x \lfloor$")

works fine (xy contains a list of points), however:

a=point2d(xy, title="$\rceil x \rfloor$")

gives an error (whatever combination I try). The same happens in the legends and axis labels.

Am I doing something wrong here or is this a bug in latex in the plot texts?

Hey, hey, hey... :-) It does work indeed. Sagemath continues to impress me! I am actually aiming to code the Riemann Xi function (s=x+iy):

s/2*(s-1)*pi^(-s/2)*gamma(s/2)*zeta(s)

using complex Ballfields only.

Thanks B r u n o !

EDIT: Here is the 'ball' code for the Riemann Xi function. It is super fast (e.g. compared to pari)

x=1/3;y=200000
z  = CBF(x,y)
z2 = z/RBF(2)
n1 = RBF(1)
K  = z2*(z-n1)*RBF(pi)^(RBF(-1)*z2)*z2.gamma()*z.zeta()
CIF(K)

Thanks Bruno. This is very helpful! I can now use functions like gamma and zeta as a CBF.method, however could not find how to e.g. add or multiply two CBF's together and get a new CBF, i.e. the basic algebra you could do in Arb doesn't seem to be supported (yet), right?

I would like to use some commands/functions from the high performance/precision Arb-library that in turn builds on the Flint-library. Searching the web, I did find references that some commands from both libraries already underpin directly or indirectly the Sagemath environment, however basic commands like for instance:

arb.pi()

don't seem to be supported (or better, I couldn't get them to work yet).

What I am after is to use Flint and Arb functions within Python like in this example, but then under the Sagemath umbrella. Is this already possible?

Thanks.

You're a star, @slelievre! It works perfectly and performs exactly as I wanted. In a dedicated pari-file the loops now function properly. I would have never figured this out myself.

@slelievre. Thanks for the detailed explanation. I have followed the instructions but unfortunately end up with an error (see screenshot attached). It seems that the limit.h file can't be found on my machine. Do you have an idea on how to resolve this? https://ibb.co/g2UZ77

Thanks for your response! I will try that out, but would appreciate some coaching. I have installed Sage/Jupyter notebook on a MAC and am not sure where exactly I have to type the command $ sage -i pari_jupyter. Is this on the command prompt (Terminal)?

I am struggling to get multiple lines of pari/gp code working in a Jupyter Sage notebook.

When I enter:

%%gp
for(x=1,10,print(x))

it all works fine, however when I for instance write:

%%gp
for(x=1,10,{
    print(x);
})

the system just 'hangs' and doesn't return any output. The same issue occurs in the Cocalc Sage cloud environment. I also tried 'pari/gp in your browser' and there it works fine.

Am I doing something wrong or isn't the multi-line pari/gp option supported in Sage?

I am using Sagemath 8.0 on my computer and also have subscribed to Cocalc (the former Sagemathcloud). I currently do quite a bit of programming in Pari/gp and noted that using %gp in the Cocalc-cloud allows me to switch to entering a sequence of Pari/gp commands, however this doesn't seem to work for the PC-version I use. Of course I can circumvent this by transforming each Pari/gp command into a string with pari(' .... '), however that can only be done line by line and produces ugly code.

Is it correct that %gp doesn't work on the PC version and if so, is there another way to write multiple lines of Pari/gp code other than producing line by line strings and/or use the terminal?

Thanks.