Pricing early-exercise and discrete barrier options by Shannon wavelet expansions

We present a pricing method based on Shannon wavelet expansions for early-exercise and discretely-monitored barrier options under exponential Lévy asset dynamics. Shannon wavelets are smooth, and thus approximate the densities that occur in finance well, resulting in exponential convergence. Application of the Fast Fourier Transform yields an efficient implementation and since wavelets give local approximations, the domain boundary errors can be naturally resolved, which is the main improvement over existing methods

Multigrid method for nonlinear poroelasticity equations

<p>In this study, a nonlinear multigrid method is applied for solving the system of incompressible poroelasticity equations considering nonlinear hydraulic conductivity. For the unsteady problem, an additional artificial term is utilized to stabilize the solutions when the equations are discretized on collocated grids. We employ two nonlinear multigrid methods, i.e. the “full approximation scheme” and “Newton multigrid” for solving the corresponding system of equations arising after discretization. For the steady case, both homogeneous and heterogeneous cases are solved and two different smoothers are examined to search for an efficient multigrid method. Numerical results show a good convergence performance for all the strategies.</p

Two-dimensional Shannon wavelet inverse Fourier technique for pricing European options

The SWIFT method for pricing European-style options on one underlying asset was recently published and presented as an accurate, robust and highly efficient technique. The purpose of this paper is to extend the method to higher dimensions by pricing exotic option contracts, called rainbow options, whose payoff depends on multiple assets. The multidimensional extension inherits the properties of the one-dimensional method, being the exponential convergence one of them. Thanks to the nature of local Shannon wavelets basis, we do not need to rely on a-priori truncation of the integration range, we have an error bound estimate and we use fast Fourier transform (FFT) algorithms to speed up computations. We test the method for similar examples with state-of-the-art methods found in the literature, and we compare our results with analytical expressions when available

Quantifying credit portfolio losses under multi-factor models

In this work, we investigate the challenging problem of estimating credit risk measures of portfolios with exposure concentration under the multi-factor Gaussian and multi-factor t-copula models. It is well-known that Monte Carlo (MC) methods are highly demanding from the computational point of view in the aforementioned situations. We present efficient and robust numerical techniques based on the Haar wavelets theory for recovering the cumulative distribution function of the loss variable from its characteristic function. To the best of our knowledge, this is the first time that multi-factor t-copula models are considered outside the MC framework. The analysis of the approximation error and the results obtained in the numerical experiments section show a reliable and useful machinery for credit risk capital measurement purposes in line with Pillar II of the Basel Accords

On the data-driven COS method

In this paper, we present the data-driven COS method, ddCOS. It is a Fourier-based financial option valuation method which assumes the availability of asset data samples: a characteristic function of the underlying asset probability density function is not required. As such, the presented technique represents a generalization of the well-known COS method [1]. The convergence of the proposed method is in line with Monte Carlo methods for pricing financial derivatives. The ddCOS method is then particularly interesting for density recovery and also for the efficient computation of the option's sensitivities Delta and Gamma. These are often used in risk management, and can be obtained at a higher accuracy with ddCOS than with plain Monte Carlo methods

Energy-stable discretization of the one-dimensional two-fluid model

In this paper we present a complete framework for the energy-stable simulation of stratified incompressible flow in channels, using the one-dimensional two-fluid model. Building on earlier energy-conserving work on the basic two-fluid model, our new framework includes diffusion, friction, and surface tension. We show that surface tension can be added in an energy-conserving manner, and that diffusion and friction have a strictly dissipative effect on the energy. We then propose spatial discretizations for these terms such that a semi-discrete model is obtained that has the same conservation properties as the continuous model. Additionally, we propose a new energy-stable advective flux scheme that is energy-conserving in smooth regions of the flow and strictly dissipative where sharp gradients appear. This is obtained by combining, using flux limiters, a previously developed energy-conserving advective flux with a novel first-order upwind scheme that is shown to be strictly dissipative. The complete framework, with diffusion, surface tension, and a bounded energy, is linearly stable to short wavelength perturbations, and exhibits nonlinear damping near shocks. The model yields smoothly converging numerical solutions, even under conditions for which the basic two-fluid model is ill-posed. With our explicit expressions for the dissipation rates, we are able to attribute the nonlinear damping to the different dissipation mechanisms, and compare their effects

Efficient Wrong-Way Risk Modelling for Funding Valuation Adjustments

Wrong-Way Risk (WWR) is an important component in Funding Valuation Adjustment (FVA) modelling. Yet, it can be challenging to compute WWR efficiently. We propose to split the relevant exposure profile into two parts: an independent part and a WWR-driven part. For the first part, already available exposures can be used where correlations between the funding spread and market risks are ignored. We express the second part of the exposure profile in terms of the stochastic drivers and approximate these by a common Gaussian stochastic factor. The proposed approximation is generic, is an add-on to the existing xVA calculations and provides an efficient and robust way to include WWR in FVA modelling. Furthermore, the approximation provides some intuition on WWR. Case studies are presented for an interest rate swap and a representative multi-currency portfolio of swaps. They illustrate that the approximation method is applicable in a practical setting due to its generic nature. We analyze the approximation error and illustrate how the approximation can be used to compute WWR sensitivities, which are needed for risk management