-
Notifications
You must be signed in to change notification settings - Fork 24
made first order actually first order and added fourth order #120
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Conversation
| [(-25, -25), (0, 0), (0, 0), (1, 0)], | ||
| [(-25, -25), (-13, -13), (0, 0), (1, 1)], | ||
| [(-25, -25), (0, 0), (0, 0), (1, 1)], | ||
| [(-25, -25), (1, 0), (0, 0), (1, 1)], |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The first order bounds only get a touch worse when we make first_order truly first order.
pynumdiff/tests/test_diff_methods.py
Outdated
| [(1, 0), (1, 0), (1, 1), (1, 1)], | ||
| [(1, 0), (1, 1), (1, 0), (1, 1)], | ||
| [(2, 2), (3, 2), (2, 2), (3, 2)], | ||
| [(2, 2), (3, 3), (2, 2), (3, 3)]], |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The "iterated_first_order" bounds get worse all over the place, because the horizontal bias introduced by a non symmetric first-order method causes drift over the iterations.
|
The |
|
Let's merge this one after #125. I think I'm going to want to make a couple follow-on commits to it. I need to:
|
Addressing #104. The resulting
first_orderis biased to one side, resulting in drift when it's iterated:This is likely why you attempted some forward-backward averaging, @florisvb.
Honestly, is it worth keeping around
first_order? My guess is it's not in wide usage, because the higher order methods are equally usable and right there for the taking. Iteratingfirst_ordermight be a thing people do, but in that case they're actually iteratingsecond_order, so we should update the code to reflect that and feed them a warning if they try to callfirst_order.