Abstract. This paper presents a computational framework that allows for a robust extraction of the extremal structure of scalar and vector fields on 2D manifolds embedded in 3D. This structure consists of critical points, separatrices, and periodic orbits. The framework is based on Forman's discrete Morse theory, which guarantees the topological consistency of the computed extremal structure. Using a graph theoretical formulation of this theory, we present an algorithmic pipeline that computes a hierarchy of extremal structures. This hierarchy is defined by an importance measure and enables the user to select an appropriate level of detail.
Keywords: Discrete Morse theory, data analysis, scalar fields, vector fields
MotivationWe propose a computational framework to extract the extremal structure of scalar and vector fields on 2D manifolds embedded in R 3 . The extremal structure of a scalar field consists of critical points and separatrices -the streamlines of the gradient field that connect the critical points. The extremal structure of a vector field additionally includes periodic orbits -the streamlines that are closed.These structures are of great interest in many applications and have a long history [2,12]. Typically, the critical points are computed by finding all zeros of the gradient or vector field. The critical points of a scalar field are classified into minima, saddles, and maxima by the eigenvalues of its Hessian, while the critical points of a vector field are classified into sinks, saddles, and sources by the eigenvalues of its Jacobian. The respective eigenvectors can be used to compute the separatrices as the solution of an autonomous ODE. For the numerical treatment of these problems and the extraction of the periodic orbits, we refer to [18,20,5].One of the biggest challenges that such numerical algorithms face is the discrete nature of the extremal structure which necessitates a lot of binary decisions. For example, the type of a critical point depends on the sign of the eigenvalues.