Improve new trap int weights convergence/behaviour
David Dickinson
Branch: feature/improve_new_trap_int_convergence
Branch: 8.1-RC
Merged
Merged pull request
Merged in feature/improve_new_trap_int_convergence (pull request #397)
Merged in feature/improve_new_trap_int_convergence (pull request #397)
Â
Add an integration test for new_trap_int
Currently this fails due to the mod vpa and related cases
Split the new_trap_int weights calculation into two sub-domains
Found to significantly reduce error in test cases
Add a high resolution tight tolerance test case
Builds on PR #395
This change is not reflected in the calculation of the trapped error using in get_verr. Not clear if this is important.
Â
Update to remove the change to the weights calculation pending further investigation. This PR is now just refactoring. See PR #523 for code reinstating this change.