Analysis of and workarounds for element reversal for a finite element-based algorithm for warping triangular and tetrahedral meshes

We consider an algorithm called FEMWARP for warping triangular and tetrahedral finite element meshes that computes the warping using the finite element method itself. The algorithm takes as input a two- or three-dimensional domain defined by a boundary mesh (segments in one dimension or triangles in two dimensions) that has a volume mesh (triangles in two dimensions or tetrahedra in three dimensions) in its interior. It also takes as input a prescribed movement of the boundary mesh. It computes as output updated positions of the vertices of the volume mesh. The first step of the algorithm is to determine from the initial mesh a set of local weights for each interior vertex that describes each interior vertex in terms of the positions of its neighbors. These weights are computed using a finite element stiffness matrix. After a boundary transformation is applied, a linear system of equations based upon the weights is solved to determine the final positions of the interior vertices. The FEMWARP algorithm has been considered in the previous literature (e.g., in a 2001 paper by Baker). FEMWARP has been succesful in computing deformed meshes for certain applications. However, sometimes FEMWARP reverses elements; this is our main concern in this paper. We analyze the causes for this undesirable behavior and propose several techniques to make the method more robust against reversals. The most successful of the proposed methods includes combining FEMWARP with an optimization-based untangler.Comment: Revision of earlier version of paper. Submitted for publication in BIT Numerical Mathematics on 27 April 2010. Accepted for publication on 7 September 2010. Published online on 9 October 2010. The final publication is available at http://www.springerlink.co

Evaluating Elevated Convection with the Downdraft Convective Inhibition

A method for evaluating the penetration of a stable layer by an elevated convective downdraft is discussed. Some controversy exists on the community’s ability to define truly elevated convection from surface-based convection. By comparing the downdraft convective inhibition (DCIN) to the downdraft convective available potential energy (DCAPE), we determine that downdraft penetration potential is progressively enabled as the DCIN is progressively smaller than the DCAPE; inversely as DCIN increases over DCAPE, so does the likelihood of purely elevated convection. Serial vertical soundings and accompanying analyses are provided to support this finding

Incoherent dynamics in the toric code subject to disorder

We numerically study the effects of two forms of quenched disorder on the anyons of the toric code. Firstly, a new class of codes based on random lattices of stabilizer operators is presented, and shown to be superior to the standard square lattice toric code for certain forms of biased noise. It is further argued that these codes are close to optimal, in that they tightly reach the upper bound of error thresholds beyond which no correctable CSS codes can exist. Additionally, we study the classical motion of anyons in toric codes with randomly distributed onsite potentials. In the presence of repulsive long-range interaction between the anyons, a surprising increase with disorder strength of the lifetime of encoded states is reported and explained by an entirely incoherent mechanism. Finally, the coherent transport of the anyons in the presence of both forms of disorder is investigated, and a significant suppression of the anyon motion is found.Comment: 13 pages, 12 figure

Exploiting asynchrony from exact forward recovery for DUE in iterative solvers

This paper presents a method to protect iterative solvers from Detected and Uncorrected Errors (DUE) relying on error detection techniques already available in commodity hardware. Detection operates at the memory page level, which enables the use of simple algorithmic redundancies to correct errors. Such redundancies would be inapplicable under coarse grain error detection, but become very powerful when the hardware is able to precisely detect errors. Relations straightforwardly extracted from the solver allow to recover lost data exactly. This method is free of the overheads of backwards recoveries like checkpointing, and does not compromise mathematical convergence properties of the solver as restarting would do. We apply this recovery to three widely used Krylov subspace methods, CG, GMRES and BiCGStab, and their preconditioned versions. We implement our resilience techniques on CG considering scenarios from small (8 cores) to large (1024 cores) scales, and demonstrate very low overheads compared to state-of-the-art solutions. We deploy our recovery techniques either by overlapping them with algorithmic computations or by forcing them to be in the critical path of the application. A trade-off exists between both approaches depending on the error rate the solver is suffering. Under realistic error rates, overlapping decreases overheads from 5.37% down to 3.59% for a non-preconditioned CG on 8 cores.This work has been partially supported by the European Research Council under the European Union's 7th FP, ERC Advanced Grant 321253, and by the Spanish Ministry of Science and Innovation under grant TIN2012-34557. L. Jaulmes has been partially supported by the Spanish Ministry of Education, Culture and Sports under grant FPU2013/06982. M. Moreto has been partially supported by the Spanish Ministry of Economy and Competitiveness under Juan de la Cierva postdoctoral fellowship JCI-2012-15047. M. Casas has been partially supported by the Secretary for Universities and Research of the Ministry of Economy and Knowledge of the Government of Catalonia and the Co-fund programme of the Marie Curie Actions of the European Union's 7th FP (contract 2013 BP B 00243).Peer ReviewedPostprint (author's final draft

Whirling Hexagons and Defect Chaos in Hexagonal Non-Boussinesq Convection

We study hexagon patterns in non-Boussinesq convection of a thin rotating layer of water. For realistic parameters and boundary conditions we identify various linear instabilities of the pattern. We focus on the dynamics arising from an oscillatory side-band instability that leads to a spatially disordered chaotic state characterized by oscillating (whirling) hexagons. Using triangulation we obtain the distribution functions for the number of pentagonal and heptagonal convection cells. In contrast to the results found for defect chaos in the complex Ginzburg-Landau equation and in inclined-layer convection, the distribution functions can show deviations from a squared Poisson distribution that suggest non-trivial correlations between the defects.Comment: 4 mpg-movies are available at http://www.esam.northwestern.edu/~riecke/lit/lit.html submitted to New J. Physic

Evidence on anti-malarial and diagnostic markets in Cambodia to guide malaria elimination strategies and policies

BACKGROUND: Understanding Cambodia's anti-malarial and diagnostic landscape in 2015 is critical for informing and monitoring strategies and policies as Cambodia moves forward with national efforts to eliminate malaria. The aim of this paper is to present timely and key findings on the public and private sector anti-malarial and diagnostic landscape in Cambodia. This evidence can serve as a baseline benchmark for guiding implementation of national strategies as well as other regional initiatives to address malaria elimination activities. METHODS: From August 17th to October 1st, 2015, a cross sectional, nationally-representative malaria outlet survey was conducted in Cambodia. A census of all public and private outlets with potential to distribute malaria testing and/or treatment was conducted among 180 communes. An audit was completed for all anti-malarials, malaria rapid diagnostic tests (RDT) and microscopy. RESULTS: A total of 26,664 outlets were screened, and 1303 outlets were eligible and interviewed. Among all screened outlets in the public sector, 75.9% of public health facilities and 67.7% of community health workers stocked both malaria diagnostic testing and a first-line artemisinin-based combination therapy (ACT). Among anti-malarial-stocking private sector outlets, 64.7% had malaria blood testing available, and 70.9% were stocking a first-line ACT. Market share data illustrate that most of the anti-malarials were sold or distributed through the private sector (58.4%), including itinerant drug vendors (23.4%). First-line ACT accounted for the majority of the market share across the public and private sectors (90.3%). Among private sector outlets stocking any anti-malarial, the proportion of outlets with a first-line ACT or RDT was higher among outlets that had reportedly received one or more forms of 'support' (e.g. reportedly received training in the previous year on malaria diagnosis [RDT and/or microscopy] and/or the national treatment guidelines for malaria) compared to outlets that did not report receiving any support (ACT: 82.1 and 60.6%, respectively; RDT: 78.2 and 64.0%, respectively). CONCLUSION: The results point to high availability and distribution of first-line ACT and widespread availability of malaria diagnosis, especially in the public sector. This suggests that there is a strong foundation for achieving elimination goals in Cambodia. However, key gaps in terms of availability of malaria commodities for case management must be addressed, particularly in the private sector where most people seek treatment. Continued engagement with the private sector will be important to ensure accelerated progress towards malaria elimination

First Neutrino Observations from the Sudbury Neutrino Observatory

The first neutrino observations from the Sudbury Neutrino Observatory are presented from preliminary analyses. Based on energy, direction and location, the data in the region of interest appear to be dominated by 8B solar neutrinos, detected by the charged current reaction on deuterium and elastic scattering from electrons, with very little background. Measurements of radioactive backgrounds indicate that the measurement of all active neutrino types via the neutral current reaction on deuterium will be possible with small systematic uncertainties. Quantitative results for the fluxes observed with these reactions will be provided when further calibrations have been completed.Comment: Latex, 7 pages, 10 figures, Invited paper at Neutrino 2000 Conference, Sudbury, Canada, June 16-21, 2000 to be published in the Proceeding