Smoothed Particle Interpolation

Smoothed particle hydrodynamics (SPH) discretization techniques are generalized to develop a method, smoothed particle interpolation (SPI), for solving initial value problems of systems of non-hydrodynamical nature. Under this approach, SPH is viewed as strictly an interpolation scheme and, as such, suitable for solving general hyperbolic and parabolic equations. The SPI method is tested on (1) the wave equation with inhomogeneous sound speed and (2) Burgers equation. The efficiency of SPI is studied by comparing SPI solutions to those obtained with standard finite difference methods. It is shown that the power of SPI arises when the smoothing particles are free to move.Comment: 13 pages (LaTeX), 9 figures (not included), [email protected]

Ellipse-preserving Hermite interpolation and subdivision

We introduce a family of piecewise-exponential functions that have the Hermite interpolation property. Our design is motivated by the search for an effective scheme for the joint interpolation of points and associated tangents on a curve with the ability to perfectly reproduce ellipses. We prove that the proposed Hermite functions form a Riesz basis and that they reproduce prescribed exponential polynomials. We present a method based on Green's functions to unravel their multi-resolution and approximation-theoretic properties. Finally, we derive the corresponding vector and scalar subdivision schemes, which lend themselves to a fast implementation. The proposed vector scheme is interpolatory and level-dependent, but its asymptotic behaviour is the same as the classical cubic Hermite spline algorithm. The same convergence properties---i.e., fourth order of approximation---are hence ensured

Efficient Numerical Methods for Pricing American Options under Lévy Models

Two new numerical methods for the valuation of American and Bermudan options are proposed, which admit a large class of asset price models for the underlying. In particular, the methods can be applied with Lévy models that admit jumps in the asset price. These models provide a more realistic description of market prices and lead to better calibration results than the well-known Black-Scholes model. The proposed methods are not based on the indirect approach via partial differential equations, but directly compute option prices as risk-neutral expectation values. The expectation values are approximated by numerical quadrature methods. While this approach is initially limited to European options, the proposed combination with interpolation methods also allows for pricing of Bermudan and American options. Two different interpolation methods are used. These are cubic splines on the one hand and a mesh-free interpolation by radial basis functions on the other hand. The resulting valuation methods allow for an adaptive space discretization and error control. Their numerical properties are analyzed and, finally, the methods are validated and tested against various single-asset and multi-asset options under different market models

Multivariate orthonormal interpolating scaling vectors

AbstractIn this paper we introduce an algorithm for the construction of interpolating scaling vectors on Rd with compact support and orthonormal integer translates. Our method is substantiated by constructing several examples of bivariate scaling vectors for quincunx and box–spline dilation matrices. As the main ingredients of our recipe we derive some implementable conditions for accuracy and orthonormality of an interpolating scaling vector in terms of its mask

Space-Varying Coefficient Models for Brain Imaging

The methodological development and the application in this paper originate from diffusion tensor imaging (DTI), a powerful nuclear magnetic resonance technique enabling diagnosis and monitoring of several diseases as well as reconstruction of neural pathways. We reformulate the current analysis framework of separate voxelwise regressions as a 3d space-varying coefficient model (VCM) for the entire set of DTI images recorded on a 3d grid of voxels. Hence by allowing to borrow strength from spatially adjacent voxels, to smooth noisy observations, and to estimate diffusion tensors at any location within the brain, the three-step cascade of standard data processing is overcome simultaneously. We conceptualize two VCM variants based on B-spline basis functions: a full tensor product approach and a sequential approximation, rendering the VCM numerically and computationally feasible even for the huge dimension of the joint model in a realistic setup. A simulation study shows that both approaches outperform the standard method of voxelwise regressions with subsequent regularization. Due to major efficacy, we apply the sequential method to a clinical DTI data set and demonstrate the inherent ability of increasing the rigid grid resolution by evaluating the incorporated basis functions at intermediate points. In conclusion, the suggested fitting methods clearly improve the current state-of-the-art, but ameloriation of local adaptivity remains desirable