CarpetX
CarpetX copied to clipboard
CarpetX: fix bug in hermite prolongation
This PR is trying to fix bug in the Hermite type of prolongation.
- Since the first order Hermite interpolation is the same with Lagrange Interpolation and quit different from higher Hermite interpolations, we removed that coeffs and calling POLY directly for first order Hermite.
- We update the coeffs for higher Hermite interpolation according to the derivation here: https://github.com/lwJi/Docs-CarpetX-Public/blob/main/notes/Prolongation/Hermite.ipynb
- We update the size the stencils required by Hermite because it's two points larger compared to the same order Lagrangian.