Search Results
Global analytic expansion of solution for a class of linear parabolic systems with coupling of first order derivatives terms
Literaturverz. We derive global analytic representations of fundamental solutions for a class of linear parabolic systems with full coupling of first order derivative terms where coefficients may depend on space and time. Pointwise convergence of the global analytic expansion is proved. This leads to analytic representations of solutions of initial-boundary problems of first and second type in terms of convolution integrals or convolution integrals and linear integral equations. The results have both analytical and numerical impact. Analytically, our representations of fundamental solutions of coupled parabolic systems may be used to define generalized stochastic processes. Moreover, some classical analytical results based on a priori estimates of elliptic equations are a simple corollary of our main result. Numerically, accurate, stable and efficient schemes for computation and error estimates in strong norms can be obtained for a considerable class of Cauchy- and initial-boundary problems of parabolic type. Furthermore, there are obvious and less obvious applications to finance and physics.
Regular polynomial interpolation and approximation of global solutions of linear partial differential equations
We consider regular polynomial interpolation algorithms on recursively defined sets of interpolation points which approximate global solutions of arbitrary well-posed systems of linear partial differential equations. Convergence of the "limit" of the recursively constructed family of polynomials to the solution and error estimates are obtained from a priori estimates for some standard classes of linear partial differential equations, i.e. elliptic and hyperbolic equations. Another variation of the algorithm allows to construct polynomial interpolations which preserve systems of linear partial differential equations at the interpolation points. We show how this can be applied in order to compute higher order terms of WKB-approximations of fundamental solutions of a large class of linear parabolic equations. The error estimates are sensitive to the regularity of the solution. Our method is compatible with recent developments for solution of higher dimensional partial differential equations, i.e. (adaptive) sparse grids, and weighted Monte-Carlo, and has obvious applications to mathematical finance and physics.
Global regularity and probabilistic schemes for free boundary surfaces of multivariate American derivatives and their Greeks
In a rather general setting of multivariate stochastic volatility market models we derive global iterative probabilistic schemes for computing the free boundary and its Greeks for a generic class of American derivative models using front-fixing methods. Establishment of convergence is closely linked to a proof of global regularity of the free boundary surface.
How to compute the length of a geodesic on a Riemannian manifold with small error in arbitrary Sobolev norms
We compute the length of geodesics on a Riemannian manifold by regular polynomial interpolation of the global solution of the eikonal equation related to the line element $ds^2=g_ijdx^idx^j$ of the manifold. Our algorithm approximates the length functional in arbitrarily strong Sobolev norms. Error estimates are obtained where the geometric information is used. It is pointed out how the algorithm can be used to get accurate approximations of solutions of linear parabolic partial differential equations leading to obvious applications in finance, physics and other sciences