Bit-index sort: a fast non-comparison integer sorting algorithm for permutations

This paper describes a fast integer sorting algorithm, herein referred as Bit-index sort, which is a non-comparison sorting algorithm for partial per-mutations, with linear complexity order in execution time. Bit-index sort uses a bit-array to classify input sequences of distinct integers, and exploits built-in bit functions in C compilers supported by machine hardware to retrieve the ordered output sequence. Results show that Bit-index sort outperforms in execution time to quicksort and counting sort algorithms. A parallel approach for Bit-index sort using two simultaneous threads is included, which obtains speedups up to 1.6.

Differences between the CME fronts tracked by an expert, an automated algorithm, and the Solar Stormwatch project

Observations from the Heliospheric Imager (HI) instruments aboard the twin STEREO spacecraft have enabled the compilation of several catalogues of coronal mass ejections (CMEs), each characterizing the propagation of CMEs through the inner heliosphere. Three such catalogues are the Rutherford Appleton Laboratory (RAL)-HI event list, the Solar Stormwatch CME catalogue, and, presented here, the J-tracker catalogue. Each catalogue uses a different method to characterize the location of CME fronts in the HI images: manual identification by an expert, the statistical reduction of the manual identifications of many citizen scientists, and an automated algorithm. We provide a quantitative comparison of the differences between these catalogues and techniques, using 51 CMEs common to each catalogue. The time-elongation profiles of these CME fronts are compared, as are the estimates of the CME kinematics derived from application of three widely used single-spacecraft-fitting techniques. The J-tracker and RAL-HI profiles are most similar, while the Solar Stormwatch profiles display a small systematic offset. Evidence is presented that these differences arise because the RAL-HI and J-tracker profiles follow the sunward edge of CME density enhancements, while Solar Stormwatch profiles track closer to the antisunward (leading) edge. We demonstrate that the method used to produce the time-elongation profile typically introduces more variability into the kinematic estimates than differences between the various single-spacecraft-fitting techniques. This has implications for the repeatability and robustness of these types of analyses, arguably especially so in the context of space weather forecasting, where it could make the results strongly dependent on the methods used by the forecaster.

A Genetic Algorithm for the selection of structural MRI features for classification of Mild Cognitive Impairment and Alzheimer's Disease

This work investigates the problem of feature selection in neuroimaging features from structural MRI brain images for the classification of subjects as healthy controls, suffering from Mild Cognitive Impairment or Alzheimer’s Disease. A Genetic Algorithm wrapper method for feature selection is adopted in conjunction with a Support Vector Machine classifier. In very large feature sets, feature selection is found to be redundant as the accuracy is often worsened when compared to an Support Vector Machine with no feature selection. However, when just the hippocampal subfields are used, feature selection shows a significant improvement of the classification accuracy. Three-class Support Vector Machines and two-class Support Vector Machines combined with weighted voting are also compared with the former and found more useful. The highest accuracy achieved at classifying the test data was 65.5% using a genetic algorithm for feature selection with a three-class Support Vector Machine classifier.

Evaluation of a polynya flux model by means of thermal infrared satellite estimates

We test the ability of a two-dimensional flux model to simulate polynya events with narrow open-water zones by comparing model results to ice-thickness and ice-production estimates derived from thermal infrared Moderate Resolution Imaging Spectroradiometer (MODIS) observations in conjunction with an atmospheric dataset. Given a polynya boundary and an atmospheric dataset, the model correctly reproduces the shape of an 11 day long event, using only a few simple conservation laws. Ice production is slightly overestimated by the model, owing to an underestimated ice thickness. We achieved best model results with the consolidation thickness parameterization developed by Biggs and others (2000). Observed regional discrepancies between model and satellite estimates might be a consequence of the missing representation of the dynamic of the thin-ice thickening (e.g. rafting). We conclude that this simplified polynya model is a valuable tool for studying polynya dynamics and estimating associated fluxes of single polynya events.

Perception of eating practices and stages of change among Brazilian adolescents

Objective. To evaluate the perception of eating practices and the stages of change among adolescents. Methods. Cross-sectional study involving a representative sample of 390 adolescents from 11 public schools in the city of Piracicaba, Brazil, in 2004. Food consumption was identified by a food frequency questionnaire and the perception of eating practices evaluation was conducted by comparing food consumption and individual classification of healthy aspects of the diet. The participants were classified within stages of change by means of a specific algorithm. A reclassification within new stages of change was proposed to identify adolescents with similar characteristics regarding food consumption and perception. Results. Low consumption of fruit and vegetables and high consumption of sweets and fats were identified. More than 44% of the adolescents had a mistaken perception of their diet. A significant relationship between the stages of change and food consumption was observed. The reclassification among stages of change, through including the pseudo-maintenance and non-reflective action stages was necessary, considering the high proportion of adolescents who erroneously classified their diets as healthy. Conclusion. Classification of the adolescents into stages of change, together with consumption and perception data, enabled identification of groups at risk, in accordance with their inadequate dietary habits and non-recognition of such habits. (C) 2009 Elsevier Inc. All rights reserved.

Time-frequency analysis of acoustic noise produced by breaking of crisp biscuits

This paper demonstrates by means of joint time-frequency analysis that the acoustic noise produced by the breaking of biscuits is dependent on relative humidity and water activity. It also shows that the time-frequency coefficients calculated using the adaptive Gabor transformation algorithm is dependent on the period of time a biscuit is exposed to humidity. This is a new methodology that can be used to assess the crispness of crisp foods. (c) 2007 Elsevier Ltd. All rights reserved.

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.

Deletion of Tumor Necrosis Factor-alpha Receptor 1 (TNFR1) Protects against Diet-induced Obesity by Means of Increased Thermogenesis

In diet-induced obesity, hypothalamic and systemic inflammatory factors trigger intracellular mechanisms that lead to resistance to the main adipostatic hormones, leptin and insulin. Tumor necrosis factor-alpha (TNF-alpha) is one of the main inflammatory factors produced during this process and its mechanistic role as an inducer of leptin and insulin resistance has been widely investigated. Most of TNF-alpha inflammatory signals are delivered by TNF receptor 1 (R1); however, the role played by this receptor in the context of obesity-associated inflammation is not completely known. Here, we show that TNFR1 knock-out (TNFR1 KO) mice are protected from diet-induced obesity due to increased thermogenesis. Under standard rodent chow or a high-fat diet, TNFR1 KO gain significantly less body mass despite increased caloric intake. Visceral adiposity and mean adipocyte diameter are reduced and blood concentrations of insulin and leptin are lower. Protection from hypothalamic leptin resistance is evidenced by increased leptin-induced suppression of food intake and preserved activation of leptin signal transduction through JAK2, STAT3, and FOXO1. Under the high-fat diet, TNFR1 KO mice present a significantly increased expression of the thermogenesis-related neurotransmitter, TRH. Further evidence of increased thermogenesis includes increased O(2) consumption in respirometry measurements, increased expressions of UCP1 and UCP3 in brown adipose tissue and skeletal muscle, respectively, and increased O(2) consumption by isolated skeletal muscle fiber mitochondria. This demonstrates that TNF-alpha signaling through TNFR1 is an important mechanism involved in obesity-associated defective thermogenesis.

The role of endothelin pathway in modulation of airway reactivity to methacholine in C57Bl/6 and BALB/c mice

Endothelin peptides have been shown to increase cholinergic neurotransmission in the airway. Genetic differences in airway responsiveness to methacholine where reported in mice. The present study compared the airway reactivity to methacholine in C57Bl/6 and BALB/c mice, the involvement of endothelin on this reactivity and endothelin levels in lung homogenates. Whole airway reactivity was analyzed by means of an isolated lung preparation where lungs were perfused through the trachea with warm gassed Krebs solution at 5 ml/min, and changes in perfusion pressure triggered by methacholine at increasing bolus doses (0.1-100 mu g) were recorded. We found that the maximal airway response to methacholine was much greater in C57Bl/6 than in BALB/c (Emax 34 +/- 2 vs 12 +/- 1 cmH(2)O, respectively). Bosentan (mixed endothelin A/B receptor antagonist; 10 mg/kg, i.p., 30 min before sacrifice) reduced lung responsiveness to methacholine in C57Bl/6 (58% at EC50 level) but had no effect in BALB/c mouse strain. This effect seems to be mediated by the endothelin ETA receptor since it was significantly reduced by the selective endothelin ETA receptor antagonist, BQ 123. Immunoreactive endothelin levels were higher in C57Bl/6 than in BALB/c lungs (43 5 vs 19 +/- 5 pg/g of tissue). In conclusion, airway reactivity to methacholine and lung endothelins content varies markedly between C57Bl/6 and BALB/c strains. Endothelins upregulate lung responsiveness to methacholine only in C57Bl/6, an effect achieved through the endothelin ETA receptor. (C) 2008 Elsevier B.V. All rights reserved.

A geometric mass-preserving redistancing scheme for the level set function

In this paper we describe and evaluate a geometric mass-preserving redistancing procedure for the level set function on general structured grids. The proposed algorithm is adapted from a recent finite element-based method and preserves the mass by means of a localized mass correction. A salient feature of the scheme is the absence of adjustable parameters. The algorithm is tested in two and three spatial dimensions and compared with the widely used partial differential equation (PDE)-based redistancing method using structured Cartesian grids. Through the use of quantitative error measures of interest in level set methods, we show that the overall performance of the proposed geometric procedure is better than PDE-based reinitialization schemes, since it is more robust with comparable accuracy. We also show that the algorithm is well-suited for the highly stretched curvilinear grids used in CFD simulations. Copyright (C) 2010 John Wiley & Sons, Ltd.

Infeasibility handling in genetic algorithm using nested domains for production planning

In this paper we present a genetic algorithm with new components to tackle capacitated lot sizing and scheduling problems with sequence dependent setups that appear in a wide range of industries, from soft drink bottling to food manufacturing. Finding a feasible solution to highly constrained problems is often a very difficult task. Various strategies have been applied to deal with infeasible solutions throughout the search. We propose a new scheme of classifying individuals based on nested domains to determine the solutions according to the level of infeasibility, which in our case represents bands of additional production hours (overtime). Within each band, individuals are just differentiated by their fitness function. As iterations are conducted, the widths of the bands are dynamically adjusted to improve the convergence of the individuals into the feasible domain. The numerical experiments on highly capacitated instances show the effectiveness of this computational tractable approach to guide the search toward the feasible domain. Our approach outperforms other state-of-the-art approaches and commercial solvers. (C) 2009 Elsevier Ltd. All rights reserved.

A Mass-Conserving Algorithm for Dynamical Lubrication Problems With Cavitation

A numerical algorithm for fully dynamical lubrication problems based on the Elrod-Adams formulation of the Reynolds equation with mass-conserving boundary conditions is described. A simple but effective relaxation scheme is used to update the solution maintaining the complementarity conditions on the variables that represent the pressure and fluid fraction. The equations of motion are discretized in time using Newmark`s scheme, and the dynamical variables are updated within the same relaxation process just mentioned. The good behavior of the proposed algorithm is illustrated in two examples: an oscillatory squeeze flow (for which the exact solution is available) and a dynamically loaded journal bearing. This article is accompanied by the ready-to-compile source code with the implementation of the proposed algorithm. [DOI: 10.1115/1.3142903]

Front tracking with moving-least-squares surfaces

The representation of interfaces by means of the algebraic moving-least-squares (AMLS) technique is addressed. This technique, in which the interface is represented by an unconnected set of points, is interesting for evolving fluid interfaces since there is]to surface connectivity. The position of the surface points can thus be updated without concerns about the quality of any surface triangulation. We introduce a novel AMLS technique especially designed for evolving-interfaces applications that we denote RAMLS (for Robust AMLS). The main advantages with respect to previous AMLS techniques are: increased robustness, computational efficiency, and being free of user-tuned parameters. Further, we propose a new front-tracking method based on the Lagrangian advection of the unconnected point set that defines the RAMLS surface. We assume that a background Eulerian grid is defined with some grid spacing h. The advection of the point set makes the surface evolve in time. The point cloud can be regenerated at any time (in particular, we regenerate it each time step) by intersecting the gridlines with the evolved surface, which guarantees that the density of points on the surface is always well balanced. The intersection algorithm is essentially a ray-tracing algorithm, well-studied in computer graphics, in which a line (ray) is traced so as to detect all intersections with a surface. Also, the tracing of each gridline is independent and can thus be performed in parallel. Several tests are reported assessing first the accuracy of the proposed RAMLS technique, and then of the front-tracking method based on it. Comparison with previous Eulerian, Lagrangian and hybrid techniques encourage further development of the proposed method for fluid mechanics applications. (C) 2008 Elsevier Inc. All rights reserved.

Fast heuristics for the Steiner tree problem with revenues, budget and hop constraints

This article describes and compares three heuristics for a variant of the Steiner tree problem with revenues, which includes budget and hop constraints. First, a greedy method which obtains good approximations in short computational times is proposed. This initial solution is then improved by means of a destroy-and-repair method or a tabu search algorithm. Computational results compare the three methods in terms of accuracy and speed. (C) 2007 Elsevier B.V. All rights reserved.