remove dead code
This commit is contained in:
@ -74,7 +74,6 @@ void PairLineLJ::compute(int eflag, int vflag)
|
|||||||
double xi[2],xj[2],fi[2],dxi,dxj,dyi,dyj;
|
double xi[2],xj[2],fi[2],dxi,dxj,dyi,dyj;
|
||||||
int *ilist,*jlist,*numneigh,**firstneigh;
|
int *ilist,*jlist,*numneigh,**firstneigh;
|
||||||
|
|
||||||
evdwl = 0.0;
|
|
||||||
ev_init(eflag,vflag);
|
ev_init(eflag,vflag);
|
||||||
|
|
||||||
double **x = atom->x;
|
double **x = atom->x;
|
||||||
|
|||||||
@ -74,7 +74,6 @@ void PairTriLJ::compute(int eflag, int vflag)
|
|||||||
double dc1[3],dc2[3],dc3[3];
|
double dc1[3],dc2[3],dc3[3];
|
||||||
int *ilist,*jlist,*numneigh,**firstneigh;
|
int *ilist,*jlist,*numneigh,**firstneigh;
|
||||||
|
|
||||||
evdwl = 0.0;
|
|
||||||
ev_init(eflag,vflag);
|
ev_init(eflag,vflag);
|
||||||
|
|
||||||
AtomVecTri::Bonus *bonus = avec->bonus;
|
AtomVecTri::Bonus *bonus = avec->bonus;
|
||||||
|
|||||||
@ -123,7 +123,7 @@ void ComputeBodyLocal::compute_local()
|
|||||||
int ncount = compute_body(0);
|
int ncount = compute_body(0);
|
||||||
if (ncount > nmax) reallocate(ncount);
|
if (ncount > nmax) reallocate(ncount);
|
||||||
size_local_rows = ncount;
|
size_local_rows = ncount;
|
||||||
ncount = compute_body(1);
|
compute_body(1);
|
||||||
}
|
}
|
||||||
|
|
||||||
/* ----------------------------------------------------------------------
|
/* ----------------------------------------------------------------------
|
||||||
|
|||||||
@ -343,7 +343,6 @@ void ComputeTempBody::compute_vector()
|
|||||||
|
|
||||||
inertia = bonus[body[i]].inertia;
|
inertia = bonus[body[i]].inertia;
|
||||||
quat = bonus[body[i]].quat;
|
quat = bonus[body[i]].quat;
|
||||||
massone = rmass[i];
|
|
||||||
|
|
||||||
// wbody = angular velocity in body frame
|
// wbody = angular velocity in body frame
|
||||||
|
|
||||||
|
|||||||
@ -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;
|
||||||
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 dtheta = acos(c) - theta0[type];
|
||||||
double dtheta2 = dtheta*dtheta;
|
double dtheta2 = dtheta*dtheta;
|
||||||
double dtheta3 = dtheta2*dtheta;
|
double dtheta3 = dtheta2*dtheta;
|
||||||
|
|||||||
Reference in New Issue
Block a user