Coalescent simulation in continuous space:Algorithms for large neighbourhood size

Many species have an essentially continuous distribution in space, in which there are no natural divisions between randomly mating subpopulations. Yet, the standard approach to modelling these populations is to impose an arbitrary grid of demes, adjusting deme sizes and migration rates in an attempt to capture the important features of the population. Such indirect methods are required because of the failure of the classical models of isolation by distance, which have been shown to have major technical flaws. A recently introduced model of extinction and recolonisation in two dimensions solves these technical problems, and provides a rigorous technical foundation for the study of populations evolving in a spatial continuum. The coalescent process for this model is simply stated, but direct simulation is very inefficient for large neighbourhood sizes. We present efficient and exact algorithms to simulate this coalescent process for arbitrary sample sizes and numbers of loci, and analyse these algorithms in detail

Volumetric diffusers : pseudorandom cylinder arrays on a periodic lattice

Most conventional diffusers take the form of a surface based treatment, and as a result can only operate in hemispherical space. Placing a diffuser in the volume of a room might provide greater efficiency by allowing scattering into the whole space. A periodic cylinder array (or sonic crystal) produces periodicity lobes and uneven scattering. Introducing defects into an array, by removing or varying the size of some of the cylinders, can enhance their diffusing abilities. This paper applies number theoretic concepts to create cylinder arrays that have more even scattering. Predictions using a Boundary Element Method are compared to measurements to verify the model, and suitable metrics are adopted to evaluate performance. Arrangements with good aperiodic autocorrelation properties tend to produce the best results. At low frequency power is controlled by object size and at high frequency diffusion is dominated by lattice spacing and structural similarity. Consequently the operational bandwidth is rather small. By using sparse arrays and varying cylinder sizes, a wider bandwidth can be achieved