13     det(J, ::CeedDim{dim})
17 @inline det(J, ::CeedDim{1}) = @inbounds J[1]  argument
18 @inline det(J, ::CeedDim{2}) = @inbounds J[1]*J[4] - J[3]*J[2]  argument
20 @inline det(J, ::CeedDim{3}) = @inbounds (  argument
21     J[1]*(J[5]*J[9] - J[6]*J[8]) -
22     J[2]*(J[4]*J[9] - J[6]*J[7]) +
23     J[3]*(J[4]*J[8] - J[5]*J[7])
28     setvoigt(J::StaticArray{Tuple{D,D},T,2})
29     setvoigt(J, ::CeedDim{dim})
31 Given a symmetric matrix `J`, return a `SVector` that encodes `J` using the [Voigt
34 The size of the symmetric matrix `J` must be known statically, either specified using
37 @inline setvoigt(J::StaticArray{Tuple{D,D}}) where {D} = setvoigt(J, CeedDim(D))
38 @inline setvoigt(J, ::CeedDim{1}) = @inbounds @SVector [J[1]]  argument
39 @inline setvoigt(J, ::CeedDim{2}) = @inbounds @SVector [J[1], J[4], J[2]]  argument
40 @inline setvoigt(J, ::CeedDim{3}) = @inbounds @SVector [J[1], J[5], J[9], J[6], J[3], J[2]]  argument
42 @inline function setvoigt!(V, J, ::CeedDim{1})  argument
43     @inbounds V[1] = J[1]
46 @inline function setvoigt!(V, J, ::CeedDim{2})  argument
48         V[1] = J[1]
49         V[2] = J[4]
50         V[3] = J[2]
54 @inline function setvoigt!(V, J, ::CeedDim{3})  argument
56         V[1] = J[1]
57         V[2] = J[5]
58         V[3] = J[9]
59         V[4] = J[6]
60         V[5] = J[3]
61         V[6] = J[2]
83 @inline function getvoigt!(J, V, ::CeedDim{1})  argument
84     @inbounds J[1, 1] = V[1]
87 @inline function getvoigt!(J, V, ::CeedDim{2})  argument
90         J[1,1] = V[1] ; J[1,2] = V[3]
91         J[2,1] = V[3] ; J[2,2] = V[2]
96 @inline function getvoigt!(J, V, ::CeedDim{3})  argument
99         J[1,1] = V[1] ; J[1,2] = V[6] ; J[1,3] = V[5]
100         J[2,1] = V[6] ; J[2,2] = V[2] ; J[2,3] = V[4]
101         J[3,1] = V[5] ; J[3,2] = V[4] ; J[3,3] = V[3]