Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Numerical integration on a triangle

  1. May 22, 2008 #1

    To solve an integral equation I have to numerically integrate the scalar Green's function on a triangle having vertices [tex](x,y) = (0,0), (1,0), (1,1)[/tex]:

    [tex]I = \int_0^1\int_0^x \frac{e^{-i 2\pi R}}{4\pi R} dy dx[/tex]



    is the distance between a fixed point outside the triangle and the integration point inside the triangle. Does anyone have any experience in this area?

    To evaluate this integral one can use precomputed weights and abcissae for triangular domains, and approximate the integral by a product rule, i.e.

    [tex]I=\sum_{n=1}^N w_i f(x_i,y_i)[/tex]

    where [tex]f(x,y)[/tex] is the Green's function. But this will only be accurate when [tex]R[/tex] is sufficiently large. When [tex]R[/tex] is small, the integrand grows increasingly singular.

    So I could use the product rule only when [tex]R[/tex] is larger than some constant, but I wonder if this is the smartest way of doing it. I guess this is a very typical problem encountered when solving any kind of real-world integral equation, but I have not been able to find any info about this problem.
  2. jcsd
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Can you offer guidance or do you also need help?
Draft saved Draft deleted