Computational Comparison of Continuous and Discontinuous Galerkin Time-Stepping Methods for Nonlinear Initial Value Problems

This article centers on the computational performance of the continuous and discontinuous Galerkin time stepping schemes for general first-order initial value problems in R n , with continuous nonlinearities. We briefly review a recent existence result for discrete solutions from [6], and provide a numerical comparison of the two time discretization methods.

High-resolution sequence stratigraphic correlation in the Upper Jurassic (Kimmeridgian)–Upper Cretaceous (Cenomanian) peritidal carbonate deposits (Western Taurides, Turkey)

Upper Jurassic (Kimmeridgian)±Upper Cretaceous (Cenomanian) inner platform carbonates in the Western Taurides are composed of metre-scale upward-shallowing cyclic deposits (parasequences) and important karstic surfaces capping some of the cycles. Peritidal cycles (shallow subtidal facies capped by tidal-¯at laminites or fenestrate limestones) are regressive- and transgressive-prone (upward-deepening followed by upward-shallowing facies trends). Subtidal cycles are of two types and indicate incomplete shallowing. Submerged subtidal cycles are composed of deeper subtidal facies overlain by shallow subtidal facies. Exposed subtidal cycles consist of deeper subtidal facies overlain by shallow subtidal facies that are capped by features indicative of prolonged subaerial exposure. Subtidal facies occur characteristically in the Jurassic, while peritidal cycles are typical for the Lower Cretaceous of the region. Within the foraminiferal and dasyclad algal biostratigraphic framework, four karst breccia levels are recognized as the boundaries of major second-order cycles, introduced for the ®rst time in this study. These levels correspond to the Kimmeridgian±Portlandian boundary, mid-Early Valanginian, mid-Early Aptian and mid-Cenomanian and represent important sea level falls which affected the distribution of foraminiferal fauna and dasyclad ¯ora of the Taurus carbonate platform. Within the Kimmeridgian±Cenomanian interval 26 third-order sequences (types 1 and 2) are recognized. These sequences are the records of eustatic sea level ¯uctuations rather than the records of local tectonic events because the boundaries of the sequences representing 1±4 Ma intervals are correlative with global sea level falls. Third-order sequences and metre-scale cyclic deposits are the major units used for long-distance, high-resolution sequence stratigraphic correlation in the Western Taurides. Metre-scale cyclic deposits (parasequences) in the Cretaceous show genetical stacking patterns within third-order sequences and correspond to fourth-order sequences representing 100±200 ka. These cycles are possibly the E2 signal (126 ka) of the orbital eccentricity cycles of the Milankovitch band. The slight deviation of values, calculated for parasequences, from the mean value of eccentricity cycles can be explained by the currently imprecise geochronology established in the Cretaceous and missed sea level oscillations when the platform lay above fluctuating sea level.

Comparison of monte carlo collimator transport methods for photon treatment planning in radiotherapy

The aim of this work was a Monte Carlo (MC) based investigation of the impact of different radiation transport methods in collimators of a linear accelerator on photon beam characteristics, dose distributions, and efficiency. Thereby it is investigated if it is possible to use different simplifications in the radiation transport for some clinical situations in order to save calculation time.

Avoiding misclassification bias with the traditional Charnley classification: rationale for a fourth Charnley class BB

Long-term follow up of patients with total hip arthroplasty (THA) revealed a marked deterioration of walking capacities in Charnley class B after postoperative year 4. We hypothesized that a specific group of patients, namely those with unilateral hip arthroplasty and an untreated but affected contralateral hip was responsible for this observation. Therefore, we conducted a study taking into consideration the two subclasses that make up Charnley class B: patients with unilateral THA and contralateral hip disease and patients with bilateral THA. A sample of 15,160 patients with 35,773 follow ups that were prospectively collected over 10 years was evaluated. The sample was categorized into four classes according to a new modified Charnley classification. Annual analyses of the proportion of patients with ambulation longer than 60 min were conducted. The traditionally labeled Charnley class B consists of two very different patient groups with respect to their walking capacities. Those with unilateral THA and contralateral hip disease have underaverage walking capacities and a deterioration of ambulation beginning 3 to 4 years after surgery. Those with bilateral THA have stable overaverage walking capacities similar to Charnley class A. An extension of the traditional Charnley classification is proposed, taking into account the two different patient groups in Charnley class B. The new fourth Charnley class consists of patients with bilateral THA and was labeled BB in order to express the presence of two artificial hip joints and to preserve the traditional classification A through C.

[Methods for investigating Trichinella infections in domestic and wild animals]

Trichinellosis is an important parasitic zoonosis that is caused by the intracellular nematode Trichinella spp.. Infection of humans occurs through consumption of raw (or undercooked) meat containing infectious larvae. In Europe, meat from pork, horse, and wild boar have been identified as most important sources of Trichinella infections in humans. In Switzerland, both the domestic pig and wild boar population are considered free of Trichinella. Conversely, Swiss foxes, lynxs and recently a wolf were found to be infected, the species identified in these animals was always referred to as Trichinella britovi. Although this species rarely infects pork and, compared to Trichinella spiralis, only causes reduced pathogenic effects in humans, the basic presence of Trichinella in Switzerland cannot be neglegted. This fact has gained increasing importance since the responsible authorities in the European Union (EU) are preparing regulations for the official Trichinella-control in meat in order to improve food safety for consumers. These regulations will be implemented as a consequence of the recent association of east European countries with the EU. This new legislation particularly takes into account, that in the past by far most cases of human trichinellosis in the EU were due to consumption of imported east European meat.Within the framework of the bilateral agreements of Switzerland with the EU, the Swiss veterinary public health authorities will have to comply with the foreseen EU regulations. Although diagnostic methods for the direct demonstation of Trichinella in pork meat are already routine practice in several Swiss abattoirs, the implementation of a meat control program for Trichinella for the entire slaughter pig population of the country would lead to an enormous increase in costs for the administration and will require an increased infrastructure in veterinary services. In order to find a reduced testing format for monitoring Trichinella infections in Swiss pork, an infection risk-oriented survey strategy is currently evaluated. In the present article, this minimized survey strategy is discussed regarding its compatibility with the EU regulations laying down rules for the official control of meat for Trichinella.

Influence of different storage methods on laser fluorescence values: a two-year study

The aim of this study was to assess the influence on the infrared laser fluorescence response of some storage methods commonly used in dental research. Forty extracted permanent teeth, selected from a pool of frozen teeth, were divided into four groups of 10. Three groups were stored at 4 degrees C in 1% chloramine, 10% formalin or 0.02% thymol solution. The fourth group was stored at -20 degrees C (no storage solution added). Fluorescence measurements were performed at 14, 77, 113, 168, 232, 486 and 737 days. After 2 years, significant decreases in fluorescence (p<0.01) for the samples in formalin (-60%), chloramine (-72%) and thymol (-54%) were observed. The frozen teeth showed a slight but non-significant increase in fluorescence of 5% (p>0.01). Storing solutions have a significant influence on the fluorescence yield. Samples used for in vitro purposes stored frozen do not significantly change their fluorescence response. Thus, cut-off values obtained under the latter conditions could be extrapolated to the in vivo situation.

Strabismus and employment: the opinion of headhunters

Abstract. Purpose: The aim of this study was to determine the influence of strabismus on an individual's ability to find employment based on the opinion of Swiss headhunters. Methods: Forty Swiss headhunters retrieved from a Swiss online telephone directory were interviewed using a validated questionnaire in order to determine if strabismus would have an impact on a person's ability to find employment. Photographs of a strabismic man and woman as well as of other computer-generated facial anomalies could be downloaded from the Internet during the interview. Results: 72.5% of headhunters judged that strabismic individuals would have more difficulties in finding a job than orthotropic persons. These difficulties were judged to be stronger in women than in men (P = 0.006), and in exotropic compared to esotropic persons (P = 0.01). Asked about seven facial anomalies, exotropia was found to have the third and esotropia the fourth (women) or fifth (men) most strongly negative impact on finding employment, after having strong acne and a visible missing tooth. Headhunters judged that persons with strabismus are significantly perceived as less attractive and less intelligent by potential employers. Conclusion: Visible strabismus negatively influences individuals' ability to find a job and therefore has an impact on their economic status. Successful strabismus surgery realigns the visual axes, producing a normal facial appearance and therefore eliminating the negative impact of strabismus on employability.

Age modeling of young non-varved lake sediments: methods and limits. Examples from two lakes in Central Chile

High-resolution and highly precise age models for recent lake sediments (last 100–150 years) are essential for quantitative paleoclimate research. These are particularly important for sedimentological and geochemical proxies, where transfer functions cannot be established and calibration must be based upon the relation of sedimentary records to instrumental data. High-precision dating for the calibration period is most critical as it determines directly the quality of the calibration statistics. Here, as an example, we compare radionuclide age models obtained on two high-elevation glacial lakes in the Central Chilean Andes (Laguna Negra: 33°38′S/70°08′W, 2,680 m a.s.l. and Laguna El Ocho: 34°02′S/70°19′W, 3,250 m a.s.l.). We show the different numerical models that produce accurate age-depth chronologies based on 210Pb profiles, and we explain how to obtain reduced age-error bars at the bottom part of the profiles, i.e., typically around the end of the 19th century. In order to constrain the age models, we propose a method with five steps: (i) sampling at irregularly-spaced intervals for 226Ra, 210Pb and 137Cs depending on the stratigraphy and microfacies, (ii) a systematic comparison of numerical models for the calculation of 210Pb-based age models: constant flux constant sedimentation (CFCS), constant initial concentration (CIC), constant rate of supply (CRS) and sediment isotope tomography (SIT), (iii) numerical constraining of the CRS and SIT models with the 137Cs chronomarker of AD 1964 and, (iv) step-wise cross-validation with independent diagnostic environmental stratigraphic markers of known age (e.g., volcanic ash layer, historical flood and earthquakes). In both examples, we also use airborne pollutants such as spheroidal carbonaceous particles (reflecting the history of fossil fuel emissions), excess atmospheric Cu deposition (reflecting the production history of a large local Cu mine), and turbidites related to historical earthquakes. Our results show that the SIT model constrained with the 137Cs AD 1964 peak performs best over the entire chronological profile (last 100–150 years) and yields the smallest standard deviations for the sediment ages. Such precision is critical for the calibration statistics, and ultimately, for the quality of the quantitative paleoclimate reconstruction. The systematic comparison of CRS and SIT models also helps to validate the robustness of the chronologies in different sections of the profile. Although surprisingly poorly known and under-explored in paleolimnological research, the SIT model has a great potential in paleoclimatological reconstructions based on lake sediments

An insulin infusion advisory system for type 1 diabetes patients based on non-linear model predictive control methods

In this paper, an Insulin Infusion Advisory System (IIAS) for Type 1 diabetes patients, which use insulin pumps for the Continuous Subcutaneous Insulin Infusion (CSII) is presented. The purpose of the system is to estimate the appropriate insulin infusion rates. The system is based on a Non-Linear Model Predictive Controller (NMPC) which uses a hybrid model. The model comprises a Compartmental Model (CM), which simulates the absorption of the glucose to the blood due to meal intakes, and a Neural Network (NN), which simulates the glucose-insulin kinetics. The NN is a Recurrent NN (RNN) trained with the Real Time Recurrent Learning (RTRL) algorithm. The output of the model consists of short term glucose predictions and provides input to the NMPC, in order for the latter to estimate the optimum insulin infusion rates. For the development and the evaluation of the IIAS, data generated from a Mathematical Model (MM) of a Type 1 diabetes patient have been used. The proposed control strategy is evaluated at multiple meal disturbances, various noise levels and additional time delays. The results indicate that the implemented IIAS is capable of handling multiple meals, which correspond to realistic meal profiles, large noise levels and time delays.

An efficient simulation algorithm for the generalized von Mises distribution of order two

In this article we propose an exact efficient simulation algorithm for the generalized von Mises circular distribution of order two. It is an acceptance-rejection algorithm with a piecewise linear envelope based on the local extrema and the inflexion points of the generalized von Mises density of order two. We show that these points can be obtained from the roots of polynomials and degrees four and eight, which can be easily obtained by the methods of Ferrari and Weierstrass. A comparative study with the von Neumann acceptance-rejection, with the ratio-of-uniforms and with a Markov chain Monte Carlo algorithms shows that this new method is generally the most efficient.

Multiscale methods for shape constraints in deconvolution: Confidence statements for qualitative features

We derive multiscale statistics for deconvolution in order to detect qualitative features of the unknown density. An important example covered within this framework is to test for local monotonicity on all scales simultaneously. We investigate the moderately ill-posed setting, where the Fourier transform of the error density in the deconvolution model is of polynomial decay. For multiscale testing, we consider a calibration, motivated by the modulus of continuity of Brownian motion. We investigate the performance of our results from both the theoretical and simulation based point of view. A major consequence of our work is that the detection of qualitative features of a density in a deconvolution problem is a doable task, although the minimax rates for pointwise estimation are very slow.