docs: def _initialize_inversion_matrix(self):
          """ the inversion """
          # The bt model is diagonal. The inversion is simply qh = -kappa**2 ph
          self.a = -(self.wv2i+self.kd2)[np.newaxis, np.newaxis, :, :]