# Cubature of Singular Integrand

**URL:** <https://discourse.julialang.org/t/cubature-of-singular-integrand/113003>\
**Category:** Numerics\
**Tags:** question\
**Created:** [April 16, 2024, 11:27am UTC](https://discourse.julialang.org/t/cubature-of-singular-integrand/113003 "2024-04-16T11:27:08Z")\
**Posts on this page:** 1\
**Showing post:** 4

<div class="post-metadata">

**Author:** ![stevengj](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/stevengj/32/71_2.png) [@stevengj](https://discourse.julialang.org/u/stevengj)\
**Post date:** [April 16, 2024, 1:33pm UTC](https://discourse.julialang.org/t/cubature-of-singular-integrand/113003/4 "2024-04-16T13:33:06Z")

</div>

See also the recent thread: [Numerical integration over 3D domain of vector-valued function](https://discourse.julialang.org/t/numerical-integration-over-3d-domain-of-vector-valued-function/112284)

You can use HCubature.jl, Cubature.jl, nested `quadgk` calls (made more efficient by IteratedIntegration.jl), MeshIntegrals.jl, Cuba.jl, or explicit tensor products of fixed-order 1d quadrature rules (which has the tradeoff off not being adaptive, but for smooth integrands can be very efficient), for example.

As discussed in e.g. the linked thread, if you have a singularity in you integrand, then if at all possible you want to do a transformation that removes the singularity (or choose a quadrature rule that builds that singularity in).

---

_[View the full topic](https://discourse.julialang.org/t/cubature-of-singular-integrand/113003)._