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