fix b_data bug

GeomScale · Aug 16, 2024 · b281839 · b281839
1 parent 5e49904
commit b281839
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/include/random_walks/uniform_accelerated_billiard_walk.hpp b/include/random_walks/uniform_accelerated_billiard_walk.hpp
@@ -235,9 +235,11 @@ struct AcceleratedBilliardWalk
             NT T;
             const NT dl = 0.995;
             int it;
+            typename Point::Coeff b;
+            NT* b_data;
             if constexpr (std::is_same<MT, Eigen::SparseMatrix<NT, Eigen::RowMajor>>::value) {
-                typename Point::Coeff b = P.get_vec();
-                NT* b_data = b.data();
+                b = P.get_vec();
+                b_data = b.data();
             }
 
             for (auto j=0u; j<walk_length; ++j)