[aspect-devel] Loading a large list of semi-unstructured point data into Aspect

[Jonathan Robey]

I am working with a Richardson Extrapolation based convergence estimator
for Aspect. Currently I need to find a better way to transfer large
amounts of quadrature point data between runs, when the triangulation
may not be consistent between the two. Do you have any advice you might
offer?

Currently, the plugin goes through the list of quadrature points and
finds the active cell around each using GridTools to do the computation,
but would you know of a more efficient method as, based on qualitative
runtime costs, I think this scales non-linearly, and perhaps not even at
log(n) per lookup.

The code for the plugin is at
https://github.com/class4kayaker/aspect/tree/Richardson_extrapolation_merge
for only the FE portions of the plugin, and
https://github.com/class4kayaker/aspect/tree/vof_richardson for the
version with some specialized code for estimating VOF convergence rates.

Thank you,

Jonathan Robey