SIAM Digital Library
 
 
 

You are not logged in Logged Out Log In

SIAM J. Sci. Comput. 28, pp. 1634-1650 (17 pages)

A Cache‐Aware Algorithm for PDEs on Hierarchical Data Structures Based on Space‐Filling Curves

Frank Günther, Miriam Mehl, Markus Pögl, and Christoph Zenger

Full Text: Download PDF | Buy PDF (US$25) | View Cart
Competitive numerical algorithms for solving partial differential equations have to work with the most efficient numerical methods like multigrid and adaptive grid refinement and thus with hierarchical data structures. Unfortunately, in most implementations, hierarchical data—typically stored in trees—cause a nonnegligible overhead in data access. To overcome this quandary—numerical efficiency versus efficient implementation—our algorithm uses space‐filling curves to build up data structures which are processed linearly. In fact, the only kind of data structure used in our implementation is stacks. Thus, data access becomes very fast—even faster than the common access to nonhierarchical data stored in matrices—and, in particular, cache misses are reduced considerably. Furthermore, the implementation of multigrid cycles and/or higher order discretizations as well as the parallelization of the whole algorithm become very easy and straightforward on these data structures.

© 2006 Society for Industrial and Applied Mathematics

RELATED DATABASES

To view database links for this article, you need to log in.

PUBLICATION DATA

ISSN

1064-8275 (print)  
1095-7197 (online)

ARTICLE DATA

History
Received February 11, 2004
Accepted February 02, 2006
Published online October 06, 2006

For access to fully linked references, you need to log in.

Close

close