Massively parallel approximate Gaussian process regression

We explore how the big-three computing paradigms -- symmetric multi-processor (SMC), graphical processing units (GPUs), and cluster computing -- can together be brought to bare on large-data Gaussian processes (GP) regression problems via a careful implementation of a newly developed local approximation scheme. Our methodological contribution focuses primarily on GPU computation, as this requires the most care and also provides the largest performance boost. However, in our empirical work we study the relative merits of all three paradigms to determine how best to combine them. The paper concludes with two case studies. One is a real data fluid-dynamics computer experiment which benefits from the local nature of our approximation; the second is a synthetic data example designed to find the largest design for which (accurate) GP emulation can performed on a commensurate predictive set under an hour.Comment: 24 pages, 6 figures, 1 tabl

A spectroscopic study of some early-type emission-line stars

Spectroscopic observations of 18 of the brighter Be and 0e stars in the Northern Hemisphere have been presented in this project. About 900 spectrograms were secured and measured using the more objective numerical technique (REDUCE and VCROSS) to obtain as many accurate radial velocities as possible for these stars. The analysis of these measurements has demonstrated clearly, that most of these early-type stars are indeed variable in radial velocity, while a real periodicity in the radial velocities of some of the programme stars is found. The long-term periodicities were attributed to the binary nature, while the short ones could be ascribed to the pulsation phenomenon. Orbital elements are determined for four newly discovered binary systems, and improved orbits are determined for four known binaries. Five other stars display short time-scale periodicity (less than 1.5 days) attributable to radial or non-radial pulsation, whilst three variables show little evidence for periodicity. Only two stars are found to have constant radial velocities. A brief summary of the properties of Be stars and all the available models to explain the Be phenomenon is presented in chapter one. Chapter two contains a brief description of the observing equipment and the techniques used. The reasons for selecting the programme stars are given in the same chapter. A complete explanation of the measuring technique used in this investigation is given in chapter three. In chapter four, all the analyses and the results for each star are given individually, while a discussion of previous work on these stars appears in the relevant sections. The final chapter contains a discussion of the results from this study together with an analysis of the energy distributions of the programme stars, compiled by the author from published fluxes and photometry from the ultraviolet to the infra-red parts of the electromagnetic spectrum. Some consideration of future work to be done on these types of stars is also given

Sub-milliarcsecond precision spectro-astrometry of Be stars

The origin of the disks around Be stars is still not known. Further progress requires a proper parametrization of their structure, both spatially and kinematically. This is challenging as the disks are very small. Here we assess whether a novel method is capable of providing these data. We obtained spectro astrometry around the Pa beta line of two bright Be stars, alpha Col and zeta Tau, to search for disk signatures. The data, with a pixel to pixel precision of the centroid position of 0.3..0.4 milliarcsecond is the most accurate such data to date. Artefacts at the 0.85 mas level are present in the data, but these are readily identified as they were non-repeatable in our redundant datasets. This does illustrate the need of taking multiple data to avoid spurious detections. The data are compared with simple model simulations of the spectro astrometric signatures due to rotating disks around Be stars. The upper limits we find for the disk radii correspond to disk sizes of a few dozen stellar radii if they rotate Keplerian. This is very close to observationally measured and theoretically expected disk sizes, and this paper therefore demonstrates that spectro-astrometry, of which we present the first such attempt, has the potential to resolve the disks around Be stars.Comment: 6 pages, A&A accepte

Comparison of Synthetic Computed Tomography Generation Methods, Incorporating Male and Female Anatomical Differences, for Magnetic Resonance Imaging-Only Definitive Pelvic Radiotherapy

Purpose: There are several means of synthetic computed tomography (sCT) generation for magnetic resonance imaging (MRI)-only planning; however, much of the research omits large pelvic treatment regions and female anatomical specific methods. This research aimed to apply four of the most popular methods of sCT creation to facilitate MRI-only radiotherapy treatment planning for male and female anorectal and gynecological neoplasms. sCT methods were validated against conventional computed tomography (CT), with regard to Hounsfield unit (HU) estimation and plan dosimetry. Methods and Materials: Paired MRI and CT scans of 40 patients were used for sCT generation and validation. Bulk density assignment, tissue class density assignment, hybrid atlas, and deep learning sCT generation methods were applied to all 40 patients. Dosimetric accuracy was assessed by dose difference at reference point, dose volume histogram (DVH) parameters, and 3D gamma dose comparison. HU estimation was assessed by mean error and mean absolute error in HU value between each sCT and CT. Results: The median percentage dose difference between the CT and sCT was &lt;1.0% for all sCT methods. The deep learning method resulted in the lowest median percentage dose difference to CT at −0.03% (IQR 0.13, −0.31) and bulk density assignment resulted in the greatest difference at −0.73% (IQR −0.10, −1.01). The mean 3D gamma dose agreement at 3%/2 mm among all sCT methods was 99.8%. The highest agreement at 1%/1 mm was 97.3% for the deep learning method and the lowest was 93.6% for the bulk density method. Deep learning and hybrid atlas techniques gave the lowest difference to CT in mean error and mean absolute error in HU estimation. Conclusions: All methods of sCT generation used in this study resulted in similarly high dosimetric agreement for MRI-only planning of male and female cancer pelvic regions. The choice of the sCT generation technique can be guided by department resources available and image guidance considerations, with minimal impact on dosimetric accuracy.</p

Properties and nature of Be stars: 29. Orbital and long-term spectral variations of \gamma\ Cassiopei\ae

A detailed analysis of more than 800 electronic high-resolution spectra of gamma Cas, which were obtained during a time interval of over 6000 days (16.84 yrs) at several observatories, documents the smooth variations in the density and/or extend of its circumstellar envelope. We found a clear anticorrelation between the peak intensity and FWHM of the H alpha emission, which seems to agree with recent models of such emission lines. The main result of this study is a confirmation of the binary nature of the object, determination of a reliable linear ephemeris T_{min.RV} = HJD (2452081.9$\pm$0.6) + (203.52$\pm$0.08)*E, and a rather definitive set of orbital elements. We clearly demonstrated that the orbit is circular within the limits of accuracy of our measurements and has a semi-amplitude of radial-velocity curve of 4.30$\pm$0.09 (km/s). No trace of the low-mass secondary was found. The time distribution of our spectra does not allow a reliable investigation of rapid spectral variations, which are undoubtedly present in the spectra. We postpone this investigation for a future study, based on series of dedicated whole-night spectral observations

Optimisation and validation of an integrated magnetic resonance imaging-only radiotherapy planning solution

Background and purpose: Magnetic resonance imaging (MRI)-only treatment planning is gaining in popularity in radiation oncology, with various methods available to generate a synthetic computed tomography (sCT) for this purpose. The aim of this study was to validate a sCT generation software for MRI-only radiotherapy planning of male and female pelvic cancers. The secondary aim of this study was to improve dose agreement by applying a derived relative electron and mass density (RED) curve to the sCT. Method and materials: Computed tomography (CT) and MRI scans of forty patients with pelvic neoplasms were used in the study. Treatment plans were copied from the CT scan to the sCT scan for dose comparison. Dose difference at reference point, 3D gamma comparison and dose volume histogram analysis was used to validate the dose impact of the sCT. The RED values were optimised to improve dose agreement by using a linear plot. Results: The average percentage dose difference at isocentre was 1.2% and the mean 3D gamma comparison with a criteria of 1%/1 mm was 84.0% ± 9.7%. The results indicate an inherent systematic difference in the dosimetry of the sCT plans, deriving from the tissue densities. With the adapted REDmod table, the average percentage dose difference was reduced to −0.1% and the mean 3D gamma analysis improved to 92.9% ± 5.7% at 1%/1 mm. Conclusions: CT generation software is a viable solution for MRI-only radiotherapy planning. The option makes it relatively easy for departments to implement a MRI-only planning workflow for cancers of male and female pelvic anatomy.</p