I see I'll have a look at that. Just curious how is the coded one different from a in place implementation? All I do is loop over mesh and average over directions.