I would like to visualize simulation results, obtained using the discontinuous Galerkin (DG) approach, within ParaView. Similarly to finite volume methods, the problem domain is divided into cube-shaped cells ("elements"). As opposed to finite volume methods, within each cell there is not just one value for the solution vector $\mathbf{u}$, but each cell contains the solution $\mathbf{u}$ at multiple Gauss integration points.

My question is whether anyone has experience with visualizing such data efficiently with ParaView/VTK, and what approach you chose to represent the data in VTK. Several possible ways come to my mind, but I do not know which one is the most promising:

(1) Use voxels
Use one voxel for each integration point.
Pro: All plugins that work with the standard VTK unstructured cell types will continue to work without changing anything.
Con: Since the integration points are not distributed evenly, it might be difficult to find the correct location of the vertices. Also, the solution can be defined twice on the cell surfaces, as the DG framework allows discontinuous solutions. Also, the hierarchical information (domain divided into elements, each element contains several points) is lost.

(2) Use polyvertices
Use one vertex per integration point.
Pro: Easiest to implement, easy to specify multiple points at the same location with different solutions.
Con: Capability to visualize data as "cells" is lost, plus the same disadvantages as above.

(3) Use VTK quadrature scheme
Use the built-in support for quadrature schemes.
Pro: Rather straightforward implementation, preserves all relations and properties of the original solution.
Con: Since this is a completely new cell type, many (most) of the existing plugins will not work anymore and will probably have to be rewritten.