Invariants function can now accept all Voigt components as separate arrays. Added test for the return shape.