Supporting 64-bit global indices in Epetra and other Trilinos packages -- Techniques used and lessons learned

The Trilinos Project is an effort to facilitate the design, development, integration and ongoing support of mathematical software libraries within an object-oriented framework. It is intended for large-scale, complex multiphysics engineering and scientific applications. Epetra is one of its basic packages. It provides serial and parallel linear algebra capabilities. Before Trilinos version 11.0, released in 2012, Epetra used the C++ int data-type for storing global and local indices for degrees of freedom (DOFs). Since int is typically 32-bit, this limited the largest problem size to be smaller than approximately two billion DOFs. This was true even if a distributed memory machine could handle larger problems. We have added optional support for C++ long long data-type, which is at least 64-bit wide, for global indices. To save memory, maintain the speed of memory-bound operations, and reduce further changes to the code, the local indices are still 32-bit. We document the changes required to achieve this feature and how the new functionality can be used. We also report on the lessons learned in modifying a mature and popular package from various perspectives -- design goals, backward compatibility, engineering decisions, C++ language features, effects on existing users and other packages, and build integration

Effects of gastrointestinal tissue structure on computed dipole vectors

<p>Abstract</p> <p>Background</p> <p>Digestive diseases are difficult to assess without using invasive measurements. Non-invasive measurements of body surface electrical and magnetic activity resulting from underlying gastro-intestinal activity are not widely used, in large due to their difficulty in interpretation. Mathematical modelling of the underlying processes may help provide additional information. When modelling myoelectrical activity, it is common for the electrical field to be represented by equivalent dipole sources. The gastrointestinal system is comprised of alternating layers of smooth muscle (SM) cells and Interstitial Cells of Cajal (ICC). In addition the small intestine has regions of high curvature as the intestine bends back upon itself. To eventually use modelling diagnostically, we must improve our understanding of the effect that intestinal structure has on dipole vector behaviour.</p> <p>Methods</p> <p>Normal intestine electrical behaviour was simulated on simple geometries using a monodomain formulation. The myoelectrical fields were then represented by their dipole vectors and an examination on the effect of structure was undertaken. The 3D intestine model was compared to a more computationally efficient 1D representation to determine the differences on the resultant dipole vectors. In addition, the conductivity values and the thickness of the different muscle layers were varied in the 3D model and the effects on the dipole vectors were investigated.</p> <p>Results</p> <p>The dipole vector orientations were largely affected by the curvature and by a transmural gradient in the electrical wavefront caused by the different properties of the SM and ICC layers. This gradient caused the dipoles to be oriented at an angle to the principal direction of electrical propagation. This angle increased when the ratio of the longitudinal and circular muscle was increased or when the the conductivity along and across the layers was increased. The 1D model was able to represent the geometry of the small intestine and successfully captured the propagation of the slow wave down the length of the mesh, however, it was unable to represent transmural diffusion within each layer, meaning the equivalent dipole sources were missing a lateral component and a reduced magnitude when compared to the full 3D models.</p> <p>Conclusion</p> <p>The structure of the intestinal wall affected the potential gradient through the wall and the orientation and magnitude of the dipole vector. We have seen that the models with a symmetrical wall structure and extreme anisotropic conductivities had similar characteristics in their dipole magnitudes and orientations to the 1D model. If efficient 1D models are used instead of 3D models, then both the differences in magnitude and orientation need to be accounted for.</p

J. L. Austin and literal meaning

Alice Crary has recently developed a radical reading of J. L. Austin's philosophy of language. The central contention of Crary's reading is that Austin gives convincing reasons to reject the idea that sentences have context-invariant literal meaning. While I am in sympathy with Crary about the continuing importance of Austin's work, and I think Crary's reading is deep and interesting, I do not think literal sentence meaning is one of Austin's targets, and the arguments that Crary attributes to Austin or finds Austinian in spirit do not provide convincing reasons to reject literal sentence meaning. In this paper, I challenge Crary's reading of Austin and defend the idea of literal sentence meaning

Validation of frequency and mode extraction calculations from time-domain simulations of accelerator cavities

The recently developed frequency extraction algorithm [G.R. Werner and J.R. Cary, J. Comp. Phys. 227, 5200 (2008)] that enables a simple FDTD algorithm to be transformed into an efficient eigenmode solver is applied to a realistic accelerator cavity modeled with embedded boundaries and Richardson extrapolation. Previously, the frequency extraction method was shown to be capable of distinguishing M degenerate modes by running M different simulations and to permit mode extraction with minimal post-processing effort that only requires solving a small eigenvalue problem. Realistic calculations for an accelerator cavity are presented in this work to establish the validity of the method for realistic modeling scenarios and to illustrate the complexities of the computational validation process. The method is found to be able to extract the frequencies with error that is less than a part in 10^5. The corrected experimental and computed values differ by about one parts in 10^$, which is accounted for (in largest part) by machining errors. The extraction of frequencies and modes from accelerator cavities provides engineers and physicists an understanding of potential cavity performance as it depends on shape without incurring manufacture and measurement costs

Electronic energy migration in Microtubules

The repeating arrangement of tubulin dimers confers great mechanical strength to microtubules, which are used as scaffolds for intracellular macromolecular transport in cells and exploited in biohybrid devices. The crystalline order in a microtubule, with lattice constants short enough to allow energy transfer between amino acid chromophores, is similar to synthetic structures designed for light harvesting. After photoexcitation, can these amino acid chromophores transfer excitation energy along the microtubule like a natural or artificial light-harvesting system? Here, we use tryptophan autofluorescence lifetimes to probe energy hopping between aromatic residues in tubulin and microtubules. By studying how the quencher concentration alters tryptophan autofluorescence lifetimes, we demonstrate that electronic energy can diffuse over 6.6 nm in microtubules. We discover that while diffusion lengths are influenced by tubulin polymerization state (free tubulin versus tubulin in the microtubule lattice), they are not significantly altered by the average number of protofilaments (13 versus 14). We also demonstrate that the presence of the anesthetics etomidate and isoflurane reduce exciton diffusion. Energy transport as explained by conventional Förster theory (accommodating for interactions between tryptophan and tyrosine residues) does not sufficiently explain our observations. Our studies indicate that microtubules are, unexpectedly, effective light harvesters

Genetic Variations Strongly Influence Phenotypic Outcome in the Mouse Retina

Variation in genetic background can significantly influence the phenotypic outcome of both disease and non-disease associated traits. Additionally, differences in temporal and strain specific gene expression can also contribute to phenotypes in the mammalian retina. This is the first report of microarray based cross-strain analysis of gene expression in the retina investigating genetic background effects. Microarray analyses were performed on retinas from the following mouse strains: C57BL6/J, AKR/J, CAST/EiJ, and NOD.NON-H2-nb1 at embryonic day 18.5 (E18.5) and postnatal day 30.5 (P30.5). Over 3000 differentially expressed genes were identified between strains and developmental stages. Differential gene expression was confirmed by qRT-PCR, Western blot, and immunohistochemistry. Three major gene networks were identified that function to regulate retinal or photoreceptor development, visual perception, cellular transport, and signal transduction. Many of the genes in these networks are implicated in retinal diseases such as bradyopsia, night-blindness, and cone-rod dystrophy. Our analysis revealed strain specific variations in cone photoreceptor cell patterning and retinal function. This study highlights the substantial impact of genetic background on both development and function of the retina and the level of gene expression differences tolerated for normal retinal function. These strain specific genetic variations may also be present in other tissues. In addition, this study will provide valuable insight for the development of more accurate models for human retinal diseases

SNAPSHOT USA 2019 : a coordinated national camera trap survey of the United States

This article is protected by copyright. All rights reserved.With the accelerating pace of global change, it is imperative that we obtain rapid inventories of the status and distribution of wildlife for ecological inferences and conservation planning. To address this challenge, we launched the SNAPSHOT USA project, a collaborative survey of terrestrial wildlife populations using camera traps across the United States. For our first annual survey, we compiled data across all 50 states during a 14-week period (17 August - 24 November of 2019). We sampled wildlife at 1509 camera trap sites from 110 camera trap arrays covering 12 different ecoregions across four development zones. This effort resulted in 166,036 unique detections of 83 species of mammals and 17 species of birds. All images were processed through the Smithsonian's eMammal camera trap data repository and included an expert review phase to ensure taxonomic accuracy of data, resulting in each picture being reviewed at least twice. The results represent a timely and standardized camera trap survey of the USA. All of the 2019 survey data are made available herein. We are currently repeating surveys in fall 2020, opening up the opportunity to other institutions and cooperators to expand coverage of all the urban-wild gradients and ecophysiographic regions of the country. Future data will be available as the database is updated at eMammal.si.edu/snapshot-usa, as well as future data paper submissions. These data will be useful for local and macroecological research including the examination of community assembly, effects of environmental and anthropogenic landscape variables, effects of fragmentation and extinction debt dynamics, as well as species-specific population dynamics and conservation action plans. There are no copyright restrictions; please cite this paper when using the data for publication.Publisher PDFPeer reviewe

Environmental effects of ozone depletion, UV radiation and interactions with climate change : UNEP Environmental Effects Assessment Panel, update 2017

Peer reviewe

Shared heritability and functional enrichment across six solid cancers

Correction: Nature Communications 10 (2019): art. 4386 DOI: 10.1038/s41467-019-12095-8Quantifying the genetic correlation between cancers can provide important insights into the mechanisms driving cancer etiology. Using genome-wide association study summary statistics across six cancer types based on a total of 296,215 cases and 301,319 controls of European ancestry, here we estimate the pair-wise genetic correlations between breast, colorectal, head/neck, lung, ovary and prostate cancer, and between cancers and 38 other diseases. We observed statistically significant genetic correlations between lung and head/neck cancer (r(g) = 0.57, p = 4.6 x 10(-8)), breast and ovarian cancer (r(g) = 0.24, p = 7 x 10(-5)), breast and lung cancer (r(g) = 0.18, p = 1.5 x 10(-6)) and breast and colorectal cancer (r(g) = 0.15, p = 1.1 x 10(-4)). We also found that multiple cancers are genetically correlated with non-cancer traits including smoking, psychiatric diseases and metabolic characteristics. Functional enrichment analysis revealed a significant excess contribution of conserved and regulatory regions to cancer heritability. Our comprehensive analysis of cross-cancer heritability suggests that solid tumors arising across tissues share in part a common germline genetic basis.Peer reviewe