Source author record

Kristian Debrabant

Kristian Debrabant appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

A Comparison of Different Approaches to Dynamic Origin-Destination Matrix Estimation in Urban Traffic

Given the counters of vehicles that traverse the roads of a traffic network, we reconstruct the travel demand that generated them expressed in terms of the number of origin-destination trips made by users. We model the problem as a bi-level optimization problem. At the inner-level, given a tentative demand, we solve a Dynamic Traffic Assignment (DTA) problem to decide the routing of the users between their origins and destinations. Finally, we adjust the number of trips and their origins and destinations at the outer-level to minimize the discrepancy between the counters generated at the inner-level and the given vehicle counts measured by sensors in the traffic network. We solve the DTA problem by employing a mesoscopic model implemented by the traffic simulator SUMO. Thus, the outer problem becomes an optimization problem that minimizes a black-box Objective Function (OF) determined by the results of the simulation, which is a costly computation. We study different approaches to the outer-level problem categorized as gradient-based and derivative-free approaches. Among the gradient-based approaches, we look at an assignment matrix-based approach and an assignment matrix-free approach that uses the Simultaneous Perturbation Stochastic Approximation (SPSA) algorithm. Among the derivative-free approaches, we investigate Machine Learning (ML) algorithms to learn a model of the simulator that can then be used as a surrogate OF in the optimization problem. We compare these approaches computationally on an artificial network. The gradient-based approaches perform the best in terms of solution quality and computational requirements. In contrast, the results obtained by the ML approach are currently less satisfactory but provide an interesting avenue for future research.

preprint2020arXiv

High order numerical integrators for single integrand Stratonovich SDEs

We show that applying any deterministic B-series method of order $p_d$ with a random step size to single integrand SDEs gives a numerical method converging in the mean-square and weak sense with order $\lfloor p_d/2\rfloor$.As an application, we derive high order energy-preserving methods for stochastic Poisson systems as well as further geometric numerical schemes for this wide class of Stratonovich SDEs.

preprint2016arXiv

Diagonally drift-implicit Runge-Kutta methods of weak order one and two for Itô SDEs and stability analysis

The class of stochastic Runge-Kutta methods for stochastic differential equations due to Rößler is considered. Coefficient families of diagonally drift-implicit stochastic Runge-Kutta (DDISRK) methods of weak order one and two are calculated. Their asymptotic stability as well as mean-square stability (MS-stability) properties are studied for a linear stochastic test equation with multiplicative noise. The stability functions for the DDISRK methods are determined and their domains of stability are compared to the corresponding domain of stability of the considered test equation. Stability regions are presented for various coefficients of the families of DDISRK methods in order to determine step size restrictions such that the numerical approximation reproduces the characteristics of the solution process.

preprint2014arXiv

On the Acceleration of the Multi-Level Monte Carlo Method

The multi-level Monte Carlo method proposed by M. Giles (2008) approximates the expectation of some functionals applied to a stochastic process with optimal order of convergence for the mean-square error. In this paper, a modified multi-level Monte Carlo estimator is proposed with significantly reduced computational costs. As the main result, it is proved that the modified estimator reduces the computational costs asymptotically by a factor $(p/α)^2$ if weak approximation methods of orders $α$ and $p$ are applied in case of computational costs growing with same order as variances decay.

preprint2014arXiv

Semi-Lagrangian schemes for linear and fully non-linear Hamilton-Jacobi-Bellman equations

We consider the numerical solution of Hamilton-Jacobi-Bellman equations arising in stochastic control theory. We introduce a class of monotone approximation schemes relying on monotone interpolation. These schemes converge under very weak assumptions, including the case of arbitrary degenerate diffusions. Besides providing a unifying framework that includes several known first order accurate schemes, stability and convergence results are given, along with two different robust error estimates. Finally, the method is applied to a super-replication problem from finance.

preprint2013arXiv

Classification of Stochastic Runge-Kutta Methods for the Weak Approximation of Stochastic Differential Equations

In the present paper, a class of stochastic Runge-Kutta methods containing the second order stochastic Runge-Kutta scheme due to E. Platen for the weak approximation of Itô stochastic differential equation systems with a multi-dimensional Wiener process is considered. Order one and order two conditions for the coefficients of explicit stochastic Runge-Kutta methods are solved and the solution space of the possible coefficients is analyzed. A full classification of the coefficients for such stochastic Runge-Kutta schemes of order one and two with minimal stage numbers is calculated. Further, within the considered class of stochastic Runge-Kutta schemes coefficients for optimal schemes in the sense that additionally some higher order conditions are fulfilled are presented.

preprint2013arXiv

Continuous Weak Approximation for Stochastic Differential Equations

A convergence theorem for the continuous weak approximation of the solution of stochastic differential equations by general one step methods is proved, which is an extension of a theorem due to Milstein. As an application, uniform second order conditions for a class of continuous stochastic Runge-Kutta methods containing the continuous extension of the second order stochastic Runge-Kutta scheme due to Platen are derived. Further, some coefficients for optimal continuous schemes applicable to Itô stochastic differential equations with respect to a multi-dimensional Wiener process are presented.

preprint2013arXiv

Convergence of Runge-Kutta Methods Applied to Linear Partial Differential-Algebraic Equations

We apply Runge-Kutta methods to linear partial differential-algebraic equations of the form $Au_t(t,x) + B(u_{xx}(t,x)+ru_x(t,x))+Cu(t,x) = f(t,x)$, where $A,B,C\in\R^{n,n}$ and the matrix $A$ is singular. We prove that under certain conditions the temporal convergence order of the fully discrete scheme depends on the time index of the partial differential-algebraic equation. In particular, fractional orders of convergence in time are encountered. Furthermore we show that the fully discrete scheme suffers an order reduction caused by the boundary conditions. Numerical examples confirm the theoretical results.

preprint2013arXiv

Families of efficient second order Runge-Kutta methods for the weak approximation of Itô stochastic differential equations

Recently, a new class of second order Runge-Kutta methods for Itô stochastic differential equations with a multidimensional Wiener process was introduced by Rößler. In contrast to second order methods earlier proposed by other authors, this class has the advantage that the number of function evaluations depends only linearly on the number of Wiener processes and not quadratically. In this paper, we give a full classification of the coefficients of all explicit methods with minimal stage number. Based on this classification, we calculate the coefficients of an extension with minimized error constant of the well-known RK32 method to the stochastic case. For three examples, this method is compared numerically with known order two methods and yields very promising results.

preprint2013arXiv

On quasi-linear PDAEs with convection: applications, indices, numerical solution

For a class of partial differential algebraic equations (PDAEs) of quasi-linear type which include nonlinear terms of convection type a possibility to determine a time and spatial index is considered. As a typical example we investigate an application from plasma physics. Especially we discuss the numerical solution of initial boundary value problems by means of a corresponding finite difference splitting procedure which is a modification of a well known fractional step method coupled with a matrix factorization. The convergence of the numerical solution towards the exact solution of the corresponding initial boundary value problem is investigated. Some results of a numerical solution of the plasma PDAE are given.

preprint2012arXiv

Semi-Lagrangian schemes for linear and fully non-linear diffusion equations

For linear and fully non-linear diffusion equations of Bellman-Isaacs type, we introduce a class of approximation schemes based on differencing and interpolation. As opposed to classical numerical methods, these schemes work for general diffusions with coefficient matrices that may be non-diagonal dominant and arbitrarily degenerate. In general such schemes have to have a wide stencil. Besides providing a unifying framework for several known first order accurate schemes, our class of schemes includes new first and higher order versions. The methods are easy to implement and more efficient than some other known schemes. We prove consistency and stability of the methods, and for the monotone first order methods, we prove convergence in the general case and robust error estimates in the convex case. The methods are extensively tested.

preprint2011arXiv

A micro/macro algorithm to accelerate Monte Carlo simulation of stochastic differential equations

We present and analyze a micro/macro acceleration technique for the Monte Carlo simulation of stochastic differential equations (SDEs) in which there is a separation between the (fast) time-scale on which individual trajectories of the SDE need to be simulated and the (slow) time-scale on which we want to observe the (macroscopic) function of interest. The method performs short bursts of microscopic simulation using an ensemble of SDE realizations, after which the ensemble is restricted to a number of macroscopic state variables. The resulting macroscopic state is then extrapolated forward in time and the ensemble is projected onto the extrapolated macroscopic state. We provide a first analysis of its convergence in terms of extrapolation time step and number of macroscopic state variables. The effects of the different approximations on the resulting error are illustrated via numerical experiments.

preprint2010arXiv

Composition of stochastic B-series with applications to implicit Taylor methods

In this article, we construct a representation formula for stochastic B-series evaluated in a B-series. This formula is used to give for the first time the order conditions of implicit Taylor methods in terms of rooted trees. Finally, as an example we apply these order conditions to derive in a simple manner a family of strong order 1.5 Taylor methods applicable to Itô SDEs.

preprint2010arXiv

Runge-Kutta methods for third order weak approximation of SDEs with multidimensional additive noise

A new class of third order Runge-Kutta methods for stochastic differential equations with additive noise is introduced. In contrast to Platen's method, which to the knowledge of the author has been up to now the only known third order Runge-Kutta scheme for weak approximation, the new class of methods affords less random variable evaluations and is also applicable to SDEs with multidimensional noise. Order conditions up to order three are calculated and coefficients of a four stage third order method are given. This method has deterministic order four and minimized error constants, and needs in addition less function evaluations than the method of Platen. Applied to some examples, the new method is compared numerically with Platen's method and some well known second order methods and yields very promising results.

preprint2010arXiv

Stochastic B-series analysis of iterated Taylor methods

For stochastic implicit Taylor methods that use an iterative scheme to compute their numerical solution, stochastic B--series and corresponding growth functions are constructed. From these, convergence results based on the order of the underlying Taylor method, the choice of the iteration method, the predictor and the number of iterations, for Itô and Stratonovich SDEs, and for weak as well as strong convergence are derived. As special case, also the application of Taylor methods to ODEs is considered. The theory is supported by numerical experiments.