Hendrik Ranocha

icon indicating a website link https://ranocha.de

icon company Johannes Gutenberg University Mainz icon location Mainz, Germany icon location Professor analyzing & developing numerical methods for differential equations, focusing on stability, mimetic properties & structure preservation

Results 595 comments of Hendrik Ranocha

Multithreaded simulations freeze sometimes

Thanks for testing it again.

User defined RHS Splitting for IMEX

You can also add a simple unit test for that.

User defined RHS Splitting for IMEX

@MarcoArtiano Could you please fix the merge conflicts and consider reducing the boilerplate code as suggested by Daniel?

Allow passing `io` to `Analysis`, `Alive`, and `Summary` `Callback`s

> We could also make the `SummaryCallback` an own `struct` similar to the `AnalysisCallback`. Not sure if this counts as breaking because `SummaryCallback()` now returns a different object, but there...

Allow passing `io` to `Analysis`, `Alive`, and `Summary` `Callback`s

> > > The Documenter.jl bug is a bit unfortunate because I don't see any good way to work around this > > > > > > Yeah looks like...

Revisit `EulerGravity` with LDG, `linear_structure` & Krylov.jl

See also https://github.com/trixi-framework/Trixi.jl/pull/419

TrixiStateVector for distributed matrix-free implicit integration

Can we also use such a `TrixiStateVector` to broadcast array operations in the ODE solver automatically, even when running on a single MPI rank but with multiple threads, respecting the...

Lower-Bound/Positivity-Preserving Limiter of Rueda-Ramirez & Gassner (1D)

What is the status of this PR - waiting for a review (by @amrueda)?

RDPK3Sp35(), RDPK3Sp49(), and RDPK3Sp510() not working with Float32

The difference comes from https://github.com/SciML/OrdinaryDiffEq.jl/blob/3c55c01f57d6cc32fba2b3c1130a779bf692a529/lib/OrdinaryDiffEqCore/src/initdt.jl#L39-L54 The working methods, such as `RDPK3SpFSAL49`, choose the first branch. The others end up in the second branch.

RDPK3Sp35(), RDPK3Sp49(), and RDPK3Sp510() not working with Float32

In any case, I would recommend using `Float32` for the time span as well. Otherwise, the intermediate results may be computed using `Float64`: ```julia julia> ode = ODEProblem(MtlArray(ones(Float32, 4)), (0.0,...