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

Fix full multi stencil

Browse Source
This commit is contained in:
Joel Clemmer
2020-11-12 15:03:44 -07:00
parent b6dfc28e38
commit 6909839ff0
2 changed files with 6 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/nstencil_full_multi2_2d.cpp
View File

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

5
src/nstencil_full_multi2_3d.cpp
View File

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