Computing heat diffusion creates weird results
I messed up I am doing a double negative and that is causing the problem:
const float laplace_coeff = -0.5f * ((1.f / tan(alpha)) + (1.f / tan(beta)));
This should be positive.