Nested parallelism for multi-core HPC systems using Java

Since its introduction in 1993, the Message Passing Interface (MPI) has become a de facto standard for writing High Performance Computing (HPC) applications on clusters and Massively Parallel Processors (MPPs). The recent emergence of multi-core processor systems presents a new challenge for established parallel programming paradigms, including those based on MPI. This paper presents a new Java messaging system called MPJ Express. Using this system, we exploit multiple levels of parallelism - messaging and threading - to improve application performance on multi-core processors. We refer to our approach as nested parallelism. This MPI-like Java library can support nested parallelism by using Java or Java OpenMP (JOMP) threads within an MPJ Express process. Practicality of this approach is assessed by porting to Java a massively parallel structure formation code from Cosmology called Gadget-2. We introduce nested parallelism in the Java version of the simulation code and report good speed-ups. To the best of our knowledge it is the first time this kind of hybrid parallelism is demonstrated in a high performance Java application. (C) 2009 Elsevier Inc. All rights reserved.

Four-dimensional variational data assimilation for high resolution nested models

Four-dimensional variational data assimilation (4D-Var) is used in environmental prediction to estimate the state of a system from measurements. When 4D-Var is applied in the context of high resolution nested models, problems may arise in the representation of spatial scales longer than the domain of the model. In this paper we study how well 4D-Var is able to estimate the whole range of spatial scales present in one-way nested models. Using a model of the one-dimensional advection–diffusion equation we show that small spatial scales that are observed can be captured by a 4D-Var assimilation, but that information in the larger scales may be degraded. We propose a modification to 4D-Var which allows a better representation of these larger scales.

Correction of atmospheric delay effects in radar interferometry using a nested mesoscale atmospheric model

The temporal variability of the atmosphere through which radio waves pass in the technique of differential radar interferometry can seriously limit the accuracy with which the method can measure surface motion. A forward, nested mesoscale model of the atmosphere can be used to simulate the variable water content along the radar path and the resultant phase delays. Using this approach we demonstrate how to correct an interferogram of Mount Etna in Sicily associated with an eruption in 2004-5. The regional mesoscale model (Unified Model) used to simulate the atmosphere at higher resolutions consists of four nested domains increasing in resolution (12, 4, 1, 0.3 km), sitting within the analysis version of a global numerical model that is used to initiate the simulation. Using the high resolution 3D model output we compute the surface pressure, temperature and the water vapour, liquid and solid water contents, enabling the dominant hydrostatic and wet delays to be calculated at specific times corresponding to the acquisition of the radar data. We can also simulate the second-order delay effects due to liquid water and ice.

Nonlinear identification using orthogonal forward regression with nested optimal regularization

An efficient data based-modeling algorithm for nonlinear system identification is introduced for radial basis function (RBF) neural networks with the aim of maximizing generalization capability based on the concept of leave-one-out (LOO) cross validation. Each of the RBF kernels has its own kernel width parameter and the basic idea is to optimize the multiple pairs of regularization parameters and kernel widths, each of which is associated with a kernel, one at a time within the orthogonal forward regression (OFR) procedure. Thus, each OFR step consists of one model term selection based on the LOO mean square error (LOOMSE), followed by the optimization of the associated kernel width and regularization parameter, also based on the LOOMSE. Since like our previous state-of-the-art local regularization assisted orthogonal least squares (LROLS) algorithm, the same LOOMSE is adopted for model selection, our proposed new OFR algorithm is also capable of producing a very sparse RBF model with excellent generalization performance. Unlike our previous LROLS algorithm which requires an additional iterative loop to optimize the regularization parameters as well as an additional procedure to optimize the kernel width, the proposed new OFR algorithm optimizes both the kernel widths and regularization parameters within the single OFR procedure, and consequently the required computational complexity is dramatically reduced. Nonlinear system identification examples are included to demonstrate the effectiveness of this new approach in comparison to the well-known approaches of support vector machine and least absolute shrinkage and selection operator as well as the LROLS algorithm.

Nested diets: a novel pattern of individual-level resource use

Many generalist populations may actually be composed of relatively specialist individuals. This `individual specialization` may have important ecological and evolutionary implications. Although this phenomenon has been documented in more than one hundred taxa, it is still unclear how individuals within a population actually partition resources. Here we applied several methods based on network theory to investigate the intrapopulation patterns of resource use in the gracile mouse opossum Gracilinanus microtarsus. We found evidence of significant individual specialization in this species and that the diets of specialists are nested within the diets of generalists. This novel pattern is consistent with a recently proposed model of optimal foraging and implies strong asymmetry in the interactions among individuals of a population.

The nested assembly of individual-resource networks

P>1. Much of the current understanding of ecological systems is based on theory that does not explicitly take into account individual variation within natural populations. However, individuals may show substantial variation in resource use. This variation in turn may be translated into topological properties of networks that depict interactions among individuals and the food resources they consume (individual-resource networks). 2. Different models derived from optimal diet theory (ODT) predict highly distinct patterns of trophic interactions at the individual level that should translate into distinct network topologies. As a consequence, individual-resource networks can be useful tools in revealing the incidence of different patterns of resource use by individuals and suggesting their mechanistic basis. 3. In the present study, using data from several dietary studies, we assembled individual-resource networks of 10 vertebrate species, previously reported to show interindividual diet variation, and used a network-based approach to investigate their structure. 4. We found significant nestedness, but no modularity, in all empirical networks, indicating that (i) these populations are composed of both opportunistic and selective individuals and (ii) the diets of the latter are ordered as predictable subsets of the diets of the more opportunistic individuals. 5. Nested patterns are a common feature of species networks, and our results extend its generality to trophic interactions at the individual level. This pattern is consistent with a recently proposed ODT model, in which individuals show similar rank preferences but differ in their acceptance rate for alternative resources. Our findings therefore suggest a common mechanism underlying interindividual variation in resource use in disparate taxa.

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.

Estabelecimento de um protocolo de Nested-PCR para detecção do vírus da anemia das galinhas e análise filogenética de cepas brasileiras.

O vírus da anemia das galinhas (CAV) pode causar imunodepressão em galinhas de todas as idades e doença nos frangos jovens, a qual é caracterizada por severa anemia, atrofia da medula óssea e hemorragias. A doença clínica é rara hoje, porém a forma subclínica é freqüentemente encontrada em criações comerciais e resulta num considerável decréscimo do desempenho. O CAV apresenta variabilidade genética, entretanto, em relação às amostras brasileiras do CAV, pouco ou quase nada desta variabilidade é conhecida. O presente trabalho descreve um protocolo de Nested-PCR para a detecção do CAV diretamente de amostras clínicas e analisa filogeneticamente as seqüências nucleotídicas das amostras positivas buscando verificar se a patologia apresentada por estas está relacionada com a variabilidade genética. Para a extração de DNA, o método baseado em tiocianato de guanidina mostrou-se mais eficiente e prático de executar que os demais testados. Foi selecionado um par de primers que amplifica uma região de 664 pb do gene vp1 e outro par que amplifica uma região interna de 539 pb para a realização da Nested-PCR. A especificidade dos primers foi avaliada utilizando amostras de lotes controlados para CAV e 30 diferentes isolados de vírus e bactérias causadoras de doenças em galinhas, as quais não geraram produto de amplificação. A sensibilidade foi determinada a partir de diluições seriadas da vacina comercial para o CAV. A Nested-PCR mostrou ser mais sensível do que a PCR e foi capaz de detectar 0,16 DICC50% da cepa vacinal. Além disso, a Nested-PCR detectou DNA viral em tecidos, soro e cama aviária de lotes com e sem sintomas clínicos.O produto de amplificação de 539 pb do gene vp1 de 44 amostras, provenientes de diferentes Estados produtores de frangos do Brasil, foi seqüenciado e foram encontradas 10 novas seqüências nucleotídicas do CAV. Estas 10 seqüências nucleotídicas foram analisadas filogeneticamente pelo método de distância neighbour joining com 1000 replicações o qual, mostrou que não houve correlação entre a patogenia apresentada nos animais e os grupos genéticos. Estas seqüências nucleotídicas também foram comparadas com 30 cepas de CAV isoladas em outros países e não foi observada correlação entre a distribuição geográfica e a variabilidade genética. Substituições de amino ácidos foram observadas em 9 posições sendo que, 65R substituindo o resíduo Q e 98F substituindo o resíduo Y ainda não haviam sido observadas. Conclui-se que, como técnica de detecção do CAV, o protocolo de Nested-PCR aqui descrito é mais sensível e menos trabalhoso do que o isolamento viral. As amostras Brasileiras de CAV possuem características filogenéticas similares às isoladas em outros países.

Infectious bronchitis virus: detection and vaccine Strain differentiation by semi-nested RT-PCR

A semi-nested reverse transcription-polymerase chain reaction (Semi-N-RT-PCR) was developed and used to detect the S glycoprotein gene of infectious bronchitis virus (IBV) strains and to discriminate H120 vaccine strain from other strains. Viral RNA was extracted from the allantoic fluid of chicken embryos and from tissues of chickens experimentally infected with different strains of IBV. Amplification and identification of the viral RNA was performed using two sets of primers complementary to a region of the S glycoprotein gene in the Semi-N-RT-PCR assay. The pair of primers used in the first PCR consisted of universal oligonucleotides flanking a more variable region of S1-S2 gene. The second primer pair was used in the Semi-N-RT-PCR and was comprised of one of the primers from the first universal pair together with either another universal internal oligolucleotide or a oligonucleotide sequence specific for the H120 strain of IBV. The universal primers detected all reference IBV strains and field isolates tested herein. The Semi-N-RT-PCR had high sensitivity and specificity, and was able to differentiate the H120 vaccine strain from other reference IBV strains; including M41 strain. All tissue samples collected from chickens experimentally infected with H120 or M41 strains were positive in the semi-nested RT-PCR using universal primers, while only the H120-infected tissue samples were amplified by the set of primers containing the H120-oligonucleotide. In conclusion, the ability of Semi-N-RT-PCR to detect distinct IBV strains and preliminarily discriminate the vaccine strain (H120) closes a diagnostic gap and offers the opportunity to use comprehensive PCR procedures for the IBV diagnosis.

In vitro culture, PCR , and nested PCR for the detection of Theileria equi in horses submitted to exercise

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Detection of bovine herpesvirus 5 (BoHV-5) in formalin-fixed, paraffin-embedded bovine brain by nested PCR in Colombian cattle

Fifteen cases of viral meningoencephalitis in Colombian cattle were tested by nested PCR analysis for the detection of bovine herpesvirus 5 (BoHV-5). All fatal cases had shown severe neurological signs and had occurred following natural outbreaks of the disease. The neurological infection was histologically characterized by mild to moderate inflammatory changes in the brain and cerebellum, including meningitis, mononuclear perivascular cuffing, gliosis, haemorrhage, and the presence of Gitter cells (macrophages) accompanying large areas of malacia. No intranuclear inclusion bodies were seen in any of the cases. Results from BoHV-5 molecular extraction analyses showed there were five positive cases thus confirming the presence of the virus in Colombia.

A nested-PCR assay for detection of Xylella fastidiosa in citrus plants and sharpshooter leafhoppers

Aims: Detection of Xylella fastidiosa in citrus plants and insect vectors.Methods and Results: Chelex 100 resin matrix was successfully standardized allowing a fast DNA extraction of X. fastidiosa. An amplicon of 500 bp was observed in samples of citrus leaf and citrus xylem extract, with and without symptoms of citrus variegated chlorosis, using PCR with a specific primer set indicating the presence of X. fastidiosa. The addition of insoluble acid-washed polyvinylpyrrolidone (PVPP) prior to DNA extraction of insect samples using Chelex 100 resin together with nested-PCR permitted the detection of X. fastidiosa within sharpshooter heads with great sensitivity. It was possible to detect up to two bacteria per reaction. From 250 sharpshooter samples comprising four species (Dilobopterus costalimai, Oncometopia facialis, Bucephalogonia xanthopis and Acrogonia sp.), 87 individuals showed positive results for X. fastidiosa in a nested-PCR assay.Conclusions: the use of Chelex 100 resin allowed a fast and efficient DNA extraction to be used in the detection of X. fastidiosa in citrus plants and insect vectors by PCR and nested-PCR assays, respectively.Significance and Impact of the study: the employment of efficient and sensitive methods to detect X. fastidiosa in citrus plants and insect vectors will greatly assist epidemiological studies.

Nested liana-tree network in three distinct neotropical vegetation formations

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)