Struct ExactNumArray 

pub struct ExactNumArray { /* private fields */ }

Batch of ExactNum values at a shared precision p, row-major.

Implementations

impl ExactNumArray

pub fn to_bytes(&self) -> Result<Vec<u8>, Error>

Encode shape, shared p, and every element.

pub fn write_bytes(&self, dest: &mut [u8]) -> Result<usize, Error>

Write the array record into dest without allocating.

pub fn from_bytes(bytes: &[u8]) -> Result<Self, Error>

Decode an array record. Shape product must match the element count.

impl ExactNumArray

pub fn to_csv_string(&self) -> Result<String, Error>

Encode each cell as Display@p=<bits>.

pub fn from_csv_str(text: &str) -> Result<Self, Error>

Parse CSV of Display / Display@p= cells. Empty → NaN.

pub fn to_csv<P: AsRef<Path>>(&self, path: P) -> Result<(), Error>

Write Self::to_csv_string to path.

pub fn from_csv<P: AsRef<Path>>(path: P) -> Result<Self, Error>

Read Self::from_csv_str from path.

impl ExactNumArray

pub fn new(p: usize) -> Self

Empty 0×0 array at precision p.

pub fn filled(p: usize, n: usize, fill: &ExactNum) -> Self

Row vector: n copies of fill (rounded to p).

pub fn filled_2d( p: usize, rows: usize, cols: usize, fill: &ExactNum, ) -> Option<Self>

rows×cols filled with fill (rounded to p).

pub fn from_values(p: usize, vals: &[ExactNum]) -> Self

From values as a row vector; each is rounded to p.

pub fn from_shape( p: usize, rows: usize, cols: usize, vals: &[ExactNum], ) -> Option<Self>

Row-major rows×cols. Length must be rows*cols.

pub fn precision(&self) -> usize

Shared precision.

pub fn shape(&self) -> (usize, usize)

(rows, cols).

pub fn reshape(&self, rows: usize, cols: usize) -> Option<Self>

Reinterpret the same buffer as rows×cols when the product matches.

pub fn len(&self) -> usize

Number of lanes.

pub fn is_empty(&self) -> bool

True if there are no lanes.

pub fn get(&self, i: usize) -> Option<&ExactNum>

Lane i in storage order.

pub fn get2(&self, i: usize, j: usize) -> Option<&ExactNum>

Entry (i, j), or None if out of range.

pub fn as_slice(&self) -> &[ExactNum]

All values in row-major order.

pub fn add(&self, rhs: &Self) -> Option<Self>

Elementwise add at p.

pub fn add_scalar(&self, s: &ExactNum) -> Self

Add a scalar to every lane.

pub fn sub(&self, rhs: &Self) -> Option<Self>

Elementwise sub.

pub fn mul(&self, rhs: &Self) -> Option<Self>

Elementwise mul.

pub fn mul_scalar(&self, s: &ExactNum) -> Self

Multiply every lane by a scalar.

pub fn div(&self, rhs: &Self) -> Option<Self>

Elementwise div.

pub fn sum(&self) -> ExactNum

Sequential sum at p.

pub fn dot(&self, rhs: &Self) -> Option<ExactNum>

Sequential dot product at p.

pub fn matmul(&self, rhs: &Self) -> Option<Self>

Software matmul: (m×k)(k×n) → (m×n). Sequential mul-then-add at p.

pub fn int(&self) -> Self

Elementwise integer part.

pub fn fract(&self) -> Self

Elementwise fractional part.

pub fn ceil(&self) -> Self

Elementwise ceil.

pub fn floor(&self) -> Self

Elementwise floor.

pub fn round(&self, n: usize, rm: RoundingMode) -> Self

Elementwise round with n binary fractional bits.

pub fn abs(&self) -> Self

Elementwise absolute value.

pub fn signum(&self) -> Self

Elementwise signum.

pub fn neg(&self) -> Self

Elementwise negation.

pub fn reciprocal(&self, p: usize, rm: RoundingMode) -> Self

Elementwise reciprocal.

pub fn nth_root(&self, n: usize, p: usize, rm: RoundingMode) -> Self

Elementwise nth_root.

pub fn powi(&self, n: usize, p: usize, rm: RoundingMode) -> Self

Elementwise powi.

pub fn powsi(&self, n: isize, p: usize, rm: RoundingMode) -> Self

Elementwise powsi.

pub fn sin(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::sin. cc is the constants cache, not a global.

pub fn cos(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::cos. cc is the constants cache, not a global.

pub fn tan(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::tan. cc is the constants cache, not a global.

pub fn asin(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::asin. cc is the constants cache, not a global.

pub fn acos(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::acos. cc is the constants cache, not a global.

pub fn atan(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::atan. cc is the constants cache, not a global.

pub fn sinh(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::sinh. cc is the constants cache, not a global.

pub fn cosh(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::cosh. cc is the constants cache, not a global.

pub fn tanh(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::tanh. cc is the constants cache, not a global.

pub fn asinh(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::asinh. cc is the constants cache, not a global.

pub fn acosh(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::acosh. cc is the constants cache, not a global.

pub fn atanh(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::atanh. cc is the constants cache, not a global.

pub fn exp(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::exp. cc is the constants cache, not a global.

pub fn exp2(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::exp2. cc is the constants cache, not a global.

pub fn exp10(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::exp10. cc is the constants cache, not a global.

pub fn expm1(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::expm1. cc is the constants cache, not a global.

pub fn ln(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::ln. cc is the constants cache, not a global.

pub fn log2(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::log2. cc is the constants cache, not a global.

pub fn log10(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::log10. cc is the constants cache, not a global.

pub fn log1p(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::log1p. cc is the constants cache, not a global.

pub fn erf(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::erf. cc is the constants cache, not a global.

pub fn erfc(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::erfc. cc is the constants cache, not a global.

pub fn gamma(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::gamma. cc is the constants cache, not a global.

pub fn ln_gamma(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::ln_gamma. cc is the constants cache, not a global.

pub fn digamma(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::digamma. cc is the constants cache, not a global.

pub fn ei(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::ei. cc is the constants cache, not a global.

pub fn si(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::si. cc is the constants cache, not a global.

pub fn ci(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::ci. cc is the constants cache, not a global.

pub fn li(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::li. cc is the constants cache, not a global.

pub fn fresnel_s(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::fresnel_s. cc is the constants cache, not a global.

pub fn fresnel_c(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::fresnel_c. cc is the constants cache, not a global.

pub fn ai(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::ai. cc is the constants cache, not a global.

pub fn bi(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::bi. cc is the constants cache, not a global.

pub fn elliptic_k(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::elliptic_k. cc is the constants cache, not a global.

pub fn elliptic_e_complete( &self, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise ExactNum::elliptic_e_complete. cc is the constants cache, not a global.

pub fn rem_pi(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Self

Elementwise ExactNum::rem_pi. cc is the constants cache, not a global.

pub fn sin_cos( &self, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> (Self, Self)

Elementwise (sin, cos) with a shared argument reduction.

pub fn sinh_cosh( &self, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> (Self, Self)

Elementwise (sinh, cosh) with a shared evaluation.

pub fn sqrt(&self, p: usize, rm: RoundingMode) -> Self

Elementwise sqrt.

pub fn cbrt(&self, p: usize, rm: RoundingMode) -> Self

Elementwise cbrt.

pub fn legendre_p(&self, n: u32, p: usize, rm: RoundingMode) -> Self

Elementwise P_n(self).

pub fn assoc_legendre_p( &self, n: u32, m: i32, p: usize, rm: RoundingMode, ) -> Self

Elementwise P_n^m(self).

pub fn hypot(&self, other: &ExactNum, p: usize, rm: RoundingMode) -> Self

Elementwise hypot(self, other).

pub fn atan2( &self, x: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise atan2(self, x). cc is the constants cache, not a global.

pub fn pow( &self, n: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise pow(self, n).

pub fn log( &self, base: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise log(self, base).

pub fn gammainc( &self, x: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise γ(self, x).

pub fn gammainc_upper( &self, x: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise Γ(self, x).

pub fn bessel_j( &self, n: usize, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise integer-order J_n(self).

pub fn bessel_j_nu( &self, nu: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise J_ν(self).

pub fn bessel_y( &self, nu: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise Y_ν(self).

pub fn bessel_i( &self, nu: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise I_ν(self).

pub fn bessel_k( &self, nu: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise K_ν(self).

pub fn elliptic_f( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise F(self | m).

pub fn jacobi_sn( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise sn(self | m).

pub fn jacobi_cn( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise cn(self | m).

pub fn jacobi_dn( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise dn(self | m).

pub fn jacobi_am( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise am(self | m).

pub fn jacobi_cd( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise cd(self | m).

pub fn jacobi_ns( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise ns(self | m).

pub fn jacobi_nc( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise nc(self | m).

pub fn jacobi_nd( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise nd(self | m).

pub fn jacobi_sc( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise sc(self | m).

pub fn jacobi_sd( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise sd(self | m).

pub fn jacobi_cs( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise cs(self | m).

pub fn jacobi_ds( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise ds(self | m).

pub fn jacobi_dc( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise dc(self | m).

pub fn elliptic_e( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise incomplete E(self | m).

pub fn elliptic_pi_complete( &self, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise complete Π(n, m) with self = n.

pub fn elliptic_pi( &self, x: &ExactNum, m: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise Π(self; x | m).

pub fn hypergeom_2f1( &self, b: &ExactNum, c: &ExactNum, z: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise {}_2F_1(self, b; c; z).

pub fn betainc( &self, b: &ExactNum, x: &ExactNum, p: usize, rm: RoundingMode, cc: &mut Consts, ) -> Self

Elementwise I_x(self, b).

pub fn lu_decomp( &self, p: usize, rm: RoundingMode, ) -> Option<(Self, Self, Vec<usize>)>

LU with partial pivoting: (L, U, P) such that row i of P·A is original row P[i], and P·A = L·U at (p, rm).

L is unit lower (n×n). U is upper (n×m). A zero pivot (singular) or a failed heap reserve (MemoryAllocation) returns None.

pub fn transpose(&self) -> Self

Row–column transpose.

pub fn qr_decomp(&self, p: usize, rm: RoundingMode) -> Option<(Self, Self)>

Modified Gram–Schmidt QR at (p, rm).

Returns (Q, R) with Q m×k having orthonormal columns, R k×n upper triangular, k = min(m, n). A rank-deficient column is a zero column of Q and a zero diagonal entry of R — not a panic.

pub fn svd_decomp( &self, p: usize, rm: RoundingMode, ) -> Option<(Self, Self, Self)>

Golub–Reinsch SVD at (p, rm).

Returns (U, Σ, V^T) with U m×k (orthonormal columns), Σ k×k diagonal (non-negative, descending), V^T k×n (orthonormal rows), k = min(m, n), so that U · Σ · V^T = A at working precision. Empty input, a NaN/Inf entry, a failed heap reserve, or failure to converge within SVD_ITER_MAX sweeps per singular value returns None.

pub fn eigen_decomp(&self, p: usize, rm: RoundingMode) -> Option<(Self, Self)>

Symmetric QR eigendecomposition at (p, rm).

Returns (Λ, V) where Λ is a 1×n row of eigenvalues (descending) and V is n×n with orthonormal columns, so A V = V diag(Λ) and V diag(Λ) V^T = A at working precision. Non-square, non-symmetric, empty, non-finite, or failure to converge within EIGEN_ITER_MAX sweeps per value returns None.

pub fn fft(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Option<Self>

Radix-2 Cooley–Tukey DFT at (p, rm, cc).

A (1, n) or (n, 1) array is real. A (2, n) array is complex (row 0 real, row 1 imaginary). n must be a power of two and at most FFT_MAX_POINTS. Returns a (2, n) spectrum (unnormalized). Empty, non-finite, or a bad shape returns None.

pub fn ifft(&self, p: usize, rm: RoundingMode, cc: &mut Consts) -> Option<Self>

Inverse radix-2 DFT at (p, rm, cc). Same layout as Self::fft. The result is divided by n (unitary inverse of the unnormalized DFT).

Trait Implementations

impl Clone for ExactNumArray

fn clone(&self) -> ExactNumArray

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for ExactNumArray

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.