Merge pull request #2770 from strlnkv/wall_region_morse_bug_fix
wall/region Morse potential derivative correction
This commit is contained in:
@ -397,7 +397,7 @@ void FixWallRegion::morse(double r)
|
|||||||
{
|
{
|
||||||
double dr = r - sigma;
|
double dr = r - sigma;
|
||||||
double dexp = exp(-alpha * dr);
|
double dexp = exp(-alpha * dr);
|
||||||
fwall = coeff1 * (dexp*dexp - dexp) / r;
|
fwall = coeff1 * (dexp*dexp - dexp);
|
||||||
eng = epsilon * (dexp*dexp - 2.0*dexp) - offset;
|
eng = epsilon * (dexp*dexp - 2.0*dexp) - offset;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|||||||
Reference in New Issue
Block a user