From b9f2d23fa768cb79d87e34692768e621c9e3ffef Mon Sep 17 00:00:00 2001 From: sjplimp Date: Fri, 13 Apr 2007 15:24:36 +0000 Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@483 f3b2605a-c512-4ea7-a41b-209d697bcdaa --- src/compute_rotate_gran.cpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/compute_rotate_gran.cpp b/src/compute_rotate_gran.cpp index 10a55bf525..db31531b9a 100644 --- a/src/compute_rotate_gran.cpp +++ b/src/compute_rotate_gran.cpp @@ -49,7 +49,7 @@ void ComputeRotateGran::init() double ComputeRotateGran::compute_scalar() { - double **phiv = atom->phiv; + double **omega = atom->omega; double *radius = atom->radius; double *rmass = atom->rmass; int *mask = atom->mask; @@ -59,8 +59,8 @@ double ComputeRotateGran::compute_scalar() for (int i = 0; i < nlocal; i++) if (mask[i] & groupbit) - erot += (phiv[i][0]*phiv[i][0] + phiv[i][1]*phiv[i][1] + - phiv[i][2]*phiv[i][2]) * radius[i]*radius[i]*rmass[i]; + erot += (omega[i][0]*omega[i][0] + omega[i][1]*omega[i][1] + + omega[i][2]*omega[i][2]) * radius[i]*radius[i]*rmass[i]; MPI_Allreduce(&erot,&scalar,1,MPI_DOUBLE,MPI_SUM,world); scalar *= pfactor;