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

The issue is that the same expression (dx0/dr0) appears in
both a user provided expression and the discretization
metrics and the code generator gets confused.

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

Test that normal DOFArrays for empty boundaries have the same number of groups as the corresponding discretization

For 687a28bf: Drop StrongWaveOperator, fix WeakWaveOperator, retool dagrt-fusion to use that directly