Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://www.adass.org/adass/proceedings/adass00/P1-13/
Дата изменения: Tue May 29 19:41:02 2001 Дата индексирования: Tue Oct 2 05:41:12 2012 Кодировка: Поисковые слова: stonehenge |
The wide-field imaging problem can be made more tractable by using parallelization. In this paper, we discuss the general wide-field imaging algorithm used in AIPS++, and the techniques used for its parallelization.
A problem occurs when imaging large fields of view with relatively long
baselines and non-coplanar arrays.
Imaging using synthesis arrays involves inverting the 3-D integral
Several algorithms exist to solve the full 3-D problem listed above (Cornwell & Perley 1992). In AIPS++ a multi-faceted transform approach has been chosen for its efficiency. This covers the image plane by a series of facets, in each of which a 2-D transform holds.
We can decompose the visibilities into a summation of re-phased faceted visibilities:
where :
The iterative multi-stage algorithm implemented in AIPS++ proceeds as follows:
Wide-field imaging is computationally expensive. The image in Figure 1 was made using the AIPS++ widefield algorithm with 225 facets. The data are a VLA observation at 74MHz in the B and C configurations. This image took close to 20days to process on a desktop workstation (SGI octane). A similar observation in the A-array of the VLA would require some ten times more computer resources to process. Along with other overheads, like better deconvolution algorithms for larger baselines, we are facing computation of 200 to 300 days on a typical desktop. This problem strongly justifies the need for parallelization of this algorithm. The problem will be more pronounced with future arrays such as the Expanded VLA (EVLA).
For the first level of parallelization we are aiming at parallelizing the nearly embarrassingly parallel sections of the widefield algorithms. There are three distinct sections which we have identified in the widefield algorithms which fall under this category:
We have made progress in parallel I/O development and evaluation of different access methods for the visibility data. This includes measuring the efficacy of parallelization with multiple processes accessing the same visibility data. We have verified that parallelization of the function to form the PSFs speeds up almost linearly for a few processors, and have parallelized the function to predict the model visibilities.
The parallel 3-D imaging approach is close to full operational use. Areas of ongoing work include: migration to larger machines or clusters and fully measuring the speed up of each algorithm computed, parallelizing the residual image formation for each facet. further work in Parallel I/O using MPI-2, and investigation of statement level parallelization using OPEN-MP.
Cornwell, T. J. & Perley, R. A. 1992 Astron. & Astrophys, 261, 353
Sault, R., Staveley-Smith, L., & Brouw, W. N. 1996 Astron. & Astrophys. Suppl., 120, 375