diff --git a/src/ASPHERE/pair_line_lj.cpp b/src/ASPHERE/pair_line_lj.cpp index ba565c5922..7e6979f0a5 100644 --- a/src/ASPHERE/pair_line_lj.cpp +++ b/src/ASPHERE/pair_line_lj.cpp @@ -74,7 +74,6 @@ void PairLineLJ::compute(int eflag, int vflag) double xi[2],xj[2],fi[2],dxi,dxj,dyi,dyj; int *ilist,*jlist,*numneigh,**firstneigh; - evdwl = 0.0; ev_init(eflag,vflag); double **x = atom->x; diff --git a/src/ASPHERE/pair_tri_lj.cpp b/src/ASPHERE/pair_tri_lj.cpp index 79253140c0..5f62d211db 100644 --- a/src/ASPHERE/pair_tri_lj.cpp +++ b/src/ASPHERE/pair_tri_lj.cpp @@ -74,7 +74,6 @@ void PairTriLJ::compute(int eflag, int vflag) double dc1[3],dc2[3],dc3[3]; int *ilist,*jlist,*numneigh,**firstneigh; - evdwl = 0.0; ev_init(eflag,vflag); AtomVecTri::Bonus *bonus = avec->bonus; diff --git a/src/BODY/compute_body_local.cpp b/src/BODY/compute_body_local.cpp index 940250d160..6443f95325 100644 --- a/src/BODY/compute_body_local.cpp +++ b/src/BODY/compute_body_local.cpp @@ -123,7 +123,7 @@ void ComputeBodyLocal::compute_local() int ncount = compute_body(0); if (ncount > nmax) reallocate(ncount); size_local_rows = ncount; - ncount = compute_body(1); + compute_body(1); } /* ---------------------------------------------------------------------- diff --git a/src/BODY/compute_temp_body.cpp b/src/BODY/compute_temp_body.cpp index 6903f7a25c..210bb07417 100644 --- a/src/BODY/compute_temp_body.cpp +++ b/src/BODY/compute_temp_body.cpp @@ -343,7 +343,6 @@ void ComputeTempBody::compute_vector() inertia = bonus[body[i]].inertia; quat = bonus[body[i]].quat; - massone = rmass[i]; // wbody = angular velocity in body frame diff --git a/src/CLASS2/angle_class2.cpp b/src/CLASS2/angle_class2.cpp index 28d081c5c9..616d191d1c 100644 --- a/src/CLASS2/angle_class2.cpp +++ b/src/CLASS2/angle_class2.cpp @@ -447,10 +447,6 @@ double AngleClass2::single(int type, int i1, int i2, int i3) if (c > 1.0) c = 1.0; if (c < -1.0) c = -1.0; - double s = sqrt(1.0 - c*c); - if (s < SMALL) s = SMALL; - s = 1.0/s; - double dtheta = acos(c) - theta0[type]; double dtheta2 = dtheta*dtheta; double dtheta3 = dtheta2*dtheta;