k1 bug in dipole-quadrupole integrator

Issue #353 resolved
Laurie Nevay created an issue

Seems related to negative values of k1.

Comments (1)

  1. William Shields

    Now fixed. The issue was only for negative K1 values. The strength parameter was calculated from literature but produced results inconsistent with PTC. It's now calculated differently for 2 cases; (-1/rho^2) < K1 < 0, and for K1 < (-1/rho^2).

  2. Log in to comment