0
$\begingroup$

I have a 3D mesh composed of tessellating triangles. I would like to perform an integral over the surface of this mesh of the form:

$$\int f(\Omega)R(\Omega)d\Omega$$

Where $R(\Omega)$ is the radius at a given angle and $f(\Omega)$ is an arbitary function. I have a method implemented but it is giving large errors.

My method is as follows, is there anything obviously wrong about it?

For each triangular face:

  1. Calculate the value of $fR$ at each vertex of the face.
  2. Calculate the average of these three values to find the average value of $fR$ across the face.
  3. Multiply this value by the solid angle covered by the face.
  4. Sum up this result for all of the faces in the mesh.

I implemented this from scratch and I thought it should work. Has anyone done anything similar and if so am I overlooking something?

Thanks,

$\endgroup$
  • $\begingroup$ If it is 3D, then each element is a tetrahedron, and your integration is with respect to $d\Omega$ (volume). $\endgroup$ – Shuhao Cao Apr 27 '13 at 22:53
  • $\begingroup$ Sorry, I think I should clarify, I would like to integrate over the surface of this mesh. I'll edit the main question. $\endgroup$ – user668074 Apr 27 '13 at 23:10
  • $\begingroup$ For numerical quadrature, the sample points should be chosen on the surface. After you done the average, you have to multiply the area of each face on each integral. Also I am sure about what "solid angle covered by the face", can you make your statement more mathematical? For example using Cartesian coordinates. $\endgroup$ – Shuhao Cao Apr 27 '13 at 23:31
  • $\begingroup$ My sample points are at the corners of triangles which stitch together to make a closed surface. Those should be appropriate right? Regarding the solid angle part, I am approximating the different $d\Omega$ with the solid angle covered by the triangle face. Mathematically and in the program, I am using the veritices to create a spherical triangle. The solid angle covered by that spherical triangle is $(A+B+C) - \pi$. Where $A$, $B$ and $C$ are the angles between the vertices as measured on the surface. link $\endgroup$ – user668074 Apr 27 '13 at 23:37
  • $\begingroup$ I don't know if you're still looking for answers to this, but if so, what is $\Omega$? $\endgroup$ – Eric Kightley Jul 12 '13 at 19:14

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.