Converting StructuredMesh to uniform grid about trixi.jl HOT 2 OPEN

You could average the solution over each cell. There may be a cleaner way, but the internal routine Trixi.wrap_array will reshape the solution vector into an array of size [num_fields, N+1, N+1, num_elements] in 2D.

The best way to average the solution is to use the quadrature nodes available in dg.basis or solver.basis to compute integral averages of the solution over each cell, but you could take the arithmetic average over each cell u[i, :, :, cell_index] as well.

from trixi.jl.

If you use trixi2vtk to convert the solution files it, by default, reinterpolates to the solution from the nnodes Legendre-Gauss-Lobatto points (which is polydeg+1) onto 2*nnodes that are uniform. This would explain your situation.

• If you want to change the number of uniform nodes onto which you reinterpolate then you can pass the additional argument nvisnodes to your trixi2vtk call, e.g., trixi2vtk("out/solution_files", nvisnodes=4).

• To deactivate this reinterpolation completely and make the solution on the original solver nodes available you can set additional argument reinterpolate to be false in your trixi2vtk call, e.g., trixi2vtk("out/solution_files", reinterpolate=false).

from trixi.jl.