LoGPC: Modeling Network Contention in Message-Passing Programs

In many real applications, for example those with frequent and irregular communication patterns or those using large messages, network contention and contention for message processing resources can be a significant part of the total execution time. This paper presents a new cost model, called LoGPC, that extends the LogP [9] and LogGP [4] models to account for the impact of network contention and network interface DMA behavior on the performance of message-passing programs. We validate LoGPC by analyzing three applications implemented with Active Messages [11, 18] on the MIT Alewife multiprocessor. Our analysis shows that network contention accounts for up to 50% of the total execution time. In addition, we show that the impact of communication locality on the communication costs is at most a factor of two on Alewife. Finally, we use the model to identify tradeoffs between synchronous and asynchronous message passing styles. 1 Introduction Users of parallel machines need good performa..

Constraining local non-Gaussianities with kSZ tomography

Kinetic Sunyaev Zel'dovich (kSZ) tomography provides a powerful probe of the radial velocity field of matter in the Universe. By cross-correlating a high resolution CMB experiment like CMB S4 and a galaxy survey like DESI or LSST, one can measure the radial velocity field with very high signal to noise over a large volume of the universe. In this paper we show how this measurement can be used to improve constraints on primordial non-Gaussianities of the local type. The velocity field provides a measurement of the unbiased matter perturbations on large scales, which can be cross-correlated with the biased large-scale galaxy density field. This results in sample variance cancellation for a measurement of scale-dependent bias due to a non-zero $f_{NL}$. Using this method we forecast that CMB S4 and LSST combined reach a sensitivity $\sigma_{f_{NL}} \sim 0.5$, which is a factor of three improvement over the sensitivity using LSST alone (without internal sample variance cancellation). We take into account critical systematics like photometric redshifts, the kSZ optical depth degeneracy, and systematics affecting the shape of the galaxy auto-power spectrum and find that these have negligible impact, thus making kSZ tomography a robust probe for primordial non-Gaussianities. We also forecast the impact of mass binning on our constraints. The techniques proposed in this paper could be an important component of achieving the theoretically important threshold of $\sigma_{f_{NL}} \lesssim 1$ with future surveys.Comment: 16 pages, 7 figure

KSZ tomography and the bispectrum

Several statistics have been proposed for measuring the kSZ effect by combining the small-scale CMB with galaxy surveys. We review five such statistics, and show that they are all mathematically equivalent to the optimal bispectrum estimator of type $\langle ggT \rangle$. Reinterpreting these kSZ statistics as special cases of bispectrum estimation makes many aspects transparent, for example optimally weighting the estimator, or incorporating photometric redshift errors. We analyze the information content of the bispectrum and show that there are two observables: the small-scale galaxy-electron power spectrum $P_{ge}(k_S)$, and the large-scale galaxy-velocity power spectrum $P_{gv}(k)$. The cosmological constraining power of the kSZ arises from its sensitivity to fluctuations on large length scales, where its effective noise level can be much better than galaxy surveys.Comment: 39 page

Wall mode dynamics and transition to chaos in magnetoconvection with a vertical magnetic field

Quasistatic magnetoconvection of a low Prandtl number fluid ($\textrm{Pr}=0.025)$ with a vertical magnetic field is considered in a unit aspect ratio box with no-slip boundaries. At high relative magnetic field strengths, given by the Hartmann number $\textrm{Ha}$, the onset of convection is known to result from a sidewall instability giving rise to the wall mode regime. Here, we carry out 3D direct numerical simulations of unprecedented length to map out the parameter space at $\textrm{Ha} = 200, 500, 1000$, varying the Rayleigh number ($\textrm{Ra}$) between $6\times10^5 \lesssim \textrm{Ra} \lesssim 5\times 10^8$. We track the development of stable equilibria produced by this primary instability, identify bifurcations leading to limit cycles, and eventually to chaotic dynamics. At {$\textrm{Ha}=200$}, the steady wall mode solution undergoes a symmetry-breaking bifurcation producing a state featuring a coexistence between wall modes and a large-scale roll in the centre of the domain which persists to higher $\textrm{Ra}$. However, under a stronger magnetic field at $\textrm{Ha}=1000$, the steady wall mode solution undergoes a Hopf bifurcation producing a limit cycle which further develops to solutions that shadow an orbit homoclinic to a saddle point. Upon a further increase in $\textrm{Ra}$, the system undergoes a subsequent symmetry break producing a coexistence between wall modes and a large-scale roll, although the large-scale roll exists only for a small range of $\textrm{Ra}$, and chaotic dynamics primarily arise due to a mixture of chaotic wall mode dynamics and arrays of cellular structures