From: Axel Vogt <mail@axelvogt.de>
To: Brian Gough <bjg@network-theory.co.uk>
Cc: gsl-discuss@sourceware.org
Subject: Re: high precision evaluation of special functions with free software
Date: Sat, 10 Feb 2007 10:06:00 -0000 [thread overview]
Message-ID: <45CD83DF.3040609@axelvogt.de> (raw)
In-Reply-To: <87lkj72vrl.wl%bjg@network-theory.co.uk>
Brian Gough wrote:
> I needed to compute some high-precision (> double precision) test
> cases for the Elliptic integrals and could not find anything to do it
> - I usually use Pari-GP (GPL'ed) but it doesn't seem to have these.
>
> If anyone has any suggestions for how to do this with free software,
> I'd be glad to hear them.
Can't you use "intnum" in Pari? If there is a singularity inbetween
one splits the integration.
For example integrating F:= x -> 1/((1-x^2)^(1/2)*(1-k^2*x^2)^(1/2))
over the unit interval is EllipticK(k) and intnum(x = 0, 1/k, f(x))
+ intnum(x = 1/k, 1, f(x) gives it ( here for 0 < 1/k < 1 ).
next prev parent reply other threads:[~2007-02-10 10:06 UTC|newest]
Thread overview: 3+ messages / expand[flat|nested] mbox.gz Atom feed top
2007-02-09 15:31 Brian Gough
2007-02-10 10:06 ` Axel Vogt [this message]
2007-02-12 9:09 ` Brian Gough
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=45CD83DF.3040609@axelvogt.de \
--to=mail@axelvogt.de \
--cc=bjg@network-theory.co.uk \
--cc=gsl-discuss@sourceware.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).