use crate::{
fem::block::element::{
ElementNodalCoordinates, ElementNodalReferenceCoordinates, ElementNodalVelocities,
FiniteElement, GradientVectors,
quadratic::pyramid::{G, M, N, Pyramid},
solid::SolidFiniteElement,
test::test_finite_element,
},
math::{ScalarList, Tensor},
mechanics::{DeformationGradient, DeformationGradientList, DeformationGradientRateList},
};
fn reference_coordinates() -> ElementNodalReferenceCoordinates<N> {
ElementNodalReferenceCoordinates::from([
[0.0, 0.0, 0.0],
[1.0, 0.0, 0.0],
[1.0, 1.0, 0.0],
[0.0, 1.0, 0.0],
[0.5, 0.5, 0.5],
[0.5, 0.0, 0.0],
[1.0, 0.5, 0.0],
[0.5, 1.0, 0.0],
[0.0, 0.5, 0.0],
[0.25, 0.25, 0.25],
[0.75, 0.25, 0.25],
[0.75, 0.75, 0.25],
[0.25, 0.75, 0.25],
])
}
test_finite_element!(Pyramid);