Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

TPE for acoustic_pulse example #351

Open
wants to merge 4 commits into
base: main
Choose a base branch
from

Conversation

MTCam
Copy link
Collaborator

@MTCam MTCam commented Jun 13, 2024

This change set updates acoustic_pulse to have a tensor product elements option.

Currently, when running with tensor product elements --tpe and overintegration --oi, only eager works properly. When run in lazy mode --lazy, then DiscretizationCollection._base_to_geoderiv_connection (called from _geometry_to_quad_if_requested) does not correctly return a quantity on the quadrature discretization.

works: python acoustic_pulse.py --tpe
works: python acoustic_pulse.py --tpe --oi
fails: python acoustic_pulse.py --tpe --oi --lazy

This problem affects the use of overintegration in mirgecom; seemingly TPE go into processing that is only valid when element groups are affine and simplices. (5b563de steers around the issue for now).

Overintegration appears to negatively affect the convergence rate for TPEs (see 00b6ca4)

CC: @inducer

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant