Heat kernel code by Yannick Kluth