Spring Meeting WSC 2009
Op vrijdag 24 april 2009 organiseert de Werkgemeenschap Scientific Computing
een voorjaarsbijeenkomst aan de TU-Eindhoven.
Op deze bijeenkomst zullen AIOs/OIOs en senior onderzoekers op het gebied van
Scientific Computing een voordracht geven.
Als locatie is de filmzaal van
Zwarte Doos gekozen van de TU-Eindhoven.
| 10.00-10.30 uur
||Registratie, koffie/thee, welkom
| 10.30-11.15 uur
|| Stefan Vandewalle (KU Leuven)
||Multigrid methods for partial differential equations with stochastic coefficients
| 11.15-11.40 uur
|| Tammo Jan Dijkema (UU)
||Approximation in high dimensional product domains
| 11.40-12.05 uur
|| Kim Volders (UA)
||Stability of finite difference schemes on nonuniform grids for the BlackScholes equation
| 12.05-12.30 uur
|| Hisham bin Zubair (TUD-EWI)
||Multigrid Preconditioners for the indefinite Helmholtz Problems on Locally Refined Grids
| 12.30-13.30 uur
||Lunch in de lounge van de Zwarte Doos
| 13.30-14.15 uur
|| Michiel Hochstenbach (TU/e)
||Eigenvalue problems and DDE stability[abstract]
| 14.15-14.40 uur
|| Ricardo Reis da Silva (UvA-KdVI)
||Computable bounds for the smallest singular value of a rank-k perturbation of a matrix
| 14.40-15.05 uur
|| Maria Ugryumova (TU/e)
||Stability and Passivity of the Super Node Algorithm for EM modelling of ICs
| 15.05-15.15 uur
| 15.15-16.00 uur
|| Kees Oosterlee (CWI,TUD-EWI)
||Numerical Mathematics Aspects in Computational Finance
Prof.dr. J.G. Verwer (CWI,UvA-KdvI), Prof.dr. W.H.A. Schilders (TU/e,NXP),
Enna van Dijk (CASA), Drs. Margreet Nool (CWI, secretaris).
Deelname aan deze voorjaarsbijeenkomst (inclusief lunch) is gratis. U wordt wel verzocht u te
registreren voor 15 april 2009.
Voor vragen: Margreet Nool en Enna van Dijk
Multigrid methods for partial differential equations with stochastic coefficients
The stochastic finite element method is an important technique for solving
certain classes of stochastic partial differential equations (PDEs). This
method approximates the solution of the PDE by a generalized polynomial
chaos expansion. By using a Galerkin projection in the stochastic dimension,
the stochastic PDE is transformed into a coupled set of deterministic PDEs.
A finite element discretization converts this deterministic PDE system into
a high dimensional algebraic system. Specialized iterative solvers are
required to solve the resulting problem.
In this talk, we shall present an overview of iterative solution approaches.
We start from iterative methods based on a block splitting of the system
matrices. Next, we extend these methods for use as preconditioner for a
Krylov method, and for use as smoother in a multilevel context. Then, the
various solvers will be compared based on their convergence properties,
computational cost and implementation effort. Our findings are illustrated
by means of two numerical problems. The first one is a steady-state
diffusion problem with a discontinuous random field as diffusion
coefficient. The second is a deterministic diffusion problem defined on a
Tammo Jan Dijkema
Approximation in high dimensional product domains
In this talk I will recall the concept of sparse grids, which can be used to approximate
functions on product domains with a rate that is almost independent of the space dimension.
So-called optimized sparse grids make this approximation rate optimal,
and truly dimension independent. However nice these results may seem,
the smoothness that is required for these methods to work, is rather high.
I will show that the solution of the Poisson equation generally is not smooth enough.
As a solution to this problem, an adaptive method can be used, yielding the same optimal
rate for a broader class of functions. I will introduce an adaptive wavelet method
that was used to solve Poisson's equation in up to 10 space dimensions.
Stability of finite difference schemes on nonuniform grids for the BlackScholes equation
We consider the well-known Black-Scholes equation from financial option pricing theory,
with given real constants
The Black-Scholes equation is a time-dependent advection-diffusion-reaction equation and is supplemented with initial
and boundary conditions.
A popular approach for the numerical solution of time-dependent partial differential equations is the
method-of-lines. It consists of two steps:
- Spatial discretization: the partial derivatives
are discretized on a finite spatial grid, yielding a (large) system of ordinary differential equations
with given fixed matrix
- Temporal discretization: the above system of ordinary differential equations is numerically integrated in time.
Our research focuses on the stability analysis of second-order finite difference methods for the spatial
discretization of the Black-Scholes equation.
We first present practical upper bounds for
denotes a scaled version of the standard spectral norm.
We subsequently present sufficient conditions for contractivity in the maximum-norm,
A virtue of our stability analysis is that it applies to spatial grids that are not uniform.
Such grids are often used in actual applications.
Numerical experiments are provided which support our theoretical results.
Finally, we briefly discuss the stability of temporal discretization schemes for (1)
Hisham bin Zubair
Multigrid Preconditioners for the indefinite Helmholtz Problems on Locally Refined Grids
In this talk we present the construction and performance of a geometric multigrid method for grids having
two different layers of refinement. The method is employed to approximately invert the Krylovpreconditioner,
i.e., the complex shifted Helmholtz operator, for the indefinite Helmholtz equation.
The usual FAC and MLAT based grid coarsening techniques only coarsen the fine layer of the grid.
The method presented here, in contrast, coarsens the whole grid simultaneously.
Combined with a simple smoother, piece-wise constant restriction and bilinear prolongation,
this gives an efficient multigrid method that works very well for the model problems,
which are 2-d Helmholtz equations with strongly varying coefficients.
Eigenvalue problems and DDE stability
We review some recent results in the field of stability of delay differential equations (DDEs) and show how
various types of eigenvalue problems (including new types) play important roles.
As these problems are often of a very high dimension, we also sketch structure-preserving
methods to numerically solve them.
Ricardo Reis da Silva
Computable bounds for the smallest singular value of a rank-k perturbation of a matrix
In this talk we look into the computation of lower bounds for the smallest singular values
of perturbations of matrices.
Our starting point is the smallest eigenvalue of a seemingly simple Hermitian rank-one
perturbation H of a Hermitian matrix A.
We will see how both Weyl's bound and a recent bound of Ipsen and Nadler are the two
first terms of a sequence of bounds.
This sequence is non-decreasing and, in general, the q-th term is the smallest
eigenvalue of a q × q matrix.
Similar bounds can be obtained for the smallest eigenvalue of rank-k
perturbations of Hermitian matrices and for the smallest singular value of perturbations
of arbitrary n×m matrices B.
Stability and Passivity of the Super Node Algorithm for EM modelling of ICs
The super node algorithm is a model order reduction technique based on physical principles.
Some of the properties of the reduced models generated by this algorithm, such as stability and
passivity have not yet been studied thoroughly.
The loss of passivity constitutes a serious problem because the reduced networks may show
artificial behavior which renders the simulations useless.
We investigate the stability and passivity properties of the algorithm.
We explain why passivity is not guaranteed and we present a way to modify the algorithm in order
to provide always passive reduced models.
Numerical Mathematics Aspects in Computational Finance
In this presentation we will discuss some topics in Finance that require
Mathematics, and, in particular, efficient numerical techniques.
We will discuss option pricing, for example, for option contracts based on
more than one underlying stock, and contracts with advanced stochastic
models for the underlying stock price dynamics.
We will focus on a mathematical framework in which we perform this research.
One of the aims is to price the options as fast as possible.