EFFICIENT CONSTRUCTION OF 2-CHAINS WITH A PRESCRIBED BOUNDARY

Let Omega be a bounded domain of R-3 whose closure (Omega) over bar is polyhedral, and let T be a triangulation of (Omega) over bar. We devise a fast algorithm for the computation of homological Seifert surfaces of any 1-boundary of T, namely, 2-chains of T whose boundary is gamma. Assuming that the boundary of Omega is sufficiently regular, we provide an explicit formula for a homological Seifert surface of any 1-boundary 7 of T. It is based on the existence of special spanning trees of the complete dual graph and on the computation of certain linking numbers associated with those spanning trees. If the triangulation T is fine, the explicit formula is too expensive to be used directly. To overcome this difficulty, we adopt an easy and very fast elimination procedure, which sometimes fails. In such a case a new unknown can be computed using the explicit formula and the elimination algorithm restarts. The numerical experiments we performed illustrate the efficiency of the resulting algorithm even when the homology of Omega is not trivial and the triangulation T of Omega consists of millions of tetrahedra.