Is outputting stress at vertices possible in Pylith? It seems that outputting at vertices is only possible for displacement and velocity.

In PyLith v2.2 stress is only available at the numerical quadrature points within a cell. Most visualization tools only understand a single point within a cell (they donâ€™t know the locations of the numerical quadrature points), so we typically use the average, which is available via `pylith.meshio.CellFilterAvg`

. See the examples and manual for details.

In the next PyLith release, it will be possible to project fields into a desire basis (for example, basis order 0 = cells, basis order 1 = vertices). Note that the field is only as accurate as the original discretization. We will try to explain this clearly in the manual.

I see, thank you very much for your help.