Trixi.jl

Trixi.jl copied to clipboard

The first steps towards differentiating through Trixi were made in #461. The next steps outlined there are - [x] Benchmark and improve performance of `jacobian_ad_forward` (#464) - [x] Allow differentiating...

ranocha

WIP: Resolve #166

1

# Work in progress - do not merge Closes #166 TODO: - [x] https://github.com/JuliaSIMD/LoopVectorization.jl/issues/262 - [x] https://github.com/JuliaSIMD/LoopVectorization.jl/issues/263 - [x] Consider changing the array type in `create_cache(mesh::TreeMesh{2}, equations, mortar_l2::LobattoLegendreMortarL2, uEltype)` when...

ranocha

WIP: Switch to P4est with CBinding v1

7

This PR prepares Trixi for working with trixi-framework/P4est.jl#42 or trixi-framework/P4est.jl#43. To use it, you need to manually `dev` a local clone of said PRs, e.g., from within a clone of...

sloede

This came up in https://github.com/trixi-framework/Trixi.jl/pull/907#discussion_r729051487: It would be nice to generalize data types used for visualization methods. For example, it would be awesome to be able to plot something like...

ranocha

Automate EC/ES tests

2

This came up one again in https://github.com/trixi-framework/Trixi.jl/pull/900#discussion_r725094522. If someone wants to earn a T-shirt, this would be a nice Hacktoberfest project. The idea would be to add an optional keyword...

ranocha

gregorgassner

Nonconservative boundary terms

3

The `split_form_kernel!` found in `solvers/dgsem_structured/dg_2d.jl` uses the matrix `derivative_split` from `dg.basis` which is zero along the main diagonal. In essence, this "peels off" the two volume contributions at the `[1,1]`...

andrewwinters5000

allow `GlmSpeedCallback` without `StepsizeCallback`

1

Is there a good way to use `RDPK3SpFSAL49()` with the existing MHD solvers? A naive swap results in errors. For example, running ```julia trixi_include("examples/unstructured_2d_dgsem/elixir_mhd_ec.jl", sol = solve(ode, RDPK3SpFSAL49(), abstol=1e-7, reltol=1e-7,...

jlchan

simplified interface

2

Hi, I am really interested in using your package for solving McKean-Vlasov equations describing the distribution of mean field of particles. All examples I have in mind are like (say...

rveltz

CI results are inconsistent for `DGMulti` MHD elixirs with `approximation_type=SBP()`

5

As noted in https://github.com/trixi-framework/Trixi.jl/pull/880, CI tests currently do not pass when running `elixir_mhd_weak_blast_wave.jl` with `approximation_type=SBP()`. To get tests to pass, a constant GLM cleaning speed was introduced in `elixir_mhd_weak_blast_wave_SBP.jl`. However,...

jlchan