From 253a17b6d06fddd945411ff70ec275dc34cfc5ae Mon Sep 17 00:00:00 2001 From: abbatux <31911482+abbatux@users.noreply.github.com> Date: Mon, 23 Oct 2017 12:45:58 +1100 Subject: [PATCH] Update pair_smd_tlsph.cpp Correction of a typo in the computation of LAME_LAMBDA. --- src/USER-SMD/pair_smd_tlsph.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/USER-SMD/pair_smd_tlsph.cpp b/src/USER-SMD/pair_smd_tlsph.cpp index 1900519935..d6ddc9c72d 100644 --- a/src/USER-SMD/pair_smd_tlsph.cpp +++ b/src/USER-SMD/pair_smd_tlsph.cpp @@ -1066,7 +1066,7 @@ void PairTlsph::coeff(int narg, char **arg) { Lookup[HEAT_CAPACITY][itype] = force->numeric(FLERR, arg[ioffset + 7]); Lookup[LAME_LAMBDA][itype] = Lookup[YOUNGS_MODULUS][itype] * Lookup[POISSON_RATIO][itype] - / ((1.0 + Lookup[POISSON_RATIO][itype] * (1.0 - 2.0 * Lookup[POISSON_RATIO][itype]))); + / ((1.0 + Lookup[POISSON_RATIO][itype]) * (1.0 - 2.0 * Lookup[POISSON_RATIO][itype])); Lookup[SHEAR_MODULUS][itype] = Lookup[YOUNGS_MODULUS][itype] / (2.0 * (1.0 + Lookup[POISSON_RATIO][itype])); Lookup[M_MODULUS][itype] = Lookup[LAME_LAMBDA][itype] + 2.0 * Lookup[SHEAR_MODULUS][itype]; Lookup[SIGNAL_VELOCITY][itype] = sqrt(