git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@3685 f3b2605a-c512-4ea7-a41b-209d697bcdaa
This commit is contained in:
@ -158,9 +158,9 @@ void PairDSMC::compute(int eflag, int vflag)
|
|||||||
convert_double_to_equivalent_int(num_of_collisions_double);
|
convert_double_to_equivalent_int(num_of_collisions_double);
|
||||||
|
|
||||||
if (num_of_collisions > number_of_A)
|
if (num_of_collisions > number_of_A)
|
||||||
error->warning("num_of_collisions > number_of_A");
|
error->warning("Pair dsmc: num_of_collisions > number_of_A");
|
||||||
if (num_of_collisions > number_of_B)
|
if (num_of_collisions > number_of_B)
|
||||||
error->warning("num_of_collisions > number_of_B");
|
error->warning("Pair dsmc: num_of_collisions > number_of_B");
|
||||||
|
|
||||||
// perform collisions on pairs of particles in icell
|
// perform collisions on pairs of particles in icell
|
||||||
|
|
||||||
|
|||||||
@ -114,7 +114,7 @@ void PairGayBerneGPU::compute(int eflag, int vflag)
|
|||||||
if (eflag || vflag) ev_setup(eflag,vflag);
|
if (eflag || vflag) ev_setup(eflag,vflag);
|
||||||
else evflag = vflag_fdotr = 0;
|
else evflag = vflag_fdotr = 0;
|
||||||
if (vflag_atom)
|
if (vflag_atom)
|
||||||
error->all("Per-atom virial not available with GPU Gay-Berne.");
|
error->all("Per-atom virial not available with GPU Gay-Berne");
|
||||||
|
|
||||||
int nlocal = atom->nlocal;
|
int nlocal = atom->nlocal;
|
||||||
int nall = nlocal + atom->nghost;
|
int nall = nlocal + atom->nghost;
|
||||||
@ -151,7 +151,7 @@ void PairGayBerneGPU::compute(int eflag, int vflag)
|
|||||||
success);
|
success);
|
||||||
#endif
|
#endif
|
||||||
if (!success)
|
if (!success)
|
||||||
error->one("Total # of atoms exceeds maximum allowed per GPGPU.\n");
|
error->one("Total # of atoms exceeds maximum allowed per GPGPU");
|
||||||
|
|
||||||
// copy atom data to GPU
|
// copy atom data to GPU
|
||||||
gb_gpu_atom(atom->x,atom->quat,atom->type,rebuild,my_thread);
|
gb_gpu_atom(atom->x,atom->quat,atom->type,rebuild,my_thread);
|
||||||
@ -324,7 +324,7 @@ void PairGayBerneGPU::init_style()
|
|||||||
lj1, lj2, lj3, lj4, offset, force->special_lj,
|
lj1, lj2, lj3, lj4, offset, force->special_lj,
|
||||||
neighbor->oneatom, my_thread, my_gpu);
|
neighbor->oneatom, my_thread, my_gpu);
|
||||||
if (!init_ok)
|
if (!init_ok)
|
||||||
error->one("AT LEAST ONE PROCESS COULD NOT ALLOCATE A CUDA-ENABLED GPU.");
|
error->one("At least one process could not allocate a CUDA-enabled gpu");
|
||||||
|
|
||||||
if (ij_new[my_thread]!=NULL)
|
if (ij_new[my_thread]!=NULL)
|
||||||
delete [] ij_new[my_thread];
|
delete [] ij_new[my_thread];
|
||||||
@ -335,7 +335,7 @@ void PairGayBerneGPU::init_style()
|
|||||||
neighbor->requests[irequest]->half = 0;
|
neighbor->requests[irequest]->half = 0;
|
||||||
neighbor->requests[irequest]->full = 1;
|
neighbor->requests[irequest]->full = 1;
|
||||||
if (force->newton_pair)
|
if (force->newton_pair)
|
||||||
error->all("Cannot use newton with GPU GayBerne pair style.");
|
error->all("Cannot use newton pair with GPU GayBerne pair style");
|
||||||
|
|
||||||
if (comm->me == 0 && screen) {
|
if (comm->me == 0 && screen) {
|
||||||
printf("\n-------------------------------------");
|
printf("\n-------------------------------------");
|
||||||
|
|||||||
@ -95,7 +95,7 @@ void PairLJCutGPU::compute(int eflag, int vflag)
|
|||||||
if (eflag || vflag) ev_setup(eflag,vflag);
|
if (eflag || vflag) ev_setup(eflag,vflag);
|
||||||
else evflag = vflag_fdotr = 0;
|
else evflag = vflag_fdotr = 0;
|
||||||
if (vflag_atom)
|
if (vflag_atom)
|
||||||
error->all("Per-atom virial not available with GPU Gay-Berne.");
|
error->all("Per-atom virial not available with GPU lj/cut");
|
||||||
|
|
||||||
int nlocal = atom->nlocal;
|
int nlocal = atom->nlocal;
|
||||||
int nall = nlocal + atom->nghost;
|
int nall = nlocal + atom->nghost;
|
||||||
@ -111,7 +111,7 @@ void PairLJCutGPU::compute(int eflag, int vflag)
|
|||||||
// copy nbors to GPU
|
// copy nbors to GPU
|
||||||
if (rebuild)
|
if (rebuild)
|
||||||
if (!lj_gpu_reset_nbors(nall, inum, ilist, list->numneigh))
|
if (!lj_gpu_reset_nbors(nall, inum, ilist, list->numneigh))
|
||||||
error->one("Total # of atoms exceed maximum allowed per GPGPU.\n");
|
error->one("Total # of atoms exceeds maximum allowed per GPGPU");
|
||||||
|
|
||||||
// copy atom data to GPU
|
// copy atom data to GPU
|
||||||
lj_gpu_atom(atom->x,atom->type,rebuild);
|
lj_gpu_atom(atom->x,atom->type,rebuild);
|
||||||
@ -225,7 +225,7 @@ void PairLJCutGPU::init_style()
|
|||||||
|
|
||||||
if (!lj_gpu_init(ij_size, atom->ntypes+1, cutsq, sigma, epsilon, lj1, lj2,lj3,
|
if (!lj_gpu_init(ij_size, atom->ntypes+1, cutsq, sigma, epsilon, lj1, lj2,lj3,
|
||||||
lj4, offset, force->special_lj, neighbor->oneatom, my_gpu))
|
lj4, offset, force->special_lj, neighbor->oneatom, my_gpu))
|
||||||
error->one("AT LEAST ONE PROCESS COULD NOT ALLOCATE A CUDA-ENABLED GPU.");
|
error->one("At least one process could not allocate a CUDA-enabled gpu");
|
||||||
|
|
||||||
if (ij_new!=NULL)
|
if (ij_new!=NULL)
|
||||||
delete [] ij_new;
|
delete [] ij_new;
|
||||||
@ -235,7 +235,7 @@ void PairLJCutGPU::init_style()
|
|||||||
neighbor->requests[irequest]->half = 0;
|
neighbor->requests[irequest]->half = 0;
|
||||||
neighbor->requests[irequest]->full = 1;
|
neighbor->requests[irequest]->full = 1;
|
||||||
if (force->newton_pair)
|
if (force->newton_pair)
|
||||||
error->all("Cannot use newton with GPU LJCut pair style.");
|
error->all("Cannot use newton pair with GPU lj/cut pair style");
|
||||||
|
|
||||||
if (comm->me == 0 && screen) {
|
if (comm->me == 0 && screen) {
|
||||||
printf("\n-------------------------------------");
|
printf("\n-------------------------------------");
|
||||||
|
|||||||
@ -2,8 +2,8 @@
|
|||||||
# this file is auto-edited when those packages are included/excluded
|
# this file is auto-edited when those packages are included/excluded
|
||||||
|
|
||||||
PKG_INC = -I../../lib/reax -I../../lib/poems -I../../lib/meam
|
PKG_INC = -I../../lib/reax -I../../lib/poems -I../../lib/meam
|
||||||
PKG_PATH = -L../../lib/reax -L../../lib/poems -L../../lib/meam
|
PKG_PATH = -L../../lib/reax -L../../lib/poems -L../../lib/meam -L../../lib/gpu
|
||||||
PKG_LIB = -lreax -lpoems -lmeam
|
PKG_LIB = -lreax -lpoems -lmeam -lgpu
|
||||||
|
|
||||||
PKG_SYSPATH = $(reax_SYSPATH) $(meam_SYSPATH)
|
PKG_SYSPATH = $(reax_SYSPATH) $(meam_SYSPATH) $(gpu_SYSPATH)
|
||||||
PKG_SYSLIB = $(reax_SYSLIB) $(meam_SYSLIB)
|
PKG_SYSLIB = $(reax_SYSLIB) $(meam_SYSLIB) $(gpu_SYSLIB)
|
||||||
|
|||||||
@ -44,7 +44,7 @@ ComputeEventDisplace::ComputeEventDisplace(LAMMPS *lmp, int narg, char **arg) :
|
|||||||
|
|
||||||
double displace_dist = atof(arg[3]);
|
double displace_dist = atof(arg[3]);
|
||||||
if (displace_dist <= 0.0)
|
if (displace_dist <= 0.0)
|
||||||
error->all("Distnace must be > 0 for compute event/displace");
|
error->all("Distance must be > 0 for compute event/displace");
|
||||||
displace_distsq = displace_dist * displace_dist;
|
displace_distsq = displace_dist * displace_dist;
|
||||||
|
|
||||||
// fix event ID will be set later by PRD
|
// fix event ID will be set later by PRD
|
||||||
|
|||||||
@ -88,7 +88,7 @@ void PRD::command(int narg, char **arg)
|
|||||||
|
|
||||||
// total # of timesteps must be multiple of t_event
|
// total # of timesteps must be multiple of t_event
|
||||||
|
|
||||||
if (t_event == 0) error->universe_all("Invalid t_event in prd command");
|
if (t_event <= 0) error->universe_all("Invalid t_event in prd command");
|
||||||
if (nsteps % t_event)
|
if (nsteps % t_event)
|
||||||
error->universe_all("PRD nsteps must be multiple of t_event");
|
error->universe_all("PRD nsteps must be multiple of t_event");
|
||||||
if (t_corr % t_event)
|
if (t_corr % t_event)
|
||||||
|
|||||||
@ -149,14 +149,11 @@ void DumpCFG::init()
|
|||||||
for (int i = 0; i < size_one; i++) {
|
for (int i = 0; i < size_one; i++) {
|
||||||
if (i == 0) ptr = strtok(format," \0");
|
if (i == 0) ptr = strtok(format," \0");
|
||||||
else ptr = strtok(NULL," \0");
|
else ptr = strtok(NULL," \0");
|
||||||
if (ptr == NULL) error->all("Dump cfg user format string error");
|
|
||||||
delete [] vformat[i];
|
delete [] vformat[i];
|
||||||
vformat[i] = new char[strlen(ptr) + 2];
|
vformat[i] = new char[strlen(ptr) + 2];
|
||||||
strcpy(vformat[i],ptr);
|
strcpy(vformat[i],ptr);
|
||||||
vformat[i] = strcat(vformat[i]," ");
|
vformat[i] = strcat(vformat[i]," ");
|
||||||
}
|
}
|
||||||
if (strtok(NULL," \0"))
|
|
||||||
error->all("Dump cfg user format string error");
|
|
||||||
|
|
||||||
// find current ptr for each compute,fix,variable
|
// find current ptr for each compute,fix,variable
|
||||||
// check that fix frequency is acceptable
|
// check that fix frequency is acceptable
|
||||||
|
|||||||
@ -202,7 +202,7 @@ FixAveHisto::FixAveHisto(LAMMPS *lmp, int narg, char **arg) :
|
|||||||
}
|
}
|
||||||
if (i == 0) kind = kindflag;
|
if (i == 0) kind = kindflag;
|
||||||
else if (kindflag != kind)
|
else if (kindflag != kind)
|
||||||
error->all("Fix ave/histo inputs are not all global, peraton, or local");
|
error->all("Fix ave/histo inputs are not all global, peratom, or local");
|
||||||
}
|
}
|
||||||
|
|
||||||
if (kind == PERATOM && mode == SCALAR)
|
if (kind == PERATOM && mode == SCALAR)
|
||||||
|
|||||||
@ -337,7 +337,7 @@ void FixBoxRelax::min_clearstore()
|
|||||||
void FixBoxRelax::min_pushstore()
|
void FixBoxRelax::min_pushstore()
|
||||||
{
|
{
|
||||||
if (current_lifo >= MAX_LIFO_DEPTH) {
|
if (current_lifo >= MAX_LIFO_DEPTH) {
|
||||||
error->all("Attempt to push beyond stack limit <FixBoxRelax>");
|
error->all("Attempt to push beyond stack limit in fix box/relax");
|
||||||
return;
|
return;
|
||||||
}
|
}
|
||||||
|
|
||||||
@ -352,7 +352,7 @@ void FixBoxRelax::min_pushstore()
|
|||||||
void FixBoxRelax::min_popstore()
|
void FixBoxRelax::min_popstore()
|
||||||
{
|
{
|
||||||
if (current_lifo <= 0) {
|
if (current_lifo <= 0) {
|
||||||
error->all("Attempt to pop empty stack <FixBoxRelax>");
|
error->all("Attempt to pop empty stack in fix box/relax");
|
||||||
return;
|
return;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|||||||
Reference in New Issue
Block a user