From 8b80d41dd5377f121a2fbd1c1628f0ff7e999166 Mon Sep 17 00:00:00 2001
From: brucefan1983 <brucenju@gmail.com>
Date: Sat, 27 Jul 2024 15:10:10 +0800
Subject: [PATCH] accumulate force

---
 src/main_gpumd/add_random_force.cu | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/main_gpumd/add_random_force.cu b/src/main_gpumd/add_random_force.cu
index 983932a3c..8a2b6e323 100644
--- a/src/main_gpumd/add_random_force.cu
+++ b/src/main_gpumd/add_random_force.cu
@@ -43,9 +43,9 @@ static __global__ void add_random_force(
   int n = blockIdx.x * blockDim.x + threadIdx.x;
   if (n < N) {
     curandState state = g_state[n];
-    g_fx[n] = force_variance * curand_normal_double(&state);
-    g_fy[n] = force_variance * curand_normal_double(&state);
-    g_fz[n] = force_variance * curand_normal_double(&state);
+    g_fx[n] += force_variance * curand_normal_double(&state);
+    g_fy[n] += force_variance * curand_normal_double(&state);
+    g_fz[n] += force_variance * curand_normal_double(&state);
     g_state[n] = state;
   }
 }