Heat kernel code by Benjamin Knorr