I have the structure now to calculate heat at time t+1 in fourier space for half of the rod
The runner file now evolves the rod as requested in the ouput, but somehow adds heat each time I reverse fourier transform.
I realize now that boundary conditions may have been the problem. I just added some bad ones and they greatly delayed the heat gain. Perhaps If I make some better ones they will completely destroy it.