Introduction

Arbitrary precision real ball arithmetic is supplied by Arb which provides a ball representation which tracks error bounds rigorously. Real numbers are represented in mid-rad interval form $[m \pm r] = [m-r, m+r]$.

The Arb real field is constructed using the ArbField constructor. This constructs the parent object for the Arb real field.

However, we define

RealField = ArbField

so that one can construct the Arb real field parent object using RealField instead of ArbField.

The types of real balls in Nemo are given in the following table, along with the libraries that provide them and the associated types of the parent objects.

All the real field types belong to the Field abstract type and the types of elements in this field, i.e. balls in this case, belong to the FieldElem abstract type.

Real field constructors

In order to construct real balls in Nemo, one must first construct the Arb real field itself. This is accomplished with the following constructor.

ArbField(prec::Int)

Return the Arb field with precision in bits prec used for operations on interval midpoints. The precision used for interval radii is a fixed implementation-defined constant (30 bits).

We define

RealField = ArbField

so that one can use RealField in place of ArbField.

Here is an example of creating an Arb real field and using the resulting parent object to coerce values into the resulting field.

RR = RealField(64)

a = RR("0.25")
b = RR("0.1")
c = RR(0.5)
d = RR(12)

Note that whilst one can coerce double precision floating point values into an Arb real field, unless those values can be represented exactly in double precision, the resulting ball can't be any more precise than the double precision supplied.

If instead, values can be represented precisely using decimal arithmetic then one can supply them to Arb using a string. In this case, Arb will store them to the precision specified when creating the Arb field.

If the values can be stored precisely as a binary floating point number, Arb will store the values exactly. See the function isexact below for more information.

Real ball constructors

Once an Arb real field is constructed, there are various ways to construct balls in that field.

Apart from coercing elements into the Arb real field as above, we offer the following functions.

#Base.zero — Method.

zero(R::ArbField)

Return exact zero in the given Arb field.

source

#Base.one — Method.

one(R::ArbField)

Return exact one in the given Arb field.

source

#Nemo.ball — Method.

ball(mid::arb, rad::arb)

Constructs an arb enclosing the range $[m-|r|, m+|r|]$, given the pair $(m, r)$.

source

Here are some examples of constructing balls.

RR = RealField(64)

a = zero(RR)
b = one(RR)
c = ball(RR(3), RR("0.0001"))

Basic functionality

The following basic functionality is provided by the default Arb real field implementation in Nemo, to support construction of generic rings over real fields. Any custom real field implementation in Nemo should provide analogues of these functions along with the usual arithmetic operations.

parent_type(::Type{arb})

Gives the type of the parent object of an Arb real field element.

elem_type(R::ArbField)

Given the parent object for a Arb field, return the type of elements of the field.

mul!(c::arb, a::arb, b::arb)

Multiply $a$ by $b$ and set the existing Arb field element $c$ to the result. This function is provided for performance reasons as it saves allocating a new object for the result and eliminates associated garbage collection.

addeq!(c::arb, a::arb)

In-place addition. Adds $a$ to $c$ and sets $c$ to the result. This function is provided for performance reasons as it saves allocating a new object for the result and eliminates associated garbage collection.

deepcopy(a::arb)

Return a copy of the Arb field element $a$, recursively copying the internal data. Arb field elements are mutable in Nemo so a shallow copy is not sufficient.

Given the parent object R for an Arb real field, the following coercion functions are provided to coerce various elements into the Arb field. Developers provide these by overloading the call operator for the real field parent objects.

R()

Coerce zero into the Arb field.

R(n::Integer)
R(f::fmpz)
R(q::fmpq)

Coerce an integer or rational value into the Arb field.

R(f::Float64)
R(f::BigFloat)

Coerce the given floating point number into the Arb field.

R(f::AbstractString)

Coerce the decimal number, given as a string, into the Arb field.

R(f::arb)

Take an Arb field element that is already in the Arb field and simply return it. A copy of the original is not made.

Here are some examples of coercing elements into the Arb field.

RR = RealField(64)

a = RR(3)
b = RR(QQ(2,3))
c = ball(RR(3), RR("0.0001"))
d = RR("3 +/- 0.0001")
f = RR("-1.24e+12345")
g = RR("nan +/- inf")

In addition to the above, developers of custom real field types must ensure that they provide the equivalent of the function base_ring(R::ArbField) which should return Union{}. In addition to this they should ensure that each real field element contains a field parent specifying the parent object of the real field element, or at least supply the equivalent of the function parent(a::arb) to return the parent object of a real field element.

Conversions

#Base.convert — Method.

convert(::Type{Float64}, x::arb)

Return the midpoint of $x$ rounded down to a machine double.

source

Basic manipulation

Numerous functions are provided to manipulate Arb field elements. Also see the section on basic functionality above.

#Nemo.base_ring — Method.

base_ring(R::ArbField)

Returns Union{} since an Arb field does not depend on any other ring.

source

#Nemo.base_ring — Method.

base_ring(x::arb)

Returns Union{} since an Arb field does not depend on any other ring.

source

#Base.parent — Method.

parent(x::arb)

Return the parent of the given Arb field element.

source

#Nemo.iszero — Method.

iszero(x::arb)

Return true if $x$ is certainly zero, otherwise return false.

source

#Nemo.isnonzero — Method.

isnonzero(x::arb)

Return true if $x$ is certainly not equal to zero, otherwise return false.

source

#Nemo.isone — Method.

isone(x::arb)

Return true if $x$ is certainly not equal to oneo, otherwise return false.

source

#Base.isfinite — Method.

isfinite(x::arb)

Return true if $x$ is finite, i.e. having finite midpoint and radius, otherwise return false.

source

#Nemo.isexact — Method.

isexact(x::arb)

Return true if $x$ is exact, i.e. has zero radius, otherwise return false.

source

#Nemo.isint — Method.

isint(x::arb)

Return true if $x$ is an exact integer, otherwise return false.

source

#Nemo.ispositive — Method.

ispositive(x::arb)

Return true if $x$ is certainly positive, otherwise return false.

source

#Nemo.isnonnegative — Method.

isnonnegative(x::arb)

Return true if $x$ is certainly nonnegative, otherwise return false.

source

#Nemo.isnegative — Method.

isnegative(x::arb)

Return true if $x$ is certainly negative, otherwise return false.

source

#Nemo.isnonpositive — Method.

isnonpositive(x::arb)

Return true if $x$ is certainly nonpositive, otherwise return false.

source

#Nemo.midpoint — Method.

midpoint(x::arb)

Return the midpoint of the ball $x$ as an Arb ball.

source

#Nemo.radius — Method.

radius(x::arb)

Return the radius of the ball $x$ as an Arb ball.

source

#Nemo.accuracy_bits — Method.

accuracy_bits(x::arb)

Return the relative accuracy of $x$ measured in bits, capped between typemax(Int) and -typemax(Int).

source

Here are some examples of basic manipulation of Arb balls.

RR = RealField(64)

a = RR("1.2 +/- 0.001")
b = RR(3)

iszero(a)
isone(b)
ispositive(a)
isfinite(b)
isint(b)
isnegative(a)
c = radius(a)
d = midpoint(b)
f = accuracy_bits(a)
S = parent(a)
T = base_ring(RR)

Arithmetic operations

Nemo provides all the standard field operations for Arb field elements, as follows. Note that division is represented by // since a field is its own fraction field and since exact division is not generally possible in an inexact field.

In addition, the following ad hoc field operations are defined.

Here are some examples of arithmetic operations on Arb balls.

RR = RealField(64)

x = RR(3)
y = RR(QQ(2,3))

a = x + y
b = x*y
d = 1//y - x
d = 3x + ZZ(100)//y
f = (x^2 + y^2) ^ QQ(1,2)

Containment

It is often necessary to determine whether a given exact value or ball is contained in a given real ball or whether two balls overlap. The following functions are provided for this purpose.

#Nemo.overlaps — Method.

overlaps(x::arb, y::arb)

Returns true if any part of the ball $x$ overlaps any part of the ball $y$, otherwise return false.

source

#Base.contains — Method.

contains(x::arb, y::arb)

Returns true if the ball $x$ contains the ball $y$, otherwise return false.

source

#Base.contains — Method.

contains(x::arb, y::Integer)

Returns true if the ball $x$ contains the given integer value, otherwise return false.

source

#Base.contains — Method.

contains(x::arb, y::fmpz)

Returns true if the ball $x$ contains the given integer value, otherwise return false.

source

#Base.contains — Method.

contains(x::arb, y::fmpq)

Returns true if the ball $x$ contains the given rational value, otherwise return false.

source

#Base.contains — Method.

contains(x::arb, y::Rational{Integer})

Returns true if the ball $x$ contains the given rational value, otherwise return false.

source

#Base.contains — Method.

contains(x::arb, y::BigFloat)

Returns true if the ball $x$ contains the given floating point value, otherwise return false.

source

The following functions are also provided for determining if a ball intersects a certain part of the real number line.

#Nemo.contains_zero — Method.

contains_zero(x::arb)

Returns true if the ball $x$ contains zero, otherwise return false.

source

#Nemo.contains_negative — Method.

contains_negative(x::arb)

Returns true if the ball $x$ contains any negative value, otherwise return false.

source

#Nemo.contains_positive — Method.

contains_positive(x::arb)

Returns true if the ball $x$ contains any positive value, otherwise return false.

source

#Nemo.contains_nonnegative — Method.

contains_nonnegative(x::arb)

Returns true if the ball $x$ contains any nonnegative value, otherwise return false.

source

#Nemo.contains_nonpositive — Method.

contains_nonpositive(x::arb)

Returns true if the ball $x$ contains any nonpositive value, otherwise return false.

source

Here are some examples of testing containment.

RR = RealField(64)
x = RR("1 +/- 0.001")
y = RR("3")

overlaps(x, y)
contains(x, y)
contains(y, 3)
contains(x, ZZ(1)//2)
contains_zero(x)
contains_positive(y)

Comparison

Nemo provides a full range of comparison operations for Arb balls. Note that a ball is considered less than another ball if every value in the first ball is less than every value in the second ball, etc.

Firstly, we introduce an exact equality which is distinct from arithmetic equality. This is distinct from arithmetic equality implemented by ==, which merely compares up to the minimum of the precisions of its operands.

#Base.isequal — Method.

isequal(x::arb, y::arb)

Return true if the balls $x$ and $y$ are precisely equal, i.e. have the same midpoints and radii.

source

A full range of functions is available for comparing balls, i.e. ==, !=, <, <=, >=, >. In fact, all these are implemented directly in C. In the table below we document these as though only == and isless had been provided to Julia.

Function

isless(x::arb, y::arb)==(x::arb, y::arb)

As well as these, we provide a full range of ad hoc comparison operators. Again, these are implemented directly in Julia, but we document them as though isless and == were provided.

Function

==(x::arb, y::Integer)==(x::Integer, y::arb)==(x::arb, y::fmpz)==(x::fmpz, y::arb)==(x::arb, y::Float64)==(x::Float64, y::arb)isless(x::arb, y::Integer)isless(x::Integer, y::arb)isless(x::arb, y::fmpz)isless(x::fmpz, y::arb)isless(x::arb, y::Float64)isless(x::Float64, y::arb)isless(x::arb, y::BigFloat)isless(x::BigFloat, y::arb)isless(x::arb, y::fmpq)isless(x::fmpq, y::arb)

Here are some examples of comparison.

RR = RealField(64)
x = RR("1 +/- 0.001")
y = RR("3")
z = RR("4")

isequal(y, z)
x <= z
x == 3
ZZ(3) < z
x != 1.23

Absolute value

#Base.abs — Method.

abs(x::arb)

Return the absolute value of $x$.

source

Here are some examples of taking the absolute value.

RR = RealField(64)
x = RR("-1 +/- 0.001")

a = abs(x)

Inverse

#Base.inv — Method.

inv(x::arb)

Return the multiplicative inverse of $x$, i.e. $1/x$.

source

Here are some examples of taking the inverse.

RR = RealField(64)
x = RR("-3 +/- 0.001")

a = inv(x)

Shifting

#Base.Math.ldexp — Method.

ldexp(x::arb, y::Int)

Return $2^yx$. Note that $y$ can be positive, zero or negative.

source

#Base.Math.ldexp — Method.

ldexp(x::arb, y::fmpz)

Return $2^yx$. Note that $y$ can be positive, zero or negative.

source

Here are some examples of shifting.

RR = RealField(64)
x = RR("-3 +/- 0.001")

a = ldexp(x, 23)
b = ldexp(x, -ZZ(15))

Miscellaneous operations

#Nemo.trim — Method.

trim(x::arb)

Return an arb interval containing $x$ but which may be more economical, by rounding off insignificant bits from the midpoint.

source

#Nemo.unique_integer — Method.

unique_integer(x::arb)

Return a pair where the first value is a boolean and the second is an fmpz integer. The boolean indicates whether the interval $x$ contains a unique integer. If this is the case, the second return value is set to this unique integer.

source

#Nemo.setunion — Method.

setunion(x::arb, y::arb)

Return an arb containing the union of the intervals represented by $x$ and $y$.

source

Here are some examples of miscellaneous operations.

RR = RealField(64)
x = RR("-3 +/- 0.001")
y = RR("2 +/- 0.5")

a = trim(x)
b, c = unique_integer(x)
d = setunion(x, y)

Constants

#Nemo.const_pi — Method.

const_pi(r::ArbField)

Return $\pi = 3.14159\ldots$ as an element of $r$.

source

#Nemo.const_e — Method.

const_e(r::ArbField)

Return $e = 2.71828\ldots$ as an element of $r$.

source

#Nemo.const_log2 — Method.

const_log2(r::ArbField)

Return $\log(2) = 0.69314\ldots$ as an element of $r$.

source

#Nemo.const_log10 — Method.

const_log10(r::ArbField)

Return $\log(10) = 2.302585\ldots$ as an element of $r$.

source

#Nemo.const_euler — Method.

const_euler(r::ArbField)

Return Euler's constant $\gamma = 0.577215\ldots$ as an element of $r$.

source

#Nemo.const_catalan — Method.

const_catalan(r::ArbField)

Return Catalan's constant $C = 0.915965\ldots$ as an element of $r$.

source

#Nemo.const_khinchin — Method.

const_khinchin(r::ArbField)

Return Khinchin's constant $K = 2.685452\ldots$ as an element of $r$.

source

#Nemo.const_glaisher — Method.

const_glaisher(r::ArbField)

Return Glaisher's constant $A = 1.282427\ldots$ as an element of $r$.

source

Here are some examples of computing real constants.

RR = RealField(200)

a = const_pi(RR)
b = const_e(RR)
c = const_euler(RR)
d = const_glaisher(RR)

Mathematical functions

#Base.floor — Method.

floor(x::arb)

Compute the floor of $x$, i.e. the greatest integer not exceeding $x$, as an Arb.

source

#Base.ceil — Method.

ceil(x::arb)

Return the ceiling of $x$, i.e. the least integer not less than $x$, as an Arb.

source

#Base.sqrt — Method.

sqrt(x::arb)

Return the square root of $x$.

source

#Nemo.rsqrt — Method.

rsqrt(x::arb)

Return the reciprocal of the square root of $x$, i.e. $1/\sqrt{x}$.

source

#Nemo.sqrt1pm1 — Method.

sqrt1pm1(x::arb)

Return $\sqrt{1+x}-1$, evaluated accurately for small $x$.

source

#Base.log — Method.

log(x::arb)

Return the principal branch of the logarithm of $x$.

source

#Base.log1p — Method.

log1p(x::arb)

Return $\log(1+x)$, evaluated accurately for small $x$.

source

#Base.exp — Method.

exp(x::arb)

Return the exponential of $x$.

source

#Base.expm1 — Method.

expm1(x::arb)

Return $\exp(x)-1$, evaluated accurately for small $x$.

source

#Base.sin — Method.

sin(x::arb)

Return the sine of $x$.

source

#Base.cos — Method.

cos(x::arb)

Return the cosine of $x$.

source

#Base.Math.sinpi — Method.

sinpi(x::arb)

Return the sine of $\pi x$.

source

#Base.Math.cospi — Method.

cospi(x::arb)

Return the cosine of $\pi x$.

source

#Base.tan — Method.

tan(x::arb)

Return the tangent of $x$.

source

#Base.Math.cot — Method.

cot(x::arb)

Return the cotangent of $x$.

source

#Nemo.tanpi — Method.

tanpi(x::arb)

Return the tangent of $\pi x$.

source

#Nemo.cotpi — Method.

cotpi(x::arb)

Return the cotangent of $\pi x$.

source

#Base.sinh — Method.

sinh(x::arb)

Return the hyperbolic sine of $x$.

source

#Base.cosh — Method.

cosh(x::arb)

Return the hyperbolic cosine of $x$.

source

#Base.tanh — Method.

tanh(x::arb)

Return the hyperbolic tangent of $x$.

source

#Base.Math.coth — Method.

coth(x::arb)

Return the hyperbolic cotangent of $x$.

source

#Base.atan — Method.

atan(x::arb)

Return the arctangent of $x$.

source

#Base.asin — Method.

asin(x::arb)

Return the arcsine of $x$.

source

#Base.acos — Method.

acos(x::arb)

Return the arccosine of $x$.

source

#Base.atanh — Method.

atanh(x::arb)

Return the hyperbolic arctangent of $x$.

source

#Base.asinh — Method.

asinh(x::arb)

Return the hyperbolic arcsine of $x$.

source

#Base.acosh — Method.

acosh(x::arb)

Return the hyperbolic arccosine of $x$.

source

#Base.Math.gamma — Method.

gamma(x::arb)

Return the Gamma function evaluated at $x$.

source

#Base.Math.lgamma — Method.

lgamma(x::arb)

Return the logarithm of the Gamma function evaluated at $x$.

source

#Nemo.rgamma — Method.

rgamma(x::arb)

Return the reciprocal of the Gamma function evaluated at $x$.

source

#Base.Math.digamma — Method.

digamma(x::arb)

Return the logarithmic derivative of the gamma function evaluated at $x$, i.e. $\psi(x)$.

source

#Base.Math.zeta — Method.

zeta(x::arb)

Return the Riemann zeta function evaluated at $x$.

source

#Nemo.sincos — Method.

sincos(x::arb)

Return a tuple $s, c$ consisting of the sine $s$ and cosine $c$ of $x$.

source

#Nemo.sincospi — Method.

sincospi(x::arb)

Return a tuple $s, c$ consisting of the sine $s$ and cosine $c$ of $\pi x$.

source

#Base.Math.sinpi — Method.

sinpi(x::fmpq, r::ArbField)

Return the sine of $\pi x$ in the given Arb field.

source

#Base.Math.cospi — Method.

cospi(x::fmpq, r::ArbField)

Return the cosine of $\pi x$ in the given Arb field.

source

#Nemo.sincospi — Method.

sincospi(x::fmpq, r::ArbField)

Return a tuple $s, c$ consisting of the sine and cosine of $\pi x$ in the given Arb field.

source

#Nemo.sinhcosh — Method.

sinhcosh(x::arb)

Return a tuple $s, c$ consisting of the hyperbolic sine and cosine of $x$.

source

#Base.Math.atan2 — Method.

atan2(x::arb, y::arb)

Return atan2$(b,a) = \arg(a+bi)$.

source

#Nemo.agm — Method.

agm(x::arb, y::arb)

Return the arithmetic-geometric mean of $x$ and $y$

source

#Base.Math.zeta — Method.

zeta(s::arb, a::arb)

Return the Hurwitz zeta function $\zeta(s,a)$.

source

#Base.Math.hypot — Method.

hypot(x::arb, y::arb)

Return $\sqrt{x^2 + y^2}$.

source

#Nemo.root — Method.

root(x::arb, n::Int)

Return the $n$-th root of $x$. We require $x \geq 0$.

source

#Nemo.fac — Method.

fac(x::arb)

Return the factorial of $x$.

source

#Nemo.fac — Method.

fac(n::Int, r::ArbField)

Return the factorial of $n$ in the given Arb field.

source

#Nemo.binom — Method.

binom(x::arb, n::UInt)

Return the binomial coefficient ${x \choose n}$.

source

#Nemo.binom — Method.

binom(n::UInt, k::UInt, r::ArbField)

Return the binomial coefficient ${n \choose k}$ in the given Arb field.

source

#Nemo.fib — Method.

fib(n::fmpz, r::ArbField)

Return the $n$-th Fibonacci number in the given Arb field.

source

#Nemo.fib — Method.

fib(n::Int, r::ArbField)

Return the $n$-th Fibonacci number in the given Arb field.

source

#Base.Math.gamma — Method.

gamma(x::fmpz, r::ArbField)

Return the Gamma function evaluated at $x$ in the given Arb field.

source

#Base.Math.gamma — Method.

gamma(x::fmpq, r::ArbField)

Return the Gamma function evaluated at $x$ in the given Arb field.

source

#Base.Math.zeta — Method.

zeta(n::Int, r::ArbField)

Return the Riemann zeta function $\zeta(n)$ as an element of the given Arb field.

source

#Nemo.bernoulli — Method.

bernoulli(n::Int, r::ArbField)

Return the $n$-th Bernoulli number as an element of the given Arb field.

source

#Nemo.risingfac — Method.

risingfac(x::arb, n::Int)

Return the rising factorial $x(x + 1)\ldots (x + n - 1)$ as an Arb.

source

#Nemo.risingfac — Method.

risingfac(x::fmpq, n::Int, r::ArbField)

Return the rising factorial $x(x + 1)\ldots (x + n - 1)$ as an element of the given Arb field.

source

#Nemo.risingfac2 — Method.

risingfac2(x::arb, n::Int)

Return a tuple containing the rising factorial $x(x + 1)\ldots (x + n - 1)$ and its derivative.

source

#Nemo.polylog — Method.

polylog(s::arb, a::arb)

Return the polylogarithm Li$_s(a)$.

source

#Nemo.polylog — Method.

polylog(s::Int, a::arb)

Return the polylogarithm Li$_s(a)$.

source

#Nemo.chebyshev_t — Method.

chebyshev_t(n::Int, x::arb)

Return the value of the Chebyshev polynomial $T_n(x)$.

source

#Nemo.chebyshev_u — Method.

chebyshev_u(n::Int, x::arb)

Return the value of the Chebyshev polynomial $U_n(x)$.

source

#Nemo.chebyshev_t2 — Method.

chebyshev_t2(n::Int, x::arb)

Return the tuple $(T_{n}(x), T_{n-1}(x))$.

source

#Nemo.chebyshev_u2 — Method.

chebyshev_u2(n::Int, x::arb)

Return the tuple $(U_{n}(x), U_{n-1}(x))$

source

#Nemo.bell — Method.

bell(n::fmpz, r::ArbField)

Return the Bell number $B_n$ as an element of $r$.

source

#Nemo.bell — Method.

bell(n::Int, r::ArbField)

Return the Bell number $B_n$ as an element of $r$.

source

#Nemo.numpart — Method.

numpart(n::fmpz, r::ArbField)

Return the number of partitions $p(n)$ as an element of $r$.

source

#Nemo.numpart — Method.

numpart(n::fmpz, r::ArbField)

Return the number of partitions $p(n)$ as an element of $r$.

source

Here are some examples of real valued mathematical functions.

RR = RealField(64)

a = floor(exp(RR(1)))
b = sinpi(QQ(5,6), RR)
c = gamma(QQ(1,3), RealField(256))
d = bernoulli(1000, RealField(53))
f = polylog(3, RR(-10))