Skip to content

Conversation

@vks
Copy link

@vks vks commented Oct 21, 2014

Can you tell me whether this looks good? (The code compiles, it is untested however.)

vks added 19 commits October 20, 2014 14:02
A wrong index was used, which resulted in nans.
Also made sure the code works for $p = 0$.
They have been calculated incorrectly. The charge-only balance function
still does not look right.
Instead output the numerator and denominator separately. This has the
advantage that the results of several separate runs can be easily
combined, because the division can be delayed until the last step, right
before plotting the balance function.
This ensures deterministic results while giving different results for
parallel runs.
The wrong signs were used, so B(pi+, pi-) was negative instead of
positive. This is fixed here.

Note that B(p, pbar) still seems to have the wrong sign for some reason.
This is a slight refactoring of Scott's and Garry's code.

Use switch statements where appropriate (in one case this should
improve performance) and return structs instead of modifying parameters.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant