Skip to content

Conversation

@DanielDoehring
Copy link
Member

Having merged #2488 one could realize #2472 and (subsequently) #2149 via this. In fact, VolumeIntegralShockCapturingHG could also be brought into this framework.

The benefit is greatly reduced computation time of the volume integral, more profound for higher polynomial degrees. For the examples employed in this (draft) PR, the volume integral computation can be roughly reduced by a factor of two.

Comment on lines +193 to +194
Possible combination: [`VolumeIntegralWeakForm`](@ref) and [`VolumeIntegralFluxDifferencing`](@ref).
TODO: Extend to [`VolumeIntegralWeakForm`](@ref) and [`VolumeIntegralShockCapturingHG`](@ref).
Copy link
Member

Choose a reason for hiding this comment

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

There are some TODO notes left. How do you want to handle them?

Copy link
Member Author

@DanielDoehring DanielDoehring Oct 20, 2025

Choose a reason for hiding this comment

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

This is probably something for the Trixi meeting. The main question is if VolumeIntegralShockCapturingHG should be (internally) re-implemented as a VolumeIntegralAdaptive.

Copy link
Member

Choose a reason for hiding this comment

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

👍

"""
VolumeIntegralAdaptive(indicator;
volume_integral_default = VolumeIntegralWeakForm(),
volume_integral_stabilized = VolumeIntegralFluxDifferencing(flux_central))
Copy link
Member

Choose a reason for hiding this comment

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

Shall we really prove this as default value? Or should we force people to make a reasonable choice?

@threaded for element in eachelement(dg, cache)
stabilized_version = decision[element]

# TODO: Generalize/Dispatch or introduce yet sub-functions of the volume integrals
Copy link
Member

Choose a reason for hiding this comment

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

Good question

@ranocha
Copy link
Member

ranocha commented Nov 20, 2025

Please ping me if you would like to have another round of feedback on this PR.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

enhancement New feature or request

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants