From e66c29fc6617955d6d2897c31c2d13d360ba4d06 Mon Sep 17 00:00:00 2001 From: sjplimp Date: Mon, 30 Apr 2007 19:26:01 +0000 Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@545 f3b2605a-c512-4ea7-a41b-209d697bcdaa --- src/fix_ave_spatial.cpp | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/src/fix_ave_spatial.cpp b/src/fix_ave_spatial.cpp index eca11bdc3f..758a0c74e7 100644 --- a/src/fix_ave_spatial.cpp +++ b/src/fix_ave_spatial.cpp @@ -307,7 +307,7 @@ void FixAveSpatial::end_of_step() // perform the computation for one sample // sum within each layer, only include atoms in fix group - // DENSITY_MASS adds mass + // DENSITY_MASS adds mass to values // DENSITY_NUM adds 1 to values // ATOM adds atom vector to values // COMPUTE adds its vector to values @@ -323,7 +323,7 @@ void FixAveSpatial::end_of_step() for (i = 0; i < nlocal; i++) { if (mask[i] & groupbit) { - ilayer = static_cast ((x[dim][i] - offset) * invdelta); + ilayer = static_cast ((x[i][dim] - offset) * invdelta); if (ilayer < 0) ilayer = 0; if (ilayer >= nlayers) ilayer = nlayers-1; count_one[ilayer] += 1.0; @@ -335,7 +335,7 @@ void FixAveSpatial::end_of_step() } else if (which == DENSITY_NUM) { for (i = 0; i < nlocal; i++) { if (mask[i] & groupbit) { - ilayer = static_cast ((x[dim][i] - offset) * invdelta); + ilayer = static_cast ((x[i][dim] - offset) * invdelta); if (ilayer < 0) ilayer = 0; if (ilayer >= nlayers) ilayer = nlayers-1; count_one[ilayer] += 1.0; @@ -356,7 +356,7 @@ void FixAveSpatial::end_of_step() m = 0; for (i = 0; i < nlocal; i++) { if (mask[i] & groupbit) { - ilayer = static_cast ((x[dim][i] - offset) * invdelta); + ilayer = static_cast ((x[i][dim] - offset) * invdelta); if (ilayer < 0) ilayer = 0; if (ilayer >= nlayers) ilayer = nlayers-1; count_one[ilayer] += 1.0; @@ -374,7 +374,7 @@ void FixAveSpatial::end_of_step() m = 0; for (i = 0; i < nlocal; i++) { if (mask[i] & groupbit) { - ilayer = static_cast ((x[dim][i] - offset) * invdelta); + ilayer = static_cast ((x[i][dim] - offset) * invdelta); if (ilayer < 0) ilayer = 0; if (ilayer >= nlayers) ilayer = nlayers-1; count_one[ilayer] += 1.0;