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

Trying to debug segfault if cudaMemcpy is limited to neighbour list update

Browse Source
This commit is contained in:
Maximilian Gaul
2021-12-25 15:36:08 +01:00
parent 0ea0587442
commit 8009b54113
3 changed files with 79 additions and 21 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
src/main.c
View File

@@ -111,13 +111,20 @@ double reneighbour(
{
double S, E;
printf("10.1\r\n");
S = getTimeStamp();
LIKWID_MARKER_START("reneighbour");
printf("10.2\r\n");
updateAtomsPbc(atom, param);
printf("10.3\r\n");
setupPbc(atom, param);
printf("10.4\r\n");
updatePbc(atom, param);
printf("10.5\r\n");
//sortAtom(atom);
buildNeighbor(atom, neighbor);
printf("10.6\r\n");
LIKWID_MARKER_STOP("reneighbour");
E = getTimeStamp();
@@ -279,6 +286,11 @@ int main(int argc, char** argv)
initialIntegrate(&param, &atom);
const bool doReneighbour = (n + 1) % param.every == 0;
const bool doesReneighbourNextRound = (n + 2) % param.every == 0;
printf("Run %d does reneighbour: %d\r\n", n, doReneighbour);
printf("10\r\n");
if(doReneighbour) {
timer[NEIGH] += reneighbour(&param, &atom, &neighbor);
@@ -286,6 +298,8 @@ int main(int argc, char** argv)
updatePbc(&atom, &param);
}
printf("11\r\n");
if(param.force_field == FF_EAM) {
timer[FORCE] += computeForceEam(&eam, &param, &atom, &neighbor, &stats, 0, n + 1);
} else {