Weird state:
Final resulting stiffness tensor matches the one from numpy, but intermediate results have one transposed dimension (ijkl vs ijlk)