Compute Q-criterion in Ascent and slice multiple times

Cases
, ,
2

I thought I will be a bit clever and a bit not, by just computing q-criterion explicitly as a field as per @WillT suggestion here for vmag:

But when passing any gradient Ascent plugin:

field-qcrit = -0.5*(grad_u_x)
; that's what I thought I need to do: grad_u_x**2 + grad_v_y**2 + grad_w_z**2 + 2*grad_u_y*grad_v_x + 2*grad_u_z*grad_w_x + 2*grad_v_z*grad_w_y)

pipeline-qslice-xy = [{'type': 'slice', 'point': {'x': 0.0, 'y': 0.0, 'z': 0.0}, 'normal': {'x': 0.0, 'y': 0.0, 'z': 1.0}}]
scene-qcrit-xy = {'type': 'pseudocolor', 'pipeline': 'qslice-xy',
                   'field': 'qcrit',
                   'render-1': {'image-width': 1024, 'image-name': 'qcrit-xy-{t:.4f}'}}

I get an error!

  File "/gpfs/mcc/HT07881/uxq12/shared/venvs/pyfr-mpich-pmix/lib/python3.13/site-packages/pyfr/__main__.py", line 533, in _process_common
    solver = get_solver(backend, mesh, soln, cfg)
  File "/gpfs/mcc/HT07881/uxq12/shared/venvs/pyfr-mpich-pmix/lib/python3.13/site-packages/pyfr/solvers/__init__.py", line 14, in get_solver
    return get_integrator(backend, systemcls, mesh, initsoln, cfg)
  File "/gpfs/mcc/HT07881/uxq12/shared/venvs/pyfr-mpich-pmix/lib/python3.13/site-packages/pyfr/integrators/__init__.py", line 34, in get_integrator
    return integrator(backend, systemcls, mesh, initsoln, cfg)
  File "/gpfs/mcc/HT07881/uxq12/shared/venvs/pyfr-mpich-pmix/lib/python3.13/site-packages/pyfr/integrators/std/controllers.py", line 24, in __init__
    self._run_plugins()
    ~~~~~~~~~~~~~~~~~^^
  File "/gpfs/mcc/HT07881/uxq12/shared/venvs/pyfr-mpich-pmix/lib/python3.13/site-packages/pyfr/integrators/base.py", line 117, in _run_plugins
    plugin(self)
    ~~~~~~^^^^^^
  File "/gpfs/mcc/HT07881/uxq12/shared/venvs/pyfr-mpich-pmix/lib/python3.13/site-packages/pyfr/plugins/ascent.py", line 554, in __call__
    self._renderer.render(_IntegratorAdapter(intg, self.cfgsect))
    ~~~~~~~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/gpfs/mcc/HT07881/uxq12/shared/venvs/pyfr-mpich-pmix/lib/python3.13/site-packages/pyfr/plugins/ascent.py", line 531, in render
    self._eval_exprs(adapter)
    ~~~~~~~~~~~~~~~~^^^^^^^^^
  File "/gpfs/mcc/HT07881/uxq12/shared/venvs/pyfr-mpich-pmix/lib/python3.13/site-packages/pyfr/plugins/ascent.py", line 493, in _eval_exprs
    grad_soln = soln_op @ grad_soln
                ~~~~~~~~^~~~~~~~~~~
ValueError: setting an array element with a sequence. The requested array has an inhomogeneous shape after 4 dimensions. The detected shape was (2, 3, 1, 5) + inhomogeneous part.

I appreciate this is a different question, but are gradients supported in current PyFR Ascent plugin? This is all on mixed mesh.