fix numeric calculation

YunDeng98 · YunDeng98 · commit 5b136981ebdf · 2024-09-10T16:24:02.000-07:00
diff --git a/SINGER/SINGER/ARG.cpp b/SINGER/SINGER/ARG.cpp
@@ -54,10 +54,10 @@ int ARG::get_index(double x) {
 void ARG::compute_rhos_thetas(double r, double m) {
     int n = (int) coordinates.size() - 1;
     for (int i = 0; i < n; i++) {
-        // rhos.push_back(r*(coordinates[i+1] - coordinates[i]));
-        // thetas.push_back(m*(coordinates[i+1] - coordinates[i]));
-        rhos.push_back(r*bin_size);
-        thetas.push_back(m*bin_size);
+        rhos.push_back(r*(coordinates[i+1] - coordinates[i]));
+        thetas.push_back(m*(coordinates[i+1] - coordinates[i]));
+        // rhos.push_back(r*bin_size);
+        // thetas.push_back(m*bin_size);
     }
 }
 

Original file line number	Diff line number	Diff line change
`@@ -54,10 +54,10 @@ int ARG::get_index(double x) {`
`54`	`54`	`void ARG::compute_rhos_thetas(double r, double m) {`
`55`	`55`	`int n = (int) coordinates.size() - 1;`
`56`	`56`	`for (int i = 0; i < n; i++) {`
`57`		`- // rhos.push_back(r*(coordinates[i+1] - coordinates[i]));`
`58`		`- // thetas.push_back(m*(coordinates[i+1] - coordinates[i]));`
`59`		`- rhos.push_back(r*bin_size);`
`60`		`- thetas.push_back(m*bin_size);`
	`57`	`+ rhos.push_back(r*(coordinates[i+1] - coordinates[i]));`
	`58`	`+ thetas.push_back(m*(coordinates[i+1] - coordinates[i]));`
	`59`	`+ // rhos.push_back(r*bin_size);`
	`60`	`+ // thetas.push_back(m*bin_size);`
`61`	`61`	`}`
`62`	`62`	`}`
`63`	`63`