/home/fmusial/MD_Benchmarks/simulation.cpp: 310 - 325
--------------------------------------------------------------------------------

310:     #pragma omp parallel for
311:     #endif
312:     for (int i = 0; i < NUM_PARTICLES; ++i) {
313:         particles.vx[i] += 0.5 * DT * particles.fx[i];
314:         particles.vy[i] += 0.5 * DT * particles.fy[i];
315:         particles.vz[i] += 0.5 * DT * particles.fz[i];
316:         
317:         // Update positions
318:         particles.x[i] += DT * particles.vx[i];
319:         particles.y[i] += DT * particles.vy[i];
320:         particles.z[i] += DT * particles.vz[i];
321:         
322:         // Apply periodic boundary conditions
323:         particles.x[i] = std::fmod(particles.x[i] + BOX_SIZE, BOX_SIZE);
324:         particles.y[i] = std::fmod(particles.y[i] + BOX_SIZE, BOX_SIZE);
325:         particles.z[i] = std::fmod(particles.z[i] + BOX_SIZE, BOX_SIZE);