def test_impossible_update_covform(self): """ Test that the _update_covform raises a LinAlgError when the precision matrix is not invertible. """ gaussian = Gaussian(prec=np.zeros([2, 2]), h_vec=[0.0, 0.0], g_val=0.0, var_names=["a", "b"]) with self.assertRaises(Exception): gaussian._update_covform()
def test_log_potential_different_forms(self): """ Test that the log_potential function returns the same value, regardless of the parameter form. """ gaussian = Gaussian(prec=[[7.0, 2.0], [2.0, 6.0]], h_vec=[4.0, 3.0], g_val=0.0, var_names=["a", "b"]) x_val = [7, 3] vrs = ["a", "b"] log_pot_canform = gaussian.log_potential(x_val=x_val, vrs=vrs) gaussian._update_covform() gaussian.canform = False log_pot_covform = gaussian.log_potential(x_val=x_val, vrs=vrs) self.assertAlmostEqual(log_pot_canform, log_pot_covform, places=1)