Source author record

Christopher J. Fluke

Christopher J. Fluke 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

22works
5topics
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

22 published item(s)

preprint2021arXiv

Shape noise and dispersion in precision weak lensing

We analyse the first measurements from precision weak lensing (PWL): a new methodology for measuring individual galaxy-galaxy weak lensing through velocity information. Our goal is to understand the observed shear distribution from PWL, which is broader than can be explained by the statistical measurement errors. We identify two possible sources of scatter to explain the observed distribution: a shape noise term associated with the underlying assumption of circular stable rotation, and an astrophysical signal consistent with a log-normal dispersion around the stellar-to-halo mass relation (SHMR). We have modelled the observed distribution as the combination of these two factors and quantified their most likely values given our data. For the current sample, we measure an effective shape noise of $σ_γ= 0.024 \pm 0.007$, highlighting the low noise impact of the method and positioning PWL as $\sim 10$ times more precise than conventional weak lensing. We also measure an average dispersion in shears of $ξ_γ= 0.53^{+0.26}_{-0.28}$\,dex over the range of $8.5 < \log M_\star < 11$. This measurement is higher than expected, which is suggestive of a relatively high dispersion in halo mass and/or profile.

preprint2020arXiv

The first shear measurements from precision weak lensing

We present an end-to-end methodology to measure the effects of weak lensing on individual galaxy-galaxy systems exploiting their kinematic information. Using this methodology, we have measured a shear signal from the velocity fields of 18 weakly-lensed galaxies. We selected a sample of systems based only on the properties of the sources, requiring them to be bright (apparent $i$-band magnitude $ < 17.4$) and in the nearby Universe ($z < 0.15$). We have observed the velocity fields of the sources with WiFeS, an optical IFU on a 2.3m telescope, and fitted them using a simple circular motion model with an external shear. We have measured an average shear of $\langle γ\rangle = 0.020 \pm 0.008$ compared to a predicted $\langle γ_{pred} \rangle = 0.005$ obtained using median stellar-to-halo relationships from the literature. While still a statistical approach, our results suggest that this new weak lensing methodology can overcome some of the limitations of traditional stacking-based techniques. We describe in detail all the steps of the methodology and make publicly available all the velocity maps for the weakly-lensed sources used in this study.

preprint2020arXiv

The influence of angular momentum and environment on the HI gas of late-type galaxies

We use high-resolution HI data from the WHISP survey to study the HI and angular momentum properties of a sample of 114 late-type galaxies. We explore the specific baryonic angular momentum -- baryonic mass ($j_b - M_b$) relation, and find that an unbroken power law of the form $j_b \propto M_b^{0.55 \pm 0.02}$ fits the data well, with an intrinsic scatter of $\sim 0.13 \pm 0.01$ dex. We revisit the relation between the atomic gas fraction, $f_{atm}$, and the integrated atomic stability parameter $q$ (the $f_{atm} - q$ relation), originally introduced by Obreschkow et al., and probe this parameter space by populating it with galaxies from different environments, in order to study the influence of the environment on their $j_b$, $f_{atm}$ and $q$ values. We find evidence that galaxies with close neighbours show a larger intrinsic scatter about the $f_{atm} - q$ relation compared to galaxies without close-neighbours. We also find enhanced SFR among the deviating galaxies with close neighbours. In addition, we use the bulge-to-total (B/T) ratio as a morphology proxy, and find a general trend of decreasing B/T values with increasing disc stability and HI fraction in the $f_{atm} - q$ plane, indicating a fundamental link between mass, specific angular momentum, gas fraction and morphology of galaxies.

preprint2020arXiv

Understanding the human in the design of cyber-human discovery systems for data-driven astronomy

High-quality, usable, and effective software is essential for supporting astronomers in the discovery-focused tasks of data analysis and visualisation. As the volume, and perhaps more crucially, the velocity of astronomical data grows, the role of the astronomer is changing. There is now an increased reliance on automated and autonomous discovery and decision-making workflows rather than visual inspection. We assert the need for an improved understanding of how astronomers (humans) currently make visual discoveries from data. This insight is a critical element for the future design, development and effective use of cyber-human discovery systems, where astronomers work in close collaboration with automated systems to gain understanding from continuous, real-time data streams. We discuss how relevant human performance data could be gathered, specifically targeting the domains of expertise and skill at visual discovery, and the identification and management of cognitive factors. By looking to other disciplines where human performance is assessed and measured, we propose four early-stage applications that would: (1) allow astronomers to evaluate, and potentially improve, their own visual discovery skills; (2) support just-in-time coaching; (3) enable talent identification; and (4) result in user interfaces that automatically respond to skill level and cognitive state. Throughout, we advocate for the importance of user studies and the incorporation of participatory design and co-design practices into the planning, implementation and evaluation of alternative user interfaces and visual discovery environments.

preprint2016arXiv

Data Compression in the Petascale Astronomy Era: a GERLUMPH case study

As the volume of data grows, astronomers are increasingly faced with choices on what data to keep -- and what to throw away. Recent work evaluating the JPEG2000 (ISO/IEC 15444) standards as a future data format standard in astronomy has shown promising results on observational data. However, there is still a need to evaluate its potential on other type of astronomical data, such as from numerical simulations. GERLUMPH (the GPU-Enabled High Resolution cosmological MicroLensing parameter survey) represents an example of a data intensive project in theoretical astrophysics. In the next phase of processing, the ~27 terabyte GERLUMPH dataset is set to grow by a factor of 100 -- well beyond the current storage capabilities of the supercomputing facility on which it resides. In order to minimise bandwidth usage, file transfer time, and storage space, this work evaluates several data compression techniques. Specifically, we investigate off-the-shelf and custom lossless compression algorithms as well as the lossy JPEG2000 compression format. Results of lossless compression algorithms on GERLUMPH data products show small compression ratios (1.35:1 to 4.69:1 of input file size) varying with the nature of the input data. Our results suggest that JPEG2000 could be suitable for other numerical datasets stored as gridded data or volumetric data. When approaching lossy data compression, one should keep in mind the intended purposes of the data to be compressed, and evaluate the effect of the loss on future analysis. In our case study, lossy compression and a high compression ratio do not significantly compromise the intended use of the data for constraining quasar source profiles from cosmological microlensing.

preprint2016arXiv

Large-scale comparative visualisation of sets of multidimensional data

We present encube $-$ a qualitative, quantitative and comparative visualisation and analysis system, with application to high-resolution, immersive three-dimensional environments and desktop displays. encube extends previous comparative visualisation systems by considering: 1) the integration of comparative visualisation and analysis into a unified system; 2) the documentation of the discovery process; and 3) an approach that enables scientists to continue the research process once back at their desktop. Our solution enables tablets, smartphones or laptops to be used as interaction units for manipulating, organising, and querying data. We highlight the modularity of encube, allowing additional functionalities to be included as required. Additionally, our approach supports a high level of collaboration within the physical environment. We show how our implementation of encube operates in a large-scale, hybrid visualisation and supercomputing environment using the CAVE2 at Monash University, and on a local desktop, making it a versatile solution. We discuss how our approach can help accelerate the discovery rate in a variety of research scenarios.

preprint2015arXiv

GERLUMPH Data Release 2: 2.5 billion simulated microlensing light curves

In the upcoming synoptic all--sky survey era of astronomy, thousands of new multiply imaged quasars are expected to be discovered and monitored regularly. Light curves from the images of gravitationally lensed quasars are further affected by superimposed variability due to microlensing. In order to disentangle the microlensing from the intrinsic variability of the light curves, the time delays between the multiple images have to be accurately measured. The resulting microlensing light curves can then be analyzed to reveal information about the background source, such as the size of the quasar accretion disc. In this paper we present the most extensive and coherent collection of simulated microlensing light curves; we have generated $>2.5$ billion light curves using the GERLUMPH high resolution microlensing magnification maps. Our simulations can be used to: train algorithms to measure lensed quasar time delays, plan future monitoring campaigns, and study light curve properties throughout parameter space. Our data are openly available to the community and are complemented by online eResearch tools, located at http://gerlumph.swin.edu.au .

preprint2015arXiv

Kinematic Modelling of Disc Galaxies using Graphics Processing Units

With large-scale Integral Field Spectroscopy (IFS) surveys of thousands of galaxies currently under-way or planned, the astronomical community is in need of methods, techniques and tools that will allow the analysis of huge amounts of data. We focus on the kinematic modelling of disc galaxies and investigate the potential use of massively parallel architectures, such as the Graphics Processing Unit (GPU), as an accelerator for the computationally expensive model-fitting procedure. We review the algorithms involved in model-fitting and evaluate their suitability for GPU implementation. We employ different optimization techniques, including the Levenberg-Marquardt and Nested Sampling algorithms, but also a naive brute-force approach based on Nested Grids. We find that the GPU can accelerate the model-fitting procedure up to a factor of ~100 when compared to a single-threaded CPU, and up to a factor of ~10 when compared to a multi-threaded dual CPU configuration. Our method's accuracy, precision and robustness are assessed by successfully recovering the kinematic properties of simulated data, and also by verifying the kinematic modelling results of galaxies from the GHASP and DYNAMO surveys as found in the literature. The resulting GBKFIT code is available for download from: http://supercomputing.swin.edu.au/gbkfit.

preprint2014arXiv

The 6dF Galaxy Survey: Peculiar Velocity Field and Cosmography

We derive peculiar velocities for the 6dF Galaxy Survey (6dFGS) and describe the velocity field of the nearby ($z<0.055$) southern hemisphere. The survey comprises 8885 galaxies for which we have previously reported Fundamental Plane data. We obtain peculiar velocity probability distributions for the redshift space positions of each of these galaxies using a Bayesian approach. Accounting for selection bias, we find that the logarithmic distance uncertainty is 0.11 dex, corresponding to $26\%$ in linear distance. We use adaptive kernel smoothing to map the observed 6dFGS velocity field out to $cz \sim 16,000$ \kms, and compare this to the predicted velocity fields from the PSCz Survey and the 2MASS Redshift Survey. We find a better fit to the PSCz prediction, although the reduced $χ^2$ for the whole sample is approximately unity for both comparisons. This means that, within the observational uncertainties due to redshift independent distance errors, observed galaxy velocities and those predicted by the linear approximation from the density field agree. However, we find peculiar velocities that are systematically more positive than model predictions in the direction of the Shapley and Vela superclusters, and systematically more negative than model predictions in the direction of the Pisces-Cetus Supercluster, suggesting contributions from volumes not covered by the models.

preprint2014arXiv

The effect of macromodel uncertainties on microlensing modelling of lensed quasars

Cosmological gravitational microlensing has been proven to be a powerful tool to constrain the structure of multiply imaged quasars, especially the accretion disc and central supermassive black-hole system. However, the derived constraints on models may be affected by large systematic errors introduced in the various stages of modelling, namely, the macromodels, the microlensing magnification maps, and the convolution with realistic disc profiles. In particular, it has been known that different macromodels of the galaxy lens that fit the observations equally well, can lead to different values of convergence, $κ$, and shear, $γ$, required to generate magnification maps. So far, 25 microlensed quasars have been studied using microlensing techniques, where each system has been modelled and analyzed individually, or in small samples. This is about to change due to the upcoming synoptic all-sky surveys, which are expected to discover thousands of quasars suitable for microlensing studies. In the present study we investigate the connection between macromodels of the galaxy lens and microlensing magnification maps throughout the parameter space in preparation for future studies of large statistical samples of systems displaying microlensing. In particular, we use 55,900 maps produced by the GERLUMPH parameter survey (available online at http://gerlumph.swin.edu.au ) and identify regions of parameter space where macromodel uncertainties ($Δκ$,$Δγ$) lead to statistically different magnification maps. Strategies for mitigating the effect of $Δκ$,$Δγ$ uncertainties are discussed in order to understand and control this potential source of systematic errors in accretion disc constraints derived from microlensing.

preprint2013arXiv

A new parameter space study of cosmological microlensing

Cosmological gravitational microlensing is a useful technique for understanding the structure of the inner parts of a quasar, especially the accretion disk and the central supermassive black hole. So far, most of the cosmological microlensing studies have focused on single objects from ~90 currently known lensed quasars. However, present and planned all-sky surveys are expected to discover thousands of new lensed systems. Using a graphics processing unit (GPU) accelerated ray-shooting code, we have generated 2550 magnification maps uniformly across the convergence (κ) and shear (γ) parameter space of interest to microlensing. We examine the effect of random realizations of the microlens positions on map properties such as the magnification probability distribution (MPD). It is shown that for most of the parameter space a single map is representative of an average behaviour. All of the simulations have been carried out on the GPU-Supercomputer for Theoretical Astrophysics Research (gSTAR).

preprint2012arXiv

Accelerating incoherent dedispersion

Incoherent dedispersion is a computationally intensive problem that appears frequently in pulsar and transient astronomy. For current and future transient pipelines, dedispersion can dominate the total execution time, meaning its computational speed acts as a constraint on the quality and quantity of science results. It is thus critical that the algorithm be able to take advantage of trends in commodity computing hardware. With this goal in mind, we present analysis of the 'direct', 'tree' and 'sub-band' dedispersion algorithms with respect to their potential for efficient execution on modern graphics processing units (GPUs). We find all three to be excellent candidates, and proceed to describe implementations in C for CUDA using insight gained from the analysis. Using recent CPU and GPU hardware, the transition to the GPU provides a speed-up of 9x for the direct algorithm when compared to an optimised quad-core CPU code. For realistic recent survey parameters, these speeds are high enough that further optimisation is unnecessary to achieve real-time processing. Where further speed-ups are desirable, we find that the tree and sub-band algorithms are able to provide 3-7x better performance at the cost of certain smearing, memory consumption and development time trade-offs. We finish with a discussion of the implications of these results for future transient surveys. Our GPU dedispersion code is publicly available as a C library at: http://dedisp.googlecode.com/

preprint2011arXiv

Accelerating the Rate of Astronomical Discovery with GPU-Powered Clusters

In recent years, the Graphics Processing Unit (GPU) has emerged as a low-cost alternative for high performance computing, enabling impressive speed-ups for a range of scientific computing applications. Early adopters in astronomy are already benefiting in adapting their codes to take advantage of the GPU's massively parallel processing paradigm. I give an introduction to, and overview of, the use of GPUs in astronomy to date, highlighting the adoption and application trends from the first ~100 GPU-related publications in astronomy. I discuss the opportunities and challenges of utilising GPU computing clusters, such as the new Australian GPU supercomputer, gSTAR, for accelerating the rate of astronomical discovery.

preprint2011arXiv

Fitting Galaxies on GPUs

Structural parameters are normally extracted from observed galaxies by fitting analytic light profiles to the observations. Obtaining accurate fits to high-resolution images is a computationally expensive task, requiring many model evaluations and convolutions with the imaging point spread function. While these algorithms contain high degrees of parallelism, current implementations do not exploit this property. With evergrowing volumes of observational data, an inability to make use of advances in computing power can act as a constraint on scientific outcomes. This is the motivation behind our work, which aims to implement the model-fitting procedure on a graphics processing unit (GPU). We begin by analysing the algorithms involved in model evaluation with respect to their suitability for modern many-core computing architectures like GPUs, finding them to be well-placed to take advantage of the high memory bandwidth offered by this hardware. Following our analysis, we briefly describe a preliminary implementation of the model fitting procedure using freely-available GPU libraries. Early results suggest a speed-up of around 10x over a CPU implementation. We discuss the opportunities such a speed-up could provide, including the ability to use more computationally expensive but better-performing fitting routines to increase the quality and robustness of fits.

preprint2011arXiv

Scientific Visualization in Astronomy: Towards the Petascale Astronomy Era

Astronomy is entering a new era of discovery, coincident with the establishment of new facilities for observation and simulation that will routinely generate petabytes of data. While an increasing reliance on automated data analysis is anticipated, a critical role will remain for visualization-based knowledge discovery. We have investigated scientific visualization applications in astronomy through an examination of the literature published during the last two decades. We identify the two most active fields for progress - visualization of large-N particle data and spectral data cubes - discuss open areas of research, and introduce a mapping between astronomical sources of data and data representations used in general purpose visualization tools. We discuss contributions using high performance computing architectures (e.g: distributed processing and GPUs), collaborative astronomy visualization, the use of workflow systems to store metadata about visualization parameters, and the use of advanced interaction devices. We examine a number of issues that may be limiting the spread of scientific visualization research in astronomy and identify six grand challenges for scientific visualization research in the Petascale Astronomy Era.

preprint2011arXiv

Shape, shear and flexion II - Quantifying the flexion formalism for extended sources with the ray-bundle method

Flexion-based weak gravitational lensing analysis is proving to be a useful adjunct to traditional shear-based techniques. As flexion arises from gradients across an image, analytic and numerical techniques are required to investigate flexion predictions for extended image/source pairs. Using the Schwarzschild lens model, we demonstrate that the ray-bundle method for gravitational lensing can be used to accurately recover second flexion, and is consistent with recovery of zero first flexion. Using lens plane to source plane bundle propagation, we find that second flexion can be recovered with an error no worse than 1% for bundle radii smaller than Δθ = 0.01 θ_E and lens plane impact pararameters greater than θ_E + Δθ, where θ_E is the angular Einstein radius. Using source plane to lens plane bundle propagation, we demonstrate the existence of a preferred flexion zone. For images at radii closer to the lens than the inner boundary of this zone, indicative of the true strong lensing regime, the flexion formalism should be used with caution (errors greater than 5% for extended image/source pairs). We also define a shear zone boundary, beyond which image shapes are essentially indistinguishable from ellipses (1% error in ellipticity). While suggestive that a traditional weak lensing analysis is satisfactory beyond this boundary, a potentially detectable non-zero flexion signal remains.

preprint2011arXiv

Spotting Radio Transients with the help of GPUs

Exploration of the time-domain radio sky has huge potential for advancing our knowledge of the dynamic universe. Past surveys have discovered large numbers of pulsars, rotating radio transients and other transient radio phenomena; however, they have typically relied upon off-line processing to cope with the high data and processing rate. This paradigm rules out the possibility of obtaining high-resolution base-band dumps of significant events or of performing immediate follow-up observations, limiting analysis power to what can be gleaned from detection data alone. To overcome this limitation, real-time processing and detection of transient radio events is required. By exploiting the significant computing power of modern graphics processing units (GPUs), we are developing a transient-detection pipeline that runs in real-time on data from the Parkes radio telescope. In this paper we discuss the algorithms used in our pipeline, the details of their implementation on the GPU and the challenges posed by the presence of radio frequency interference.

preprint2010arXiv

Advanced Architectures for Astrophysical Supercomputing

Astronomers have come to rely on the increasing performance of computers to reduce, analyze, simulate and visualize their data. In this environment, faster computation can mean more science outcomes or the opening up of new parameter spaces for investigation. If we are to avoid major issues when implementing codes on advanced architectures, it is important that we have a solid understanding of our algorithms. A recent addition to the high-performance computing scene that highlights this point is the graphics processing unit (GPU). The hardware originally designed for speeding-up graphics rendering in video games is now achieving speed-ups of $O(100\times)$ in general-purpose computation -- performance that cannot be ignored. We are using a generalized approach, based on the analysis of astronomy algorithms, to identify the optimal problem-types and techniques for taking advantage of both current GPU hardware and future developments in computing architectures.

preprint2010arXiv

Astrophysical Supercomputing with GPUs: Critical Decisions for Early Adopters

General purpose computing on graphics processing units (GPGPU) is dramatically changing the landscape of high performance computing in astronomy. In this paper, we identify and investigate several key decision areas, with a goal of simplyfing the early adoption of GPGPU in astronomy. We consider the merits of OpenCL as an open standard in order to reduce risks associated with coding in a native, vendor-specific programming environment, and present a GPU programming philosophy based on using brute force solutions. We assert that effective use of new GPU-based supercomputing facilities will require a change in approach from astronomers. This will likely include improved programming training, an increased need for software development best-practice through the use of profiling and related optimisation tools, and a greater reliance on third-party code libraries. As with any new technology, those willing to take the risks, and make the investment of time and effort to become early adopters of GPGPU in astronomy, stand to reap great benefits.

preprint2010arXiv

GPU-Based Volume Rendering of Noisy Multi-Spectral Astronomical Data

Traditional analysis techniques may not be sufficient for astronomers to make the best use of the data sets that current and future instruments, such as the Square Kilometre Array and its Pathfinders, will produce. By utilizing the incredible pattern-recognition ability of the human mind, scientific visualization provides an excellent opportunity for astronomers to gain valuable new insight and understanding of their data, particularly when used interactively in 3D. The goal of our work is to establish the feasibility of a real-time 3D monitoring system for data going into the Australian SKA Pathfinder archive. Based on CUDA, an increasingly popular development tool, our work utilizes the massively parallel architecture of modern graphics processing units (GPUs) to provide astronomers with an interactive 3D volume rendering for multi-spectral data sets. Unlike other approaches, we are targeting real time interactive visualization of datasets larger than GPU memory while giving special attention to data with low signal to noise ratio - two critical aspects for astronomy that are missing from most existing scientific visualization software packages. Our framework enables the astronomer to interact with the geometrical representation of the data, and to control the volume rendering process to generate a better representation of their datasets.

preprint2010arXiv

Visualisation and Analysis Challenges for WALLABY

Visualisation and analysis of terabyte-scale datacubes, as will be produced with the Australian Square Kilometre Array Pathfinder (ASKAP), will pose challenges for existing astronomy software and the work practices of astronomers. Focusing on the proposed outcomes of WALLABY (Widefield ASKAP L-Band Legacy All-Sky Blind Survey), and using lessons learnt from HIPASS (HI Parkes All Sky Survey), we identify issues that astronomers will face with WALLABY data cubes. We comment on potential research directions and possible solutions to these challenges.

preprint2010arXiv

Visualization-Directed Interactive Model-Fitting to Spectral Data Cubes

Spectral datasets obtained at radio frequencies and optical/IR wavelengths are increasing in complexity as new facilities and instruments come online, resulting in an increased need to visualize and quantitatively analyze the velocity structures. As the visible structure in spectral data cubes is not purely spatial, additional insight is required to relate structures in 2D space plus line-of-sight velocity to their true three-dimensional (3D) structures. This can be achieved through the use of models that are converted to velocity-space representations. We have used the S2PLOT programming library to enable intuitive, interactive comparison between 3D models and spectral data, with potential for improved understanding of the spatial configurations. We also report on the use of 3D Cartesian shapelets to support quantitative analysis.