accumulate force

brucefan1983 · Jul 27, 2024 · 8b80d41 · 8b80d41
1 parent 8dcc497
commit 8b80d41
Showing 1 changed file with 3 additions and 3 deletions.
diff --git 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;
   }
 }