[aspect-devel] Visualizing actual solution

Wolfgang Bangerth bangerth at tamu.edu
Thu Oct 29 12:40:23 PDT 2015

On 10/29/2015 01:22 PM, Jonathan Perry-Houts wrote:
> I guess it's not super important, but is there an easy way to
> interrogate an FEValues object for the solution at an arbitrary point?

Yes, you just need to create a quadrature formula that includes this one 

> Even just outputting a regularly spaced grid of point values would make
> a nice, mesh-independent file type that could be plotted with gnuplot,
> or matplotlib, or something.

That's basically what the subdivide option to DataOut::build_patches 
does -- it subdivides each cell into a regular set of sub-cells, each of 
which can then be visualized.


Wolfgang Bangerth               email:            bangerth at math.tamu.edu
                                 www: http://www.math.tamu.edu/~bangerth/

More information about the Aspect-devel mailing list