pith. sign in

arxiv: 1907.07224 · v1 · pith:KHRFEWWZnew · submitted 2019-07-16 · 💻 cs.CE · cs.GR

The Effect of Data Transformations on Scalar Field Topological Analysis of High-Order FEM Solutions

Pith reviewed 2026-05-24 20:23 UTC · model grok-4.3

classification 💻 cs.CE cs.GR
keywords high-order FEMtopological analysisdata transformationscalar fieldsL-SIAC filterfinite element simulationsdiscontinuous dataanalysis pipelines
0
0 comments X

The pith

The choice of transformation method for high-order FEM data changes the topological features extracted from scalar fields.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tests how converting high-order finite element simulation outputs into continuous low-order data affects later topological analysis. It applies two standard transformation methods plus the L-SIAC filter to simulation results that contain element-boundary discontinuities and measures the resulting differences in extracted features. A reader would care because topological tools are now used to inspect complex flow simulations, yet those tools require transformed inputs whose side effects had not been compared side by side. The study finds that the same underlying data can yield visibly different topological structures depending on the chosen conversion step, so analysis pipelines must select the transformation deliberately rather than by default.

Core claim

The central claim is that transforming high-order FEM simulation data using different methodologies produces diverse behaviors in topological analysis. The empirical comparison of two commonly used transformation methodologies together with the L-SIAC filter shows that each produces distinct topological features from the same high-order scalar fields, and therefore pipelines built with current topological tools require explicit consideration of the transformation choice.

What carries the argument

Data transformation methodologies that convert discontinuous high-order polynomial data into continuous low-order representations for input to topological analysis tools.

If this is right

  • Topological features extracted from HO-FEM scalar fields vary with the transformation method applied before analysis.
  • The L-SIAC filter produces results that differ from the two standard transformation approaches in at least some test cases.
  • Analysis pipelines must evaluate the effect of the chosen transformation on the specific features of interest.
  • Inconsistent topological structures can appear from the same simulation data when different preparation steps are used.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Directly extending topological tools to accept high-order discontinuous inputs would remove the transformation step and its observed variability.
  • Quantifying how each transformation preserves or distorts specific topological invariants could guide method selection for particular analysis goals.
  • Applying the same comparison to vector or tensor fields from the same simulations would test whether the sensitivity is limited to scalar cases.

Load-bearing premise

The two commonly used transformation methodologies together with the L-SIAC filter are representative of possible transformations and the diverse behaviors will generalize to other high-order FEM datasets and topological tools.

What would settle it

Running the same topological extraction on a high-order FEM dataset whose analytic topology is known in advance and finding identical features after every transformation method would falsify the claim of diverse behaviors.

Figures

Figures reproduced from arXiv: 1907.07224 by Ashok Jallepalli, Joshua A. Levine, Robert M. Kirby.

Figure 1
Figure 1. Figure 1: Topological segmentation of counter-rotating vortex sampled using different methodologies discussed in the paper and by [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: The persistence curve, the persistence diagram, and the segmen [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Persistence curves of the sampled data (at three resolutions) as [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 2
Figure 2. Figure 2: The function f(x, y), simulation mesh, and the projected data. To create a higher order representation of f(x, y), we project it on an unstructured triangular element mesh (Figure 2b) over the domain [0,1] × [0,1] with degree 2 polynomial expansions on each element (hereafter referred to as the simulation mesh). With this simple example, we can mimic the elementwise discontinuities we get when computing de… view at source ↗
Figure 6
Figure 6. Figure 6: Persistence curves of the subdivided data (at three resolutions) [PITH_FULL_IMAGE:figures/full_fig_p005_6.png] view at source ↗
Figure 5
Figure 5. Figure 5: Visualizations of the persistence diagrams ((a), (b), and (c)) and [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: Visualizations of the persistence diagrams ((a), (b), and (c)) and [PITH_FULL_IMAGE:figures/full_fig_p005_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: Visualizations of the persistence diagrams ((a), (b), and (c)) and [PITH_FULL_IMAGE:figures/full_fig_p006_9.png] view at source ↗
Figure 8
Figure 8. Figure 8: Persistence curves of the L-SIAC data (at three resolutions) as [PITH_FULL_IMAGE:figures/full_fig_p006_8.png] view at source ↗
Figure 10
Figure 10. Figure 10: (a) Simulation mesh of the flow over a 2D cylinder. (b) Vorticity [PITH_FULL_IMAGE:figures/full_fig_p006_10.png] view at source ↗
Figure 13
Figure 13. Figure 13: The persistence diagrams of the vorticity for the fluid flow past [PITH_FULL_IMAGE:figures/full_fig_p007_13.png] view at source ↗
Figure 11
Figure 11. Figure 11: Persistence curves for the sampled, subdivided, and L-SIAC [PITH_FULL_IMAGE:figures/full_fig_p007_11.png] view at source ↗
Figure 14
Figure 14. Figure 14: Segmentation of the vorticity over a flow past a cylinder described [PITH_FULL_IMAGE:figures/full_fig_p007_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Segmentation of the vorticity for the flow past a cylinder de [PITH_FULL_IMAGE:figures/full_fig_p008_15.png] view at source ↗
Figure 18
Figure 18. Figure 18: Zoom in on the region associated with Figure 17. For each of [PITH_FULL_IMAGE:figures/full_fig_p008_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: Persistence diagrams for the vorticity of the counter-rotating [PITH_FULL_IMAGE:figures/full_fig_p008_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: Segmentation corresponding to the topologically simplified vor [PITH_FULL_IMAGE:figures/full_fig_p009_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: Persistence diagrams for the vorticity of the counter-rotating [PITH_FULL_IMAGE:figures/full_fig_p009_21.png] view at source ↗
read the original abstract

High-order finite element methods (HO-FEM) are gaining popularity in the simulation community due to their success in solving complex flow dynamics. There is an increasing need to analyze the data produced as output by these simulations. Simultaneously, topological analysis tools are emerging as powerful methods for investigating simulation data. However, most of the current approaches to topological analysis have had limited application to HO-FEM simulation data for two reasons. First, the current topological tools are designed for linear data (polynomial degree one), but the polynomial degree of the data output by these simulations is typically higher (routinely up to polynomial degree six). Second, the simulation data and derived quantities of the simulation data have discontinuities at element boundaries, and these discontinuities do not match the input requirements for the topological tools. One solution to both issues is to transform the high-order data to achieve low-order, continuous inputs for topological analysis. Nevertheless, there has been little work evaluating the possible transformation choices and their downstream effect on the topological analysis. We perform an empirical study to evaluate two commonly used data transformation methodologies along with the recently introduced L-SIAC filter for processing high-order simulation data. Our results show diverse behaviors are possible. We offer some guidance about how best to consider a pipeline of topological analysis of HO-FEM simulations with the currently available implementations of topological analysis.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper presents an empirical study evaluating the effects of two commonly used data transformation methodologies and the L-SIAC filter when preparing high-order FEM simulation outputs (typically polynomial degree up to 6, with element-boundary discontinuities) for input to topological analysis tools designed for linear continuous data. The central claim is that these transformations produce visibly diverse behaviors in the resulting topological outputs, and the authors provide guidance on pipeline construction using currently available tools.

Significance. If the reported differences are reproducible and not artifacts of the specific test cases, the work is significant for the growing community using topological analysis on HO-FEM data. It supplies concrete evidence that preprocessing choices are not neutral and must be documented, which is a practical contribution even if the study does not claim the three methods exhaust the design space.

major comments (2)
  1. [Results] Results section: the claim of 'diverse behaviors' rests on visual comparison of topological features; without quantitative metrics (e.g., persistence diagram distances, Betti number differences, or statistical tests across multiple datasets), it is difficult to judge whether the observed differences are systematic or case-specific.
  2. [Methods] Methods: the description of the two 'commonly used' transformation methodologies lacks sufficient implementation detail (exact interpolation orders, handling of quadrature points, and boundary treatment) to permit independent reproduction or extension of the comparison.
minor comments (2)
  1. [Introduction] Abstract and introduction: the statement that 'there has been little work evaluating the possible transformation choices' would benefit from one or two additional citations to related preprocessing studies in visualization or CFD.
  2. [Figures] Figure captions: several figures showing topological outputs would be clearer if they explicitly labeled the transformation method, polynomial degree, and mesh resolution used for each panel.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We appreciate the referee's feedback on our manuscript. The comments highlight areas where we can enhance the presentation of our empirical study. We respond to each major comment below.

read point-by-point responses
  1. Referee: [Results] Results section: the claim of 'diverse behaviors' rests on visual comparison of topological features; without quantitative metrics (e.g., persistence diagram distances, Betti number differences, or statistical tests across multiple datasets), it is difficult to judge whether the observed differences are systematic or case-specific.

    Authors: We agree that quantitative metrics would strengthen the evaluation. In the revised manuscript we will add Betti number differences computed on the presented test cases to provide a quantitative complement to the visual comparisons. revision: yes

  2. Referee: [Methods] Methods: the description of the two 'commonly used' transformation methodologies lacks sufficient implementation detail (exact interpolation orders, handling of quadrature points, and boundary treatment) to permit independent reproduction or extension of the comparison.

    Authors: We thank the referee for noting this gap. The revised version will expand the Methods section with the requested implementation details on interpolation orders, quadrature point handling, and boundary treatment for both transformation methodologies. revision: yes

Circularity Check

0 steps flagged

No significant circularity: purely empirical comparison

full rationale

The paper conducts an empirical study comparing the effects of two common data transformation methodologies and the L-SIAC filter on topological analysis outputs for high-order FEM scalar fields. No mathematical derivations, parameter fitting, predictions from first principles, or uniqueness theorems are present. The central claim rests on direct observation of differing topological features (e.g., critical points, persistence diagrams) across the tested transformations on provided datasets. Self-citations, if any, are incidental and not load-bearing for the reported behaviors. The work is self-contained as an observational benchmark against external tools and data.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

This is an empirical comparison study; the abstract identifies no free parameters, mathematical axioms, or new postulated entities.

pith-pipeline@v0.9.0 · 5777 in / 1010 out tokens · 28034 ms · 2026-05-24T20:23:41.341252+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

73 extracted references · 73 canonical work pages

  1. [1]

    Acharya and V

    A. Acharya and V . Natarajan. A parallel and memory efficient algorithm for constructing the contour tree. In 2015 IEEE Pacific Visualization Symposium (PacificVis), pp. 271–278. IEEE, 2015. 9

  2. [2]

    Ahrens, B

    J. Ahrens, B. Geveci, and C. Law. Paraview: An end-user tool for large- data visualization. The Visualization Handbook, pp. 717–731, 2005

  3. [3]

    Al Sulti and K

    F. Al Sulti and K. Ohkitani. V ortex merger and topological changes in two-dimensional turbulence. Physical Review E, 86(1):016309, 2012

  4. [4]

    T. F. Banchoff. Critical points and curvature for embedded polyhedral surfaces. The American Mathematical Monthly, 1970

  5. [5]

    Biasotti, D

    S. Biasotti, D. Giorgio, M. Spagnuolo, and B. Falcidieno. Reeb graphs for shape analysis and applications. Theoretical Computer Science, 2008

  6. [6]

    J. H. Bramble and A. H. Schatz. Higher order local accuracy by averaging in the finite element method. Mathematics of Computation, 31(137):94– 111, 1977

  7. [7]

    Bremer, G

    P.-T. Bremer, G. Weber, V . Pascucci, M. Day, and J. Bell. Analyzing and tracking burning structures in lean premixed hydrogen flames. IEEE Transactions on Visualization and Computer Graphics, 16(2):248–260, 2010

  8. [8]

    C. D. Cantwell, D. Moxey, A. Comerford, A. Bolis, G. Rocco, G. Men- galdo, D. De Grazia, S. Yakovlev, J.-E. Lombard, D. Ekelschot, B. Jordi, H. Xu, Y . Mohamied, C. Eskilsson, B. Nelson, P. V os, C. Biotto, R. M. Kirby, and S. J. Sherwin. Nektar++ : An open-source spectral/hp element framework. Computer Physics Communications, 192:205–219, Jul 2015. doi: ...

  9. [9]

    Carr and J

    H. Carr and J. Snoeyink. Representing interpolant topology for contour tree computation. In Topology-Based Methods in Visualization II , pp. 59–73. Springer, 2009

  10. [10]

    H. Carr, J. Snoeyink, and U. Axen. Computing contour trees in all dimen- sions. Computational Geometry, 24(2):75–94, 2003

  11. [11]

    H. Carr, J. Snoeyink, and M. van de Panne. Flexible isosurfaces: Simplify- ing and displaying scalar topology using the contour tree. Computational Geometry, 43(1):42 – 58, 2010

  12. [12]

    Cazals, F

    F. Cazals, F. Chazal, and T. Lewiner. Molecular shape analysis based upon the Morse-Smale complex and the Connolly function. In Symp. on Comp. Geom., 2003

  13. [13]

    G. Chen, K. Mischaikow, R. S. Laramee, and E. Zhang. Efficient morse decompositions of vector fields. IEEE Trans. Visualization and Computer Graphics, 14(4), 2008

  14. [14]

    Cockburn, M

    B. Cockburn, M. Luskin, C.-W. Shu, and E. S ¨uli. Post-processing of Galerkin methods for hyperbolic problems. In Discontinuous Galerkin Methods, pp. 291–300. Springer, 2000

  15. [15]

    Cockburn, M

    B. Cockburn, M. Luskin, C.-W. Shu, and E. S ´uli. Enhanced accuracy by post-processing for finite element methods for hyperbolic equations. Mathematics of Computation, 72(242):577–606, 2003

  16. [16]

    Cohen-Steiner, H

    D. Cohen-Steiner, H. Edelsbrunner, and J. Harer. Stability of persistence diagrams. In Symp. on Comp. Geom., 2005

  17. [17]

    Cole-McLaughlin, H

    K. Cole-McLaughlin, H. Edelsbrunner, J. Harer, V . Natarajan, and V . Pas- cucci. Loops in Reeb graphs of 2-manifolds. In Symp. on Comp. Geom., pp. 344–350, 2003

  18. [18]

    M. Day, J. Bell, P.-T. Bremer, V . Pascucci, V . Beckner, and M. Lijew- ski. Turbulence effects on cellular burning structures in lean premixed hydrogen flames. Combustion and Flame, 156:1035–1045, 2009

  19. [19]

    De Floriani, U

    L. De Floriani, U. Fugacci, F. Iuricich, and P. Magillo. Morse complexes for shape segmentation and homological analysis: discrete models and algorithms. Comp. Grap. For., 2015

  20. [20]

    S. E. Dillard, V . Natarajan, G. H. Weber, V . Pascucci, and B. Hamann. Topology-guided tessellation of quadratic elements. International Journal of Computational Geometry & Applications, 19(02):195–211, 2009

  21. [21]

    Docampo-S´anchez, J

    J. Docampo-S´anchez, J. K. Ryan, M. Mirzargar, and R. M. Kirby. Multi- dimensional filtering: Reducing the dimension through rotation. SIAM Journal on Scientific Computing, 39(5):A2179–A2200, 2017

  22. [22]

    Doraiswamy and V

    H. Doraiswamy and V . Natarajan. Output-sensitive construction of reeb graphs. IEEE Trans. on Vis. and Comp. Graph., 2012

  23. [23]

    Doraiswamy and V

    H. Doraiswamy and V . Natarajan. Computing reeb graphs as a union of contour trees. IEEE Trans. on Vis. and Comp. Graph., 2013

  24. [24]

    Edelsbrunner and J

    H. Edelsbrunner and J. Harer. Computational Topology: An Introduction. American Mathematical Society, 2009

  25. [25]

    Edelsbrunner, D

    H. Edelsbrunner, D. Letscher, and A. Zomorodian. Topological persistence and simplification. Disc. Compu. Geom., 2002

  26. [26]

    Etiene, L

    T. Etiene, L. Nonato, C. Scheidegger, J. Tienry, T. Peters, V . Pascucci, R. Kirby, and C. Silva. Topology verification for isosurface extraction. IEEE Transactions on Visualization and Computer Graphics, 18(6):952– 965, June 2012. doi: 10.1109/TVCG.2011.109

  27. [27]

    N. Fehn, W. A. Wall, and M. Kronbichler. Robust and efficient discontinu- ous Galerkin methods for under-resolved turbulent incompressible flows. Journal of Computational Physics, 372:667–693, 2018

  28. [28]

    R. Forman. A user’s guide to discrete Morse theory. Adv. in Math., 1998

  29. [29]

    Guenther, R

    D. Guenther, R. Alvarez-Boto, J. Contreras-Garcia, J.-P. Piquemal, and J. Tierny. Characterizing molecular interactions in chemical systems.IEEE Trans. on Vis. and Comp. Graph., 2014

  30. [30]

    Gueunet, P

    C. Gueunet, P. Fortin, J. Jomier, and J. Tierny. Contour forests: Fast multi-threaded augmented contour trees. In IEEE LDAV, 2016

  31. [31]

    Gyulassy, P

    A. Gyulassy, P. Bremer, R. Grout, H. Kolla, J. Chen, and V . Pascucci. Stability of dissipation elements: A case study in combustion. Comp. Graph. For., 2014

  32. [32]

    Gyulassy, P

    A. Gyulassy, P. T. Bremer, B. Hamann, and V . Pascucci. A practical approach to morse-smale complex computation: Scalability and generality. IEEE Trans. on Vis. and Comp. Graph., 2008

  33. [33]

    Gyulassy, A

    A. Gyulassy, A. Knoll, K. Lau, B. Wang, P. Bremer, M. Papka, L. A. Curtiss, and V . Pascucci. Interstitial and interlayer ion diffusion geometry extraction in graphitic nanosphere battery materials. IEEE Trans. on Vis. and Comp. Graph., 2015

  34. [34]

    Gyulassy, V

    A. Gyulassy, V . Natarajan, M. Duchaineau, V . Pascucci, E. Bringa, A. Hig- ginbotham, and B. Hamann. Topologically Clean Distance Fields. IEEE Trans. on Vis. and Comp. Graph., 2007

  35. [35]

    Heine, H

    C. Heine, H. Leitte, M. Hlawitschka, F. Iuricich, L. De Floriani, G. Scheuermann, H. Hagen, and C. Garth. A survey of topology-based methods in visualization. Comp. Grap. For., 2016

  36. [36]

    Z. W. H.T. Huynh and P. Vincent. High-order methods for computational fluid dynamics: A brief review of compact differential formulations on unstructured grids. Computers & Fluids, 98:209–220, 2014

  37. [37]

    Jallepalli, J

    A. Jallepalli, J. Docampo-S ´anchez, J. K. Ryan, R. Haimes, and R. M. Kirby. On the treatment of field quantities and elemental continuity in FEM solutions. IEEE Transactions on Visualization and Computer Graphics, 24(1):903–912, 2017

  38. [38]

    Jallepalli, R

    A. Jallepalli, R. Haimes, and R. M. Kirby. Adaptive characteristic length for L-SIAC filtering of fem data. Journal of Scientific Computing , pp. 1–22, 2018

  39. [39]

    Jallepalli and R

    A. Jallepalli and R. M. Kirby. Efficient algorithms for the line-siac filter. Journal of Scientific Computing, 2019

  40. [40]

    G. E. Karniadakis and S. J. Sherwin. Spectral/hp element methods for Computational Fluid Dynamics (Second Edition). Oxford University Press, 2005

  41. [41]

    Kasten, J

    J. Kasten, J. Reininghaus, I. Hotz, and H. Hege. Two-dimensional time- dependent vortex regions based on the acceleration magnitude. IEEE Trans. on Vis. and Comp. Graph., 2011

  42. [42]

    D. E. Laney, P. Bremer, A. Mascarenhas, P. Miller, and V . Pascucci. Understanding the structure of the turbulent mixing layer in hydrodynamic instabilities. IEEE Trans. on Vis. and Comp. Graph., 2006

  43. [43]

    X. Li, J. K. Ryan, R. M. Kirby, and C. Vuik. Smoothness-Increasing Accuracy-Conserving (SIAC) filters for derivative approximations of dis- continuous Galerkin (DG) solutions over nonuniform meshes and near boundaries. Journal of Computational and Applied Mathematics, 294:275– 296, 2016

  44. [44]

    J.-E. W. Lombard, D. Moxey, S. J. Sherwin, J. F. A. Hoessler, S. Dhanda- pani, and M. J. Taylor. Implicit large-eddy simulation of a wingtip vortex. AIAA Journal, 54(2):506–518, 2016

  45. [45]

    Mirzaee, L

    H. Mirzaee, L. Ji, J. K. Ryan, and R. M. Kirby. Smoothness-increasing accuracy-conserving (SIAC) postprocessing for discontinuous Galerkin solutions over structured triangular meshes. SIAM Journal on Numerical Analysis, 49(5):1899–1920, 2011

  46. [46]

    Mirzaee, J

    H. Mirzaee, J. King, J. K. Ryan, and R. M. Kirby. Smoothness-increasing accuracy-conserving filters for discontinuous Galerkin solutions over un- structured triangular meshes. SIAM Journal on Scientific Computing , 35(1):A212–A230, 2013

  47. [47]

    Mirzaee, J

    H. Mirzaee, J. K. Ryan, and R. M. Kirby. Efficient implementation of smoothness-increasing accuracy-conserving (SIAC) filters for discontinu- ous Galerkin solutions. Journal of Scientific Computing, 52(1):85–112, 2012

  48. [48]

    Mirzaee, J

    H. Mirzaee, J. K. Ryan, and R. M. Kirby. Smoothness-increasing accuracy- conserving (SIAC) filters for discontinuous Galerkin solutions: Applica- tion to structured tetrahedral meshes. Journal of Scientific Computing , 58(3):690–704, 2014

  49. [49]

    Mirzargar, A

    M. Mirzargar, A. Jallepalli, J. K. Ryan, and R. M. Kirby. Hexagonal smoothness-increasing accuracy-conserving filtering. Journal of Scientific Computing, 73(2-3):1072–1093, 2017

  50. [50]

    Moxey, R

    D. Moxey, R. Amici, and R. M. Kirby. Efficient matrix-free high-order fi- nite element evaluation for simplicial elements.SIAM Journal on Scientific Computing, Under Review, 2019. 10 To appear in IEEE Transactions on Visualization and Computer Graphics

  51. [51]

    Nguyen and J

    D.-M. Nguyen and J. Peters. Nonuniform discontinuous Galerkin filters via shift and scale. SIAM Journal on Numerical Analysis, 54(3):1401–1422, 2016

  52. [52]

    Norgard and P.-T

    G. Norgard and P.-T. Bremer. Robust computation of morse–smale complexes of bilinear functions. Computer Aided Geometric Design , 30(6):577–587, 2013

  53. [53]

    Nucha, G.-P

    G. Nucha, G.-P. Bonneau, S. Hahmann, and V . Natarajan. Computing con- tour trees for 2d piecewise polynomial functions. In Computer Graphics Forum, vol. 36, pp. 23–33. Wiley Online Library, 2017

  54. [54]

    S. Parsa. A deterministic o(m log m) time algorithm for the reeb graph. In Symp. on Comp. Geom., 2012

  55. [55]

    Pascucci and K

    V . Pascucci and K. Cole-McLaughlin. Parallel computation of the topology of level sets. Algorithmica, 38(1):249–268, 2004

  56. [56]

    Pascucci, G

    V . Pascucci, G. Scorzelli, P. T. Bremer, and A. Mascarenhas. Robust on-line computation of Reeb graphs: simplicity and speed. ACM Trans. on Graph., 2007

  57. [57]

    Patane, M

    G. Patane, M. Spagnuolo, and B. Falcidieno. Reeb graph computation based on a minimal contouring. In Proc. of IEEE Shape Modeling Inter- national, 2008

  58. [58]

    J. K. Ryan, X. Li, R. M. Kirby, and K. Vuik. One-sided position-dependent smoothness-increasing accuracy-conserving (SIAC) filtering over uniform and non-uniform meshes. Journal of Scientific Computing, 64(3):773–817, 2015

  59. [59]

    Shivashankar, P

    N. Shivashankar, P. Pranav, V . Natarajan, R. van de Weygaert, E. P. Bos, and S. Rieder. Felix: A topology based framework for visual exploration of cosmic filaments. IEEE Trans. on Vis. and Comp. Graph., 2016. http: //vgl.serc.iisc.ernet.in/felix/index.html

  60. [60]

    A. Sidot. Spectral/hp Element Methods Applied To Vortex Structures In A Low Incidence Delta Wing Wake Compared To Experiments. M.S. Thesis in Aeronautics, Imperial College London, September 2017

  61. [61]

    B. S. Sohn and C. L. Bajaj. Time varying contour topology. IEEE Trans. on Vis. and Comp. Graph., 2006

  62. [62]

    T. Sousbie. The persistent cosmic web and its filamentary structure: Theory and implementations. Royal Astronomical Society, 2011. http: //www2.iap.fr/users/sousbie/web/html/indexd41d.html

  63. [63]

    Steffen, S

    M. Steffen, S. Curtis, R. M. Kirby, and J. K. Ryan. Investigation of smoothness-increasing accuracy-conserving filters for improving stream- line integration through discontinuous fields. IEEE Transactions on Visu- alization and Computer Graphics, 14(3):680–692, 2008

  64. [64]

    Tarasov and M

    S. Tarasov and M. Vyali. Construction of contour trees in 3d in o(n log n) steps. In Symp. on Comp. Geom., 1998

  65. [65]

    D. M. Thomas and V . Natarajan. Multiscale symmetry detection in scalar fields by clustering contours. IEEE Trans. on Vis. and Comp. Graph. , 2014

  66. [66]

    Tierny, G

    J. Tierny, G. Favelier, J. A. Levine, C. Gueunet, and M. Michaux. The topology toolkit. IEEE Trans. Vis. Comput. Graph., 24(1):832–842, 2018

  67. [67]

    Tierny, A

    J. Tierny, A. Gyulassy, E. Simon, and V . Pascucci. Loop surgery for volumetric meshes: Reeb graphs reduced to contour trees. IEEE Trans. on Vis. and Comp. Graph., 2009

  68. [68]

    Tierny and V

    J. Tierny and V . Pascucci. Generalized topological simplification of scalar fields on surfaces. IEEE Trans. on Vis. and Comp. Graph., 2012

  69. [69]

    Tricoche, C

    X. Tricoche, C. Garth, G. Kindlmann, E. Deines, G. Scheuermann, M. Ruetten, and C. Hansen. Visualization of intricate flow structures for vortex breakdown analysis. In IEEE Visualization 2004, pp. 187–194. IEEE, 2004

  70. [70]

    van Kreveld, R

    M. van Kreveld, R. van Oostrum, C. Bajaj, V . Pasucci, and D. Schikore. Contour trees and small seed sets for isosurface traversal. In Symp. on Comp. Geom., 1997

  71. [71]

    van Slingerland, J

    P. van Slingerland, J. K. Ryan, and C. Vuik. Position-dependent smoothness-increasing accuracy-conserving (SIAC) filtering for improv- ing discontinuous Galerkin solutions. SIAM Journal on Scientific Comput- ing, 33(2):802–825, 2011

  72. [72]

    Walfisch, J

    D. Walfisch, J. K. Ryan, R. M. Kirby, and R. Haimes. One-sided smoothness-increasing accuracy-conserving filtering for enhanced stream- line integration through discontinuous fields.Journal of Scientific Com- puting, 38(2):164–184, 2009

  73. [73]

    Weber, S

    G. Weber, S. E. Dillard, H. Carr, V . Pascucci, and B. Hamann. Topology- controlled volume rendering. IEEE Trans. on Vis. and Comp. Graph. , 2007. 11