remove dead code
This commit is contained in:
@ -474,7 +474,7 @@ double PairSpinExchangeBiquadratic::compute_energy(int i, int j, double rsq,
|
||||
int *type = atom->type;
|
||||
int itype,jtype;
|
||||
double Jex,Kex,ra,sdots;
|
||||
double rj,rk,r2j,r2k,ir3j,ir3k;
|
||||
double rj,rk,r2j,r2k;
|
||||
double energy = 0.0;
|
||||
itype = type[i];
|
||||
jtype = type[j];
|
||||
@ -482,10 +482,8 @@ double PairSpinExchangeBiquadratic::compute_energy(int i, int j, double rsq,
|
||||
ra = sqrt(rsq);
|
||||
rj = ra/J3[itype][jtype];
|
||||
r2j = rsq/J3[itype][jtype]/J3[itype][jtype];
|
||||
ir3j = 1.0/(rj*rj*rj);
|
||||
rk = ra/K3[itype][jtype];
|
||||
r2k = rsq/K3[itype][jtype]/K3[itype][jtype];
|
||||
ir3k = 1.0/(rk*rk*rk);
|
||||
|
||||
Jex = 4.0*J1_mech[itype][jtype]*r2j;
|
||||
Jex *= (1.0-J2[itype][jtype]*r2j);
|
||||
|
||||
Reference in New Issue
Block a user