On the Ludwig integration algorithm for triangular subregions

By using area coordinates one obtains algebraic expressions for the computation of two-dimensional integrals with oscillatory integrands over triangular subregions that are simpler than those published recently [1]. All the advantages of previous methods [1], [2] are retained. Moreover, the new equation is evaluated only once for each triangular subregion.