* Add overintegration option for wave-op-mpi

* Add OI+lazy to examples CI testing

* Draft: Use P^0 discr on affine groups even if quad grid requested

* Refactor geoderiv connection/interpolation steps

* Fix quadrature interpolation

* Add pytato pytest array context factory

* Update metric unit tests

* Clarifying comment

Co-authored-by: Andreas Klöckner <inform@tiker.net>

* Add nonaffine option for wave op example

* Run lazy + nonaffine + overintegration in examples CI testing

* Add quadrature case for metric unit tests

* Add optional import for debugging outside of pytest

* Add docstrings for _use_geoderiv_connection

Co-authored-by: Andreas Kloeckner <inform@tiker.net>

Alexandru Fikl authored 3 years ago and Andreas Klöckner committed 3 years ago

gmsh generation was failing in v4.9 and couldn't find a setup
that worked in v4.8 and v4.9. Possible bug?

* reimplement norm() without recursion

* disallow norm of a scalar; fix tests

* Update function: dt_non_geometric_factor -> dt_non_geometric_factors

* Add timestep estimation helper function

* Add clarifying remark about timestep estimation

Co-authored-by: Andreas Klöckner <inform@tiker.net>

* Make wavespeed argument option

* Drop dt estimator; provide characteristic lengthscales

* Expand and update documentation

Co-authored-by: Andreas Klöckner <inform@tiker.net>

* Allow HyperbolicOperator.max_characteristic_velocity to return local data

* Add warning to dt_geometric_factors about non-volume discretizations

* Add fudge factor to surface advection time step

* Time step estimation: clean up freeze/thaw semantics

* Wave models: clean up freeze/thaw

* Variable coefficient wave model: do not nodal-reduce characteristic velocity

* Add missing freeze(...) for dim=1 case

* Add missing thaw in dt utils unit test

* Fix variable coefficient wave example

* Memoize characteristic lengths rather than geometric factors

* Use actx from wave speed DOFArray in variable coefficient operator

* Pass array context to estimate_rk4_timestep

* Ensure application actx is passed through to timestep estimation

* Update Operator method interfaces for timestepping

* Always return thawed characteristic lengthscales

Co-authored-by: Andreas Klöckner <inform@tiker.net>

* Simplify WADG mass inverse

* Test mass inverse on warped geometry

* Always WADG all the time also for symbolic

* Fold test_nonaffine_mass_operator_inverse into test_mass_operator_inverse

Co-authored-by: Thomas Gibson <gibsonthomas1120@hotmail.com>

- Rename, document the supporting functions for improved clarity
- Drop dim argument from (former) surface_normal

The function registry bits are only relevant for symbolic operators,
and the broadcasting test duplicates what's already being done in
arraycontext.

This ensures that the symbolic stuff remains tested until we remove it.