Fast Multipole Method as a Matrix-Free Hierarchical Low-Rank Approximation
نویسندگان
چکیده
There has been a large increase in the amount of work on hierarchical lowrank approximation methods, where the interest is shared by multiple communities that previously did not intersect. This objective of this article is two-fold; to provide a thorough review of the recent advancements in this field from both analytical and algebraic perspectives, and to present a comparative benchmark of two highly optimized implementations of contrasting methods for some simple yet representative test cases. We categorize the recent advances in this field from the perspective of compute-memory tradeoff, which has not been considered in much detail in this area. Benchmark tests reveal that there is a large difference in the memory consumption and performance between the different methods.
منابع مشابه
Communication Complexity of the Fast Multipole Method and its Algebraic Variants
A combination of hierarchical tree-like data structures and data access patterns from fast multipole methods and hierarchical low-rank approximation of linear operators from H-matrix methods appears to form an algorithmic path forward for efficient implementation of many linear algebraic operations of scientific computing at the exascale. The combination provides asymptotically optimal computat...
متن کاملIntroduction to Fast Multipole Methods
This is a simple introduction to fast multipole methods for the N-body summation problems. Low rank approximation plus hierarchical decomposition leads to fast O(N) or O(N logN) algorithms for the summation problem or equivalently the computation of a matrix-vector product. Code for 1-D problem is presented for a better illustration. We consider a summation problem in the form:
متن کاملA fast direct solver for the advection-diffusion equation using low-rank approximation
Today’s most powerful supercomputers make use of many thousands or millions of parallel processors to achieve high speeds. However, some of the core solution algorithms in computational physics codes are not well suited to this environment. Linear scaling is generally lost as the number of processors becomes very large. New solution algorithms specifically designed for massively parallel execut...
متن کاملPreprint ANL/MCS-P5112-0314 DATA STRUCTURE AND ALGORITHMS FOR RECURSIVELY LOW-RANK COMPRESSED MATRICES
We present a data structure and several operations it supports for recursively lowrank compressed matrices; that is, the diagonal blocks of the matrix are recursively partitioned, and the off-diagonal blocks in each partition level admit a low-rank approximation. Such a compression is embraced in many linearor near-linear-time methods for kernel matrices, including the fast multipole method, th...
متن کاملA Direct Solver for the Advection-diffusion Equation Using Green’s Functions and Low-rank Approximation
A new direct solution method for the advection-diffusion equation is presented. By employing a semi-implicit time discretisation, the equation is rewritten as a heat equation with source terms. The solution is obtained by discretely approximating the integral convolution of the associated Green’s function with advective source terms. The heat equation has an exponentially decaying Green’s funct...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- CoRR
دوره abs/1602.02244 شماره
صفحات -
تاریخ انتشار 2016