Search results for "numerical"

Mathematical Methods for Research Excellence : Book of Abstracts

2022

Reverse Catmull-Clark Subdivision

2006

Reverse subdivision consists in constructing a coarse mesh of a model from a finer mesh of this same model. In this paper, we give formulas for reverse Catmull-Clark subdivision. These formulas allow the constructing of a coarse mesh for almost all meshes. The condition for being able to apply these formulas is that the mesh to be reversed must be generated by the subdivision of a coarse mesh. Except for this condition, the mesh can be arbitrary. Vertices can be regular or extraordinary and the mesh itself can be arbitrary (triangular, quadrilateral…).

A Continuous Approach to FETI-DP Mortar Methods: Application to Dirichlet and Stokes Problem

2013

In this contribution we extend the FETI-DP mortar method for elliptic problems introduced by Bernardi et al. [2] and Chacon Vera [3] to the case of the incompressible Stokes equations showing that the same results hold in the two dimensional setting. These ideas extend easily to three dimensional problems. Finally some numerical tests are shown as a conclusion. This contribution is a condensed version of a more detailed forthcoming paper. We use standard notation, see for instance [1].

Liftings and extensions of operators in Brownian setting

2020

We investigate the operators T on a Hilbert space H which have 2-isometric liftings S with the property S ∗ S H ⊂ H . We show that such liftings are closely related to some extensions of T, which h...

Surrogate-assisted evolutionary multiobjective shape optimization of an air intake ventilation system

2017

We tackle three different challenges in solving a real-world industrial problem: formulating the optimization problem, connecting different simulation tools and dealing with computationally expensive objective functions. The problem to be optimized is an air intake ventilation system of a tractor and consists of three computationally expensive objective functions. We describe the modeling of the system and its numerical evaluation with a commercial software. To obtain solutions in few function evaluations, a recently proposed surrogate-assisted evolutionary algorithm K-RVEA is applied. The diameters of four different outlets of the ventilation system are considered as decision variables. Fr…

Numerical Recovery of Source Singularities via the Radiative Transfer Equation with Partial Data

2013

The inverse source problem for the radiative transfer equation is considered, with partial data. Here we demonstrate numerical computation of the normal operator $X_{V}^{*}X_{V}$ where $X_{V}$ is the partial data solution operator to the radiative transfer equation. The numerical scheme is based in part on a forward solver designed by F. Monard and G. Bal. We will see that one can detect quite well the visible singularities of an internal optical source $f$ for generic anisotropic $k$ and $\sigma$, with or without noise added to the accessible data $X_{V}f$. In particular, we use a truncated Neumann series to estimate $X_{V}$ and $X_{V}^{*}$, which provides a good approximation of $X_{V}^{*…

A posteriori error estimates for time-dependent reaction-diffusion problems based on the Payne-Weinberger inequality

2015

We consider evolutionary reaction-diffusion problem with mixed Dirichlet--Robin boundary conditions. For this class of problems, we derive two-sided estimates of the distance between any function in the admissible energy space and exact solution of the problem. The estimates (majorants and minorants) are explicitly computable and do not contain unknown functions or constants. Moreover, it is proved that the estimates are equivalent to the energy norm of the deviation from the exact solution.

Reduced Order Models for Pricing European and American Options under Stochastic Volatility and Jump-Diffusion Models

2017

Abstract European options can be priced by solving parabolic partial(-integro) differential equations under stochastic volatility and jump-diffusion models like the Heston, Merton, and Bates models. American option prices can be obtained by solving linear complementary problems (LCPs) with the same operators. A finite difference discretization leads to a so-called full order model (FOM). Reduced order models (ROMs) are derived employing proper orthogonal decomposition (POD). The early exercise constraint of American options is enforced by a penalty on subset of grid points. The presented numerical experiments demonstrate that pricing with ROMs can be orders of magnitude faster within a give…

Reduced Order Models for Pricing American Options under Stochastic Volatility and Jump-diffusion Models

2016

American options can be priced by solving linear complementary problems (LCPs) with parabolic partial(-integro) differential operators under stochastic volatility and jump-diffusion models like Heston, Merton, and Bates models. These operators are discretized using finite difference methods leading to a so-called full order model (FOM). Here reduced order models (ROMs) are derived employing proper orthogonal decomposition (POD) and non negative matrix factorization (NNMF) in order to make pricing much faster within a given model parameter variation range. The numerical experiments demonstrate orders of magnitude faster pricing with ROMs. peerReviewed

IMEX schemes for pricing options under jump–diffusion models

2014

We propose families of IMEX time discretization schemes for the partial integro-differential equation derived for the pricing of options under a jump-diffusion process. The schemes include the families of IMEX-midpoint, IMEX-CNAB and IMEX-BDF2 schemes. Each family is defined by a convex combination parameter [email protected]?[0,1], which divides the zeroth-order term due to the jumps between the implicit and explicit parts in the time discretization. These IMEX schemes lead to tridiagonal systems, which can be solved extremely efficiently. The schemes are studied through Fourier stability analysis and numerical experiments. It is found that, under suitable assumptions and time step restric…