Line search multilevel optimization as computational methods for dense optical flow

We evaluate the performance of different optimization techniques developed in the context of optical flowcomputation with different variational models. In particular, based on truncated Newton methods (TN) that have been an effective approach for large-scale unconstrained optimization, we develop the use of efficient multilevel schemes for computing the optical flow. More precisely, we evaluate the performance of a standard unidirectional multilevel algorithm - called multiresolution optimization (MR/OPT), to a bidrectional multilevel algorithm - called full multigrid optimization (FMG/OPT). The FMG/OPT algorithm treats the coarse grid correction as an optimization search direction and eventually scales it using a line search. Experimental results on different image sequences using four models of optical flow computation show that the FMG/OPT algorithm outperforms both the TN and MR/OPT algorithms in terms of the computational work and the quality of the optical flow estimation.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

An algorithm for semi-infinite polynomial optimization

"Vegeu el resum a l'inici del document del fitxer adjunt."

A linear optimization technique for graph pebbling

Graph pebbling is a network model for studying whether or not a given supply of discrete pebbles can satisfy a given demand via pebbling moves. A pebbling move across an edge of a graph takes two pebbles from one endpoint and places one pebble at the other endpoint; the other pebble is lost in transit as a toll. It has been shown that deciding whether a supply can meet a demand on a graph is NP-complete. The pebbling number of a graph is the smallest t such that every supply of t pebbles can satisfy every demand of one pebble. Deciding if the pebbling number is at most k is NP 2 -complete. In this paper we develop a tool, called theWeight Function Lemma, for computing upper bounds and sometimes exact values for pebbling numbers with the assistance of linear optimization. With this tool we are able to calculate the pebbling numbers of much larger graphs than in previous algorithms, and much more quickly as well. We also obtain results for many families of graphs, in many cases by hand, with much simpler and remarkably shorter proofs than given in previously existing arguments (certificates typically of size at most the number of vertices times the maximum degree), especially for highly symmetric graphs. Here we apply theWeight Function Lemma to several specific graphs, including the Petersen, Lemke, 4th weak Bruhat, Lemke squared, and two random graphs, as well as to a number of infinite families of graphs, such as trees, cycles, graph powers of cycles, cubes, and some generalized Petersen and Coxeter graphs. This partly answers a question of Pachter, et al., by computing the pebbling exponent of cycles to within an asymptotically small range. It is conceivable that this method yields an approximation algorithm for graph pebbling.

A parameter-free approach for solving combinatorial optimization problems through biased randomization of efficient heuristics

This paper discusses the use of probabilistic or randomized algorithms for solving combinatorial optimization problems. Our approach employs non-uniform probability distributions to add a biased random behavior to classical heuristics so a large set of alternative good solutions can be quickly obtained in a natural way and without complex conguration processes. This procedure is especially useful in problems where properties such as non-smoothness or non-convexity lead to a highly irregular solution space, for which the traditional optimization methods, both of exact and approximate nature, may fail to reach their full potential. The results obtained are promising enough to suggest that randomizing classical heuristics is a powerful method that can be successfully applied in a variety of cases.

On the stability of the optimal value and the optimal set in optimization problems

The paper develops a stability theory for the optimal value and the optimal set mapping of optimization problems posed in a Banach space. The problems considered in this paper have an arbitrary number of inequality constraints involving lower semicontinuous (not necessarily convex) functions and one closed abstract constraint set. The considered perturbations lead to problems of the same type as the nominal one (with the same space of variables and the same number of constraints), where the abstract constraint set can also be perturbed. The spaces of functions involved in the problems (objective and constraints) are equipped with the metric of the uniform convergence on the bounded sets, meanwhile in the space of closed sets we consider, coherently, the Attouch-Wets topology. The paper examines, in a unified way, the lower and upper semicontinuity of the optimal value function, and the closedness, lower and upper semicontinuity (in the sense of Berge) of the optimal set mapping. This paper can be seen as a second part of the stability theory presented in [17], where we studied the stability of the feasible set mapping (completed here with the analysis of the Lipschitz-like property).

Artificial Snow Optimization in Winter Sport Destinations Using a Multi-agent Simulation

This paper presents the Juste-Neige system for predicting the snow height on the ski runs of a resort using a multi-agent simulation software. Its aim is to facilitate snow cover management in order to i) reduce the production cost of artificial snow and to improve the profit margin for the companies managing the ski resorts; and ii) to reduce the water and energy consumption, and thus to reduce the environmental impact, by producing only the snow needed for a good skiing experience. The software provides maps with the predicted snow heights for up to 13 days. On these maps, the areas most exposed to snow erosion are highlighted. The software proceeds in three steps: i) interpolation of snow height measurements with a neural network; ii) local meteorological forecasts for every ski resort; iii) simulation of the impact caused by skiers using a multi-agent system. The software has been evaluated in the Swiss ski resort of Verbier and provides useful predictions.

Copper, selenium, zinc, and thiamine balances during continuous venovenous hemodiafiltration in critically ill patients.

BACKGROUND: Acute renal failure is a serious complication in critically ill patients and frequently requires renal replacement therapy, which alters trace element and vitamin metabolism. OBJECTIVE: The objective was to study trace element balances during continuous renal replacement therapy (CRRT) in intensive care patients. DESIGN: In a prospective randomized crossover trial, patients with acute renal failure received CRRT with either sodium bicarbonate (Bic) or sodium lactate (Lac) as a buffering agent over 2 consecutive 24-h periods. Copper, selenium, zinc, and thiamine were measured with highly sensitive analytic methods in plasma, replacement solutions, and effluent during 8-h periods. Balances were calculated as the difference between fluids administered and effluent losses and were compared with the recommended intakes (RI) from parenteral nutrition. RESULTS: Nineteen sessions were conducted in 11 patients aged 65 +/- 10 y. Baseline plasma concentrations of copper were normal, whereas those of selenium and zinc were below reference ranges; glutathione peroxidase was in the lower range of normal. The replacement solutions contained no detectable copper, 0.01 micromol Se/L (Bic and Lac), and 1.42 (Bic) and 0.85 (Lac) micromol Zn/L. Micronutrients were detectable in all effluents, and losses were stable in each patient; no significant differences were found between the Bic and Lac groups. The 24-h balances were negative for selenium (-0.97 micromol, or 2 times the daily RI), copper (-6.54 micromol, or 0.3 times the daily RI), and thiamine (-4.12 mg, or 1.5 times the RI) and modestly positive for zinc (20.7 micromol, or 0.2 times the RI). CONCLUSIONS: CRRT results in significant losses and negative balances of selenium, copper, and thiamine, which contribute to low plasma concentrations. Prolonged CRRT is likely to result in selenium and thiamine depletion despite supplementation at recommended amounts.

Establishment and characterization of a new continuous cell line from Lutzomyia longipalpis (Diptera: Psychodidae) and its susceptibility to infections with arboviruses and Leishmania chagasi

Embryonic tissue explants of the sand fly Lutzomyia longipalpis (Lutz & Neiva 1912) the main vector of Leishmania chagasi (Cunha and Chagas), were used to obtain a continuous cell line (Lulo). The tissues were seeded in MM/VP12 medium and these were incubated at 28ºC. The first subculture was obtained 45 days after explanting and 96 passages have been made to date. Lulo is composed of epithelioid cells, showed a 0.04 generations/hour exponential growth rate and population doubling time at 24.7 h. The cell line isoenzymatic profiles were determined by using PGI, PGM, MPI and 6-PGDH systems, coinciding with patterns obtained from the same species and colony's pupae and adults. The species karyotype characteristics were recognized (2n = 8), in which pair 1 is subtelocentric and pairs 2, 3 and 4 are metacentric. Lulo was free from bacterial, fungal, mycoplasmic and viral infection. Susceptibility to five arbovirus was determined, the same as Lulo interaction with Leishmania promastigotes.

Continuous monitoring of liver oxygenation with near infrared spectroscopy during naso-gastric tube feeding in neonates.

Near infrared spectroscopy (NIRS) is a non-invasive method of estimating the haemoglobin concentration changes in certain tissues. It is frequently used to monitor oxygenation of the brain in neonates. At present it is not clear whether near infrared spectroscopy of other organs (e.g. the liver as a corresponding site in the splanchnic region, which reacts very sensitively to haemodynamic instability) provides reliable values on their tissue oxygenation. The aim of the study was to test near infrared spectroscopy by measuring known physiologic changes in tissue oxygenation of the liver in newborn infants during and after feeding via a naso-gastric tube. The test-retest variability of such measurements was also determined. On 28 occasions in 25 infants we measured the tissue oxygenation index (TOI) of the liver and the brain continuously before, during and 30 minutes after feeding via a gastric tube. Simultaneously we measured arterial oxygen saturation (SaO2), heart rate (HR) and mean arterial blood pressure (MAP). In 10 other newborn infants we performed a test-retest analysis of the liver tissue oxygenation index to estimate the variability in repeated intra-individual measurements. The tissue oxygenation index of the liver increased significantly from 56.7 +/- 7.5% before to 60.3 +/- 5.6% after feeding (p < 0.005), and remained unchanged for the next 30 minutes. The tissue oxygenation index of the brain (62.1 +/- 9.7%), SaO2 (94.4 +/- 7.1%), heart rate (145 +/- 17.3 min-1) and mean arterial blood pressure (52.8 +/- 10.2 mm Hg) did not change significantly. The test-retest variability for intra-individual measurements was 2.7 +/- 2.1%. After bolus feeding the tissue oxygenation index of the liver increased as expected. This indicates that near infrared spectroscopy is suitable for monitoring changes in tissue oxygenation of the liver in newborn infants.

A new continuous cell line from the mosquito Psorophora confinnis (Diptera: Culicidae) and its susceptibility to infections with some arboviruses

A new cell line, PC-0199-BR, was established from embryonated eggs of the mosquito Psorophora confinnis. To date (September 2000) it has had 62 continuous passages. This is the first report of a cell line of mosquitoes belonging to the genus Psorophora. Cell growth initially was achieved in the MM/VP12 medium, supplemented with 20% fetal bovine serum; however, the subcultures were later adapted to Grace's medium with 10% fetal bovine serum. Cell morphology in the primary cultures was heterogeneous; but later in the established cell line, the predominant cell type was epithelioid. Cultured cells were predominantly diploid (2n=6); however, chromosome abnormalities were observed in a small proportion of the cells in later passages. C and G band patterns were also determined in the karyotype. The cell line isozyme profiles coincided with pupae and adult samples of the species taken from the same colony. A preliminary arbovirus susceptibility study for the cell line was undertaken. No evidence was observed of contamination of the cell line with bacteria, fungi or mycoplasma.

Optimization of spatial resolution for peripheral magnetic resonance angiography.

RATIONALE AND OBJECTIVES: To determine optimum spatial resolution when imaging peripheral arteries with magnetic resonance angiography (MRA). MATERIALS AND METHODS: Eight vessel diameters ranging from 1.0 to 8.0 mm were simulated in a vascular phantom. A total of 40 three-dimensional flash MRA sequences were acquired with incremental variations of fields of view, matrix size, and slice thickness. The accurately known eight diameters were combined pairwise to generate 22 "exact" degrees of stenosis ranging from 42% to 87%. Then, the diameters were measured in the MRA images by three independent observers and with quantitative angiography (QA) software and used to compute the degrees of stenosis corresponding to the 22 "exact" ones. The accuracy and reproducibility of vessel diameter measurements and stenosis calculations were assessed for vessel size ranging from 6 to 8 mm (iliac artery), 4 to 5 mm (femoro-popliteal arteries), and 1 to 3 mm (infrapopliteal arteries). Maximum pixel dimension and slice thickness to obtain a mean error in stenosis evaluation of less than 10% were determined by linear regression analysis. RESULTS: Mean errors on stenosis quantification were 8.8% +/- 6.3% for 6- to 8-mm vessels, 15.5% +/- 8.2% for 4- to 5-mm vessels, and 18.9% +/- 7.5% for 1- to 3-mm vessels. Mean errors on stenosis calculation were 12.3% +/- 8.2% for observers and 11.4% +/- 15.1% for QA software (P = .0342). To evaluate stenosis with a mean error of less than 10%, maximum pixel surface, the pixel size in the phase direction, and the slice thickness should be less than 1.56 mm2, 1.34 mm, 1.70 mm, respectively (voxel size 2.65 mm3) for 6- to 8-mm vessels; 1.31 mm2, 1.10 mm, 1.34 mm (voxel size 1.76 mm3), for 4- to 5-mm vessels; and 1.17 mm2, 0.90 mm, 0.9 mm (voxel size 1.05 mm3) for 1- to 3-mm vessels. CONCLUSION: Higher spatial resolution than currently used should be selected for imaging peripheral vessels.

A new automated method versus continuous positive airway pressure method for measuring pressure-volume curves in patients with acute lung injury

Objective: To compare pressure–volume (P–V) curves obtained with the Galileo ventilator with those obtained with the CPAP method in patients with ALI or ARDS receiving mechanical ventilation. P–V curves were fitted to a sigmoidal equation with a mean R2 of 0.994 ± 0.003. Lower (LIP) and upper inflection (UIP), and deflation maximum curvature (PMC) points calculated from the fitted variables showed a good correlation between methods with high intraclass correlation coefficients. Bias and limits of agreement for LIP, UIP and PMC obtained with the two methods in the same patient were clinically acceptable.