I've played around with the heat example provided in the tutorial.
However, when I increased the domain size from 8 by 16 to
const unsigned int vec_height = 1024;
const unsigned int vec_width = 1024;
const unsigned int vec_size = vec_height * vec_width;
and remove the output at the end, I get segmentation faults.