Automatic decomposition of unstructured meshes employing genetic algorithms for parallel FEM computations

Parallel execution of computational mechanics codes requires efficient mesh-partitioning techniques. These mesh-partitioning techniques divide the mesh into specified number of submeshes of approximately the same size and at the same time, minimise the interface nodes of the submeshes. This paper describes a new mesh partitioning technique, employing Genetic Algorithms. The proposed algorithm operates on the deduced graph (dual or nodal graph) of the given finite element mesh rather than directly on the mesh itself. The algorithm works by first constructing a coarse graph approximation using an automatic graph coarsening method. The coarse graph is partitioned and the results are interpolated onto the original graph to initialise an optimisation of the graph partition problem. In practice, hierarchy of (usually more than two) graphs are used to obtain the final graph partition. The proposed partitioning algorithm is applied to graphs derived from unstructured finite element meshes describing practical engineering problems and also several example graphs related to finite element meshes given in the literature. The test results indicate that the proposed GA based graph partitioning algorithm generates high quality partitions and are superior to spectral and multilevel graph partitioning algorithms.

Evolution of metastable phases during crystallization of pulsed laser deposited amorphous Al–Fe films

Amorphous thin films of different Al–Fe compositions were produced by plasma/vapor quenching during pulsed laser deposition. The chosen compositions Al72Fe28, Al40Fe60, and Al18Fe82 correspond to Al5Fe2 and B2-ordered AlFe intermetallic compounds and α–Fe solid solution, respectively. The films contained fine clusters that increased with iron content. The sequences of phase evolution observed in the heating stage transmission electron microscopy studies of the pulsed laser ablation deposited films of Al72Fe28, Al40Fe60, and Al18Fe82 compositions showed evidence of composition partitioning during crystallization for films of all three compositions. This composition partitioning, in turn, resulted in the evolution of phases of compositions richer in Fe, as well as richer in Al, compared to the overall film composition in each case. The evidence of Fe-rich phases was the B2 phase in Al72Fe28 film, the L12- and DO3-ordered phases in Al40Fe60 film, and the hexagonal ε–Fe in the case of the Al18Fe82 film. On the other hand, the Al-rich phases were Al13Fe4 for both Al72Fe28 and Al40Fe60 films and DO3 and Al5Fe2 phases in the case of Al18Fe82 film. We believe that this tendency of composition partitioning during crystallization from amorphous phase is a consequence of the tendency of clustering of the Fe atoms in the amorphous phase during nucleation. The body-centered cubic phase has a nucleation advantage over other metastable phases for all three compositions. The amorphization of Al18Fe82 composition and the evolution of L12 and ε–Fe phases in the Al–Fe system were new observations of this work.

Numerical Modelling of the Observed Diurnal Cycle of Indian Monsoon Rainfall

We have analysed the diurnal cycle of rainfall over the Indian region (10S-35N, 60E-100E) using both satellite and in-situ data, and found many interesting features associated with this fundamental, yet under-explored, mode of variability. Since there is a distinct and strong diurnal mode of variability associated with the Indian summer monsoon rainfall, we evaluate the ability of the Weather Research and Forecasting Model (WRF) to simulate the observed diurnal rainfall characteristics. The model (at 54km grid-spacing) is integrated for the month of July, 2006, since this period was particularly favourable for the study of diurnal cycle. We first evaluate the sensitivity of the model to the prescribed sea surface temperature (SST), by using two different SST datasets, namely, Final Analyses (FNL) and Real-time Global (RTG). It was found that with RTG SST the rainfall simulation over central India (CI) was significantly better than that with FNL. On the other hand, over the Bay of Bengal (BoB), rainfall simulated with FNL was marginally better than with RTG. However, the overall performance of RTG SST was found to be better than FNL, and hence it was used for further model simulations. Next, we investigated the role of the convective parameterization scheme on the simulation of diurnal cycle of rainfall. We found that the Kain-Fritsch (KF) scheme performs significantly better than Betts-Miller-Janjić (BMJ) and Grell-Devenyi schemes. We also studied the impact of other physical parameterizations, namely, microphysics, boundary layer, land surface, and the radiation parameterization, on the simulation of diurnal cycle of rainfall, and identified the “best” model configuration. We used this configuration of the “best” model to perform a sensitivity study on the role of various convective components used in the KF scheme. In particular, we studied the role of convective downdrafts, convective timescale, and feedback fraction, on the simulated diurnal cycle of rainfall. The “best” model simulations, in general, show a good agreement with observations. Specifically, (i) Over CI, the simulated diurnal rainfall peak is at 1430 IST, in comparison to the observed 1430-1730 IST peak; (ii) Over Western Ghats and Burmese mountains, the model simulates a diurnal rainfall peak at 1430 IST, as opposed to the observed peak of 1430-1730 IST; (iii) Over Sumatra, both model and observations show a diurnal peak at 1730 IST; (iv) The observed southward propagating diurnal rainfall bands over BoB are weakly simulated by WRF. Besides the diurnal cycle of rainfall, the mean spatial pattern of total rainfall and its partitioning between the convective and stratiform components, are also well simulated. The “best” model configuration was used to conduct two nested simulations with one-way, three-level nesting (54-18-6km) over CI and BoB. While, the 54km and 18km simulations were conducted for the whole of July, 2006, the 6km simulation was carried out for the period 18 - 24 July, 2006. The results of our coarse- and fine-scale numerical simulations of the diurnal cycle of monsoon rainfall will be discussed.

Ligand-modulated Parallel Mechanical Unfolding Pathways of Maltose-binding Proteins

Protein folding and unfolding are complex phenomena, and it is accepted that multidomain proteins generally follow multiple pathways. Maltose-binding protein (MBP) is a large (a two-domain, 370-amino acid residue) bacterial periplasmic protein involved in maltose uptake. Despite the large size, it has been shown to exhibit an apparent two-state equilibrium unfolding in bulk experiments. Single-molecule studies can uncover rare events that are masked by averaging in bulk studies. Here, we use single-molecule force spectroscopy to study the mechanical unfolding pathways of MBP and its precursor protein (preMBP) in the presence and absence of ligands. Our results show that MBP exhibits kinetic partitioning on mechanical stretching and unfolds via two parallel pathways: one of them involves a mechanically stable intermediate (path I) whereas the other is devoid of it (path II). The apoMBP unfolds via path I in 62% of the mechanical unfolding events, and the remaining 38% follow path II. In the case of maltose-bound MBP, the protein unfolds via the intermediate in 79% of the cases, the remaining 21% via path II. Similarly, on binding to maltotriose, a ligand whose binding strength with the polyprotein is similar to that of maltose, the occurrence of the intermediate is comparable (82% via path I) with that of maltose. The precursor protein preMBP also shows a similar behavior upon mechanical unfolding. The percentages of molecules unfolding via path I are 53% in the apo form and 68% and 72% upon binding to maltose and maltotriose, respectively, for preMBP. These observations demonstrate that ligand binding can modulate the mechanical unfolding pathways of proteins by a kinetic partitioning mechanism. This could be a general mechanism in the unfolding of other large two-domain ligand-binding proteins of the bacterial periplasmic space.

Memory Architecture Exploration Framework for Cache Based Embedded SoC

Today's feature-rich multimedia products require embedded system solution with complex System-on-Chip (SoC) to meet market expectations of high performance at a low cost and lower energy consumption. The memory architecture of the embedded system strongly influences critical system design objectives like area, power and performance. Hence the embedded system designer performs a complete memory architecture exploration to custom design a memory architecture for a given set of applications. Further, the designer would be interested in multiple optimal design points to address various market segments. However, tight time-to-market constraints enforces short design cycle time. In this paper we address the multi-level multi-objective memory architecture exploration problem through a combination of exhaustive-search based memory exploration at the outer level and a two step based integrated data layout for SPRAM-Cache based architectures at the inner level. We present a two step integrated approach for data layout for SPRAM-Cache based hybrid architectures with the first step as data-partitioning that partitions data between SPRAM and Cache, and the second step is the cache conscious data layout. We formulate the cache-conscious data layout as a graph partitioning problem and show that our approach gives up to 34% improvement over an existing approach and also optimizes the off-chip memory address space. We experimented our approach with 3 embedded multimedia applications and our approach explores several hundred memory configurations for each application, yielding several optimal design points in a few hours of computation on a standard desktop.

Nature's Swiss Army Knives: Ovipositor Structure Mirrors Ecology in a Multitrophic Fig Wasp Community

Background: Resource partitioning is facilitated by adaptations along niche dimensions that range from morphology to behaviour. The exploitation of hidden resources may require specially adapted morphological or sensory tools for resource location and utilisation. Differences in tool diversity and complexity can determine not only how many species can utilize these hidden resources but also how they do so. Methodology and Principal Findings: The sclerotisation, gross morphology and ultrastructure of the ovipositors of a seven-member community of parasitic wasps comprising of gallers and parasitoids developing within the globular syconia (closed inflorescences) of Ficus racemosa (Moraceae) was investigated. These wasps also differ in their parasitism mode (external versus internal oviposition) and their timing of oviposition into the expanding syconium during its development. The number and diversity of sensilla, as well as ovipositor teeth, increased from internally ovipositing to externally ovipositing species and from gallers to parasitoids. The extent of sclerotisation of the ovipositor tip matched the force required to penetrate the syconium at the time of oviposition of each species. The internally ovipositing pollinator had only one type of sensillum and a single notch on the ovipositor tip. Externally ovipositing species had multiple sensilla types and teeth on their ovipositors. Chemosensilla were most concentrated at ovipositor tips while mechanoreceptors were more widely distributed, facilitating the precise location of hidden hosts in these wasps which lack larval host-seeking behaviour. Ovipositor traits of one parasitoid differed from those of its syntopic galler congeners and clustered with those of parasitoids within a different wasp subfamily. Thus ovipositor tools can show lability based on adaptive necessity, and are not constrained by phylogeny. Conclusions/Significance: Ovipositor structure mirrored the increasingly complex trophic ecology and requirements for host accessibility in this parasite community. Ovipositor structure could be a useful surrogate for predicting the biology of parasites in other communities.

Rough Core Vector Clustering

Support Vector Clustering has gained reasonable attention from the researchers in exploratory data analysis due to firm theoretical foundation in statistical learning theory. Hard Partitioning of the data set achieved by support vector clustering may not be acceptable in real world scenarios. Rough Support Vector Clustering is an extension of Support Vector Clustering to attain a soft partitioning of the data set. But the Quadratic Programming Problem involved in Rough Support Vector Clustering makes it computationally expensive to handle large datasets. In this paper, we propose Rough Core Vector Clustering algorithm which is a computationally efficient realization of Rough Support Vector Clustering. Here Rough Support Vector Clustering problem is formulated using an approximate Minimum Enclosing Ball problem and is solved using an approximate Minimum Enclosing Ball finding algorithm. Experiments done with several Large Multi class datasets such as Forest cover type, and other Multi class datasets taken from LIBSVM page shows that the proposed strategy is efficient, finds meaningful soft cluster abstractions which provide a superior generalization performance than the SVM classifier.

Climate forcing and response to idealized changes in surface latent and sensible heat

Land use and land cover changes affect the partitioning of latent and sensible heat, which impacts the broader climate system. Increased latent heat flux to the atmosphere has a local cooling influence known as `evaporative cooling', but this energy will be released back to the atmosphere wherever the water condenses. However, the extent to which local evaporative cooling provides a global cooling influence has not been well characterized. Here, we perform a highly idealized set of climate model simulations aimed at understanding the effects that changes in the balance between surface sensible and latent heating have on the global climate system. We find that globally adding a uniform 1 W m(-2) source of latent heat flux along with a uniform 1 W m(-2) sink of sensible heat leads to a decrease in global mean surface air temperature of 0.54 +/- 0.04 K. This occurs largely as a consequence of planetary albedo increases associated with an increase in low elevation cloudiness caused by increased evaporation. Thus, our model results indicate that, on average, when latent heating replaces sensible heating, global, and not merely local, surface temperatures decrease.

Tie lines and activities in the system CoO MgO SiO2 at 1373 K

The tie lines between (CoXMg1−X)O solid solution with rock salt structure and orthosilicate solid solution (CoYMg1−Y)-Si0.5O2, and between orthosilicate and metasilicate (CoZMg1-Z)SiO3 crystalline solutions, have been determined experimentally at 1373 K. The compositions of coexisting phases have been determined by electron probe microanalysis (EPMA) and lattice parameter measurement on equilibrated samples. The metasilicate solid solution exists only for 0 > Z > 0.213. The activity of CoO in the rock salt solid solution was determined as a function of composition and temperature in the range of 1023 to 1373 K using a solid-state galvanic cell: Pt, (CoXMg1−X)O+Co|(Y2O3)ZrO2|Co+CoO, Pt The free energy of mixing of (CoXMg1−X)O crystalline solution can be expressed by the equation ΔGE=X(1 −X)[(6048 − 2.146T)X+ (8745 − 3.09T)(1 −X)] J·mol−1 The thermodynamic data for the rock salt phase is combined with information on interphase partitioning of Co and Mg to generate the mixing properties for the ortho- and metasilicate solid solutions. For the orthosilicate solution (CoYMg1 −Y)Si0.5O2 at 1373 K, the excess Gibbs free energy of mixing is given by the relation ΔGE=Y(1 −Y)[2805Y+ 3261(1 −Y)] J·mol−1 For the metasilicate solution (CoZMg1 −Z)SiO3 at the same temperature, the excess free energy can be expressed by the relation ΔGE=Z(1 −Z)[2570Z+ 3627(1 −Z)] J·mol−1

On the cut‐on frequencies of a large round duct with azimuthal as well as radial partitions

This letter deals with a three‐dimensional analysis of circular sectors and annular segments resulting from the partitioning of a round (cylindrical) duct for use in an active noise control system. The relevant frequency equations are derived for stationary medium and solved numerically to arrive at the cut‐on frequencies of the first few modes. The resultant table indicates among other things that azimuthal partitioning does not raise the cutoff frequency (the smallest cut‐on frequency) beyond a particular value, and that radial partitioning is counterproductive in that respect.

Improving Superscalar Instruction Dispatch And Issue By Exploiting Dynamic Code Sequences

Superscalar processors currently have the potential to fetch multiple basic blocks per cycle by employing one of several recently proposed instruction fetch mechanisms. However, this increased fetch bandwidth cannot be exploited unless pipeline stages further downstream correspondingly improve. In particular,register renaming a large number of instructions per cycle is diDcult. A large instruction window, needed to receive multiple basic blocks per cycle, will slow down dependence resolution and instruction issue. This paper addresses these and related issues by proposing (i) partitioning of the instruction window into multiple blocks, each holding a dynamic code sequence; (ii) logical partitioning of the registerjle into a global file and several local jles, the latter holding registers local to a dynamic code sequence; (iii) the dynamic recording and reuse of register renaming information for registers local to a dynamic code sequence. Performance studies show these mechanisms improve performance over traditional superscalar processors by factors ranging from 1.5 to a little over 3 for the SPEC Integer programs. Next, it is observed that several of the loops in the benchmarks display vector-like behavior during execution, even if the static loop bodies are likely complex for compile-time vectorization. A dynamic loop vectorization mechanism that builds on top of the above mechanisms is briefly outlined. The mechanism vectorizes up to 60% of the dynamic instructions for some programs, albeit the average number of iterations per loop is quite small.

Location Area Planning Using Simulated Annealing With A New Solution Representation

Location area planning problem is to partition the cellular/mobile network into location areas with the objective of minimizing the total cost. This partitioning problem is a difficult combinatorial optimization problem. In this paper, we use the simulated annealing with a new solution representation. In our method, we can automatically generate different number of location areas using Compact Index (CI) to obtain the optimal/best partitions. We compare the results obtained in our method with the earlier results available in literature. We show that our methodology is able to perform better than earlier methods.

On the Presumed Kinetic Consequences of Pre-equilibrium. Implications for the Michaelis-Menten Equation

The overall rate equation for a reaction sequence consisting of a pre-equilibrium and rate-determining steps should not be derived on the basis of the concentration of the intermediate product (X). This is apparently indicated by transition state theory (as the path followed to reach the highest energy transition state is irrelevant), but also proved by a straight-forward mathematical approach. The thesis is further supported by the equations of concurrent reactions as applied to the partitioning of X between the two competing routes (reversal of the pre-equilibrium and formation of product). The rate equation may only be derived rigorously on the basis of the law of mass action. It is proposed that the reactants acquire the overall activation energy prior to the pre-equilibrium, thus forming X in a high-energy state en route to the rate-determining transition state. (It is argued that conventional energy profile diagrams are misleading and need to be reinterpreted.) Also, these arguments invalidate the Michaelis-Menten equation of enzyme kinetics, and necessitate a fundamental revision of our present understanding of enzyme catalysis. (The observed ``saturation kinetics'' possibly arises from weak binding of a second molecule of substrate at the active site; analogous conclusions apply to reactions at surfaces).

Tight co-degree condition for perfect matchings in 4-graphs

We will give a tight minimum co-degree condition for a 4-uniform hypergraph to contain a perfect matching.

Cross-Guided Clustering: Transfer of Relevant Supervision across Tasks

Lack of supervision in clustering algorithms often leads to clusters that are not useful or interesting to human reviewers. We investigate if supervision can be automatically transferred for clustering a target task, by providing a relevant supervised partitioning of a dataset from a different source task. The target clustering is made more meaningful for the human user by trading-off intrinsic clustering goodness on the target task for alignment with relevant supervised partitions in the source task, wherever possible. We propose a cross-guided clustering algorithm that builds on traditional k-means by aligning the target clusters with source partitions. The alignment process makes use of a cross-task similarity measure that discovers hidden relationships across tasks. When the source and target tasks correspond to different domains with potentially different vocabularies, we propose a projection approach using pivot vocabularies for the cross-domain similarity measure. Using multiple real-world and synthetic datasets, we show that our approach improves clustering accuracy significantly over traditional k-means and state-of-the-art semi-supervised clustering baselines, over a wide range of data characteristics and parameter settings.