code: def get_strain_from_stress(self, stress):
          
          compl_exp = self.get_compliance_expansion()
          strain = 0
          for n, compl in enumerate(compl_exp):
              strain += compl.einsum_sequence([stress]*(n+1)) / factorial(n+1)
          return strain
docs: Gets the strain from a stress state according
          to the compliance expansion corresponding to the
          tensor expansion.