TY - THES A1 - Gathungu, Duncan Kioi T1 - On Multigrid and H-Matrix Methods for Partial Integro-Differential Equations T1 - Mehrgitter-Verfahren und hierarchische Matrix-Lösungsverfahren für Klassen von partiellen Integro-Differential-Problemen N2 - The main theme of this thesis is the development of multigrid and hierarchical matrix solution procedures with almost linear computational complexity for classes of partial integro-differential problems. An elliptic partial integro-differential equation, a convection-diffusion partial integro-differential equation and a convection-diffusion partial integro-differential optimality system are investigated. In the first part of this work, an efficient multigrid finite-differences scheme for solving an elliptic Fredholm partial integro-differential equation (PIDE) is discussed. This scheme combines a second-order accurate finite difference discretization and a Simpson's quadrature rule to approximate the PIDE problem and a multigrid scheme and a fast multilevel integration method of the Fredholm operator allowing the fast solution of the PIDE problem. Theoretical estimates of second-order accuracy and results of local Fourier analysis of convergence of the proposed multigrid scheme are presented. Results of numerical experiments validate these estimates and demonstrate optimal computational complexity of the proposed framework that includes numerical experiments for elliptic PIDE problems with singular kernels. The experience gained in this part of the work is used for the investigation of convection diffusion partial-integro differential equations in the second part of this thesis. Convection-diffusion PIDE problems are discretized using a finite volume scheme referred to as the Chang and Cooper (CC) scheme and a quadrature rule. Also for this class of PIDE problems and this numerical setting, a stability and accuracy analysis of the CC scheme combined with a Simpson's quadrature rule is presented proving second-order accuracy of the numerical solution. To extend and investigate the proposed approximation and solution strategy to the case of systems of convection-diffusion PIDE, an optimal control problem governed by this model is considered. In this case the research focus is the CC-Simpson's discretization of the optimality system and its solution by the proposed multigrid strategy. Second-order accuracy of the optimization solution is proved and results of local Fourier analysis are presented that provide sharp convergence estimates of the optimal computational complexity of the multigrid-fast integration technique. While (geometric) multigrid techniques require ad-hoc implementation depending on the structure of the PIDE problem and on the dimensionality of the domain where the problem is considered, the hierarchical matrix framework allows a more general treatment that exploits the algebraic structure of the problem at hand. In this thesis, this framework is extended to the case of combined differential and integral problems considering the case of a convection-diffusion PIDE. In this case, the starting point is the CC discretization of the convection-diffusion operator combined with the trapezoidal quadrature rule. The hierarchical matrix approach exploits the algebraic nature of the hierarchical matrices for blockwise approximations by low-rank matrices of the sparse convection-diffusion approximation and enables data sparse representation of the fully populated matrix where all essential matrix operations are performed with at most logarithmic optimal complexity. The factorization of part of or the whole coefficient matrix is used as a preconditioner to the solution of the PIDE problem using a generalized minimum residual (GMRes) procedure as a solver. Numerical analysis estimates of the accuracy of the finite-volume and trapezoidal rule approximation are presented and combined with estimates of the hierarchical matrix approximation and with the accuracy of the GMRes iterates. Results of numerical experiments are reported that successfully validate the theoretical estimates and the optimal computational complexity of the proposed hierarchical matrix solution procedure. These results include an extension to higher dimensions and an application to the time evolution of the probability density function of a jump diffusion process. N2 - Das Hauptthema dieser Arbeit ist die Entwicklung von Mehrgitter-Verfahren und hierarchischer Matrix-Lösungsverfahren mit nahezu linearer Rechenkomplexität für Klassen von partiellen Integro-Differential-Problemen. Es werden eine elliptische partielle Integro-Differentialgleichung, eine partielle Konvektions-Diffusions-Integro-Differentialgleichung und ein partielles Konvektions-Diffusions-Integro-Differential-Optimalitätssystem untersucht. Im ersten Teil dieser Arbeit wurde ein effizientes Mehrgitter-Finite-Differenzen-Schema zur Lösung einer elliptischen Fredholm partiellen Integro-Differentialgleichungen (PIDE) diskutiert. Dieses Schema kombiniert eine exakte finite Differenzen-Diskretisierung zweiter Ordnung mit einer Quadraturregel von Simpson, um das PIDE-Problem mit einem Mehrgitter-Schema und einer schnellen Multilevel-Integrationsmethode des Fredholm-Operators zu lösen, was eine schnelle Lösung des PIDE-Probleme ermöglicht. Theoretische Abschätzungen der Genauigkeit zweiter Ordnung und Ergebnisse der lokalen Fourier-Analyse der Konvergenz des vorgeschlagenen Mehrgitter-Systems werden präsentiert. Ergebnisse von numerischen Experimenten validieren diese Schätzungen und demonstrieren die optimale rechnerische Komplexität des vorgeschlagenen Frameworks, das numerische Experimente für elliptische PIDE mit singulären Kernen beinhaltet. Die in diesem Teil der Arbeit gewonnenen Erfahrungen werden zur Untersuchung einer partielle Konvektions-Diffusions-Integro-Differentialgleichungen im zweiten Teil verwendet. Konvektions-Diffusions-PIDE-Probleme werden unter Verwendung eines Finite-Volumen-Schemas, das als das Chang- Cooper- (CC-) Schema bezeichnet wird, und einer Quadraturregel diskretisiert. Auch für diese Klasse von PIDE-Problemen und diese numerische Einstellung wird eine Stabilitäts- und Genauigkeitsanalyse des CC-Schemas in Kombination mit einer Quadraturregel von Simpson vorgestellt, die die Genauigkeit der numerischen Lösung zweiter Ordnung beweist. Um die vorgeschlagene Approximations- und Lösungsstrategie auf den Fall von Konvektions-Diffusions-PIDE-Systemen auszudehnen und zu untersuchen, wird ein Optimalsteuerungsproblem mit diesem Modell als Nebenbedingung untersucht. Der Forschungsschwerpunkt liegt dabei auf der Diskretisierung des Optimalitätssystems durch die CC-Simpson-Lösung und dessen Lösung durch die vorgeschlagene Mehrgitter-Strategie. Die Genauigkeit der optimalen Lösung zweiter Ordnung wird bewiesen und es werden Ergebnisse der lokalen Fourier-Analyse präsentiert, die scharfe Konvergenz-Schätzungen der optimalen Berechnungskomplexität der schnellen Mehrgitter Integrationstechnik liefern. Während (geometrische) Mehrgitterverfahren je nach Struktur des PIDE-Problems und der Dimensionalität des Gebietes, in dem das Problem berücksichtigt wird, eine Ad-hoc-Implementierung erfordern, ermöglicht das hierarchische Matrix-Framework eine allgemeinere Behandlung, die die algebraische Struktur des Problems nutzt. In dieser Arbeit wird dieses Verfahren auf den Fall kombinierter Differential- und Integralprobleme im Fall einer Konvektions-Diffusions-PIDE erweitert. In diesem Fall ist der Startpunkt die CC-Diskretisierung des Konvektions-Diffusions-Operators in Kombination mit der Trapez-Quadratur-Regel. Der hierarchische Matrixansatz nutzt die algebraische Natur der hierarchischen Matrizen für blockweise Approximationen durch niedrigrangige Matrizen der dünn besetzten Konvektions-Diffusionsmatrix und ermöglicht eine datenarme Darstellung der vollständig besetzten Matrix, bei der alle wesentlichen Matrixoperationen mit höchstens logarithmisch optimaler Komplexität durchgeführt werden. Die Faktorisierung eines Teils oder der gesamten Koeffizientenmatrix wird als Vorbedingung für die Lösung der PIDE-Probleme unter Verwendung eines verallgemeinerten minimalen Restwert-Verfahrens (GMRes) als Löser verwendet. Eine numerische Analyse der Abschätzungen der Genauigkeit der Finite-Volumen- und Trapezregel-Approximation werden präsentiert und kombiniert mit Abschätzungen der hierarchischen Matrix-Näherung und mit der Genauigkeit der GMRes iterationen kombiniert. Ergebnisse numerischer Experimente werden vorgestellt, die theoretischen Abschätzungen und die optimale rechnerische Komplexität der vorgeschlagenen hierarchischen Matrix Lösungsverfahren erfolgreich validieren. Diese Ergebnisse beinhalten eine Erweiterung auf höhere Dimensionen und eine Anwendung auf die zeitliche Entwicklung der Wahrscheinlichkeitsdichtefunktion des Sprungdiffusionsprozesses. KW - Mehrgitterverfahren KW - Hierarchische Matrix KW - Integrodifferentialgleichung KW - Numerical analysis KW - multigrid KW - hierarchical matrix KW - partial integro-differential equations Y1 - 2018 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-156430 ER - TY - THES A1 - Gaviraghi, Beatrice T1 - Theoretical and numerical analysis of Fokker-Planck optimal control problems for jump-diffusion processes T1 - Theoretische und numerische Analyse von Fokker-Planck Optimalsteuerungsproblemen von Sprung-Diffusions-Prozessen N2 - The topic of this thesis is the theoretical and numerical analysis of optimal control problems, whose differential constraints are given by Fokker-Planck models related to jump-diffusion processes. We tackle the issue of controlling a stochastic process by formulating a deterministic optimization problem. The key idea of our approach is to focus on the probability density function of the process, whose time evolution is modeled by the Fokker-Planck equation. Our control framework is advantageous since it allows to model the action of the control over the entire range of the process, whose statistics are characterized by the shape of its probability density function. We first investigate jump-diffusion processes, illustrating their main properties. We define stochastic initial-value problems and present results on the existence and uniqueness of their solutions. We then discuss how numerical solutions of stochastic problems are computed, focusing on the Euler-Maruyama method. We put our attention to jump-diffusion models with time- and space-dependent coefficients and jumps given by a compound Poisson process. We derive the related Fokker-Planck equations, which take the form of partial integro-differential equations. Their differential term is governed by a parabolic operator, while the nonlocal integral operator is due to the presence of the jumps. The derivation is carried out in two cases. On the one hand, we consider a process with unbounded range. On the other hand, we confine the dynamic of the sample paths to a bounded domain, and thus the behavior of the process in proximity of the boundaries has to be specified. Throughout this thesis, we set the barriers of the domain to be reflecting. The Fokker-Planck equation, endowed with initial and boundary conditions, gives rise to Fokker-Planck problems. Their solvability is discussed in suitable functional spaces. The properties of their solutions are examined, namely their regularity, positivity and probability mass conservation. Since closed-form solutions to Fokker-Planck problems are usually not available, one has to resort to numerical methods. The first main achievement of this thesis is the definition and analysis of conservative and positive-preserving numerical methods for Fokker-Planck problems. Our SIMEX1 and SIMEX2 (Splitting-Implicit-Explicit) schemes are defined within the framework given by the method of lines. The differential operator is discretized by a finite volume scheme given by the Chang-Cooper method, while the integral operator is approximated by a mid-point rule. This leads to a large system of ordinary differential equations, that we approximate with the Strang-Marchuk splitting method. This technique decomposes the original problem in a sequence of different subproblems with simpler structure, which are separately solved and linked to each other through initial conditions and final solutions. After performing the splitting step, we carry out the time integration with first- and second-order time-differencing methods. These steps give rise to the SIMEX1 and SIMEX2 methods, respectively. A full convergence and stability analysis of our schemes is included. Moreover, we are able to prove that the positivity and the mass conservation of the solution to Fokker-Planck problems are satisfied at the discrete level by the numerical solutions computed with the SIMEX schemes. The second main achievement of this thesis is the theoretical analysis and the numerical solution of optimal control problems governed by Fokker-Planck models. The field of optimal control deals with finding control functions in such a way that given cost functionals are minimized. Our framework aims at the minimization of the difference between a known sequence of values and the first moment of a jump-diffusion process; therefore, this formulation can also be considered as a parameter estimation problem for stochastic processes. Two cases are discussed, in which the form of the cost functional is continuous-in-time and discrete-in-time, respectively. The control variable enters the state equation as a coefficient of the Fokker-Planck partial integro-differential operator. We also include in the cost functional a $L^1$-penalization term, which enhances the sparsity of the solution. Therefore, the resulting optimization problem is nonconvex and nonsmooth. We derive the first-order optimality systems satisfied by the optimal solution. The computation of the optimal solution is carried out by means of proximal iterative schemes in an infinite-dimensional framework. N2 - Die vorliegende Arbeit beschäftigt sich mit der theoretischen und numerischen Analyse von Optimalsteuerungsproblemen, deren Nebenbedingungen die Fokker-Planck-Gleichungen von Sprung-Diffusions-Prozessen sind. Unsere Strategie baut auf der Formulierung eines deterministischen Problems auf, um einen stochastischen Prozess zu steuern. Der Ausgangspunkt ist, die Wahrscheinlichkeitsdichtefunktion des Prozesses zu betrachten, deren zeitliche Entwicklung durch die Fokker-Planck-Gleichung modelliert wird. Dieser Ansatz ist vorteilhaft, da er es ermöglicht, den gesamten Bereich des Prozesses durch die Wirkung der Steuerung zu beeinflussen. Zuerst beschäftigen wir uns mit Sprung-Diffusions-Prozessen. Wir definieren Ausgangswertprobleme, die durch stochastische Differentialgleichungen beschrieben werden, und präsentieren Ergebnisse zur Existenz und Eindeutigkeit ihrer Lösungen. Danach diskutieren wir, wie numerische Lösungen stochastischer Probleme berechnet werden, wobei wir uns auf die Euler-Maruyama-Methode konzentrieren. Wir wenden unsere Aufmerksamkeit auf Sprung-Diffusions-Modelle mit zeit- und raumabhängigen Koeffizienten und Sprüngen, die durch einen zusammengesetzten Poisson-Prozess modelliert sind. Wir leiten die zugehörigen Fokker-Planck-Glei-chungen her, die die Form von partiellen Integro-Differentialgleichungen haben. Ihr Differentialterm wird durch einen parabolischen Operator beschrieben, während der nichtlokale Integraloperator Spr\"{u}nge modelliert. Die Ableitung wird auf zwei unterschiedlichen Arten ausgef\"{u}hrt, je nachdem, ob wir einen Prozess mit unbegrenztem oder beschränktem Bereich betrachten. In dem zweiten Fall muss das Verhalten des Prozesses in der Nähe der Grenzen spezifiziert werden; in dieser Arbeit setzen wir reflektierende Grenzen. Die Fokker-Planck-Gleichung, zusammen mit einem Anfangswert und geeigneten Randbedingungen, erzeugt das Fokker-Planck-Problem. Die Lösbarkeit dieses Pro-blems in geeigneten Funktionenräumen und die Eigenschaften dessen Lösung werden diskutiert, nämlich die Positivität und die Wahrscheinlichkeitsmassenerhaltung. Da analytische Lösungen von Fokker-Planck-Problemen oft nicht verfügbar sind, m\"{u}ssen numerische Methoden verwendet werden. Die erste bemerkenswerte Leistung dieser Arbeit ist die Definition und Analyse von konservativen numerischen Verfahren, die Fokker-Planck-Probleme lösen. Unsere SIMEX1 und SIMEX2 (Splitting-Implizit-Explizit) Schemen basieren auf der Linienmethode. Der Differentialoperator wird durch das Finite-Volumen-Schema von Chang und Cooper diskretisiert, während der Integraloperator durch eine Mittelpunktregel angenähert wird. Dies führt zu einem großen System von gewöhnlichen Differentialgleichungen, das mit der Strang-Marchuk-Splitting-Methode gelöst wird. Diese Technik teilt das ursprüngliche Problem in eine Folge verschiedener Teilprobleme mit einer einfachen Struktur, die getrennt gelöst werden und danach durch deren Anfangswerte miteinander verbunden werden. Dank der Splitting-Methode kann jedes Teilproblem implizit oder explizit gelöst werden. Schließlich wird die numerische Integration des Anfangswertsproblems mit zwei Verfahren durchgeführt, n\"{a}mlich dem Euler-Verfahren und dem Predictor-Corrector-Verfahren. Eine umfassende Konvergenz- und Stabilitätsanalyse unserer Systeme ist enthalten. Darüber hinaus können wir beweisen, dass die Positivität und die Massenerhaltung der Lösung von Fokker-Planck-Problemen auf diskreter Ebene durch die numerischen Lösungen erfüllt werden, die mit den SIMEX-Schemen berechnet wurden. Die zweite bemerkenswerte Leistung dieser Arbeit ist die theoretische Analyse und die numerische Behandlung von Optimalsteuerungsproblemen, deren Nebenbedingungen die Fokker-Planck-Probleme von Sprung-Diffusions-Prozessen sind. Der Bereich der optimalen Steuerung befasst sich mit der Suche nach einer optimalen Funktion, die eine gegebene Zielfunktion minimiert. Wir zielen auf die Minimierung des Unterschieds zwischen einer bekannten Folge von Werten und dem ersten Moment eines Sprung-Diffusions-Prozesses. Auf diese Weise kann unsere Formulierung auch als ein Parameterschätzungsproblem für stochastische Prozesse angesehen werden. Zwei Fälle sind erläutert, in denen die Zielfunktion zeitstetig beziehungsweise zeitdiskret ist. Da die Steuerung ein Koeffizient des Integro-Differentialoperators der Zustandsglei-chung ist und die Zielfunktion einen $ L^1 $-Term beinhaltet, der die dünne Besetzung der Lösung erhöht, ist das Optimierungsproblem nichtkonvex und nichtglatt. Die von der optimalen L\"{o}sung erf\"{u}llten notwendigen Bedingungen werden hergeleitet, die man mit einem System beschreiben kann. Die Berechnung optimaler Lösungen wird mithilfe von Proximal-Methoden durchgeführt, die entsprechend um den unendlichdimensionalen Fall erweitert wurden. KW - Numerical analysis KW - Fokker-Planck KW - optimal control problems KW - jump-diffusion processes KW - Fokker-Planck-Gleichung KW - Optimale Kontrolle Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-145645 ER -