This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
gibson:juliablog:kuramoto_sivashinksy [2017/07/09 06:23] gibson [The benchmark algorithm] |
gibson:juliablog:kuramoto_sivashinksy [2017/07/09 07:30] (current) gibson |
||
---|---|---|---|
Line 21: | Line 21: | ||
==== The results ==== | ==== The results ==== | ||
- | First, the results. The left plot shows execution time of 3200 time steps of the algorithm as a function of $N_x$, the number of gridpoints in the Fourier decomposition. The dominant cost of the algorithm should be the FFTs | + | |
+ | **The left plot** shows execution time of 3200 time steps of the algorithm as a function of $N_x$, the number of gridpoints in the Fourier decomposition. The dominant cost of the algorithm should be the FFTs, which should scale as $N_x \log N_x$. All the codes use the same FFTW libraries, so ideally, they should all collapse onto the same $N_x \log N_x$ line as linear and fixed-size overheads costs decrease relative to that. | ||
{{:gibson:juliablog:cputime.png?400|}} {{:gibson:juliablog:timeloc.png?400|}} | {{:gibson:juliablog:cputime.png?400|}} {{:gibson:juliablog:timeloc.png?400|}} | ||
The right plot shows | The right plot shows |