dg1d: Add Bsplines1 (!215)
Continuation of !210 This adds a dispatch type `Bspline1` similar to `Bspline2`. Based on this we implement the necessary method to be able to use it in `DGElement`. I also fixed a subtle bug in the quadrature rules I used in some of the `Bspline` methods. That is, I accidentally used a LG rule with `Npts=2`, which corresponds to a `N=Npts-1=1` order rule that is exact for polynomials up to `2*N=2` order. But for `Bspline2` we encounter fourth order polynomials, so we should be using `Npts=3`. Also adds a new option to the `Mesh` parameters section: * `kind = "modal_bspline1"`: Utilize a 1st order modal Bspline Ansatz for the DG evolution.
parent
0ff6f2b0
No related branches found
No related tags found
Showing
- plot/bspline.jl 72 additions, 20 deletionsplot/bspline.jl
- src/ScalarEq/setup.jl 1 addition, 1 deletionsrc/ScalarEq/setup.jl
- src/bspline.jl 188 additions, 4 deletionssrc/bspline.jl
- src/dgelement.jl 46 additions, 15 deletionssrc/dgelement.jl
- src/parameters.jl 10 additions, 1 deletionsrc/parameters.jl
- test/IntegrationTests/refs/modal_bspline2_advection_sine/output1d.h5 0 additions, 0 deletions...ationTests/refs/modal_bspline2_advection_sine/output1d.h5
- test/IntegrationTests/refs/modal_bspline2_advection_sine/substep_output.h5 0 additions, 0 deletions...ests/refs/modal_bspline2_advection_sine/substep_output.h5
- test/UnitTests/src/test_dgelement.jl 36 additions, 2 deletionstest/UnitTests/src/test_dgelement.jl
Loading
Please register or sign in to comment