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

Fixing bug in max rate threshold

Browse Source
This commit is contained in:
jtclemm
2023-07-15 16:09:53 -06:00
parent d6532d3550
commit 2ef326273d
1 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/fix_deform.cpp
View File

@ -1214,8 +1214,9 @@ void FixDeform::set_pressure()
} }
if (max_h_rate != 0) if (max_h_rate != 0)
if (fabs(set[i].ptarget) > max_h_rate) if (fabs(h_rate[i]) > max_h_rate)
h_rate[i] = max_h_rate * h_rate[i] / fabs(h_rate[i]); h_rate[i] = max_h_rate * h_rate[i] / fabs(h_rate[i]);
h_ratelo[i] = -0.5 * h_rate[i]; h_ratelo[i] = -0.5 * h_rate[i];
double offset = 0.5 * (domain->boxhi[i] - domain->boxlo[i]) * (1.0 + update->dt * h_rate[i]); double offset = 0.5 * (domain->boxhi[i] - domain->boxlo[i]) * (1.0 + update->dt * h_rate[i]);
@ -1378,7 +1379,7 @@ void FixDeform::set_iso()
} }
if (max_h_rate != 0) if (max_h_rate != 0)
if (fabs(set[6].ptarget) > max_h_rate) if (fabs(v_rate) > max_h_rate)
v_rate = max_h_rate * v_rate / fabs(v_rate); v_rate = max_h_rate * v_rate / fabs(v_rate);
set[6].cumulative_strain += update->dt * v_rate; set[6].cumulative_strain += update->dt * v_rate;