[sage-support] Re: QQbar operations extremely slow ; is this normal ?

2021-06-11 Thread Fredrik Johansson
What is reasonable depends on what you expect to be able to do with the solutions. Numerical evaluation should be easy, but if you want canonical forms (minimal polynomials) or the ability to check equalities, that's going to be far more costly. It looks like the eigenvalues of this matrix will

[sage-support] Re: SEGV in mpmath.findroot

2012-10-02 Thread Fredrik Johansson
On Monday, October 1, 2012 4:19:47 PM UTC+2, Georgi Guninski wrote: > > import mpmath > mpmath.mp.pretty=True > mpmath.mp.dps=30 > def F(x): > return mpmath.zeta(x)+mpmath.zeta(x,derivative=1) > > r=mpmath.findroot(F,[0.1+mpmath.j],solver="muller") > > > Unhandled SIGSEGV: A seg

[sage-support] Re: get values of log integral without branch cut

2011-05-12 Thread Fredrik Johansson
On May 13, 2:50 am, kcrisman wrote: > On May 12, 8:39 pm, Fredrik Johansson > wrote: > > > On May 12, 11:19 pm, kcrisman wrote: > > > > This should, in theory, give a plot of li(20^z) along the critical > > > line of the Riemann zeta function.  Unfortunatel

[sage-support] Re: get values of log integral without branch cut

2011-05-12 Thread Fredrik Johansson
On May 12, 11:19 pm, kcrisman wrote: > This should, in theory, give a plot of li(20^z) along the critical > line of the Riemann zeta function.  Unfortunately, as you will see if > you plot this, it succeeds until it hits a branch cut (I assume), and > does not look so nice, not to mention missing

[sage-support] Re: Using scipy special functions

2011-04-13 Thread Fredrik Johansson
On Apr 13, 7:48 pm, ObsessiveMathsFreak wrote: > On Apr 12, 8:52 pm, ObsessiveMathsFreak > > wrote: > > On Apr 12, 3:51 am, kcrisman wrote: > > > > > > Or simply legendre_P, legendre_Q > > > > > Unfortunately, these functions do not support non integer values of n, > > > > i.e. they don't suppor

[sage-support] Re: mpmath's incomplete gamma function evaluates much faster for small arguments

2010-12-01 Thread Fredrik Johansson
On Nov 27, 2:15 pm, KvS wrote: > Dear all, > > just a quick question/remark, today I was working (plotting etc.) with > some quantities that involve (upper) incomplete gamma functions. Some > operations took very long, as it turned out due to the incomplete > gamma function evaluating very slowly

[sage-support] Re: Bessel Function Zeroes

2010-07-10 Thread Fredrik Johansson
On Jul 2, 1:02 am, Pat LeSmithe wrote: > Possibly naive questions: > >   * Are there analogous arbitrary-precision routines in or wrapped by Sage? Arbitrary-precision computation of Bessel function zeros will be provided in the next version of mpmath. See http://fredrik-j.blogspot.com/2010/07/sa