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

Port pbc to new macros. Enable debugging flags.

Browse Source
This commit is contained in:
Jan Eitzinger 2021-03-23 10:03:55 +01:00
parent 6679b6c8aa
commit fc1fc9fd45
3 changed files with 22 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Makefile
View File

@ -1,4 +1,4 @@
TAG = CLANG
TAG = GCC
#CONFIGURE BUILD SYSTEM
TARGET = MDBench-$(TAG)

8
include_GCC.mk
View File

@ -8,14 +8,12 @@ ANSI_CFLAGS += -std=c99
ANSI_CFLAGS += -pedantic
ANSI_CFLAGS += -Wextra
# CFLAGS = -O0 -g -std=c99 -fargument-noalias
CFLAGS = -O3 -march=znver1 -ffast-math -funroll-loops -fopenmp
CFLAGS = -O0 -g -std=c99 -fargument-noalias
# CFLAGS = -O3 -march=znver1 -ffast-math -funroll-loops -fopenmp
CXXFLAGS = $(CFLAGS)
ASFLAGS = -masm=intel
FCFLAGS =
LFLAGS =
DEFINES = -D_GNU_SOURCE
DEFINES = -D_GNU_SOURCE -DALIGNMENT=64 -DPRECISION=2
INCLUDES =
LIBS =

34
src/pbc.c
View File

@ -68,26 +68,23 @@ void updateAtomsPbc(Atom *atom, Parameter *param)
MD_FLOAT zprd = param->zprd;
for(int i = 0; i < atom->Nlocal; i++) {
MD_FLOAT x = atom_x(i);
MD_FLOAT y = atom_y(i);
MD_FLOAT z = atom_z(i);
if(x < 0.0) {
atom_x(i) = x + xprd;
} else if(x >= xprd) {
atom_x(i) = x - xprd;
if(atom_x(i) < 0.0) {
atom_x(i) += xprd;
} else if(atom_x(i) >= xprd) {
atom_x(i) -= xprd;
}
if(y < 0.0) {
atom_y(i) = y + yprd;
} else if(y >= yprd) {
atom_y(i) = y - yprd;
if(atom_y(i) < 0.0) {
atom_y(i) += yprd;
} else if(atom_y(i) >= yprd) {
atom_y(i) -= yprd;
}
if(z < 0.0) {
atom_z(i) = z + zprd;
} else if(z >= zprd) {
atom_z(i) = z - zprd;
if(atom_z(i) < 0.0) {
atom_z(i) += zprd;
} else if(atom_z(i) >= zprd) {
atom_z(i) -= zprd;
}
}
}
@ -96,7 +93,12 @@ void updateAtomsPbc(Atom *atom, Parameter *param)
* defining ghost atoms around domain
* only creates mapping and coordinate corrections
* that are then enforced in updatePbc */
#define ADDGHOST(dx,dy,dz) Nghost++; BorderMap[Nghost] = i; PBCx[Nghost] = dx; PBCy[Nghost] = dy; PBCz[Nghost] = dz;
#define ADDGHOST(dx,dy,dz) \
Nghost++; \
BorderMap[Nghost] = i; \
PBCx[Nghost] = dx; \
PBCy[Nghost] = dy; \
PBCz[Nghost] = dz;
void setupPbc(Atom *atom, Parameter *param)
{
MD_FLOAT xprd = param->xprd;