Gaussian HMC Reflection Assertions - is_in fix

GeomScale · Jun 5, 2024 · 11d9662 · 11d9662
1 parent 8afe4d0
commit 11d9662
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/include/random_walks/gaussian_hamiltonian_monte_carlo_exact_walk.hpp b/include/random_walks/gaussian_hamiltonian_monte_carlo_exact_walk.hpp
@@ -100,15 +100,15 @@ struct Walk
                 auto pbpair = P.trigonometric_positive_intersect(_p, _v, _omega, _facet_prev);
                 if (T <= pbpair.first) {
                     update_position(_p, _v, T, _omega);
-                    assert(P.is_in(_p) != 0 && "Point is outside the polytope after update_position");
+                    assert(P.is_in(_p) == 0);
                     break;
                 }
                 _lambda_prev = pbpair.first;
                 T -= _lambda_prev;
                 update_position(_p, _v, _lambda_prev, _omega);
-                assert(P.is_in(_p) != 0 && "Point is outside the polytope after update_position");
+                assert(P.is_in(_p) == 0);
                 P.compute_reflection(_v, _p, pbpair.second);
-                assert(P.is_in(_p) != 0 && "Point is outside the polytope after compute_reflection");
+                assert(P.is_in(_p) == 0);
                 it++;
             }
             if (it == _rho) {