Multilevel evaluation of integral transforms with asymptotically smooth kernels

A. Brandt, C.H. Venner

    Research output: Contribution to journalArticleAcademicpeer-review

    21 Citations (Scopus)
    16 Downloads (Pure)

    Abstract

    Multilevel algorithms developed for the fast evaluation of integral transforms and the solution of the corresponding integral and integrodifferential equations rely on smoothness properties of the discrete kernel (matrix) and thereby on grid uniformity (see [A. Brandt and A.A. Lubrecht, J. Comput. Phys., 90 (1990), pp. 348--370], [C.H. Venner and A.A. Lubrecht, Multigrid Methods IV: Proc. 4th European Multigrid Conference, Amsterdam1993, P. Hemker and P. Wesseling, eds., Birkhäuser, Basel, 1994]). However, in actual applications, e.g., in contact mechanics, in many cases a substantial increase of efficiency can be obtained using nonuniform grids, since the solution is smooth in large parts of the domain with large gradients that occur only locally. In this paper a new algorithm is presented which relies on the smoothness of the continuum kernel only, independent of the grid configuration. This will facilitate the introduction of local refinements, wherever needed. Also, the evaluations will generally be faster; for a d-dimensional problem only O(sd + 1 ) operations per gridpoint are needed if s is the order of discretization. The algorithm is tested using a one-dimensional model problem with logarithmic kernel. Results are presented using both a second- and a fourth-order discretization. For testing purposes and to compare with results presented in [A. Brandt and A.A. Lubrecht, J. Comput. Phys., 90 (1990), pp. 348--370], uniform grids covering the entire domain were considered first.
    Original languageEnglish
    Pages (from-to)468-492
    JournalSIAM journal on scientific computing
    Volume19
    Issue number2
    DOIs
    Publication statusPublished - 1998

    Keywords

    • Multigrid
    • Integral transform
    • Singular smooth kernel
    • Fast evaluation
    • Local grid refinement

    Fingerprint Dive into the research topics of 'Multilevel evaluation of integral transforms with asymptotically smooth kernels'. Together they form a unique fingerprint.

  • Cite this