Commit 6b289b07 authored by Axel Kohlmeyer's avatar Axel Kohlmeyer
Browse files

change incorrect EINERTIA constant in rigid body integrators from 4.0 to 2.0... 

change incorrect EINERTIA constant in rigid body integrators from 4.0 to 2.0 (same as in other integrators)
parent 0ecdb998

src/RIGID/fix_rigid.cpp

+1 −1
Original line number Diff line number Diff line
@@ -50,7 +50,7 @@ enum{ISO,ANISO,TRICLINIC}; 
#define EPSILON 1.0e-7



#define SINERTIA 0.4            // moment of inertia prefactor for sphere

#define EINERTIA 0.4            // moment of inertia prefactor for ellipsoid

#define EINERTIA 0.2            // moment of inertia prefactor for ellipsoid

#define LINERTIA (1.0/12.0)     // moment of inertia prefactor for line segment



/* ---------------------------------------------------------------------- */

src/RIGID/fix_rigid_small.cpp

+1 −1
Original line number Diff line number Diff line
@@ -54,7 +54,7 @@ FixRigidSmall *FixRigidSmall::frsptr; 
#define BIG 1.0e20



#define SINERTIA 0.4            // moment of inertia prefactor for sphere

#define EINERTIA 0.4            // moment of inertia prefactor for ellipsoid

#define EINERTIA 0.2            // moment of inertia prefactor for ellipsoid

#define LINERTIA (1.0/12.0)     // moment of inertia prefactor for line segment



#define DELTA_BODY 10000

src/USER-OMP/fix_rigid_omp.cpp

+1 −1
Original line number Diff line number Diff line
@@ -39,7 +39,7 @@ using namespace MathConst; 


enum{SINGLE,MOLECULE,GROUP};	// same as in FixRigid



#define EINERTIA 0.4            // moment of inertia prefactor for ellipsoid

#define EINERTIA 0.2            // moment of inertia prefactor for ellipsoid



typedef struct { double x,y,z; } dbl3_t;

src/USER-OMP/fix_rigid_small_omp.cpp

+1 −1
Original line number Diff line number Diff line
@@ -37,7 +37,7 @@ using namespace LAMMPS_NS; 
using namespace FixConst;

using namespace MathConst;



#define EINERTIA 0.4            // moment of inertia prefactor for ellipsoid

#define EINERTIA 0.2            // moment of inertia prefactor for ellipsoid



enum{FULL_BODY,INITIAL,FINAL,FORCE_TORQUE,VCM_ANGMOM,XCM_MASS,ITENSOR,DOF};

CRA Git | Maintained and supported by SUSTech CRA and CCSE