  The presented paper concentrates on the boundary element method (BEM) for the
heat equation in three spatial dimensions. In particular, we deal with tensor
product space-time meshes allowing for quadrature schemes analytic in time and
numerical in space. The spatial integrals can be treated by standard BEM
techniques known from three dimensional stationary problems. The contribution
of the paper is twofold. First, we provide temporal antiderivatives of the heat
kernel necessary for the assembly of BEM matrices and the evaluation of the
representation formula. Secondly, the presented approach has been implemented
in a publicly available library besthea allowing researchers to reuse the
formulae and BEM routines straightaway. The results are validated by numerical
experiments in an HPC environment.
