moebiusband/MD-Bench
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

♻️ further refactoring step fixing

Browse Source
This commit is contained in:
Martin Bauernfeind 2022-06-23 19:34:16 +02:00
parent 7f068a6959
commit f4313f64e5
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/force.cu
View File

@ -91,14 +91,14 @@ __global__ void calc_force(
atom_fz(i) = fiz; atom_fz(i) = fiz;
} }
__global__ void kernel_initial_integrate(MD_FLOAT dtforce, MD_FLOAT dt, int Nlocal, Atom a) { __global__ void kernel_initial_integrate(MD_FLOAT dtforce, MD_FLOAT dt, int Nlocal, Atom *a) {
const int i = blockIdx.x * blockDim.x + threadIdx.x; const int i = blockIdx.x * blockDim.x + threadIdx.x;
if( i >= Nlocal ) { if( i >= Nlocal ) {
return; return;
} }
Atom *atom = &a; Atom *atom = a;
atom_vx(i) += dtforce * atom_fx(i); atom_vx(i) += dtforce * atom_fx(i);
atom_vy(i) += dtforce * atom_fy(i); atom_vy(i) += dtforce * atom_fy(i);
@ -108,14 +108,14 @@ __global__ void kernel_initial_integrate(MD_FLOAT dtforce, MD_FLOAT dt, int Nloc
atom_z(i) = atom_z(i) + dt * atom_vz(i); atom_z(i) = atom_z(i) + dt * atom_vz(i);
} }
__global__ void kernel_final_integrate(MD_FLOAT dtforce, int Nlocal, Atom a) { __global__ void kernel_final_integrate(MD_FLOAT dtforce, int Nlocal, Atom *a) {
const int i = blockIdx.x * blockDim.x + threadIdx.x; const int i = blockIdx.x * blockDim.x + threadIdx.x;
if( i >= Nlocal ) { if( i >= Nlocal ) {
return; return;
} }
Atom *atom = &a; Atom *atom = a;
atom_vx(i) += dtforce * atom_fx(i); atom_vx(i) += dtforce * atom_fx(i);
atom_vy(i) += dtforce * atom_fy(i); atom_vy(i) += dtforce * atom_fy(i);