fenris 0.0.1

A library for advanced finite element computations in Rust
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

//! Quadrature rules parametrized by polynomial total-order accuracy.
//!
//! TODO: Docs
//!
//! TODO: Tests? Can test that we have equivalence with `fenris-quadrature` maybe

use fenris_quadrature::polyquad;

use crate::nalgebra::RealField;
use crate::quadrature;
use crate::quadrature::{QuadratureError, QuadraturePair2d, QuadraturePair3d};

pub fn triangle<T: RealField>(strength: usize) -> Result<QuadraturePair2d<T>, QuadratureError> {
    let (weights, points) = polyquad::triangle(strength)?;
    Ok(quadrature::convert_quadrature_rule_from_2d_f64((weights, points)))
}

pub fn quadrilateral<T: RealField>(strength: usize) -> Result<QuadraturePair2d<T>, QuadratureError> {
    let (weights, points) = polyquad::quadrilateral(strength)?;
    Ok(quadrature::convert_quadrature_rule_from_2d_f64((weights, points)))
}

pub fn tetrahedron<T: RealField>(strength: usize) -> Result<QuadraturePair3d<T>, QuadratureError> {
    let (weights, points) = polyquad::tetrahedron(strength)?;
    Ok(quadrature::convert_quadrature_rule_from_3d_f64((weights, points)))
}

pub fn hexahedron<T: RealField>(strength: usize) -> Result<QuadraturePair3d<T>, QuadratureError> {
    let (weights, points) = polyquad::hexahedron(strength)?;
    Ok(quadrature::convert_quadrature_rule_from_3d_f64((weights, points)))
}

pub fn prism<T: RealField>(strength: usize) -> Result<QuadraturePair3d<T>, QuadratureError> {
    let (weights, points) = polyquad::prism(strength)?;
    Ok(quadrature::convert_quadrature_rule_from_3d_f64((weights, points)))
}

pub fn pyramid<T: RealField>(strength: usize) -> Result<QuadraturePair3d<T>, QuadratureError> {
    let (weights, points) = polyquad::pyramid(strength)?;
    Ok(quadrature::convert_quadrature_rule_from_3d_f64((weights, points)))
}