Commit 86810c2d authored by julient31's avatar julient31
Browse files

Commit2 JT 031919 

- correct error in fix_precession_spin
- only the sign of the force needed a correction
parent a508f1de

src/SPIN/fix_precession_spin.cpp

+2 −2
Original line number Diff line number Diff line
@@ -230,8 +230,8 @@ void FixPrecessionSpin::compute_single_precession(int i, double spi[3], double f 
void FixPrecessionSpin::compute_zeeman(int i, double fmi[3])

{

  double **sp = atom->sp;

  fmi[0] += sp[i][0]*hx;

  fmi[1] += sp[i][1]*hy;

  fmi[0] += sp[i][3]*hx;

  fmi[1] += sp[i][3]*hy;

  fmi[2] += sp[i][3]*hz;

}

CRA Git | Maintained and supported by SUSTech CRA and CCSE