Additional comment
Bug fix: Had to revert to full Voigt matrices for stiffness and compliance.
Stiffness is really compliance. Name change in elasticity module need to be reflected here.
Return, no raise the result.