GPU accelerated algorithms for computing matrix function vector products with applications to exponential integrators and fractional diffusion

The efficient computation of matrix function vector products has become an important area of research in recent times, driven in particular by two important applications: the numerical solution of fractional partial differential equations and the integration of large systems of ordinary differential equations. In this work we consider a problem that combines these two applications, in the form of a numerical solution algorithm for fractional reaction diffusion equations that after spatial discretisation, is advanced in time using the exponential Euler method. We focus on the efficient implementation of the algorithm on Graphics Processing Units (GPU), as we wish to make use of the increased computational power available with this hardware. We compute the matrix function vector products using the contour integration method in [N. Hale, N. Higham, and L. Trefethen. Computing Aα, log(A), and related matrix functions by contour integrals. SIAM J. Numer. Anal., 46(5):2505–2523, 2008]. Multiple levels of preconditioning are applied to reduce the GPU memory footprint and to further accelerate convergence. We also derive an error bound for the convergence of the contour integral method that allows us to pre-determine the appropriate number of quadrature points. Results are presented that demonstrate the effectiveness of the method for large two-dimensional problems, showing a speedup of more than an order of magnitude compared to a CPU-only implementation.

The short-term accommodation response to aniso-accommodative stimuli in isometropia

Purpose There have been only a limited number of studies examining the accommodative response that occurs when the two eyes are provided with disparate accommodative stimuli, and the results from these studies to date have been equivocal. In this study, we therefore aimed to examine the capacity of the visual system to aniso-accommodate by objectively measuring the interocular difference in the accommodation response between fellow dominant and non-dominant eyes under controlled monocular and binocular viewing conditions during short-term exposure to aniso-accommodative stimuli. Methods The accommodative response of each eye of sixteen young isometropic adults (mean age 22 ± 2 years) with normal binocular vision was measured using an open-field autorefractor during a range of testing conditions; monocularly (accommodative demands ranging from 1.32 to 4.55 D) and binocularly while altering the accommodation demand for each eye (aniso-accommodative stimuli ranging from 0.24 to 2.05 D). Results Under monocular viewing conditions, the dominant and non-dominant eyes displayed a highly symmetric accommodative response; mean interocular difference in spherical equivalent 0.01 ± 0.06 D (relative) and 0.22 ± 0.06 D (absolute) (p>0.05). During binocular viewing, the dominant eye displayed a greater accommodative response (0.11 ± 0.34 D relative and 0.24 ± 0.26 D absolute) irrespective of whether the demand of the dominant or non-dominant eye was altered (p = 0.01). Astigmatic power vectors J0 and J45 did not vary between eyes or with increasing accommodation demands under monocular or binocular viewing conditions (p>0.05). Conclusion The dominant and non-dominant eyes of young isometropic individuals display a similar consensual lag of accommodation under both monocular and binocular viewing conditions, with the dominant eye showing a small but significantly greater (by 0.12 to 0.25 D) accommodative response. Evidence of short-term aniso-accommodation in response to asymmetric accommodation demands was not observed.

Bitboard methods for games

Bitboards allow the efficient encoding of games for computer play and the application of fast bitwiseparallel algorithms for common game-related operations. This article describes: (1) a selection of bitboard techniques including an introduction to bitboards and bitwise operations; (2) a classification scheme that distinguishes filter, query and update methods, and; (3) a sampling of bitboard algorithms for a range of games other than chess, with notes on their performance and practical application.

Positioning as a source of competitive advantage: Benchmarking Rotorua’s position as a domestic short break holiday destination

Travellers are spoilt by holiday choice, and yet will usually only seriously consider a few destinations during the decision process. With thousands of destination marketing organisations (DMOs) competing for attention, places are becoming increasingly substitutable. The study of destination competitiveness is an emerging field, and thesis contributes to an enhanced understanding by addressing three topics that have received relatively little attention in the tourism literature: destination positioning, the context of short break holidays, and domestic travel in New Zealand. A descriptive model of positioning as a source of competitive advantage is developed, and tested through 12 propositions. The destination of interest is Rotorua, which was arguably New Zealand’s first tourist destination. The market of interest is Auckland, which is Rotorua’s largest visitor market. Rotorua’s history is explored to identify factors that may have contributed to the destination’s current image in the Auckland market. A mix of qualitative and quantitative procedures is then utilised to determine Rotorua’s position, relative to a competing set of destinations. Based on an applied research problem, the thesis attempts to bridge the gap between academia and industry by providing useable results and benchmarks for five regional tourism organisations (RTOs). It is proposed that, in New Zealand, the domestic short break market represents a valuable opportunity not explicitly targeted by the competitive set of destinations. Conceptually, the thesis demonstrates the importance of analysing a destination’s competitive position, from the demand perspective, in a travel context; and then the value of comparing this ‘ideal’ position with that projected by the RTO. The thesis concludes Rotorua’s market position in the Auckland short break segment represents a source of comparative advantage, but is not congruent with the current promotional theme, which is being used in all markets. The findings also have implications for destinations beyond the context of the thesis. In particular, a new definition for ‘destination attractiveness’ is proposed, which warrants consideration in the design of future destination positioning analyses.

A survey: Particle Swarm Optimization based algorithms to solve premature convergence problem

Particle Swarm Optimization (PSO) is a biologically inspired computational search and optimization method based on the social behaviors of birds flocking or fish schooling. Although, PSO is represented in solving many well-known numerical test problems, but it suffers from the premature convergence. A number of basic variations have been developed due to solve the premature convergence problem and improve quality of solution founded by the PSO. This study presents a comprehensive survey of the various PSO-based algorithms. As part of this survey, the authors have included a classification of the approaches and they have identify the main features of each proposal. In the last part of the study, some of the topics within this field that are considered as promising areas of future research are listed.

Comparison of IQS and verbal-performance IQ discrepancies estimated from two seven-subtest short forms of the WAIS-R

The present study compared IQs and Verbal-Performance IQ discrepancies estimated from two seven-subtest short forms of the Wechsler Adult Intelligence Scale-Revised (WAIS-R) in a sample of 100 subjects referred for neuropsychological assessment. The short forms of Warrington, James, and Maciejewski (1986) and Ward (1990) yielded similar correlation coefficients and absolute error rates with respect to WAIS-R IQs, although the Warrington short form requires more time to administer and score. Both short forms were able to detect significant Verbal-Performance IQ discrepancies 70% of the time. However, they incorrectly yielded significant discrepancies for approximately 25% of the sample who did not have significant differences on the full WAIS-R. The results do not support reporting and interpreting significant Verbal-Performance IQ discrepancies estimated from these short forms.

Brain network efficiency and topology depend on the fiber tracking method: 11 tractography algorithms compared in 536 subjects

As connectivity analyses become more popular, claims are often made about how the brain's anatomical networks depend on age, sex, or disease. It is unclear how results depend on tractography methods used to compute fiber networks. We applied 11 tractography methods to high angular resolution diffusion images of the brain (4-Tesla 105-gradient HARDI) from 536 healthy young adults. We parcellated 70 cortical regions, yielding 70×70 connectivity matrices, encoding fiber density. We computed popular graph theory metrics, including network efficiency, and characteristic path lengths. Both metrics were robust to the number of spherical harmonics used to model diffusion (4th-8th order). Age effects were detected only for networks computed with the probabilistic Hough transform method, which excludes smaller fibers. Sex and total brain volume affected networks measured with deterministic, tensor-based fiber tracking but not with the Hough method. Each tractography method includes different fibers, which affects inferences made about the reconstructed networks.

A case study of short-sale constraints and limits to arbitrage

This paper uses a unique data set of trades in a unique pair of securities that enables the precise identification of individual broker activity and the trade direction of that activity. We find direct evidence that the imposition (removal) of short-sale constraints limits (generates) trading activity consistent with brokers exploiting apparent mispricing.

Will Obama fast-track the Trans-Pacific Partnership?

This week there has been discussions between leaders from the Pacific Rim over the Trans-Pacific Partnership in Bali, Indonesia at APEC...

New graph-based algorithms to efficiently solve large scale open pit mining optimisation problems

In the mining optimisation literature, most researchers focused on two strategic-level and tactical-level open-pit mine optimisation problems, which are respectively termed ultimate pit limit (UPIT) or constrained pit limit (CPIT). However, many researchers indicate that the substantial numbers of variables and constraints in real-world instances (e.g., with 50-1000 thousand blocks) make the CPIT’s mixed integer programming (MIP) model intractable for use. Thus, it becomes a considerable challenge to solve the large scale CPIT instances without relying on exact MIP optimiser as well as the complicated MIP relaxation/decomposition methods. To take this challenge, two new graph-based algorithms based on network flow graph and conjunctive graph theory are developed by taking advantage of problem properties. The performance of our proposed algorithms is validated by testing recent large scale benchmark UPIT and CPIT instances’ datasets of MineLib in 2013. In comparison to best known results from MineLib, it is shown that the proposed algorithms outperform other CPIT solution approaches existing in the literature. The proposed graph-based algorithms leads to a more competent mine scheduling optimisation expert system because the third-party MIP optimiser is no longer indispensable and random neighbourhood search is not necessary.

Fast growth of polycrystalline graphene by chemical vapor deposition of ethanol on copper

High conductive graphene films can be grown on metal foils by chemical vapor deposition (CVD). We here analyzed the use of ethanol, an economic precursor, which results also safer than commonly-used methane. A comprehensive range of process parameters were explored in order to obtain graphene films with optimal characteristics in view of their use in optoelectronics and photovoltaics. Commercially-available and electro-polished copper foils were used as substrates. By finely tuning the CVD conditions, we obtained few-layer (2-4) graphene films with good conductivity (-500 Ohm/sq) and optical transmittance around 92-94% at 550 nm on unpolished copper foils. The growth on electro-polished copper provides instead predominantly mono-layer films with lower conductivity (>1000 Ohm/sq) and with a transmittance of 97.4% at 550 nm. As for the device properties, graphene with optimal properties as transparent conductive film were produced by CVD on standard copper with specific process conditions.

Channel selection in the short-time modulation domain for distant speech recognition

Automatic speech recognition from multiple distant micro- phones poses significant challenges because of noise and reverberations. The quality of speech acquisition may vary between microphones because of movements of speakers and channel distortions. This paper proposes a channel selection approach for selecting reliable channels based on selection criterion operating in the short-term modulation spectrum domain. The proposed approach quantifies the relative strength of speech from each microphone and speech obtained from beamforming modulations. The new technique is compared experimentally in the real reverb conditions in terms of perceptual evaluation of speech quality (PESQ) measures and word error rate (WER). Overall improvement in recognition rate is observed using delay-sum and superdirective beamformers compared to the case when the channel is selected randomly using circular microphone arrays.

Short-rib polydactyly and jeune syndromes are caused by mutations in WDR60

Short-rib polydactyly syndromes (SRPS I-V) are a group of lethal congenital disorders characterized by shortening of the ribs and long bones, polydactyly, and a range of extraskeletal phenotypes. A number of other disorders in this grouping, including Jeune and Ellis-van Creveld syndromes, have an overlapping but generally milder phenotype. Collectively, these short-rib dysplasias (with or without polydactyly) share a common underlying defect in primary cilium function and form a subset of the ciliopathy disease spectrum. By using whole-exome capture and massive parallel sequencing of DNA from an affected Australian individual with SRPS type III, we detected two novel heterozygous mutations in WDR60, a relatively uncharacterized gene. These mutations segregated appropriately in the unaffected parents and another affected family member, confirming compound heterozygosity, and both were predicted to have a damaging effect on the protein. Analysis of an additional 54 skeletal ciliopathy exomes identified compound heterozygous mutations in WDR60 in a Spanish individual with Jeune syndrome of relatively mild presentation. Of note, these two families share one novel WDR60 missense mutation, although haplotype analysis suggested no shared ancestry. We further show that WDR60 localizes at the base of the primary cilium in wild-type human chondrocytes, and analysis of fibroblasts from affected individuals revealed a defect in ciliogenesis and aberrant accumulation of the GLI2 transcription factor at the centrosome or basal body in the absence of an obvious axoneme. These findings show that WDR60 mutations can cause skeletal ciliopathies and suggest a role for WDR60 in ciliogenesis.

Short-term effects of electronic and tobacco cigarettes on exhaled nitric oxide

The objective of this study was to compare the short-term respiratory effects due to the inhalation of electronic and conventional tobacco cigarette-generated mainstream aerosols through the measurement of the exhaled nitric oxide (eNO). To this purpose, twenty-five smokers were asked to smoke a conventional cigarette and to vape an electronic cigarette (with and without nicotine), and an electronic cigarette without liquid (control session). Electronic and tobacco cigarette mainstream aerosols were characterized in terms of total particle number concentrations and size distributions. On the basis of the measured total particle number concentrations and size distributions, the average particle doses deposited in alveolar and tracheobronchial regions of the lungs for a single 2-s puff were also estimated considering a subject performing resting (sitting) activity. Total particle number concentrations in the mainstream resulted equal to 3.5 ± 0.4 × 109, 5.1 ± 0.1 × 109, and 3.1 ± 0.6 × 109 part. cm− 3 for electronic cigarettes without nicotine, with nicotine, and for conventional cigarettes, respectively. The corresponding alveolar doses for a resting subject were estimated equal to 3.8 × 1010, 5.2 × 1010 and 2.3 × 1010 particles. The mean eNO variations measured after each smoking/vaping session were equal to 3.2 ppb, 2.7 ppb and 2.8 ppb for electronic cigarettes without nicotine, with nicotine, and for conventional cigarettes, respectively; whereas, negligible eNO changes were measured in the control session. Statistical tests performed on eNO data showed statistically significant differences between smoking/vaping sessions and the control session, thus confirming a similar effect on human airways whatever the cigarette smoked/vaped, the nicotine content, and the particle dose received.

High-resolution SNP microarray investigation of copy number variations on chromosome 18 in a control cohort

Copy number variations (CNVs) as described in the healthy population are purported to contribute significantly to genetic heterogeneity. Recent studies have described CNVs using lymphoblastoid cell lines or by application of specifically developed algorithms to interrogate previously described data. However, the full extent of CNVs remains unclear. Using high-density SNP array, we have undertaken a comprehensive investigation of chromosome 18 for CNV discovery and characterisation of distribution and association with chromosome architecture. We identified 399 CNVs, of which loss represents 98%, 58% are less than 2.5 kb in size and 71% are intergenic. Intronic deletions account for the majority of copy number changes with gene involvement. Furthermore, one-third of CNVs do not have putative breakpoints within repetitive sequences. We conclude that replicative processes, mediated either by repetitive elements or microhomology, account for the majority of CNVs in the healthy population. Genomic instability involving the formation of a non-B structure is demonstrated in one region.