Fix trigonal system

Looks like we have two problems here. First, we need to
make sure we make use of the symmetry relations for
trigonal system to calculate C24 etc. We also need to
use different strains for high-trigonal compared to
hexagonal.
Fix bug for trigonal - cMatrix matrices were wrong

Also avoid double calculating stuff.