Skip to content

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