Loading src/compute_angmom_chunk.cpp +3 −3 Original line number Diff line number Diff line Loading @@ -155,9 +155,9 @@ void ComputeAngmomChunk::compute_array() dz = unwrap[2] - comall[index][2]; if (rmass) massone = rmass[i]; else massone = mass[type[i]]; angmom[i][0] += massone * (dy*v[i][2] - dz*v[i][1]); angmom[i][1] += massone * (dz*v[i][0] - dx*v[i][2]); angmom[i][2] += massone * (dx*v[i][1] - dy*v[i][0]); angmom[index][0] += massone * (dy*v[i][2] - dz*v[i][1]); angmom[index][1] += massone * (dz*v[i][0] - dx*v[i][2]); angmom[index][2] += massone * (dx*v[i][1] - dy*v[i][0]); } MPI_Allreduce(&angmom[0][0],&angmomall[0][0],3*nchunk, Loading src/compute_omega_chunk.cpp +3 −3 Original line number Diff line number Diff line Loading @@ -180,9 +180,9 @@ void ComputeOmegaChunk::compute_array() dz = unwrap[2] - comall[index][2]; if (rmass) massone = rmass[i]; else massone = mass[type[i]]; angmom[i][0] += massone * (dy*v[i][2] - dz*v[i][1]); angmom[i][1] += massone * (dz*v[i][0] - dx*v[i][2]); angmom[i][2] += massone * (dx*v[i][1] - dy*v[i][0]); angmom[index][0] += massone * (dy*v[i][2] - dz*v[i][1]); angmom[index][1] += massone * (dz*v[i][0] - dx*v[i][2]); angmom[index][2] += massone * (dx*v[i][1] - dy*v[i][0]); } MPI_Allreduce(&angmom[0][0],&angmomall[0][0],3*nchunk, Loading src/compute_torque_chunk.cpp +3 −3 Original line number Diff line number Diff line Loading @@ -153,9 +153,9 @@ void ComputeTorqueChunk::compute_array() dx = unwrap[0] - comall[index][0]; dy = unwrap[1] - comall[index][1]; dz = unwrap[2] - comall[index][2]; torque[i][0] += dy*f[i][2] - dz*f[i][1]; torque[i][1] += dz*f[i][0] - dx*f[i][2]; torque[i][2] += dx*f[i][1] - dy*f[i][0]; torque[index][0] += dy*f[i][2] - dz*f[i][1]; torque[index][1] += dz*f[i][0] - dx*f[i][2]; torque[index][2] += dx*f[i][1] - dy*f[i][0]; } MPI_Allreduce(&torque[0][0],&torqueall[0][0],3*nchunk, Loading Loading
src/compute_angmom_chunk.cpp +3 −3 Original line number Diff line number Diff line Loading @@ -155,9 +155,9 @@ void ComputeAngmomChunk::compute_array() dz = unwrap[2] - comall[index][2]; if (rmass) massone = rmass[i]; else massone = mass[type[i]]; angmom[i][0] += massone * (dy*v[i][2] - dz*v[i][1]); angmom[i][1] += massone * (dz*v[i][0] - dx*v[i][2]); angmom[i][2] += massone * (dx*v[i][1] - dy*v[i][0]); angmom[index][0] += massone * (dy*v[i][2] - dz*v[i][1]); angmom[index][1] += massone * (dz*v[i][0] - dx*v[i][2]); angmom[index][2] += massone * (dx*v[i][1] - dy*v[i][0]); } MPI_Allreduce(&angmom[0][0],&angmomall[0][0],3*nchunk, Loading
src/compute_omega_chunk.cpp +3 −3 Original line number Diff line number Diff line Loading @@ -180,9 +180,9 @@ void ComputeOmegaChunk::compute_array() dz = unwrap[2] - comall[index][2]; if (rmass) massone = rmass[i]; else massone = mass[type[i]]; angmom[i][0] += massone * (dy*v[i][2] - dz*v[i][1]); angmom[i][1] += massone * (dz*v[i][0] - dx*v[i][2]); angmom[i][2] += massone * (dx*v[i][1] - dy*v[i][0]); angmom[index][0] += massone * (dy*v[i][2] - dz*v[i][1]); angmom[index][1] += massone * (dz*v[i][0] - dx*v[i][2]); angmom[index][2] += massone * (dx*v[i][1] - dy*v[i][0]); } MPI_Allreduce(&angmom[0][0],&angmomall[0][0],3*nchunk, Loading
src/compute_torque_chunk.cpp +3 −3 Original line number Diff line number Diff line Loading @@ -153,9 +153,9 @@ void ComputeTorqueChunk::compute_array() dx = unwrap[0] - comall[index][0]; dy = unwrap[1] - comall[index][1]; dz = unwrap[2] - comall[index][2]; torque[i][0] += dy*f[i][2] - dz*f[i][1]; torque[i][1] += dz*f[i][0] - dx*f[i][2]; torque[i][2] += dx*f[i][1] - dy*f[i][0]; torque[index][0] += dy*f[i][2] - dz*f[i][1]; torque[index][1] += dz*f[i][0] - dx*f[i][2]; torque[index][2] += dx*f[i][1] - dy*f[i][0]; } MPI_Allreduce(&torque[0][0],&torqueall[0][0],3*nchunk, Loading
