CPU times:
The mesh is made here of one single zones varying from 643 to 5123.
FastS runs at ~ 0.22 μs/pts/iter
i.e. for a 5123 cube, an adimensioned run time t*U0/L0=16 requires about 43.8 hours on 12 cores (x86 Xeon processors)
Results:
Animated isoQ surfaces:
 
Enstrophy, Kinetic energy and Kinetic energy dissipation rate compared with a very accurate spectral code.