An iterative Multiscale Finite-Volume algorithm converging to the exact solution

The multiscale finite volume (MsFV) method has been developed to efficiently solve large heterogeneous problems (elliptic or parabolic); it is usually employed for pressure equations and delivers conservative flux fields to be used in transport problems. The method essentially relies on the hypothesis that the (fine-scale) problem can be reasonably described by a set of local solutions coupled by a conservative global (coarse-scale) problem. In most cases, the boundary conditions assigned for the local problems are satisfactory and the approximate conservative fluxes provided by the method are accurate. In numerically challenging cases, however, a more accurate localization is required to obtain a good approximation of the fine-scale solution. In this paper we develop a procedure to iteratively improve the boundary conditions of the local problems. The algorithm relies on the data structure of the MsFV method and employs a Krylov-subspace projection method to obtain an unconditionally stable scheme and accelerate convergence. Two variants are considered: in the first, only the MsFV operator is used; in the second, the MsFV operator is combined in a two-step method with an operator derived from the problem solved to construct the conservative flux field. The resulting iterative MsFV algorithms allow arbitrary reduction of the solution error without compromising the construction of a conservative flux field, which is guaranteed at any iteration. Since it converges to the exact solution, the method can be regarded as a linear solver. In this context, the schemes proposed here can be viewed as preconditioned versions of the Generalized Minimal Residual method (GMRES), with a very peculiar characteristic that the residual on the coarse grid is zero at any iteration (thus conservative fluxes can be obtained).

Proposal of a parallel architecture for a motion detection algorithm

This paper proposes a parallel architecture for estimation of the motion of an underwater robot. It is well known that image processing requires a huge amount of computation, mainly at low-level processing where the algorithms are dealing with a great number of data. In a motion estimation algorithm, correspondences between two images have to be solved at the low level. In the underwater imaging, normalised correlation can be a solution in the presence of non-uniform illumination. Due to its regular processing scheme, parallel implementation of the correspondence problem can be an adequate approach to reduce the computation time. Taking into consideration the complexity of the normalised correlation criteria, a new approach using parallel organisation of every processor from the architecture is proposed

Speculative parallelization of multipath radiosity algorithm

In computer graphics, global illumination algorithms take into account not only the light that comes directly from the sources, but also the light interreflections. This kind of algorithms produce very realistic images, but at a high computational cost, especially when dealing with complex environments. Parallel computation has been successfully applied to such algorithms in order to make it possible to compute highly-realistic images in a reasonable time. We introduce here a speculation-based parallel solution for a global illumination algorithm in the context of radiosity, in which we have taken advantage of the hierarchical nature of such an algorithm

A Monte Carlo-Based Fiber Tracking Algorithm using Diffusion Tensor MRI

Diffusion tensor magnetic resonance imaging, which measures directional information of water diffusion in the brain, has emerged as a powerful tool for human brain studies. In this paper, we introduce a new Monte Carlo-based fiber tracking approach to estimate brain connectivity. One of the main characteristics of this approach is that all parameters of the algorithm are automatically determined at each point using the entropy of the eigenvalues of the diffusion tensor. Experimental results show the good performance of the proposed approach

Aplicació del càlcul d’esforços en 2n ordre en estructures metàl·liques segons el CTE i NBE-EA-95

Aquest estudi se centra en el document del CTE que es refereix a les estructuresmetàl·liques d’acer, concretament en els càlculs a tenir en compte a l’hora de realitzar pòrtics d’estructura d’edificació amb perfils metàl·lics laminats. Aquest document és el DB-SE-A, que són les sigles de “Document Bàsic Seguretat Estructural Acer”. Aquest document introdueix algunes modificacions en el mètode de càlcul de les estructures metàl·liques respecte l’antiga norma.La idea d’aquest estudi és fer una comparació entre les dues normes – la vella i ja derogada, i la nova i vigent – i com tracten el càlcul d’esforços en segon ordre i el tractament de les imperfeccions.Com ja s’ha comentat aquest Document Bàsic de la Seguretat Estructural en Acer del CTE (DB-SE-A) substitueix l’antiga norma que es referia al càlcul d’estructura metàl·lica, la NBE-EA-95. En l’aplicació de les accions i el càlcul d’esforços s’utilitza del CTE el Document Bàsic dedicat a les accions en l’edificació (DB-SE-AE) i que substitueix l’antiga NBE-AE-88

Evolutionary histories of expanded peptidase families in Schistosoma mansoni

Schistosoma mansoni is one of the three main causative agents of human schistosomiasis, a major health problem with a vast socio-economic impact. Recent advances in the proteomic analysis of schistosomes have revealed that peptidases are the main virulence factors involved in the pathogenesis of this disease. In this context, evolutionary studies can be applied to identify peptidase families that have been expanded in genomes over time in response to different selection pressures. Using a phylogenomic approach, we searched for expanded endopeptidase families in the S. mansoni predicted proteome with the aim of contributing to the knowledge of such enzymes as potential therapeutic targets. We found three endopeptidase families that comprise leishmanolysins (metallopeptidase M8 family), cercarial elastases (serine peptidase S1 family) and cathepsin D proteins (aspartic peptidase A1 family). Our results suggest that the Schistosoma members of these families originated from successive gene duplication events in the parasite lineage after its diversification from other metazoans. Overall, critical residues are conserved among the duplicated genes/proteins. Furthermore, each protein family displays a distinct evolutionary history. Altogether, this work provides an evolutionary view of three S. mansoni peptidase families, which allows for a deeper understanding of the genomic complexity and lineage-specific adaptations potentially related to the parasitic lifestyle.

Pfcrt haplotypes and the evolutionary history of chloroquine-resistant Plasmodium falciparum

Mutations in the Pfcrt gene that change the resulting amino acids and form different haplotypes are common and correlate with the prevalence of chloroquine resistant (CQR) field isolates of the malaria parasite, Plasmodium falciparum. This correlation provides opportunities to infer the global evolutionary history of CQ resistance by analysing CQR Pfcrt haplotype data. We collated data on the Pfcrt haplotypes from different global studies and performed evolutionary genetic analysis to present comprehensive and comparative information on the global distribution of five major CQR-Pfcrt haplotypes and evolutionary inter-relationships among 38 different countries. Using the haplotype diversity data, inter-continental genetic differentiation was also ascertained.

1.° Joannis Boccacii de Certaldo liber de mulieribus claris. — 2.° Scholium de Sibyllis, imprimìs de Tiburtina et de ejus Romam adventu, ubi ea de Imperio Romano, de Christo, de Francorum Regibus vaticinari fingitur, quae postea evenerunt. — 3.° M. T. Ciceronis oratio pro lege Manilia.

Colbertinus

Constraint algorithm for k-presymplectic Hamiltonian systems. Application to singular field theories

The k-symplectic formulation of field theories is especially simple, since only tangent and cotangent bundles are needed in its description. Its defining elements show a close relationship with those in the symplectic formulation of mechanics. It will be shown that this relationship also stands in the presymplectic case. In a natural way,one can mimick the presymplectic constraint algorithm to obtain a constraint algorithmthat can be applied to k-presymplectic field theory, and more particularly to the Lagrangian and Hamiltonian formulations offield theories defined by a singular Lagrangian, as well as to the unified Lagrangian-Hamiltonian formalism (Skinner--Rusk formalism) for k-presymplectic field theory. Two examples of application of the algorithm are also analyzed.

Plastic and evolutionary responses of cell size and number to larval malnutrition in Drosophila melanogaster.

Both development and evolution under chronic malnutrition lead to reduced adult size in Drosophila. We studied the contribution of changes in size vs. number of epidermal cells to plastic and evolutionary reduction of wing size in response to poor larval food. We used flies from six populations selected for tolerance to larval malnutrition and from six unselected control populations, raised either under standard conditions or under larval malnutrition. In the control populations, phenotypic plasticity of wing size was mediated by both cell size and cell number. In contrast, evolutionary change in wing size, which was only observed as a correlated response expressed on standard food, was mediated entirely by reduction in cell number. Plasticity of cell number had been lost in the selected populations, and cell number did not differ between the sexes despite males having smaller wings. Results of this and other experimental evolution studies are consistent with the hypothesis that alleles which increase body size through prolonged growth affect wing size mostly via cell number, whereas alleles which increase size through higher growth rate do so via cell size.

Chromosomal divergence and evolutionary inferences in Rhodniini based on the chromosomal location of ribosomal genes

In this study, we used fluorescence in situ hybridisation to determine the chromosomal location of 45S rDNA clusters in 10 species of the tribe Rhodniini (Hemiptera: Reduviidae: Triatominae). The results showed striking inter and intraspecific variability, with the location of the rDNA clusters restricted to sex chromosomes with two patterns: either on one (X chromosome) or both sex chromosomes (X and Y chromosomes). This variation occurs within a genus that has an unchanging diploid chromosome number (2n = 22, including 20 autosomes and 2 sex chromosomes) and a similar chromosome size and genomic DNA content, reflecting a genome dynamic not revealed by these chromosome traits. The rDNA variation in closely related species and the intraspecific polymorphism in Rhodnius ecuadoriensis suggested that the chromosomal position of rDNA clusters might be a useful marker to identify recently diverged species or populations. We discuss the ancestral position of ribosomal genes in the tribe Rhodniini and the possible mechanisms involved in the variation of the rDNA clusters, including the loss of rDNA loci on the Y chromosome, transposition and ectopic pairing. The last two processes involve chromosomal exchanges between both sex chromosomes, in contrast to the widely accepted idea that the achiasmatic sex chromosomes of Heteroptera do not interchange sequences.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.