It looks like when all_eigs=False, add_rn_eig_jump randomly changes either the gamma or amplitude of each psr's RN:
|
which_eig = np.random.choice(2, size=Npsr) |
I think it is intended to pick a pulsar to perturb randomly.
With that said, we may be only using it with all_eigs=True right now.