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

All degree six local unitary invariants of k qudits

We give explicit index-free formulae for all the degree six (and also degree four and two) algebraically independent local unitary invariant polynomials for finite dimensional k-partite pure and mixed quantum states. We carry out this by the use of graph-technical methods, which provides illustrations for this abstract topic.Comment: 18 pages, 6 figures, extended version. Comments are welcom

Mid-Miocene cooling and the extinction of tundra in continental Antarctica

A major obstacle in understanding the evolution of Cenozoic climate has been the lack of well dated terrestrial evidence from high-latitude, glaciated regions. Here, we report the discovery of exceptionally well preserved fossils of lacustrine and terrestrial organisms from the McMurdo Dry Valleys sector of the Transantarctic Mountains for which we have established a precise radiometric chronology. The fossils, which include diatoms, palynomorphs, mosses, ostracodes, and insects, represent the last vestige of a tundra community that inhabited the mountains before stepped cooling that first brought a full polar climate to Antarctica. Paleoecological analyses, 40Ar/39Ar analyses of associated ash fall, and climate inferences from glaciological modeling together suggest that mean summer temperatures in the region cooled by at least 8°C between 14.07 ± 0.05 Ma and 13.85 ± 0.03 Ma. These results provide novel constraints for the timing and amplitude of middle-Miocene cooling in Antarctica and reveal the ecological legacy of this global climate transition

Trilinos users guide.

The Trilinos Project is an effort to facilitate the design, development, integration and ongoing support of mathematical software libraries. A new software capability is introduced into Trilinos as a package. A Trilinos package is an integral unit usually developed by a small team of experts in a particular algorithms area such as algebraic preconditioners, nonlinear solvers, etc. The Trilinos Users Guide is a resource for new and existing Trilinos users. Topics covered include how to configure and build Trilinos, what is required to integrate an existing package into Trilinos and examples of how those requirements can be met, as well as what tools and services are available to Trilinos packages. Also discussed are some common practices that are followed by many Trilinos package developers. Finally, a snapshot of current Trilinos packages and their interoperability status is provided, along with a list of supported computer platforms

On the algebra of local unitary invariants of pure and mixed quantum states

We study the structure of the inverse limit of the graded algebras of local unitary invariant polynomials using its Hilbert series. For k subsystems, we conjecture that the inverse limit is a free algebra and the number of algebraically independent generators with homogenous degree 2m equals the number of conjugacy classes of index m subgroups in a free group on k-1 generators. Similarly, we conjecture that the inverse limit in the case of k-partite mixed state invariants is free and the number of algebraically independent generators with homogenous degree m equals the number of conjugacy classes of index m subgroups in a free group on k generators. The two conjectures are shown to be equivalent. To illustrate the equivalence, using the representation theory of the unitary groups, we obtain all invariants in the m=2 graded parts and express them in a simple form both in the case of mixed and pure states. The transformation between the two forms is also derived. Analogous invariants of higher degree are also introduced.Comment: 14 pages, no figure

An actionable anti-racism plan for geoscience organizations.

Geoscience organizations shape the discipline. They influence attitudes and expectations, set standards, and provide benefits to their members. Today, racism and discrimination limit the participation of, and promote hostility towards, members of minoritized groups within these critical geoscience spaces. This is particularly harmful for Black, Indigenous, and other people of color in geoscience and is further exacerbated along other axes of marginalization, including disability status and gender identity. Here we present a twenty-point anti-racism plan that organizations can implement to build an inclusive, equitable and accessible geoscience community. Enacting it will combat racism, discrimination, and the harassment of all members

Overview of the VA Quality Enhancement Research Initiative (QUERI) and QUERI theme articles: QUERI Series

<p>Abstract</p> <p>Background</p> <p>Continuing challenges to timely adoption of evidence-based clinical practices in healthcare have generated intense interest in the development and application of new implementation methods and frameworks. These challenges led the United States (U.S.) Department of Veterans Affairs (VA) to create the Quality Enhancement Research Initiative (QUERI) in the late 1990s. QUERI's purpose was to harness VA's health services research expertise and resources in an ongoing system-wide effort to improve the performance of the VA healthcare system and, thus, quality of care for veterans. QUERI in turn created a systematic means of involving VA researchers both in enhancing VA healthcare quality, by implementing evidence-based practices, and in contributing to the continuing development of implementation science.</p> <p>The efforts of VA researchers to improve healthcare delivery practices through QUERI and related initiatives are documented in a growing body of literature. The scientific frameworks and methodological approaches developed and employed by QUERI are less well described. A QUERI Series of articles in <it>Implementation Science </it>will illustrate many of these QUERI tools. This <it>Overview </it>article introduces both QUERI and the Series.</p> <p>Methods</p> <p>The <it>Overview </it>briefly explains the purpose and context of the QUERI Program. It then describes the following: the key operational structure of QUERI Centers, guiding frameworks designed to enhance implementation and related research, QUERI's progress and promise to date, and the Series' general content. QUERI's frameworks include a core set of steps for diagnosing and closing quality gaps and, simultaneously, advancing implementation science. Throughout the paper, the envisioned involvement and activities of VA researchers within QUERI Centers also are highlighted. The Series is then described, illustrating the use of QUERI frameworks and other tools designed to respond to implementation challenges.</p> <p>Conclusion</p> <p>QUERI's simultaneous pursuit of improvement and research goals within a large healthcare system may be unique. However, descriptions of this still-evolving effort, including its conceptual frameworks, methodological approaches, and enabling processes, should have applicability to implementation researchers in a range of health care settings. Thus, the <it>Series </it>is offered as a resource for other implementation research programs and researchers pursuing common goals in improving care and developing the field of implementation science.</p

General Anesthetics Predicted to Block the GLIC Pore with Micromolar Affinity

Although general anesthetics are known to modulate the activity of ligand-gated ion channels in the Cys-loop superfamily, there is at present neither consensus on the underlying mechanisms, nor predictive models of this modulation. Viable models need to offer quantitative assessment of the relative importance of several identified anesthetic binding sites. However, to date, precise affinity data for individual sites has been challenging to obtain by biophysical means. Here, the likely role of pore block inhibition by the general anesthetics isoflurane and propofol of the prokaryotic pentameric channel GLIC is investigated by molecular simulations. Microscopic affinities are calculated for both single and double occupancy binding of isoflurane and propofol to the GLIC pore. Computations are carried out for an open-pore conformation in which the pore is restrained to crystallographic radius, and a closed-pore conformation that results from unrestrained molecular dynamics equilibration of the structure. The GLIC pore is predicted to be blocked at the micromolar concentrations for which inhibition by isofluorane and propofol is observed experimentally. Calculated affinities suggest that pore block by propofol occurs at signifcantly lower concentrations than those for which inhibition is observed: we argue that this discrepancy may result from binding of propofol to an allosteric site recently identified by X-ray crystallography, which may cause a competing gain-of-function effect. Affinities of isoflurane and propofol to the allosteric site are also calculated, and shown to be 3 mM for isoflurane and for propofol; both anesthetics have a lower affinity for the allosteric site than for the unoccupied pore

Single-Cell Analysis of Ploidy and Centrosomes Underscores the Peculiarity of Normal Hepatocytes

Polyploidization is the most well recognized feature of the liver. Yet, a quantitative and behavioral analysis of centrosomes and DNA content in normal hepatocytes has been limited by the technical challenges of methods available. By using a novel approach employing FISH for chromosomes 18, X and Y we provide, for the first time, a detailed analysis of DNA copies during physiological development in the liver at single cell level. We demonstrate that aneuploidy and unbalanced DNA content in binucleated hepatocytes are common features in normal adult liver. Despite the common belief that hepatocytes contain 1, 2 or no more than 4 centrosomes, our double staining for centrosome associated proteins reveals extranumerary centrosomes in a high percentage of cells as early as 15 days of age. We show that in murine liver the period between 15 days and 1.5 months marks the transition from a prevalence of mononucleated cells to up to 75% of binucleated cells. Our data demonstrate that this timing correlates with a switch in centrosomes number. At 15 days the expected 1 or 2 centrosomes converge with several hepatocytes that contain 3 centrosomes; at 1.5 months the percentage of cells with 3 centrosomes decreases concomitantly with the increase of cells with more than 4 centrosomes. Our analysis shows that the extranumerary centrosomes emerge in concomitance with the process of binucleation and polyploidization and maintain α-tubulin nucleation activity. Finally, by integrating interphase FISH and immunofluorescent approaches, we detected an imbalance between centrosome number and DNA content in liver cells that deviates from the equilibrium expected in normal cells. We speculate that these unique features are relevant to the peculiar biological function of liver cells which are continuously challenged by stress, a condition that could predispose to genomic instability