Re: [patch, fortran] Implement blocked eoshift for eoshift0

2017-07-02 Thread Paul Richard Thomas
Hi Thomas, The timings are impressive! OK for trunk. Thanks Paul On 1 July 2017 at 14:48, Thomas Koenig wrote: > Hello world, > > the attached patch implements the blocked algorithm for > constant shift for dim > 1 for eoshift0 (which handles > the case of constant shift and constant fill valu

[patch, fortran] Implement blocked eoshift for eoshift0

2017-07-01 Thread Thomas Koenig
Hello world, the attached patch implements the blocked algorithm for constant shift for dim > 1 for eoshift0 (which handles the case of constant shift and constant fill value). Speedup, as for cshift, is large. Moving a 500*500*500 array by -3 with eo_bench.f90 (also attached): $ gfortran -O3