×

Using distributed computers to deterministically approximate higher dimensional convection diffusion equations. (English) Zbl 1067.68021

Summary: A new approach to solving \(D > 3\) spatial dimensional convection-diffusion equations on clusters of workstations is derived by exploiting the stability and scalability of the combination of a generalized \(D\)-dimensional High-Order Compact (HOC) implicit finite difference scheme and parallelized GMRES(m). We then consider its application to multifactor option pricing using the Black-Scholes equation and further show that an isotropic fourth order compact difference scheme is numerically stable and determine conditions under which its coefficient matrix is positive definite. The performance of GMRES(m) on distributed computers is limited by the inter-processor communication required by the matrix-vector multiplication. It is shown that the compact scheme requires approximately half the number of communications as a non-compact difference scheme of the same order of truncation error. As the dimensionality is increased, the ratio of computation that can be overlapped with communication also increases. CPU times and parallel efficiency graphs for single time step approximation of up to a 7D HOC scheme on 16 processors confirm the numerical stability constraint and demonstrate improved parallel scalability over non-compact difference schemes.

MSC:

68M14 Distributed systems
68W15 Distributed algorithms
65Y10 Numerical algorithms for specific classes of architectures
91B28 Finance etc. (MSC2000)

Software:

PETSc
PDFBibTeX XMLCite
Full Text: DOI