Investigating In Situ Reduction via Lagrangian Representations for Cosmology and Seismology Applications

Although many types of computational simulations produce timevarying vector fields, subsequent analysis is often limited to single time slices due to excessive costs. Fortunately, a new approach using a Lagrangian representation can enable time-varying vector field analysis while mitigating these costs. With this approach, a Lagrangian representation is calculated while the simulation code is running, and the result is explored after the simulation. Importantly, the effectiveness of this approach varies based on the nature of the vector field, requiring in-depth investigation for each application area. With this study, we evaluate the effectiveness for previously unexplored cosmology and seismology applications. We do this by considering encumbrance (on the simulation) and accuracy (of the reconstructed result). To inform encumbrance, we integrated in situ infrastructure with two simulation codes, and evaluated on representative HPC environments, performing Lagrangian in situ reduction using GPUs as well as CPUs. To inform accuracy, our study conducted a statistical analysis across a range of spatiotemporal configurations as well as a qualitative evaluation. In all, we demonstrate effectiveness for both cosmology and seismology — time-varying vector fields from these domains can be reduced to less than 1% of the total data via Lagrangian representations, while maintaining accurate reconstruction and requiring under 10% of total execution time in over 80% of our experiments.

[1]  Alexander Sen Gupta,et al.  Quantification of errors induced by temporal resolution on Lagrangian particles in an eddy-resolving model , 2014 .

[2]  Kenneth I. Joy,et al.  Improved post hoc flow analysis via Lagrangian representations , 2014, 2014 IEEE 4th Symposium on Large Data Analysis and Visualization (LDAV).

[3]  David Pugmire,et al.  Performance-Portable Particle Advection with VTK-m , 2018, EGPGV@EuroVis.

[4]  Hank Childs,et al.  VisIt: An End-User Tool for Visualizing and Analyzing Very Large Data , 2011 .

[5]  Daniel Weiskopf,et al.  Hierarchical Line Integration , 2011, IEEE Transactions on Visualization and Computer Graphics.

[6]  K. Döös,et al.  Calculating Lagrangian Trajectories Using Time-Dependent Velocity Fields , 2001 .

[7]  Kenneth I. Joy,et al.  Interpolation-Based Pathline Tracing in Particle-Based Flow Visualization , 2015, IEEE Transactions on Visualization and Computer Graphics.

[8]  Kenneth I. Joy,et al.  Analysis of Error in Interpolation-Based Pathline Tracing , 2016, EuroVis.

[9]  Roxana Bujack,et al.  Revisiting the Evaluation of In Situ Lagrangian Analysis , 2018, EGPGV@EuroVis.

[10]  Roxana Bujack,et al.  An Interpolation Scheme for VDVP Lagrangian Basis Flows , 2019, EGPGV@EuroVis.

[11]  Martin Isenburg,et al.  Fast and Efficient Compression of Floating-Point Data , 2006, IEEE Transactions on Visualization and Computer Graphics.

[12]  Jean M. Sexton,et al.  Nyx: A MASSIVELY PARALLEL AMR CODE FOR COMPUTATIONAL COSMOLOGY , 2013, J. Open Source Softw..

[13]  Leigh Orf A Violently Tornadic Supercell Thunderstorm Simulation Spanning a Quarter-Trillion Grid Volumes: Computational Challenges, I/O Framework, and Visualizations of Tornadogenesis , 2019, Atmosphere.

[14]  Kenneth I. Joy,et al.  Subsampling-based compression and flow visualization , 2015, Electronic Imaging.

[15]  B. Blanke,et al.  Lagrangian methods for flow climatologies and trajectory error assessment , 2004 .

[16]  Kenneth I. Joy,et al.  Error Estimates for Lagrangian Flow Field Representations , 2016, EuroVis.

[17]  Markus Gross,et al.  A Fluid Flow Data Set for Machine Learning and its Application to Neural Flow Map Interpolation , 2020, IEEE Transactions on Visualization and Computer Graphics.

[18]  Christoph Garth,et al.  Scalable In Situ Lagrangian Flow Map Extraction: Demonstrating the Viability of a Communication-Free Model , 2020, EGPGV@EuroVis.

[19]  Suresh K. Lodha,et al.  Topology Preserving Top-Down Compression of 2D Vector Fields Using Bintree and Triangular Quadtrees , 2003, IEEE Trans. Vis. Comput. Graph..

[20]  Kenneth I. Joy,et al.  Lagrangian representations of flow fields with parameter curves , 2015, 2015 IEEE 5th Symposium on Large Data Analysis and Visualization (LDAV).

[21]  Yu Cheng,et al.  Lagrangian ocean analysis: Fundamentals and practices , 2018 .

[22]  Xin Tong,et al.  Salient time steps selection from large scale time-varying data sets with dynamic time warping , 2012, IEEE Symposium on Large Data Analysis and Visualization (LDAV).

[23]  Kwan-Liu Ma,et al.  VTK-m: Accelerating the Visualization Toolkit for Massively Threaded Architectures , 2016, IEEE Computer Graphics and Applications.

[24]  Christian Rössl,et al.  Combining topological simplification and topology preserving compression for 2D vector fields , 2003, 11th Pacific Conference onComputer Graphics and Applications, 2003. Proceedings..

[25]  M. Böttinger,et al.  The tropical-subtropical coupling in the Southeast Atlantic from the perspective of the northern Benguela upwelling system , 2019, PloS one.

[26]  Christoph Peters,et al.  Void-and-Cluster Sampling of Large Scattered Data and Trajectories , 2019, IEEE Transactions on Visualization and Computer Graphics.

[27]  N. Anders Petersson,et al.  Wave propagation in anisotropic elastic materials and curvilinear coordinates using a summation-by-parts finite difference method , 2015, J. Comput. Phys..