* Bump tolerance for test_mass_mat_trig in test_grudge_sym_old

* Set tolerances in test_mass_trig (for old sym code) equal to ones in new-style tests

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

https://github.com/inducer/meshmode/pull/225

* 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>

* Add jump property to TracePair.

* Change jump --> diff

* Add simple smoke test for trace pair

* placate flake8

* Undo spurious chage to primitives.py

* Tweak TracePair docstrings

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

* Make nodal reductions mpi-aware

* Add clarifying remark for nodal reductions

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

* Use timestep estimation in mpi wave test

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.

* vectorize local_grad and weak_local_grad

* added stack argument to local_grad/weak_local_grad

* add test for gradient using op

* propagate stack argument

* tweak comment to avoid flake8 warning

* stack -> (not) nested

* add visualization to gradient test

* disable strong form gradient test for now

* minor cosmetic changes

* tweak local_div/weak_local_div vectorization to work with both nested and non-nested arrays

* suppress pylint false positives

* remove inverse_mass around grad/div in strong form tests

* add check for valid 'form' argument

* remove unneeded flux terms from strong form tests

* Tweak phrasing: "no jumps" instead of "smooth"

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