Design Principles for Sparse Matrix Multiplication on the GPU

We implement two novel algorithms for sparse-matrix dense-matrix multiplication (SpMM) on the GPU. Our algorithms expect the sparse input in the popular compressed-sparse-row (CSR) format and thus do not require expensive format conversion. While previous SpMM work concentrates on thread-level parallelism, we additionally focus on latency hiding with instruction-level parallelism and load-balancing. We show, both theoretically and experimentally, that the proposed SpMM is a better fit for the GPU than previous approaches. We identify a key memory access pattern that allows efficient access into both input and output matrices that is crucial to getting excellent performance on SpMM. By combining these two ingredients---(i) merge-based load-balancing and (ii) row-major coalesced memory access---we demonstrate a 4.1x peak speedup and a 31.7% geomean speedup over state-of-the-art SpMM implementations on real-world datasets.Comment: 16 pages, 7 figures, International European Conference on Parallel and Distributed Computing (Euro-Par) 201

Aortic dissection at the University hospital of the West Indies: A 20-year clinicopathological study of autopsy cases

<p>Abstract</p> <p>Background</p> <p>An autopsy study of aortic dissection (AD) at our institution was previously reported. In the approximately 20 years since then, however, many aspects of diagnosis and treatment of this disease have changed, with a fall in mortality reported in many centers around the world. An impression amongst our pathologists that, there might be an increase in the prevalence of AD in the autopsy service at our hospital, since that earlier report, led to this repeated study, in an attempt to validate that notion. We also sought to identify any changes in clinicopathological features between the two series or any occurring during this study period itself.</p> <p>Findings</p> <p>All cases of AD identified at autopsy, during the 20-year period since the conclusion of the last study, were collected and pertinent clinical and pathological data were analyzed and compared, both within the two decades of this study period and against the results of the last study.</p> <p>Fifty-six cases comprised this study group including 36 males and 20 females, with a mean age of 63.9 years. There were, more patients in the second decade (n = 33; 59%) compared with the first decade (n = 23; 41%). Hypertension as a risk factor was identified in 52 (93%) cases and rupture occurred in 49 (88%) cases. A clinical diagnosis of AD was considered prior to surgery or autopsy in 25 (45%) cases overall, more during the second decade. Surgery was attempted in 25% of all cases with an increase in the second decade compared with the first.</p> <p>Conclusions</p> <p>Compared with the earlier review, a variety of changes in the profile of patients with AD in the autopsy service has been noted, including a reversal in the female predominance seen previously. Other observations include an increase in cases where the correct clinical diagnosis was considered and in which surgical treatment was attempted, changes also evident when the second decade of the present study was compared with the earlier decade. Overall, there were many positive trends. However, areas that could still be improved include an increased index of suspicion for the diagnosis of AD and perhaps in the initiation of treatment, earlier, in those cases where the correct diagnosis was considered.</p

Identification and characterization of a novel non-structural protein of bluetongue virus

Bluetongue virus (BTV) is the causative agent of a major disease of livestock (bluetongue). For over two decades, it has been widely accepted that the 10 segments of the dsRNA genome of BTV encode for 7 structural and 3 non-structural proteins. The non-structural proteins (NS1, NS2, NS3/NS3a) play different key roles during the viral replication cycle. In this study we show that BTV expresses a fourth non-structural protein (that we designated NS4) encoded by an open reading frame in segment 9 overlapping the open reading frame encoding VP6. NS4 is 77–79 amino acid residues in length and highly conserved among several BTV serotypes/strains. NS4 was expressed early post-infection and localized in the nucleoli of BTV infected cells. By reverse genetics, we showed that NS4 is dispensable for BTV replication in vitro, both in mammalian and insect cells, and does not affect viral virulence in murine models of bluetongue infection. Interestingly, NS4 conferred a replication advantage to BTV-8, but not to BTV-1, in cells in an interferon (IFN)-induced antiviral state. However, the BTV-1 NS4 conferred a replication advantage both to a BTV-8 reassortant containing the entire segment 9 of BTV-1 and to a BTV-8 mutant with the NS4 identical to the homologous BTV-1 protein. Collectively, this study suggests that NS4 plays an important role in virus-host interaction and is one of the mechanisms played, at least by BTV-8, to counteract the antiviral response of the host. In addition, the distinct nucleolar localization of NS4, being expressed by a virus that replicates exclusively in the cytoplasm, offers new avenues to investigate the multiple roles played by the nucleolus in the biology of the cell

Mesoscopic organization reveals the constraints governing C. elegans nervous system

One of the biggest challenges in biology is to understand how activity at the cellular level of neurons, as a result of their mutual interactions, leads to the observed behavior of an organism responding to a variety of environmental stimuli. Investigating the intermediate or mesoscopic level of organization in the nervous system is a vital step towards understanding how the integration of micro-level dynamics results in macro-level functioning. In this paper, we have considered the somatic nervous system of the nematode Caenorhabditis elegans, for which the entire neuronal connectivity diagram is known. We focus on the organization of the system into modules, i.e., neuronal groups having relatively higher connection density compared to that of the overall network. We show that this mesoscopic feature cannot be explained exclusively in terms of considerations, such as optimizing for resource constraints (viz., total wiring cost) and communication efficiency (i.e., network path length). Comparison with other complex networks designed for efficient transport (of signals or resources) implies that neuronal networks form a distinct class. This suggests that the principal function of the network, viz., processing of sensory information resulting in appropriate motor response, may be playing a vital role in determining the connection topology. Using modular spectral analysis, we make explicit the intimate relation between function and structure in the nervous system. This is further brought out by identifying functionally critical neurons purely on the basis of patterns of intra- and inter-modular connections. Our study reveals how the design of the nervous system reflects several constraints, including its key functional role as a processor of information.Comment: Published version, Minor modifications, 16 pages, 9 figure

Body Mass Index is Associated with USF1 Haplotype in Korean Premenopausal Women

The upstream stimulatory factor 1 (USF1) gene has been shown to play an essential role as the cause of familial combined hyperlipidemia, and there are several association studies on the relationship between USF1 and metabolic disorders. In this study, we analyzed two single nucleotide polymorphisms in USF1 rs2073653 (306A>G) and rs2516840 (1748C>T) between the case (dyslipidemia or obesity) group and the control group in premenopausal females, postmenopausal females, and males among 275 Korean subjects. We observed a statistically significant difference in the GC haplotype between body mass index (BMI) ≥25 kg/m2 and BMI <25 kg/m2 groups in premenopausal females (χ2=4.23, p=0.04). It seems that the USF1 GC haplotype is associated with BMI in premenopausal Korean females

Search for CP violation in D0 and D+ decays

A high statistics sample of photoproduced charm particles from the FOCUS (E831) experiment at Fermilab has been used to search for CP violation in the Cabibbo suppressed decay modes D+ to K-K+pi+, D0 to K-K+ and D0 to pi-pi+. We have measured the following CP asymmetry parameters: A_CP(K-K+pi+) = +0.006 +/- 0.011 +/- 0.005, A_CP(K-K+) = -0.001 +/- 0.022 +/- 0.015 and A_CP(pi-pi+) = +0.048 +/- 0.039 +/- 0.025 where the first error is statistical and the second error is systematic. These asymmetries are consistent with zero with smaller errors than previous measurements.Comment: 12 pages, 4 figure

A Study of D0 --> K0(S) K0(S) X Decay Channels

Using data from the FOCUS experiment (FNAL-E831), we report on the decay of $D^0$ mesons into final states containing more than one $K^0_S$. We present evidence for two Cabibbo favored decay modes, $D^0\to K^0_SK^0_S K^- \pi^+$ and $D^0\to K^0_SK^0_S K^+ \pi^-$, and measure their combined branching fraction relative to $D^0\to \bar{K} ^0\pi^+\pi^-$ to be $\frac{\Gamma(D^0\to K^0_SK^0_SK^{\pm}\pi^{\mp})}{\Gamma(D^0\to \bar{K} ^0\pi^+\pi^-)}$ = 0.0106 $\pm$ 0.0019 $\pm$ 0.0010. Further, we report new measurements of $\frac{\Gamma(D^0\to K^0_SK^0_SK^0_S)}{\Gamma(D^0\to \bar{K} ^0\pi^+\pi^-)}$ = 0.0179 $\pm$ 0.0027 $\pm$ 0.0026, $\frac{\Gamma(D^0\to K^0\bar{K} ^0)}{\Gamma(D^0\to \bar{K} ^0\pi^+\pi^-)}$ = 0.0144 $\pm$ 0.0032 $\pm$ 0.0016, and $\frac{\Gamma(D^0\to K^0_SK^0_S\pi^+\pi^-)}{\Gamma(D^0\to \bar{K} ^0\pi^+\pi^-)}$ = 0.0208 $\pm$ 0.0035 $\pm$ 0.0021 where the first error is statistical and the second is systematic.Comment: 11 pages, 3 figures, typos correcte

Architecture of an Antagonistic Tree/Fungus Network: The Asymmetric Influence of Past Evolutionary History

Compartmentalization and nestedness are common patterns in ecological networks. The aim of this study was to elucidate some of the processes shaping these patterns in a well resolved network of host/pathogen interactions.Based on a long-term (1972-2005) survey of forest health at the regional scale (all French forests; 15 million ha), we uncovered an almost fully connected network of 51 tree taxa and 157 parasitic fungal species. Our analyses revealed that the compartmentalization of the network maps out the ancient evolutionary history of seed plants, but not the ancient evolutionary history of fungal species. The very early divergence of the major fungal phyla may account for this asymmetric influence of past evolutionary history. Unlike compartmentalization, nestedness did not reflect any consistent phylogenetic signal. Instead, it seemed to reflect the ecological features of the current species, such as the relative abundance of tree species and the life-history strategies of fungal pathogens. We discussed how the evolution of host range in fungal species may account for the observed nested patterns.Overall, our analyses emphasized how the current complexity of ecological networks results from the diversification of the species and their interactions over evolutionary times. They confirmed that the current architecture of ecological networks is not only dependent on recent ecological processes

Evolving network structure of academic institutions

Today’s colleges and universities consist of highly complex structures that dictate interactions between the administration, faculty, and student body. These structures can play a role in dictating the efficiency of policy enacted by the administration and determine the effect that curriculum changes in one department have on other departments. Despite the fact that the features of these complex structures have a strong impact on the institutions, they remain by-and-large unknown in many cases. In this paper we study the academic structure of our home institution of Trinity College in Hartford, CT using the major and minor patterns between graduating students to build a temporal multiplex network describing the interactions between different departments. Using recent network science techniques developed for such temporal networks we identify the evolving community structures that organize departments’ interactions, as well as quantify the interdisciplinary centrality of each department. We implement this framework for Trinity College, finding practical insights and applications, but also present it as a general framework for colleges and universities to better understand their own structural makeup in order to better inform academic and administrative policy

The dependence of dijet production on photon virtuality in ep collisions at HERA

The dependence of dijet production on the virtuality of the exchanged photon, Q^2, has been studied by measuring dijet cross sections in the range 0 < Q^2 < 2000 GeV^2 with the ZEUS detector at HERA using an integrated luminosity of 38.6 pb^-1. Dijet cross sections were measured for jets with transverse energy E_T^jet > 7.5 and 6.5 GeV and pseudorapidities in the photon-proton centre-of-mass frame in the range -3 < eta^jet <0. The variable xg^obs, a measure of the photon momentum entering the hard process, was used to enhance the sensitivity of the measurement to the photon structure. The Q^2 dependence of the ratio of low- to high-xg^obs events was measured. Next-to-leading-order QCD predictions were found to generally underestimate the low-xg^obs contribution relative to that at high xg^obs. Monte Carlo models based on leading-logarithmic parton-showers, using a partonic structure for the photon which falls smoothly with increasing Q^2, provide a qualitative description of the data.Comment: 35 pages, 6 eps figures, submitted to Eur.Phys.J.