Merge pull request #276 from brucefan1983/sdc_bug

Fixed a bug in sdc when num_atom*run_steps is too large
brucefan1983 · Sep 20, 2022 · 6df93ed · 6df93ed
2 parents 087c411 + aa9cca4
commit 6df93ed
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/measure/sdc.cu b/src/measure/sdc.cu
@@ -194,7 +194,7 @@ void SDC::postprocess(const char* input_dir)
   CHECK(cudaDeviceSynchronize()); // needed for pre-Pascal GPU
 
   // normalize by the number of atoms and number of time origins
-  const double vac_scaler = 1.0 / (num_atoms_ * num_time_origins_);
+  const double vac_scaler = 1.0 / ((double)num_atoms_ * (double)num_time_origins_);
   for (int nc = 0; nc < num_correlation_steps_; nc++) {
     vacx_[nc] *= vac_scaler;
     vacy_[nc] *= vac_scaler;