A new parallel algorithm for parsing arithmetic infix expressions

A new parallel algorithm for transforming an arithmetic infix expression into a par se tree is presented. The technique is based on a result due to Fischer (1980) which enables the construction of the parse tree, by appropriately scanning the vector of precedence values associated with the elements of the expression. The algorithm presented here is suitable for execution on a shared memory model of an SIMD machine with no read/write conflicts permitted. It uses O(n) processors and has a time complexity of O(log2n) where n is the expression length. Parallel algorithms for generating code for an SIMD machine are also presented.

DECT evaluation of noncalc fied coronary artery plaque

Purpose: Composition of the coronary artery plaque is known to have critical role in heart attack. While calcified plaque can easily be diagnosed by conventional CT, it fails to distinguish between fibrous and lipid rich plaques. In the present paper, the authors discuss the experimental techniques and obtain a numerical algorithm by which the electron density (rho(e)) and the effective atomic number (Z(eff)) can be obtained from the dual energy computed tomography (DECT) data. The idea is to use this inversion method to characterize and distinguish between the lipid and fibrous coronary artery plaques. Methods: For the purpose of calibration of the CT machine, the authors prepare aqueous samples whose calculated values of (rho(e), Z(eff)) lie in the range of (2.65 x 10(23) <= rho(e) <= 3.64 x 10(23)/cm(3)) and (6.80 <= Z(eff) <= 8.90). The authors fill the phantom with these known samples and experimentally determine HU(V-1) and HU(V-2), with V-1,V-2 = 100 and 140 kVp, for the same pixels and thus determine the coefficients of inversion that allow us to determine (rho(e), Z(eff)) from the DECT data. The HU(100) and HU(140) for the coronary artery plaque are obtained by filling the channel of the coronary artery with a viscous solution of methyl cellulose in water, containing 2% contrast. These (rho(e), Z(eff)) values of the coronary artery plaque are used for their characterization on the basis of theoretical models of atomic compositions of the plaque materials. These results are compared with histopathological report. Results: The authors find that the calibration gives Pc with an accuracy of 3.5% while Z(eff) is found within 1% of the actual value, the confidence being 95%. The HU(100) and HU(140) are found to be considerably different for the same plaque at the same position and there is a linear trend between these two HU values. It is noted that pure lipid type plaques are practically nonexistent, and microcalcification, as observed in histopathology, has to be taken into account to explain the nature of the observed (rho(e), Z(eff)) data. This also enables us to judge the composition of the plaque in terms of basic model which considers the plaque to be composed of fibres, lipids, and microcalcification. Conclusions: This simple and reliable method has the potential as an effective modality to investigate the composition of noncalcified coronary artery plaques and thus help in their characterization. In this inversion method, (rho(e), Z(eff)) of the scanned sample can be found by eliminating the effects of the CT machine and also by ensuring that the determination of the two unknowns (rho(e), Z(eff)) does not interfere with each other and the nature of the plaque can be identified in terms of a three component model. (C) 2015 American Association of Physicists in Medicine.

Using tau polarization to probe the stau co-annihilation region of the minimal supergravity model at the LHC

The minimal supergravity model predicts the polarization of the tau coming from the stau to bino decay in the co-annihilation region to +1. This can be exploited to extract this soft tau signal at LHC and also to measure the tiny mass differences between the stau and the bi lightest superparticle. Moreover, this strategy will be applicable for a wider class of bino lightest superparticle models, where the lighter stau has a right component at least of similar size as the left.

Measurement of aerodynamic forces for missile shaped body in hypersonic shock tunnel using 6-component accelerometer based balance system

This paper reports the basic design of a new six component force balance system using miniature piezoelectric accelerometers to measure all aerodynamic forces and moments for a test model in hypersonic shock tunnel (HST2). Since the flow duration in a hypersonic shock tunnel is of the order of $1$ ms, the balance system [1] uses fast response accelerometers (PCB Piezotronics; frequency range of 1-10 kHz) for obtaining the aerodynamic data. The alance system has been used to measure the basic aerodynamic forces and moments on a missile shaped body at Mach $8$ in the IISc hypersonic shock tunnel. The experimentally measured values match well with theoretical predictions.

Grids with multiple batch systems for performance enhancement of multi-component and parameter sweep parallel applications

In this work, we evaluate the benefits of using Grids with multiple batch systems to improve the performance of multi-component and parameter sweep parallel applications by reduction in queue waiting times. Using different job traces of different loads, job distributions and queue waiting times corresponding to three different queuing policies(FCFS, conservative and EASY backfilling), we conducted a large number of experiments using simulators of two important classes of applications. The first simulator models Community Climate System Model (CCSM), a prominent multi-component application and the second simulator models parameter sweep applications. We compare the performance of the applications when executed on multiple batch systems and on a single batch system for different system and application configurations. We show that there are a large number of configurations for which application execution using multiple batch systems can give improved performance over execution on a single system.

Two-Phase Thermodynamic Model for Efficient and Accurate Absolute Entropy of Water from Molecular Dynamics Simulations

Presented here is the two-phase thermodynamic (2PT) model for the calculation of energy and entropy of molecular fluids from the trajectory of molecular dynamics (MD) simulations. In this method, the density of state (DoS) functions (including the normal modes of translation, rotation, and intramolecular vibration motions) are determined from the Fourier transform of the corresponding velocity autocorrelation functions. A fluidicity parameter (f), extracted from the thermodynamic state of the system derived from the same MD, is used to partition the translation and rotation modes into a diffusive, gas-like component (with 3Nf degrees of freedom) and a nondiffusive, solid-like component. The thermodynamic properties, including the absolute value of entropy, are then obtained by applying quantum statistics to the solid component and applying hard sphere/rigid rotor thermodynamics to the gas component. The 2PT method produces exact thermodynamic properties of the system in two limiting states: the nondiffusive solid state (where the fluidicity is zero) and the ideal gas state (where the fluidicity becomes unity). We examine the 2PT entropy for various water models (F3C, SPC, SPC/E, TIP3P, and TIP4P-Ew) at ambient conditions and find good agreement with literature results obtained based on other simulation techniques. We also validate the entropy of water in the liquid and vapor phases along the vapor-liquid equilibrium curve from the triple point to the critical point. We show that this method produces converged liquid phase entropy in tens of picoseconds, making it an efficient means for extracting thermodynamic properties from MD simulations.

Induction reheating of A356.2 aluminum alloy and thixocasting as automobile component

The development work for producing an automobile component by thixocasting using A356.2 alloy was introduced. As the first step, the alloy was electromagnetically stirred and solidified to produce a billet with non-dendritic microstructure. The microstructure depended on several process parameters such as stirring intensity, stirring frequency, cooling rate, and melt initial superheat. Through a series of computational studies and controlled experiments, a set of process parameters were identified to produce the best microstructures. Reheating of a billet with non-dendritic microstructure to a semisolid temperature was the next step for thixo-casting of the components. The reheating process was characterized for various reheating cycles using a vertical-type reheating machine. The induction heating cycle was optimized to obtain a near-uniform temperature distribution in radial as well as axial direction of the billet, and the heating was continued until the liquid fraction reached about 50%. These parameters were determined with the help of a computational fluid dynamics (CFD) model of die filling and solidification of the semisolid alloy. The heated billets were subsequently thixo-cast into automobile components using a real-time controlled die casting machine. The results show that the castings are near net shape, free from porosity, good surface finish and have superior mechanical properties compared to those produced by conventional die casting processes using the same alloy.

A rationalisation of fatigue thresholds in pearlitic steels using a theoretical model

From a detailed re-examination of results in the literature, the effects of microstructure sizes, namely interlamellar spacing, pearlitic colony size and the prior austentitic grain size on the thresholds for fatigue crack growth (ΔKth) and crack closure (Kcl, th) have been illustrated. It is shown that while interlamellar spacing explicitly controls yield strength, a similar effect on ΔKth cannot be expected. On the other hand, the pearlitic colony size is shown to strongly influence ΔKth and Kcl, th through the deflection and retardation of cracks at colony boundaries. Consequently, an increase in ΔKth and Kcl, th with colony size has been found. The development of a theoretical model to illustrate the effects of colony size, shear flow stress in the slip band and macroscopic yield strength on Kcl, th and ΔKth is presented. the model assumes colony boundaries as potential sites for slip band pile-up formation and subsequent crack deflection finally leading to zig-zag crack growth. Using the concepts of roughness induced crack closure, the magnitude of Kcl, th is quantified as a function of colony size. In deriving the model, the flow stress in the slip band has been considered to represent the work hardened state in pearlite. Comparison of the theoretically predicted trend with the experimental data demonstrates very good agreement. Further, the intrinsic or closure free component of the fatigue threshold, ΔKeff, th is found to be insensitive to colony size and interlamellar spacing. Using a criterion for intrinsic fatigue threshold which considers the attainment of a critical fracture stress over a characteristic distance corresponding to interlamellar spacing, ΔKth values at high R values can be estimated with reasonable accuracy. The magnitude of ΔKth as a function of colony size is then obtained by summing up the average value of experimentally obtained ΔKeff, th values and the predicted Kcl, th values as a function of colony size. Again, very good agreement of the theoretically predicted ΔKth values with those experimentally obtained has been demonstrated.

A Bilinear Constitutive Model for Isotropic Bimodulus Materials

Proper formulation of stress-strain relations, particularly in tension-compression situations for isotropic biomodulus materials, is an unresolved problem. Ambartsumyan's model [8] and Jones' weighted compliance matrix model [9] do not satisfy the principle of coordinate invariance. Shapiro's first stress invariant model [10] is too simple a model to describe the behavior of real materials. In fact, Rigbi [13] has raised a question about the compatibility of bimodularity with isotropy in a solid. Medri [2] has opined that linear principal strain-principal stress relations are fictitious, and warned that the bilinear approximation of uniaxial stress-strain behavior leads to ill-working bimodulus material model under combined loading. In the present work, a general bilinear constitutive model has been presented and described in biaxial principal stress plane with zonewise linear principal strain-principal stress relations. Elastic coefficients in the model are characterized based on the signs of (i) principal stresses, (ii) principal strains, and (iii) on the value of strain energy component ratio ER greater than or less than unity. The last criterion is used in tension-compression and compression-tension situations to account for different shear moduli in pure shear stress and pure shear strain states as well as unequal cross compliances.

Two-Fermi liquid model for high Tc superconductivity involving suppression of tunnelling by decoherence

We consider a model system of two interacting Fermi-liquids, one of which is light and the other much heavier. In the normal state the lighter component provides a quantum mechanical bath coupled 'ohmically' to the heavier component in the sense of Caldeira and Leggett, suppressing thereby the band (tunnelling) matrix elements of the heavier component. Thus we lose the energy of delocalization. On the other hand, a superconducting ordering stiffens the bath spectral function at low energies and so restores the tunnelling. The resulting regain of the delocalization energy bootstraps so as to stabilize the superconducting order that caused it. It is conceivable that the motions parallel to the easy ab-plane and along the hard c-axis may also effectively correspond to the light and the heavy Fermi-liquids, respectively.

A model for static foam drainage

A model for static foam drainage, based on the pentagonal dodecahedral shape of bubbles, that takes into account the surface mobility of both films and Plateau border walls has been developed. The model divides the Plateau borders into nearly horizontal and nearly vertical categories and assigns different roles to them. The films are assumed to drain into all the adjacent Plateau borders equally. The horizontal Plateau borders are assumed to receive liquid from films and drain into vertical Plateau borders, which in turn form the main component for gravity drainage. The model yields the liquid holdup values for films, horizontal Plateau borders and vertical Plateau borders as functions of height and time. The model has been tested on static foams whose cumulative drainage was measured as a function of time. The experimental data on the effect of foam height, initial holdup, surface viscosity, etc. can be explained by the model quantitatively.

ECG component delineation by Prony's method

A simple, non-iterative method for component wave delineation from the electrocardiogram (ECG) is derived by modelling its discrete cosine transform (DCT) as a sum of damped cosinusoids. Amplitude, phase, damping factor and frequency parameters of each of the cosinusoids are estimated by the extended Prony method. Different component waves are represented by non-overlapping clusters of model poles in the z plane and thus a component wave is derived by the addition of the inverse transformed (IDCT) impulse responses of the poles in the cluster. Akaike's information criterion (AIC) is used to determine the model order. The method performed satisfactory even in the presence of artifacts. The efficacy of the method is illustrated by analysis of continuous strips of ECG data.

Prediction of the Indian summer monsoon rainfall using a state-of-the-art coupled ocean-atmosphere model

A state-of-the-art model of the coupled ocean-atmosphere system, the climate forecast system (CFS), from the National Centres for Environmental Prediction (NCEP), USA, has been ported onto the PARAM Padma parallel computing system at the Centre for Development of Advanced Computing (CDAC), Bangalore and retrospective predictions for the summer monsoon (June-September) season of 2009 have been generated, using five initial conditions for the atmosphere and one initial condition for the ocean for May 2009. Whereas a large deficit in the Indian summer monsoon rainfall (ISMR; June-September) was experienced over the Indian region (with the all-India rainfall deficit by 22% of the average), the ensemble average prediction was for above-average rainfall during the summer monsoon. The retrospective predictions of ISMR with CFS from NCEP for 1981-2008 have been analysed. The retrospective predictions from NCEP for the summer monsoon of 1994 and that from CDAC for 2009 have been compared with the simulations for each of the seasons with the stand-alone atmospheric component of the model, the global forecast system (GFS), and observations. It has been shown that the simulation with GFS for 2009 showed deficit rainfall as observed. The large error in the prediction for the monsoon of 2009 can be attributed to a positive Indian Ocean Dipole event seen in the prediction from July onwards, which was not present in the observations. This suggests that the error could be reduced with improvement of the ocean model over the equatorial Indian Ocean.

Non-Fermi liquid behavior in a mean-field study of the t-J model: Role of zero-point spin fluctuations

We present analytic results to show that the Schwinger-boson hole-fermion mean-field state exhibits non-Fermi liquid behavior due to spin-charge separation. The physical electron Green's function consists of three additive components. (a) A Fermi-liquid component associated with the bose condensate. (b) A non-Fermi liquid component which has a logarithmic peak and a long tail that gives rise to a linear density of states that is symmetric about the Fermi level and a momentum distribution function with a logarithmic discontinuity at the Fermi surface. (c) A second non-Fermi liquid component associated with the thermal bosons which leads to a constant density of states. It is shown that zero-point fluctuations associated with the spin-degrees of freedom are responsible for the logarithmic instabilities and the restoration of particle-hole symmetry close to the Fermi surface.

WOMM: A Weak Operational Memory Model

Memory models of shared memory concurrent programs define the values a read of a shared memory location is allowed to see. Such memory models are typically weaker than the intuitive sequential consistency semantics to allow efficient execution. In this paper, we present WOMM (abbreviation for Weak Operational Memory Model) that formally unifies two sources of weak behavior in hardware memory models: reordering of instructions and weakly consistent memory. We show that a large number of optimizations are allowed by WOMM. We also show that WOMM is weaker than a number of hardware memory models. Consequently, if a program behaves correctly under WOMM, it will be correct with respect to those hardware memory models. Hence, WOMM can be used as a formally specified abstraction of the hardware memory models. Moreover; unlike most weak memory models, WOMM is described using operational semantics, making it easy to integrate into a model checker for concurrent programs. We further show that WOMM has an important property - it has sequential consistency semantics for datarace-free programs.