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

fixed case where no axial components are active

Browse Source 
git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@15160 f3b2605a-c512-4ea7-a41b-209d697bcdaa
This commit is contained in:
athomps
2016-06-09 18:31:54 +00:00
parent 50ac419d3f
commit 84de575cc7
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/fix_nh.cpp
View File

@ -2163,7 +2163,7 @@ void FixNH::compute_press_target()
p_target[i] = p_start[i] + delta * (p_stop[i]-p_start[i]); p_target[i] = p_start[i] + delta * (p_stop[i]-p_start[i]);
p_hydro += p_target[i]; p_hydro += p_target[i];
} }
p_hydro /= pdim; if (pdim > 0) p_hydro /= pdim;
if (pstyle == TRICLINIC) if (pstyle == TRICLINIC)
for (int i = 3; i < 6; i++) for (int i = 3; i < 6; i++)
@ -2215,7 +2215,7 @@ void FixNH::nh_omega_dot()
for (int i = 0; i < 3; i++) for (int i = 0; i < 3; i++)
if (p_flag[i]) if (p_flag[i])
mtk_term2 += omega_dot[i]; mtk_term2 += omega_dot[i];
mtk_term2 /= pdim * atom->natoms; if (pdim > 0) mtk_term2 /= pdim * atom->natoms;
} }
if (pstyle == TRICLINIC) { if (pstyle == TRICLINIC) {