FourierFlows / FourierFlows.jl

Tools for building fast, hackable, pseudospectral partial differential equation solvers on periodic domains
https://bit.ly/FourierFlows
MIT License
207 stars 29 forks source link

Use FastBroadcast for multithreaded broadcasting? #337

Open glwagner opened 1 year ago

glwagner commented 1 year ago

For some time we've lamented that while FFTW is multithreaded, broadcasting operations such as

https://github.com/FourierFlows/FourierFlows.jl/blob/e059ba96260f483b493659fcefb343a08f0db9ae/src/timesteppers.jl#L109

are not.

FastBroadcast.jl provides an answer:

https://github.com/YingboMa/FastBroadcast.jl

It seems that even without multithreading, it can speed things up.