From bd29b4adb28520e27009ee4cb0655b5d8f3dd22e Mon Sep 17 00:00:00 2001
From: LeiGeng <lei.geng1997@gmail.com>
Date: Fri, 7 Feb 2025 16:35:38 +0100
Subject: [PATCH] add "round" function in "coord_to_id"

---
 include/xfac/grid.h | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/include/xfac/grid.h b/include/xfac/grid.h
index 7b6e74f..e81c9ae 100644
--- a/include/xfac/grid.h
+++ b/include/xfac/grid.h
@@ -89,7 +89,7 @@ struct Quantics {
         assert(dim==us.size());
         vector<int> id(tensorLen,0);
         for(auto i=0u; i<us.size(); i++) {
-            std::bitset<64> bi=(us[i]-a)/deltaX;
+            std::bitset<64> bi=round((us[i]-a)/deltaX);
             for(auto d=0; d<nBit; d++)
                 if (fused) id[d] |= (bi[d]<<i);
                 else id[i+d*dim]=bi[d];