Hi All,

I am trying to calculate the "S_xy" component of stress tensor in two

different ways using the following commands.

variable natoms equal count(flow)

compute peratom flow stress/atom

compute sxy flow reduce sum c_peratom[4] -------------------> I

variable sxy equal c_sxy/${natoms} ---------------------------> II

I expect the following commands to produce the same results. But I get

different results.

I am doing anything wrong here.

Srikanth