Commit af1e7972 authored by Pierre de Buyl's avatar Pierre de Buyl
Browse files

fix mass and nve/asphere 

1. Set mass to 1. The radius is 1/2, so the density must be
   1/(4/3 pi 1/2^3)
2. Use nve/asphere instead of nve. Else, the particle's orientation are
   not integrated.
parent 074e18d9

examples/USER/misc/propel_self/3d_quaternion/in.3d_quaternion

+2 −2
Original line number Diff line number Diff line
@@ -18,11 +18,11 @@ pair_coeff * * 1.0 1.0 
pair_modify shift yes



# mass * 1.0

set type * density 1.0

set type * shape 1.0 1.0 1.0

set type * density 1.9098593171027443

set type * quat 0 0 1 0



fix step all nve

fix step all nve/asphere

fix temp all langevin 1.0 1.0 1.0 13 angmom 3.333333333



neighbor 0.6 bin

CRA Git | Maintained and supported by SUSTech CRA and CCSE