Hochdimensionale effiziente Approximationsalgorithmen basierend auf Abtastmengen mit Rang-1 Strukturen und Anwendungen
Zusammenfassung der Projektergebnisse
Im Rahmen des Projektes gelang es, ein hocheffizientes Verfahren zur approximativen Bestimmung von Lösungen partieller Differentialgleichungen mit zufälligen Koeffizienten zu erstellen. Größere Anzahlen zufälliger Parameter stellen dabei besondere Herausforderungen dar, so dass nur hocheffiziente Verfahren zum Einsatz kommen können. Die numerisch ermittelten Lösungen dienen dann der vertiefenden Charakterisierung von Parameterabhängigkeiten, welche im Anwendungsfall oft bereits durch theoretische Überlegungen grundlegend charakterisiert sind. Im konkreten Approximationsproblem treten vor allem auch durch die in hochdimensionalen Problemstellungen hohe Anzahl verschiedenartiger Kopplungen der Parameter jedoch Feinheiten auf, die mittels theoretischer Ergebnisse meist nur unzureichend beschrieben werden und damit zu mangelhaften Ergebnissen nichtadaptiver Verfahren führen können. Der von uns entwickelte Algorithmus ist in der Lage genau diese Feinheiten in den ermittelten Approximationen abzubilden und bereits während des Approximationsprozesses z.B. in Form reduzierter Abtastungen davon zu profitieren. Um die Effizienz der eingesetzten dimensionsinkrementellen Methode deutlich zu steigern, wurden effizientere Abtaststrategien zur Rekonstruktion von und Approximation mit trigonometrischen Polynomen entwickelt. Besonderer Augenmerk lag hierbei auf dem allgegenwärtigen Ziel der Effizienz bzgl. Konstruktion der Abtastschemata, Anzahl benötigter Abtastknoten und Rekonstruktionaufwand. Als Folge lassen sich sehr komplexe Probleme, insb. partielle Differentialgleichungen mit zufälligen Koeffizienten, deutlich ausführlicher numerisch analysieren als das bislang möglich war. Verschiedene Implementierungen der entwickelten Algorithmen wurden in unterschiedlichen Software-Paketen frei zugänglich veröffentlicht. Ein umfassenderes Software-Paket zur dimensionsinkrementellen Methode ist aktuell in Arbeit. Dies schließt eine universelle Schnittstelle zur Einbindung weiterer Rekonstruktionsverfahren, wie z.B. den im Projekt entwickelten Verfahren zur Rekonstruktion mittels multivariater Chebyshev-Polynome, ein. Diese Verfahren sind tatsächlich eine Weiterentwicklung der für multivariate trigonometrische Polynome entwickelten Verfahren. Dabei nutzen wir Chebyshev-Polynome als Basisfunktionen für nichtperiodische Problemstellungen. Entsprechende Algorithmen zur Bestimmung geeigneter Abtastschemata und zu deren Anwendung in dimensionsinkrementellen Approximationsalgorithmen stehen mit Abschluss des Projektes zur Verfügung. Diese wurden auch bereits ausgiebig erfolgreich getestet. Die bloße Anwendung der Algorithmen auf einige Differentialgleichungen mit zufälligen Koeffizienten stellt natürlich nur das Ende des vorliegenden Projektes dar. Selbstverständlich werfen die entwickelten, innovativen Algorithmen neue Fragestellungen auf. Auf der einen Seite geben numerische Tests Anlass, die dimensionsinkrementelle Strategie selbst weiter zu entwicklen, wie z.B. weitere Verbesserungen des verwendeten Thresholdings vorzunehmen. Auf der anderen Seite ist die Anwendbarkeit auf komplexere Differentialgleichungen mit zufülligen Koeffizienten nun ein zentrales Thema. So stellt sich die Frage nach den umfassenden Möglichkeiten aber auch nach den Limitierungen der existierenden Algorithmen. Darüber hinaus benötigt eine adäquate Interpretation vorliegender Approximanten weiterer Behandlung. Wir nutzten bereits den Ansatz der Varianzanalyse (ANOVA), um Abhängigkeiten von verschiedenen Variablenkombinationen zu quantifizieren. Mit einem solchen Ansatz wurde in einer von mir betreuten Dissertation für Approximationen aus gestreuten Daten (scattered data approximation) gearbeitet und auch sehr erfolgreich auf verschiedene Problemstellungen aus der Praxis angewandt. Die Überschneidungen zur Approximation unter Nutzung des dimensionsinkrementellen Algorithmus sind offensichtlich, die strukturierte Anwendung der ANOVA zur Approximationsverbesserung muss Thema zukünftiger Forschung sein. Zu den Forschungsergebnissen sind mehrere Publikationen in internationalen Zeitschriften erschienen.
Projektbezogene Publikationen (Auswahl)
-
Approximation of multivariate periodic functions based on sampling along multiple rank-1 lattices. Journal of Approximation Theory, 246, 1-27.
Kämmerer, Lutz & Volkmer, Toni
-
Constructing spatial discretizations for sparse multivariate trigonometric polynomials that allow for a fast discrete Fourier transform. Applied and Computational Harmonic Analysis, 47(3), 702-729.
Kämmerer, Lutz
-
A sparse FFT approach for ODE with random coefficients. Advances in Computational Mathematics, 46(5).
Bochmann, Maximilian; Kämmerer, Lutz & Potts, Daniel
-
A deterministic algorithm for constructing multiple rank-1 lattices of near-optimal size. Advances in Computational Mathematics, 47(6).
Gross, Craig; Iwen, Mark A.; Kämmerer, Lutz & Volkmer, Toni
-
Sparse Fourier transforms on rank-1 lattices for the rapid and low-memory approximation of functions of many variables. Sampling Theory, Signal Processing, and Data Analysis, 20(1).
Gross, Craig; Iwen, Mark; Kämmerer, Lutz & Volkmer, Toni
-
Worst-case Recovery Guarantees for Least Squares Approximation Using Random Samples. Constructive Approximation, 54(2), 295-352.
Kämmerer, Lutz; Ullrich, Tino & Volkmer, Toni
-
The uniform sparse FFT with application to PDEs with random coefficients. Sampling Theory, Signal Processing, and Data Analysis, 20(2).
Kämmerer, Lutz; Potts, Daniel & Taubert, Fabian
-
Nonlinear approximation in bounded orthonormal product bases. Sampling Theory, Signal Processing, and Data Analysis, 21(1).
Kämmerer, Lutz; Potts, Daniel & Taubert, Fabian
-
On the reconstruction of functions from values at subsampled quadrature points. Mathematics of Computation, 93(346), 785-809.
Bartel, Felix; Kämmerer, Lutz; Potts, Daniel & Ullrich, Tino
-
Constructing efficient spatial discretizations of spans of multivariate Chebyshev polynomials
L. Kämmerer
-
Exact discretization, tight frames and recovery via D-optimal designs. The SMAI Journal of computational mathematics, 11, 607-636.
Bartel, Felix; Kämmerer, Lutz; Pozharska, Kateryna; Schäfer, Martin & Ullrich, Tino
