4,492 research outputs found
Fast and Efficient Numerical Methods for an Extended Black-Scholes Model
An efficient linear solver plays an important role while solving partial
differential equations (PDEs) and partial integro-differential equations
(PIDEs) type mathematical models. In most cases, the efficiency depends on the
stability and accuracy of the numerical scheme considered. In this article we
consider a PIDE that arises in option pricing theory (financial problems) as
well as in various scientific modeling and deal with two different topics. In
the first part of the article, we study several iterative techniques
(preconditioned) for the PIDE model. A wavelet basis and a Fourier sine basis
have been used to design various preconditioners to improve the convergence
criteria of iterative solvers. We implement a multigrid (MG) iterative method.
In fact, we approximate the problem using a finite difference scheme, then
implement a few preconditioned Krylov subspace methods as well as a MG method
to speed up the computation. Then, in the second part in this study, we analyze
the stability and the accuracy of two different one step schemes to approximate
the model.Comment: 29 pages; 10 figure
Time Domain Computation of a Nonlinear Nonlocal Cochlear Model with Applications to Multitone Interaction in Hearing
A nonlinear nonlocal cochlear model of the transmission line type is studied
in order to capture the multitone interactions and resulting tonal suppression
effects. The model can serve as a module for voice signal processing, it is a
one dimensional (in space) damped dispersive nonlinear PDE based on mechanics
and phenomenology of hearing. It describes the motion of basilar membrane (BM)
in the cochlea driven by input pressure waves. Both elastic damping and
selective longitudinal fluid damping are present. The former is nonlinear and
nonlocal in BM displacement, and plays a key role in capturing tonal
interactions. The latter is active only near the exit boundary (helicotrema),
and is built in to damp out the remaining long waves. The initial boundary
value problem is numerically solved with a semi-implicit second order finite
difference method. Solutions reach a multi-frequency quasi-steady state.
Numerical results are shown on two tone suppression from both high-frequency
and low-frequency sides, consistent with known behavior of two tone
suppression. Suppression effects among three tones are demonstrated by showing
how the response magnitudes of the fixed two tones are reduced as we vary the
third tone in frequency and amplitude. We observe qualitative agreement of our
model solutions with existing cat auditory neural data. The model is thus
simple and efficient as a processing tool for voice signals.Comment: 23 pages,7 figures; added reference
Electronic structure and phase stability of oxide semiconductors: Performance of dielectric-dependent hybrid functional DFT, benchmarked against band structure calculations and experiments
We investigate band gaps, equilibrium structures, and phase stabilities of
several bulk polymorphs of wide-gap oxide semiconductors ZnO, TiO2,ZrO2, and
WO3. We are particularly concerned with assessing the performance of hybrid
functionals built with the fraction of Hartree-Fock exact exchange obtained
from the computed electronic dielectric constant of the material. We provide
comparison with more standard density-functional theory and GW methods. We
finally analyze the chemical reduction of TiO2 into Ti2O3, involving a change
in oxide stoichiometry. We show that the dielectric-dependent hybrid functional
is generally good at reproducing both ground-state (lattice constants, phase
stability sequences, and reaction energies) and excited-state (photoemission
gaps) properties within a single, fully ab initio framework.Comment: Minor changes in the final published versio
Optimising Spatial and Tonal Data for PDE-based Inpainting
Some recent methods for lossy signal and image compression store only a few
selected pixels and fill in the missing structures by inpainting with a partial
differential equation (PDE). Suitable operators include the Laplacian, the
biharmonic operator, and edge-enhancing anisotropic diffusion (EED). The
quality of such approaches depends substantially on the selection of the data
that is kept. Optimising this data in the domain and codomain gives rise to
challenging mathematical problems that shall be addressed in our work.
In the 1D case, we prove results that provide insights into the difficulty of
this problem, and we give evidence that a splitting into spatial and tonal
(i.e. function value) optimisation does hardly deteriorate the results. In the
2D setting, we present generic algorithms that achieve a high reconstruction
quality even if the specified data is very sparse. To optimise the spatial
data, we use a probabilistic sparsification, followed by a nonlocal pixel
exchange that avoids getting trapped in bad local optima. After this spatial
optimisation we perform a tonal optimisation that modifies the function values
in order to reduce the global reconstruction error. For homogeneous diffusion
inpainting, this comes down to a least squares problem for which we prove that
it has a unique solution. We demonstrate that it can be found efficiently with
a gradient descent approach that is accelerated with fast explicit diffusion
(FED) cycles. Our framework allows to specify the desired density of the
inpainting mask a priori. Moreover, is more generic than other data
optimisation approaches for the sparse inpainting problem, since it can also be
extended to nonlinear inpainting operators such as EED. This is exploited to
achieve reconstructions with state-of-the-art quality.
We also give an extensive literature survey on PDE-based image compression
methods
- …