code: def to_xyz100(self, data, description):
          
          rgb_c = compute_to(data, description, self)
          # Step 6: Calculate R, G and B
          # rgb = (rgb_c.T / self.D_RGB).T
          # Step 7: Calculate X, Y and Z
          # xyz = self.solve_M16(rgb)
          return dot(self.invM_, rgb_c)
docs: Input: J or Q; C, M or s; H or h