Struct MaterialEllipticOperator 

pub struct MaterialEllipticOperator<'a, Material>(/* private fields */);

A wrapper that turns any hyper elastic material into an elliptic operator for use with fenris assembly operations.

The wrapper assumes a displacement-based formulation, i.e. that the solution field is the displacement $\vec u (\vec X) = \vec x(\vec X) - \vec X$. In other words, the nodal weights should correspond to displacements, not deformed positions. Alternatively, you may transform deformed positions to displacements as a preprocessing step before handing off the resulting displacements to assembly functionality relying on this operator wrapper.

This implies the following relations:

$$ \begin{aligned} \vec F &= \vec I + (\nabla \vec u)^T, \\ \vec P (\vec F) &= g^T (\nabla \vec u), \\ \mathcal{C}_{\vec P}(\vec F, \vec a, \vec b) &= \mathcal{C}_g(\nabla \vec u, \vec a, \vec b). \\ \end{aligned} $$

Implementations

impl<'a, Material> MaterialEllipticOperator<'a, Material>

pub fn new(material: &'a Material) -> Self

Trait Implementations

impl<'a, T, GeometryDim, Material> EllipticContraction<T, GeometryDim> for MaterialEllipticOperator<'a, Material>

where T: Real, GeometryDim: SmallDim, Material: HyperelasticMaterial<T, GeometryDim>, DefaultAllocator: DimAllocator<T, GeometryDim>,

fn contract( &self, u_grad: &OMatrix<T, GeometryDim, GeometryDim>, a: &OVector<T, GeometryDim>, b: &OVector<T, GeometryDim>, parameters: &Self::Parameters, ) -> OMatrix<T, Self::SolutionDim, Self::SolutionDim>

Compute $ \mathcal{C}_g(\nabla u, a, b)$ with the given parameters.

fn symmetry(&self) -> Symmetry

Whether the contraction operator is symmetric.

fn accumulate_contractions_into( &self, output: DMatrixViewMut<'_, T>, alpha: T, u_grad: &OMatrix<T, GeometryDim, Self::SolutionDim>, a: DVectorView<'_, T>, b: DVectorView<'_, T>, parameters: &Self::Parameters, )

Compute the contraction for a number of vectors at the same time, with the given parameters.

impl<'a, T, GeometryDim, Material> EllipticEnergy<T, GeometryDim> for MaterialEllipticOperator<'a, Material>

where T: Real, GeometryDim: SmallDim, Material: HyperelasticMaterial<T, GeometryDim>, DefaultAllocator: DimAllocator<T, GeometryDim>,

fn compute_energy( &self, u_grad: &OMatrix<T, GeometryDim, GeometryDim>, parameters: &Self::Parameters, ) -> T

impl<'a, T, GeometryDim, Material> EllipticOperator<T, GeometryDim> for MaterialEllipticOperator<'a, Material>

where T: Real, GeometryDim: SmallDim, Material: HyperelasticMaterial<T, GeometryDim>, DefaultAllocator: DimAllocator<T, GeometryDim>,

fn compute_elliptic_operator( &self, u_grad: &OMatrix<T, GeometryDim, GeometryDim>, parameters: &Self::Parameters, ) -> OMatrix<T, GeometryDim, Self::SolutionDim>

Compute the elliptic operator $g = g(\nabla u)$ with the provided operator parameters.

fn compute_elliptic_operator_transpose( &self, u_grad: &OMatrix<T, GeometryDim, Self::SolutionDim>, parameters: &Self::Parameters, ) -> OMatrix<T, Self::SolutionDim, GeometryDim>

Compute the transpose $g^T$ of the elliptic operator $g$.

impl<'a, T, GeometryDim, Material> Operator<T, GeometryDim> for MaterialEllipticOperator<'a, Material>

where T: Real, GeometryDim: SmallDim, Material: HyperelasticMaterial<T, GeometryDim>, DefaultAllocator: DimAllocator<T, GeometryDim>,

type SolutionDim = GeometryDim

type Parameters = <Material as HyperelasticMaterial<T, GeometryDim>>::Parameters

The parameters associated with the operator.