Bernstein-Bézier bases for tetrahedral finite elements

Procedures are given for the efficient computation of the mass and stiffness matrices with these basis functions without using quadrature rules for (piecewise) constant coefficients on affine tetrahedra. Numerical results are presented to illustrate the use of the basis to approximate representative problems.