Skip to content

Commit

Permalink
added complete elliptic integrals of the second kind E
Browse files Browse the repository at this point in the history
  • Loading branch information
bmatschke committed Mar 12, 2021
1 parent ab53d18 commit 27002c2
Show file tree
Hide file tree
Showing 3 changed files with 857 additions and 0 deletions.
@@ -0,0 +1,52 @@
ID: INPUT{id.yaml}

Title: >
Complete elliptic integral of the second kind $E(m)$
Definition: >
The list contains the evaluations of the
complete elliptic integral of the second kind
$E(m) = \int_0^1 \frac{\sqrt{1-mt^2}}{\sqrt{1-t^2}} dt$
for some parameters $m=k^2$,
where $k$ is the ellpitic modulus.
Parameters:
m:
display: $m$
type: R
constraints: $0\leq m \leq 1$

Comments:
comment-k: >
$E(m)$ is also denoted as $E(k)$ despite the relation
$m=k^2$ between the parameter $m$ and the elliptic modulus $k$.
Formulas:
formula-sin: >
$E(m) = \int_0^{\pi/2} \sqrt{1-m\sin^2\theta}} d\theta$.
Programs:

References:

Links:
Wiki:
title: "Wikipedia: Complete elliptic integral of the second kind"
url: https://en.wikipedia.org/wiki/Elliptic_integral#Complete_elliptic_integral_of_the_second_kind

Similar tables:

Keywords:

Tags:
- period
- ellipic curves

Data properties:
type: R
complete: no

Display properties:
number-header: $E(m)$

Numbers: INPUT{numbers.yaml}
@@ -0,0 +1,30 @@
import yaml
import os
from utils.utils import numbers_to_yaml
from utils.utils import real_interval_to_sage_string


path = 'data/Integrals/Periods/Complete_elliptic_integral_of_the_second_kind_E/'

prec10 = 100 #relative precision in base 10

m_range = [a/b for b, a in cartesian_product(([1..50],[0..50]))
if gcd(a,b) == 1 and a/b <= 1]

RIFprec = RealIntervalField(prec10 * 3.4 * 2)

numbers = {}
for m in m_range:
print("m:",m)

number = elliptic_ec(m)

number_str = real_interval_to_sage_string(
RIFprec(number),
max_digits = prec10,
).replace('?','')

numbers[str(m)] = number_str

filename = os.path.join(path, 'numbers.yaml')
yaml.dump(numbers, stream = open(filename, 'w'), sort_keys = False)

0 comments on commit 27002c2

Please sign in to comment.