In radio astronomy obtaining a high dynamic range in synthesis imaging of wide fields requires a correction for time and direction-dependent effects. Applying direction-dependent correction can be done by either partitioning the image in facets and applying a direction-independent correction per facet, or by including the correction in the gridding kernel (AW-projection). An advantage of AW-projection over faceting is that the effectively applied beam is a sinc interpolation of the sampled beam, where the correction applied in the faceting approach is a discontinuous piece wise constant beam. However, AW-projection quickly becomes prohibitively expensive when the corrections vary over short time scales. This occurs for example when ionospheric effects are included in the correction. The cost of the frequent recomputation of the oversampled convolution kernels then dominates the total cost of gridding. Image domain gridding is a new approach that avoids the costly step of computing oversampled convolution kernels. Instead low-resolution images are made directly for small groups of visibilities which are then transformed and added to the large $uv$ grid. The computations have a simple, highly parallel structure that maps very well onto massively parallel hardware such as graphical processing units (GPUs). Despite being more expensive in pure computation count, the throughput is comparable to classical W-projection. The accuracy is close to classical gridding with a continuous convolution kernel. Compared to gridding methods that use a sampled convolution function, the new method is more accurate. Hence the new method is at least as fast and accurate as classical W-projection, while allowing for the correction for quickly varying direction-dependent effects.

We analyse the accuracy of radio interferometric gridding of visibilities with the aim to quantify the Epoch of Reionization (EoR) 21-cm power spectrum bias caused by gridding, ultimately to determine the suitability of different imaging algorithms and gridding settings for 21-cm power spectrum analysis. We simulate realistic LOFAR data, and construct power spectra with convolutional gridding and w-stacking, w-projection, image domain gridding and without w-correction. These are compared against directly Fourier transformed data. The influence of oversampling, kernel size, w-quantization, kernel windowing function and image padding are quantified. The gridding excess power is measured with a foreground subtraction strategy, for which foregrounds have been subtracted using Gaussian progress regression, as well as with a foreground avoidance strategy. Constructing a power spectrum that has a bias significantly lower compared to the expected EoR signals is possible with the tested methods, but requires a kernel oversampling factor > 4000 and, when using w-correction, > 500 w-quantization levels. These values are higher than typical values used for imaging, but are computationally feasible. The kernel size and padding factor parameters are less crucial. Among the tested methods, image domain gridding shows the highest accuracy with the lowest imaging time. LOFAR 21-cm power spectrum results are not affected by gridding. Image domain gridding is overall the most suitable algorithm for 21-cm EoR experiments, including for future SKA EoR analyses. Nevertheless, convolutional gridding with tuned parameters results in sufficient accuracy. This holds also for w-stacking for wide-field imaging. The w-projection algorithm is less suitable because of the kernel oversampling requirements, and a faceting approach is unsuitable due to the resulting spatial discontinuities.

We present a modular framework, the Workload Characterisation Framework (WCF), that is developed to reproducibly obtain, store and compare key characteristics of radio astronomy processing software. As a demonstration, we discuss the experiences using the framework to characterise a LOFAR calibration and imaging pipeline.