Source author record

David Harvey

David Harvey appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2026arXiv

Deterministic methods for finding elements of large multiplicative order

We revisit the problem of rigorously and deterministically finding elements of large order in the multiplicative group of integers modulo a natural number $N$. Solving this problem is an essential step in several recent deterministic algorithms for factoring $N$, including the currently fastest ones. In 2018, the second author gave an algorithm that for a given target order $D \geq N^{2/5}$, finds either an element of order exceeding $D$, or a nontrivial divisor of $N$, or proves that $N$ is prime. The running time was \[ O\left(\frac{D^{1/2}}{(\log \log D)^{1/2}} \log^2 N \right) \] bit operations, asymptotically the same as the cost of computing the order of a single element using Sutherland's optimisation of the classical babystep-giantstep method. Subsequent work by several authors weakened the hypothesis $D \geq N^{2/5}$ to $D \geq N^{1/6}$. In this paper, we show that the hypothesis may be dropped altogether. Moreover, if $N$ is prime, we can guarantee returning an element of order exceeding $D$, rather than a proof that $N$ is prime.

preprint2021arXiv

From dwarf galaxies to galaxy clusters: Self-Interacting Dark Matter over 7 orders of magnitude in halo mass

In this paper we study the density profiles of self-interacting dark matter (SIDM) haloes spanning the full observable mass range, from dwarf galaxies to galaxy clusters. Using realistic simulations that model the baryonic physics relevant for galaxy formation, we compare the density profiles of haloes simulated with either SIDM or cold and collisionless dark matter (CDM) to those inferred from observations of stellar velocity dispersion, gas rotation curves, weak and strong gravitational lensing, and/or X-ray maps. We make our comparison in terms of the maximal surface density of haloes, circumventing the need for semi-analytic or parametric models for dark matter density profiles. We find that the maximal surface density as a function of halo mass is well reproduced by CDM simulations that include baryons, while for SIDM with a velocity-independent cross-section of 1 cm$^2/$g, the simulated galaxy clusters have mean maximal surface densities that are below those of observed systems by an amount greater than the standard deviation of the observed maximal surface density at fixed mass. For less massive systems both CDM and SIDM agree with the observation equally well.

preprint2020arXiv

Distributed Possibilistic Learning in Multi-Agent Systems

Possibility theory is proposed as an uncertainty representation framework for distributed learning in multi-agent systems and robot swarms. In particular, we investigate its application to the best-of-n problem where the aim is for a population of agents to identify the highest quality out of n options through local interactions between individuals and limited direct feedback from the environment. In this context we claim that possibility theory provides efficient mechanisms by which an agent can learn about the state of the world, and which can allow them to handle inconsistencies between what they and others believe by varying the level of imprecision of their own beliefs. We introduce a discrete time model of a population of agents applying possibility theory to the best-of-n problem. Simulation experiments are then used to investigate the accuracy of possibility theory in this context as well as its robustness to noise under varying amounts of direct evidence. Finally, we compare possibility theory in this context with a similar probabilistic approach.

preprint2020arXiv

Observable tests of self-interacting dark matter in galaxy clusters: cosmological simulations with SIDM and baryons

We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing.

preprint2020arXiv

The BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

preprint2020arXiv

The distribution of dark matter and gas spanning six megaparsecs around the post-merger galaxy cluster MS0451-03

Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions.

preprint2019arXiv

Exploiting flux ratio anomalies to probe warm dark matter in future large scale surveys

Flux ratio anomalies in strong gravitationally lensed quasars constitute a unique way to probe the abundance of non-luminous dark matter haloes, and hence the nature of dark matter. In this paper we identify double imaged quasars as a statistically efficient probe of dark matter, since they are 20 times more abundant than quadruply imaged quasars. Using N-body simulations that include realistic baryonic feedback, we measure the full distribution of flux ratios in doubly imaged quasars for cold (CDM) and warm dark matter (WDM) cosmologies. Through this method, we fold in two key systematics - quasar variability and line-of-sight structures. We find that WDM cosmologies predict a ~6 per cent difference in the cumulative distribution functions of flux ratios relative to CDM, with CDM predicting many more small ratios. Finally, we estimate that ~600 doubly imaged quasars will need to be observed in order to be able to unambiguously discern between CDM and the two WDM models studied here. Such sample sizes will be easily within reach of future large scale surveys such as Euclid. In preparation for this survey data we require discerning the scale of the uncertainties in modelling lens galaxies and their substructure in simulations, plus a strong understanding of the selection function of observed lensed quasars.

preprint2019arXiv

Setting the scene for BUFFALO: A study of the matter distribution in the HFF galaxy cluster MACS J0416.1-2403 and its parallel field

In the context of the BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations) survey, we present a new analysis of the merging galaxy cluster MACS\,J0416.1-2403 ($z = 0.397$) and its parallel field using the data collected by the Hubble Frontier Fields (HFF) campaign. In this work, we measure the surface mass density from a weak-lensing analysis, and characterise the overall matter distribution in both the cluster and parallel fields. The surface mass distribution derived for the parallel field shows clumpy overdensities connected by filament-like structures elongated in the direction of the cluster core. We also characterise the X-ray emission of the cluster, and compare it with the lensing mass distribution. We identify five substructures at the $>5σ$ level over the two fields, four of them being in the cluster one. Furthermore, three of them are located close to the edges of the field of view, and border issues can significantly hamper the determination of their physical parameters. Finally, we compare our results with the predicted subhalo distribution of one of the Hydrangea/C-EAGLE simulated cluster. Significant differences are obtained suggesting the simulated cluster is at a more advanced evolutionary state than MACS\,J0416.1-2403. Our results anticipate the upcoming BUFFALO observations that will link the two HFF fields, extending further the \emph{HST} coverage, and thus allowing a better characterisation of the reported substructures.

preprint2016arXiv

Computing L-series of geometrically hyperelliptic curves of genus three

Let C/Q be a curve of genus three, given as a double cover of a plane conic. Such a curve is hyperelliptic over the algebraic closure of Q, but may not have a hyperelliptic model of the usual form over Q. We describe an algorithm that computes the local zeta functions of C at all odd primes of good reduction up to a prescribed bound N. The algorithm relies on an adaptation of the "accumulating remainder tree" to matrices with entries in a quadratic field. We report on an implementation, and compare its performance to previous algorithms for the ordinary hyperelliptic case.

preprint2016arXiv

Irregular primes to two billion

We compute all irregular primes less than 2^31 = 2 147 483 648. We verify the Kummer-Vandiver conjecture for each of these primes, and we check that the p-part of the class group of Q(zeta_p) has the simplest possible structure consistent with the index of irregularity of p. Our method for computing the irregular indices saves a constant factor in time relative to previous methods, by adapting Rader's algorithm for evaluating discrete Fourier transforms.

preprint2016arXiv

Systematic or Signal? How dark matter misalignments can bias strong lensing models of galaxy clusters

We explore how assuming that mass traces light in strong gravitational lensing models can lead to systematic errors in the predicted position of multiple images. Using a model based on the galaxy cluster MACSJ0416 (z = 0.397) from the Hubble Frontier Fields, we split each galactic halo into a baryonic and dark matter component. We then shift the dark matter halo such that it no longer aligns with the baryonic halo and investigate how this affects the resulting position of multiple images. We find for physically motivated misalignments in dark halo position, ellipticity, position angle and density profile, that multiple images can move on average by more than 0.2" with individual images moving greater than 1". We finally estimate the full error induced by assuming that light traces mass and find that this assumption leads to an expected RMS error of 0.5", almost the entire error budget observed in the Frontier Fields. Given the large potential contribution from the assumption that light traces mass to the error budget in mass reconstructions, we predict that it should be possible to make a first significant detection and characterisation of dark halo misalignments in the Hubble Frontier Fields with strong lensing. Finally, we find that it may be possible to detect ~1kpc offsets between dark matter and baryons, the smoking gun for self-interacting dark matter, should the correct alignment of multiple images be observed.

preprint2015arXiv

A weak lensing comparability study of galaxy mergers that host AGNs

We compared the total mass density profiles of three different types of galaxies using weak gravitational lensing: (i) 29 galaxies that host quasars at z~0.32 that are in a post-starburst (PSQ) phase with high star formation indicating recent merger activity, (ii) 22 large elliptical galaxies from the SLACS sample that do not host a quasar at z~0.23, and (iii) 17 galaxies that host moderately luminous quasars at z~0.36 powered by disk instabilities, but with no intense star formation. On an initial test we found no evidence for a connection between the merger state of a galaxy and the profile of the halo, with the PSQ profile comparable to that of the other two samples and consistent with the Leauthaud et al. (2014) study of moderately luminous quasars in COSMOS. Given the compatibility of the two quasar samples, we combined these and found no evidence for any connection between black hole activity and the dark matter halo. All three mass profiles remained compatible with isothermality given the present data.

preprint2015arXiv

Computing zeta functions of arithmetic schemes

We present new algorithms for computing zeta functions of algebraic varieties over finite fields. In particular, let X be an arithmetic scheme (scheme of finite type over Z), and for a prime p let zeta_{X_p}(s) be the local factor of its zeta function. We present an algorithm that computes zeta_{X_p}(s) for a single prime p in time p^(1/2+o(1)), and another algorithm that computes zeta_{X_p}(s) for all primes p < N in time N (log N)^(3+o(1)). These generalise previous results of the author from hyperelliptic curves to completely arbitrary varieties.

preprint2015arXiv

The behaviour of dark matter associated with 4 bright cluster galaxies in the 10kpc core of Abell 3827

Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright elliptical galaxies within a core of radius 10kpc. Such corrugation of the stellar distribution is very rare, and suggests recent formation by several simultaneous mergers. We map the distribution of associated dark matter, using new Hubble Space Telescope imaging and VLT/MUSE integral field spectroscopy of a gravitationally lensed system threaded through the cluster core. We find that each of the central galaxies retains a dark matter halo, but that (at least) one of these is spatially offset from its stars. The best-constrained offset is 1.62+/-0.48kpc, where the 68% confidence limit includes both statistical error and systematic biases in mass modelling. Such offsets are not seen in field galaxies, but are predicted during the long infall to a cluster, if dark matter self-interactions generate an extra drag force. With such a small physical separation, it is difficult to definitively rule out astrophysical effects operating exclusively in dense cluster core environments - but if interpreted solely as evidence for self-interacting dark matter, this offset implies a cross-section sigma/m=(1.7+/-0.7)x10^{-4}cm^2/g x (t/10^9yrs)^{-2}, where t is the infall duration.

preprint2015arXiv

The non-gravitational interactions of dark matter in colliding galaxy clusters

Collisions between galaxy clusters provide a test of the non-gravitational forces acting on dark matter. Dark matter's lack of deceleration in the `bullet cluster collision' constrained its self-interaction cross-section σ_DM/m < 1.25cm2/g (68% confidence limit) for long-ranged forces. Using the Chandra and Hubble Space Telescopes we have now observed 72 collisions, including both `major' and `minor' mergers. Combining these measurements statistically, we detect the existence of dark mass at 7.6σsignificance. The position of the dark mass has remained closely aligned within 5.8+/-8.2 kpc of associated stars: implying a self-interaction cross-section σ_DM/m < 0.47 cm2/g (95% CL) and disfavoring some proposed extensions to the standard model.

preprint2014arXiv

Computing Hasse-Witt matrices of hyperelliptic curves in average polynomial time

We present an efficient algorithm to compute the Hasse-Witt matrix of a hyperelliptic curve C/Q modulo all primes of good reduction up to a given bound N, based on the average polynomial-time algorithm recently introduced by Harvey. An implementation for hyperelliptic curves of genus 2 and 3 is more than an order of magnitude faster than alternative methods for N = 2^26.

preprint2014arXiv

Even faster integer multiplication

We give a new proof of Fürer's bound for the cost of multiplying n-bit integers in the bit complexity model. Unlike Fürer, our method does not require constructing special coefficient rings with "fast" roots of unity. Moreover, we prove the more explicit bound O(n log n K^(log^* n))$ with K = 8. We show that an optimised variant of Fürer's algorithm achieves only K = 16, suggesting that the new algorithm is faster than Fürer's by a factor of 2^(log^* n). Assuming standard conjectures about the distribution of Mersenne primes, we give yet another algorithm that achieves K = 4.

preprint2014arXiv

Faster polynomial multiplication over finite fields

Let p be a prime, and let M_p(n) denote the bit complexity of multiplying two polynomials in F_p[X] of degree less than n. For n large compared to p, we establish the bound M_p(n) = O(n log n 8^(log^* n) log p), where log^* is the iterated logarithm. This is the first known Fürer-type complexity bound for F_p[X], and improves on the previously best known bound M_p(n) = O(n log n log log n log p).

preprint2014arXiv

Hubble Frontier Fields: The Geometry and Dynamics of the Massive Galaxy Cluster Merger MACSJ0416.1-2403

We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dynamics, we complement this lensing analysis with a study of the intra-cluster gas using shallow Chandra data, and a three-dimensional model of the distribution and motions of cluster galaxies derived from over 100 spectroscopic redshifts. The multi-scale grid model obtained from our combined lensing analysis extends the high-precision strong-lensing mass reconstruction recently performed to cluster-centric distances of almost 1 Mpc. Our analysis detects the two well known mass concentrations in the cluster core. A pronounced offset between collisional and collisionless matter is only observed for the SW cluster component, while excellent alignment is found for the NE cluster. Both the lensing analysis and the distribution of cluster light strongly suggest the presence of a third massive structure, almost 2 arcmin SW of the cluster centre. Since no X-ray emission is detected in this region, we conclude that this structure is non-virialised and speculate that it might be part of a large-scale filament almost aligned with our line of sight. Combining all evidence from the distribution of dark and luminous matter, we propose two alternative scenarios for the trajectories of the components of MACSJ0416.1-2403. Upcoming deep X-ray observations that allow the detection of shock fronts, cold cores, and sloshing gas (all key diagnostics for studies of cluster collisions) will allow us to test, and distinguish between these two scenarios.

preprint2013arXiv

Counting points on hyperelliptic curves in average polynomial time

Let g >= 1 and let Q be a monic, squarefree polynomial of degree 2g + 1 in Z[x]. For an odd prime p not dividing the discriminant of Q, let Z_p(T) denote the zeta function of the hyperelliptic curve of genus g over the finite field F_p obtained by reducing the coefficients of the equation y^2 = Q(x) modulo p. We present an explicit deterministic algorithm that given as input Q and a positive integer N, computes Z_p(T) simultaneously for all such primes p < N, whose average complexity per prime is polynomial in g, log N, and the number of bits required to represent Q.

preprint2013arXiv

Dark Matter Astrometry: Accuracy of sub-halo positions for the measurement of self-interaction cross sections

Direct evidence for the existence of dark matter and measurements of its interaction cross-section have been provided by the physical offset between dark matter and intra- cluster gas in merging systems like the Bullet Cluster. Although a smaller signal, this effect is more abundant in minor mergers where infalling substructure dark matter and gas are segregated. In such low-mass systems the gravitational lensing signal comes primarily from weak lensing. A fundamental step in determining such an offset in sub- structure is the ability to accurately measure the positions of dark matter sub-peaks. Using simulated Hubble Space Telescope observations, we make a first assessment of the precision and accuracy with which we can measure infalling groups using weak gravitational lensing. We demonstrate that using an existing and well-used mass re- construction algorithm can measure the positions of 1.5x1013 M substructures that have parent halos ten times more massive with a bias of less than 0.3 . In this regime, our analysis suggests the precision is sufficient to detect (at 3 σ statistical significance) the expected mean offset between dark matter and baryonic gas in infalling groups from a sample of 50 massive clusters.

preprint2013arXiv

Faster arithmetic for number-theoretic transforms

We show how to improve the efficiency of the computation of fast Fourier transforms over F_p where p is a word-sized prime. Our main technique is optimisation of the basic arithmetic, in effect decreasing the total number of reductions modulo p, by making use of a redundant representation for integers modulo p. We give performance results showing a significant improvement over Shoup's NTL library.

preprint2013arXiv

Observing Dark Worlds: A crowdsourcing experiment for dark matter mapping

We present the results and conclusions from the citizen science competition `Observing Dark Worlds', where we asked participants to calculate the positions of dark matter halos from 120 catalogues of simulated weak lensing galaxy data, using computational methods. In partnership with Kaggle (http://www.kaggle.com), 357 users participated in the competition which saw 2278 downloads of the data and 3358 submissions. We found that the best algorithms improved on the benchmark code, LENSTOOL by > 30% and could measure the positions of > 3x10^14MSun halos to less than 5'' and < 10^14MSun to within 1'. In this paper, we present a brief overview of the winning algorithms with links to available code. We also discuss the implications of the experiment for future citizen science competitions.

preprint2013arXiv

On the cross-section of Dark Matter using substructure infall into galaxy clusters

We develop a statistical method to measure the interaction cross-section of Dark Matter, exploiting the continuous minor merger events in which small substructures fall into galaxy clusters. We find that by taking the ratio of the distances between the galaxies and Dark Matter, and galaxies and gas in accreting sub-halos, we form a quantity that can be statistically averaged over a large sample of systems whilst removing any inherent line-of-sight projections. In order to interpret this ratio as a cross-section of Dark Matter we derive an analytical description of sub-halo infall which encompasses; the force of the main cluster potential, the drag on a gas sub-halo, a model for Dark Matter self-interactions and the resulting sub-halo drag, the force on the gas and galaxies due to the Dark Matter sub-halo potential, and finally the buoyancy on the gas and Dark Matter. We create mock observations from cosmological simulations of structure formation and find that collisionless Dark Matter becomes physically separated from X-ray gas by up to 20h^-1 kpc. Adding realistic levels of noise, we are able to predict achievable constraints from observational data. Current archival data should be able to detect a difference in the dynamical behaviour of Dark Matter and standard model particles at 6 sigma, and measure the total interaction cross-section sigma/m with 68% confidence limits of +/- 1cm2g^-1. We note that this method is not restricted by the limited number of major merging events and is easily extended to large samples of clusters from future surveys which could potentially push statistical errors to 0.1cm^2g^-1.

preprint2012arXiv

BK Lyncis: The Oldest Old Nova?... And a Bellwether for Cataclysmic-Variable Evolution

We summarize the results of a 20-year campaign to study the light curves of BK Lyncis, a nova-like star strangely located below the 2-3 hour orbital period gap in the family of cataclysmic variables. Two apparent "superhumps" dominate the nightly light curves - with periods 4.6% longer, and 3.0% shorter, than P_orb. The first appears to be associated with the star's brighter states (V~14), while the second appears to be present throughout and becomes very dominant in the low state (V~15.7). Starting in the year 2005, the star's light curve became indistinguishable from that of a dwarf nova - in particular, that of the ER UMa subclass. Reviewing all the star's oddities, we speculate: (a) BK Lyn is the remnant of the probable nova on 30 December 101, and (b) it has been fading ever since, but has taken ~2000 years for the accretion rate to drop sufficiently to permit dwarf-nova eruptions. If such behavior is common, it can explain other puzzles of CV evolution. One: why the ER UMa class even exists (because all members can be remnants of recent novae). Two: why ER UMa stars and short-period novalikes are rare (because their lifetimes, which are essentially cooling times, are short). Three: why short-period novae all decline to luminosity states far above their true quiescence (because they're just getting started in their postnova cooling). Four: why the orbital periods, accretion rates, and white-dwarf temperatures of short-period CVs are somewhat too large to arise purely from the effects of gravitational radiation (because the unexpectedly long interval of enhanced postnova brightness boosts the mean mass-transfer rate). These are substantial rewards in return for one investment of hypothesis: that the second parameter in CV evolution, besides P_orb, is time since the last classical-nova eruption.

preprint2012arXiv

Origins of weak lensing systematics, and requirements on future instrumentation (or knowledge of instrumentation)

The first half of this paper explores the origin of systematic biases in the measurement of weak gravitational lensing. Compared to previous work, we expand the investigation of PSF instability and fold in for the first time the effects of non-idealities in electronic imaging detectors and imperfect galaxy shape measurement algorithms. Together, these now explain the additive A(l) and multiplicative M(l) systematics typically reported in current lensing measurements. We find that overall performance is driven by a product of a telescope/camera's *absolute performance*, and our *knowledge about its performance*. The second half of this paper propagates any residual shear measurement biases through to their effect on cosmological parameter constraints. Fully exploiting the statistical power of Stage IV weak lensing surveys will require additive biases A<1.8e-12 and multiplicative biases M<4.0e-3. These can be allocated between individual budgets in hardware, calibration data and software, using results from the first half of the paper. If instrumentation is stable and well-calibrated, we find extant shear measurement software from GREAT10 already meet requirements on galaxies detected at S/N=40. Averaging over a population of galaxies with a realistic distribution of sizes, it also meets requirements for a 2D cosmic shear analysis from space. If used on fainter galaxies or for 3D cosmic shear tomography, existing algorithms would need calibration on simulations to avoid introducing bias at a level similar to the statistical error. Requirements on hardware and calibration data are discussed in more detail in a companion paper. Our analysis is intentionally general, but is specifically being used to drive the hardware and ground segment performance budget for the design of the European Space Agency's recently-selected Euclid mission.

preprint2010arXiv

An in-place truncated Fourier transform and applications to polynomial multiplication

The truncated Fourier transform (TFT) was introduced by van der Hoeven in 2004 as a means of smoothing the "jumps" in running time of the ordinary FFT algorithm that occur at power-of-two input sizes. However, the TFT still introduces these jumps in memory usage. We describe in-place variants of the forward and inverse TFT algorithms, achieving time complexity O(n log n) with only O(1) auxiliary space. As an application, we extend the second author's results on space-restricted FFT-based polynomial multiplication to polynomials of arbitrary degree.

preprint2010arXiv

Characterizing projective spaces on deformations of Hilbert schemes of K3 surfaces

We seek to characterize homology classes of Lagrangian projective spaces embedded in irreducible holomorphic-symplectic manifolds, up to the action of the monodromy group. This paper addresses the case of manifolds deformation-equivalent to the Hilbert scheme of length-three subschemes of a K3 surface. The class of the projective space in the cohomology ring has prescribed intersection properties, which translate into Diophantine equations. Possible homology classes correspond to integral points on an explicit elliptic curve; our proof entails showing that the only such point is two-torsion.