Struct GLFixedTable

pub struct GLFixedTable { /* private fields */ }

The fixed-order Gauss-Legendre integration routines are provided for fast integration of smooth functions with known polynomial order. The n-point Gauss-Legendre rule is exact for polynomials of order 2*n-1 or less. For example, these rules are useful when integrating basis functions to form mass matrices for the Galerkin method. Unlike other numerical integration routines within the library, these routines do not accept absolute or relative error bounds.

Implementations

impl GLFixedTable

pub fn new(n: usize) -> Option<GLFixedTable>

This function determines the Gauss-Legendre abscissae and weights necessary for an n-point fixed order integration scheme. If possible, high precision precomputed coefficients are used. If precomputed weights are not available, lower precision coefficients are computed on the fly.

pub fn glfixed_point( &self, a: f64, b: f64, i: usize, ) -> Result<(f64, f64), Value>

For i in [0, …, t->n - 1], this function obtains the i-th Gauss-Legendre point xi and weight wi on the interval [a,b]. The points and weights are ordered by increasing point value. A function f may be integrated on [a,b] by summing wi * f(xi) over i.

Returns (xi, wi) if it succeeded.

pub fn glfixed<F: Fn(f64) -> f64>(&self, f: F, a: f64, b: f64) -> f64

This function applies the Gauss-Legendre integration rule contained in table self and returns the result.

Trait Implementations

impl Drop for GLFixedTable

fn drop(&mut self)

Executes the destructor for this type.