code: def g_voigt(self):
          
          return (2. * self.voigt[:3, :3].trace() -
                  np.triu(self.voigt[:3, :3]).sum() +
                  3 * self.voigt[3:, 3:].trace()) / 15.
docs: returns the G_v shear modulus