ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

change radthresh initialization

Browse Source
This commit is contained in:
Axel Kohlmeyer
2024-04-10 09:52:58 -04:00
parent 83cd20dfeb
commit 7ca4bf05cf
1 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
src/create_atoms.cpp
View File

@ -157,9 +157,9 @@ void CreateAtoms::command(int narg, char **arg)
maxtry = DEFAULT_MAXTRY; maxtry = DEFAULT_MAXTRY;
radscale = 1.0; radscale = 1.0;
mesh_style = BISECTION; mesh_style = BISECTION;
radthresh = domain->lattice->xlattice; // NOTE to Axel - I think this should be 1.0 by default radthresh = 1.0;
mesh_density = 1.0; // similar to how this setting is 1.0 mesh_density = 1.0;
// see rescaling of both below if units = lattice
nbasis = domain->lattice->nbasis; nbasis = domain->lattice->nbasis;
basistype = new int[nbasis]; basistype = new int[nbasis];
for (int i = 0; i < nbasis; i++) basistype[i] = ntype; for (int i = 0; i < nbasis; i++) basistype[i] = ntype;
@ -375,8 +375,8 @@ void CreateAtoms::command(int narg, char **arg)
xone[2] *= domain->lattice->zlattice; xone[2] *= domain->lattice->zlattice;
} else if (style == RANDOM) { } else if (style == RANDOM) {
if (overlapflag) overlap *= domain->lattice->xlattice; if (overlapflag) overlap *= domain->lattice->xlattice;
} else if (style == MESH) { // NOTE to Axel - here is the rescaling of both params } else if (style == MESH) {
if (mesh_style == BISECTION) { // by lattice spacings if units = lattice, similar to xone,overlap if (mesh_style == BISECTION) {
radthresh *= domain->lattice->xlattice; radthresh *= domain->lattice->xlattice;
} else if (mesh_style == QUASIRANDOM) { } else if (mesh_style == QUASIRANDOM) {
mesh_density /= (domain->lattice->xlattice * domain->lattice->xlattice); mesh_density /= (domain->lattice->xlattice * domain->lattice->xlattice);