I have carefully read the codes of fix wall/lj126 of LAMMPS Jan08 . In function of void FixWallLJ126pore::post_force(int vflag), the force between wall and particle reads that fwall = r6inv*(coeff1*r6inv - coeff2) * side, that is, fwall =epsilon*[48*(sigma/r)^12-24*(sigma/r)^6]*side. But I think it is Fwall =epsilon*[48*(sigma/r)^12-24*(sigma/r)^6]* (1/r) *side. I am confused. Maybe I can’t completely understand the code.
If I code the Steele potential in a cylinder wall. Need I only modify the length of matrix wall and wall_all in FixWallLJ93::compute_scalar() and FixWallLJ93::compute_vector(int n)?
What’s the function of MPI_Allreduce(wall,wall_all,4,MPI_DOUBLE,MPI_SUM,world).
Thanks a million.