Abstract. Arbitrage-free prices u of European contracts on risky assets whose log-returns are modelled by Lévy processes satisfy a parabolic partial integro-differential equationThis PIDE is localized to bounded domains and the error due to this localization is estimated. The localized PIDE is discretized by the θ-scheme in time and a wavelet Galerkin method with N degrees of freedom in log-price space. The dense matrix for A can be replaced by a sparse matrix in the wavelet basis, and the linear systems in each implicit time step are solved approximatively with GMRES in linear complexity. The total work of the algorithm for M time steps is bounded by O (MN(log(N )) 2 ) operations and O(N log(N )) memory. The deterministic algorithm gives optimal convergence rates (up to logarithmic terms) for the computed solution in the same complexity as finite difference approximations of the standard Black-Scholes equation. Computational examples for various Lévy price processes are presented.Mathematics Subject Classification. 65N30, 60J75.

No abstract

Abstract. The convergence of a two-scale FEM for elliptic problems in divergence form with coefficients and geometries oscillating at length scale ε 1 is analyzed. Full elliptic regularity independent of ε is shown when the solution is viewed as mapping from the slow into the fast scale. Two-scale FE spaces which are able to resolve the ε scale of the solution with work independent of ε and without analytical homogenization are introduced. Robust in ε error estimates for the two-scale FE spaces are proved. Numerical experiments confirm the theoretical analysis.Mathematics Subject Classification. 65N30.

No abstract

A new nite element method for elliptic problems with locally periodic microstructure of length " > 0 is developed and analyzed. It is shown that the method converges, as " ! 0, to the solution of the homogenized problem with optimal order in " and exponentially in the number of degrees of freedom independent o f " > 0. The computational work of the method is bounded independently of ". Numerical experiments demonstrate the feasibility and con rm the theoretical results.