Measurement equation
A lot of this content was taken from Jack Line's
WODEN
.
The measurement equation (also known as the Radio Interferometer Measurement
Equation; RIME) used in hyperdrive
's calculations is:
\[ V_{s,f}(u_f,v_f,w_f) = \int\int S_{s,f}(l,m) e^{2 \pi i \phi} \frac{dl dm}{n} \]
where
- \( V_{s,f}(u_f,v_f,w_f) \) is the measured visibility in some Stokes polarisation \( s \) for some frequency \( f \) at baseline coordinates \( u_f, v_f, w_f \);
- \( S_{s,f} \) is the apparent brightness in the direction \( l, m \) at the same frequency \( f \);
- \( i \) is the imaginary unit;
- \( \phi = \left(u_fl + v_fm + w_f(n-1) \right) \); and
- \( n = \sqrt{1 - l^2 - m^2} \).
As we cannot ever have the true \( S_{s,f} \) function, we approximate with a sky-model source list, which details the expected positions and brightnesses of sources. This effectively turns the above continuous equation into a discrete one:
\[ V_{s,f}(u_f,v_f,w_f) = \sum S_{s,f}(l,m) e^{2 \pi i \left(u_fl + v_fm + w_f(n-1) \right)} \]
hyperdrive
implements this equation as code, either on the CPU or GPU
(preferred), and it is a good example of an embarrassingly parallel problem.