File Exchange

image thumbnail

Gauss Quadrature for Tetrahedra

version 1.0 (1.78 KB) by

Compute Gauss weights and nodes for a specied tetrahedron

10 Downloads

Updated

View License

This script computes the (x,y,z) Gauss quadrature nodes and weights for a tetrahedron with user-specied vertices.

Comments and Ratings (4)

yj2011

yj2011 (view profile)

Work wonderfully. One question. Not quite understand the Gauss Quadrature, but I noticed that the code by default outputs 64 note points. Can I increase the number of note points in order to improve the integration results?

saurabh gupta

John D'Errico

Its fixed. Thanks Greg.

John D'Errico

Wait until Greg fixes this. The following code fragment should compute the volume of a simple tetrahedron.
At the very least, since the function is constant over the
simplex domain, its integral should be positive.

V = [0 0 0;0 0 1;0 1 0;1 0 0];
[X,Y,Z,W]=tetraquad(4,V);
F = ones(size(X));
Q = W'*F

Q =
   -0.1667

The magnitude is correct, just the sign is wrong. Q should be 1/6.

MATLAB Release
MATLAB 7 (R14)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video