Efficient adaptive integration of functions with sharp gradients and cusps in n-dimensional parallelepipeds

Details Efficient adaptive integration of functions with sharp gradients and cusps in n-dimensional parallelepipeds Efficient adaptive integration of functions with sharp gradients and cusps in n-dimensional parallelepipeds

2012-02-23

arxiv.org/abs/1202.5341v1

Mathematics

Numerical Analysis

22 pages

Scientific paper

In this paper, we study the efficient numerical integration of functions with sharp gradients and cusps. An adaptive integration algorithm is presented that systematically improves the accuracy of the integration of a set of functions. The algorithm is based on a divide and conquer strategy and is independent of the location of the sharp gradient or cusp. The error analysis reveals that for a $C^0$ function (derivative-discontinuity at a point), a rate of convergence of $n+1$ is obtained in $R^n$. Two applications of the adaptive integration scheme are studied. First, we use the adaptive quadratures for the integration of the regularized Heaviside function---a strongly localized function that is used for modeling sharp gradients. Then, the adaptive quadratures are employed in the enriched finite element solution of the all-electron Coulomb problem in crystalline diamond. The source term and enrichment functions of this problem have sharp gradients and cusps at the nuclei. We show that the optimal rate of convergence is obtained with only a marginal increase in the number of integration points with respect to the pure finite element solution with the same number of elements. The adaptive integration scheme is simple, robust, and directly applicable to any generalized finite element method employing enrichments with sharp local variations or cusps in $n$-dimensional parallelepiped elements.

Affiliated with

Also associated with

No associations

Rating

Efficient adaptive integration of functions with sharp gradients and cusps in n-dimensional parallelepipeds 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 Efficient adaptive integration of functions with sharp gradients and cusps in n-dimensional parallelepipeds, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Efficient adaptive integration of functions with sharp gradients and cusps in n-dimensional parallelepipeds will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-78027