Skip to content

investigated #116 and added some comments to help clarify the code#123

Merged
pavelkomarov merged 1 commit intomasterfrom
higher-order-tvr-fd
Jul 2, 2025
Merged

investigated #116 and added some comments to help clarify the code#123
pavelkomarov merged 1 commit intomasterfrom
higher-order-tvr-fd

Conversation

@pavelkomarov
Copy link
Collaborator

No description provided.

dxdt_hat_f = dxdt_hat[-1] + ddxdt_hat_f # What is this doing? Could we use a 2nd order FD above natively?
dxdt_hat = np.hstack((dxdt_hat, dxdt_hat_f))

# fix first point
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Having to find a special last value makes sense, but I experimented for a while and determined the first point doesn't really need to be fixed. In fact, having these lines increased the error for a few of the test functions (the ones with curvature). Zooming in, you could see it in the plots.

@pavelkomarov pavelkomarov merged commit 89f3109 into master Jul 2, 2025
1 check passed
@pavelkomarov pavelkomarov deleted the higher-order-tvr-fd branch July 2, 2025 01:24
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