Max Kontak

Max Kontak, Jules Vidal, Julien Tierny

In urgent decision making applications, ensemble simulations are an important way to determine different outcome scenarios based on currently available data. In this paper, we will analyze the output of ensemble simulations by considering so-called persistence diagrams, which are reduced representations of the original data, motivated by the extraction of topological features. Based on a recently published progressive algorithm for the clustering of persistence diagrams, we determine the optimal number of clusters, and therefore the number of significantly different outcome scenarios, by the minimization of established statistical score functions. Furthermore, we present a proof-of-concept prototype implementation of the statistical selection of the number of clusters and provide the results of an experimental study, where this implementation has been applied to real-world ensemble data sets.

Simone Gramsch, Max Kontak, Volker Michel

An elementary algorithm is used to simulate the industrial production of a fiber of a 3-dimensional nonwoven fabric. The algorithm simulates the fiber as a polyline where the direction of each segment is stochastically drawn based on a given probability density function (PDF) on the unit sphere. This PDF is obtained from data of directions of fiber fragments which originate from computer tomography scans of a real non-woven fabric. However, the simulation algorithm requires numerous evaluations of the PDF. Since the established technique of a kernel density estimator leads to very high computational costs, a novel greedy algorithm for estimating a sparse representation of the PDF is introduced. Numerical tests for a synthetic and a real example are presented. In a realistic scenario, the introduced sparsity ansatz leads to a reduction of the computation time for 100 fibers from nearly 40 days to 41 minutes.