tx258/lammps-gran-kokkos
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Changing particle size metric to neighcutsq

Browse Source
This commit is contained in:
Joel Clemmer
2020-11-22 20:40:47 -07:00
parent 00a74558ff
commit de5df539c9
6 changed files with 10 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/nstencil_full_multi2_2d.cpp
View File

@ -51,7 +51,7 @@ void NStencilFullMulti22d::set_stencil_properties()
stencil_half[i][j] = 0;
stencil_skip[i][j] = 0;
if(cuttypesq[i] <= cuttypesq[j]){
if(cutneighsq[i][i] <= cutneighsq[j][j]){
stencil_cut[i][j] = sqrt(cutneighsq[j][j]);
stencil_bin_type[i][j] = j;
} else {

2
src/nstencil_full_multi2_3d.cpp
View File

@ -51,7 +51,7 @@ void NStencilFullMulti23d::set_stencil_properties()
stencil_half[i][j] = 0;
stencil_skip[i][j] = 0;
if(cuttypesq[i] <= cuttypesq[j]){
if(cutneighsq[i][i] <= cutneighsq[j][j]){
stencil_cut[i][j] = sqrt(cutneighsq[j][j]);
stencil_bin_type[i][j] = j;
} else {

4
src/nstencil_half_multi2_2d.cpp
View File

@ -49,12 +49,12 @@ void NStencilHalfMulti22d::set_stencil_properties()
for (i = 1; i <= n; i++) {
for (j = 1; j <= n; j++) {
if(i == j) continue;
if(cuttypesq[i] > cuttypesq[j]) continue;
if(cutneighsq[i][i] > cutneighsq[j][j]) continue;
stencil_skip[i][j] = 0;
stencil_cut[i][j] = sqrt(cutneighsq[i][j]);
if(cuttypesq[i] == cuttypesq[j]){
if(cutneighsq[i][i] == cutneighsq[j][j]){
stencil_half[i][j] = 1;
stencil_bin_type[i][j] = i;
} else {

4
src/nstencil_half_multi2_2d_tri.cpp
View File

@ -49,12 +49,12 @@ void NStencilHalfMulti22dTri::set_stencil_properties()
for (i = 1; i <= n; i++) {
for (j = 1; j <= n; j++) {
if(i == j) continue;
if(cuttypesq[i] > cuttypesq[j]) continue;
if(cutneighsq[i][i] > cutneighsq[j][j]) continue;
stencil_skip[i][j] = 0;
stencil_cut[i][j] = sqrt(cutneighsq[i][j]);
if(cuttypesq[i] == cuttypesq[j]){
if(cutneighsq[i][i] == cutneighsq[j][j]){
stencil_half[i][j] = 1;
stencil_bin_type[i][j] = i;
} else {

4
src/nstencil_half_multi2_3d.cpp
View File

@ -49,12 +49,12 @@ void NStencilHalfMulti23d::set_stencil_properties()
for (i = 1; i <= n; i++) {
for (j = 1; j <= n; j++) {
if(i == j) continue;
if(cuttypesq[i] > cuttypesq[j]) continue;
if(cutneighsq[i][i] > cutneighsq[j][j]) continue;
stencil_skip[i][j] = 0;
stencil_cut[i][j] = sqrt(cutneighsq[i][j]);
if(cuttypesq[i] == cuttypesq[j]){
if(cutneighsq[i][i] == cutneighsq[j][j]){
stencil_half[i][j] = 1;
stencil_bin_type[i][j] = i;
} else {

4
src/nstencil_half_multi2_3d_tri.cpp
View File

@ -49,12 +49,12 @@ void NStencilHalfMulti23dTri::set_stencil_properties()
for (i = 1; i <= n; i++) {
for (j = 1; j <= n; j++) {
if(i == j) continue;
if(cuttypesq[i] > cuttypesq[j]) continue;
if(cutneighsq[i][i] > cutneighsq[j][j]) continue;
stencil_skip[i][j] = 0;
stencil_cut[i][j] = sqrt(cutneighsq[i][j]);
if(cuttypesq[i] == cuttypesq[j]){
if(cutneighsq[i][i] == cutneighsq[j][j]){
stencil_half[i][j] = 1;
stencil_bin_type[i][j] = i;
} else {