Mathematics – Numerical Analysis
Scientific paper
2011-05-26
Mathematics
Numerical Analysis
Scientific paper
An algorithm for the direct inversion of the linear systems arising from Nystrom discretization of integral equations on one-dimensional domains is described. The method typically has O(N) complexity when applied to boundary integral equations (BIEs) in the plane with non-oscillatory kernels such as those associated with the Laplace and Stokes' equations. The scaling coefficient suppressed by the "big-O" notation depends logarithmically on the requested accuracy. The method can also be applied to BIEs with oscillatory kernels such as those associated with the Helmholtz and Maxwell equations; it is efficient at long and intermediate wave-lengths, but will eventually become prohibitively slow as the wave-length decreases. To achieve linear complexity, rank deficiencies in the off-diagonal blocks of the coefficient matrix are exploited. The technique is conceptually related to the H and H^2 matrix arithmetic of Hackbusch and co-workers, and is closely related to previous work on Hierarchically Semi-Separable matrices.
Gillman Adrianna
Martinsson Per-Gunnar
Young Patrick
No associations
LandOfFree
A direct solver with O(N) complexity for integral equations on one-dimensional domains does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with A direct solver with O(N) complexity for integral equations on one-dimensional domains, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and A direct solver with O(N) complexity for integral equations on one-dimensional domains will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-218429