Runtime dependence computation and execution of loops on heterogeneous systems

GPUs have been used for parallel execution of DOALL loops. However, loops with indirect array references can potentially cause cross iteration dependences which are hard to detect using existing compilation techniques. Applications with such loops cannot easily use the GPU and hence do not benefit from the tremendous compute capabilities of GPUs. In this paper, we present an algorithm to compute at runtime the cross iteration dependences in such loops. The algorithm uses both the CPU and the GPU to compute the dependences. Specifically, it effectively uses the compute capabilities of the GPU to quickly collect the memory accesses performed by the iterations by executing the slice functions generated for the indirect array accesses. Using the dependence information, the loop iterations are levelized such that each level contains independent iterations which can be executed in parallel. Another interesting aspect of the proposed solution is that it pipelines the dependence computation of the future level with the actual computation of the current level to effectively utilize the resources available in the GPU. We use NVIDIA Tesla C2070 to evaluate our implementation using benchmarks from Polybench suite and some synthetic benchmarks. Our experiments show that the proposed technique can achieve an average speedup of 6.4x on loops with a reasonable number of cross iteration dependences.

Nonlinear dielectric response of BaBi4Ti4O15 ceramics

Nonlinear dielectric response of BaBi4Ti4O15 ceramics synthesized via the conventional solid-state reaction route has been monitored over a wide range of electric field strengths (E-0 = 0.5 - 5 kV/cm). Dielectric permittivity was found to increase linearly within the range of applied field. Rayleigh relations were employed to interpret the nonlinear dielectric response and the contribution of irreversible domain wall motion to the macroscopic permittivity was separated. The values of room temperature Rayleigh dielectric coefficient (alpha) and relative initial permittivity (epsilon'(init)) were found to be 2.28 +/- 0.02 cm/kV and 146.10 +/- 0.07, respectively. A reasonable agreement between the simulated and measured polarization-electric field (P-E) hysteresis loops was observed at an applied electric field of 5 kV/cm.

Constraining large earthquakes along the Andaman trench using deepwater turbidites: prospects and challenges

The 2004 Sumatra-Andaman earthquake was unprecedented in terms of its magnitude (M-w 9.2), rupture length along the plate boundary (1300 km) and size of the resultant tsunami. Since 2004, efforts are being made to improve the understanding of the seismic hazard in the Sumatra-Andaman subduction zone in terms of recurrence patterns of major earthquakes and tsunamis. It is reasonable to assume that previous earthquake events in the Myanmar Andaman segment must be preserved in the geological record in the form of seismo-turbidite sequences. Here we present the prospects of conducting deep ocean palaeoseismicity investigations in order to refine the quantification of the recurrence pattern of large subduction-zone earthquakes along the Andaman-Myanmar arc. Our participation in the Sagar Kanya cruise SK-273 (in June 2010) was to test the efficacy of such a survey. The primary mission of the cruise, along a short length (300 km) of the Sumatra Andaman subduction front was to collect bathymetric data of the ocean floor trenchward of the Andaman Islands. The agenda of our piggyback survey was to fix potential coring sites that might preserve seismo-turbidite deposits. In this article we present the possibilities and challenges of such an exercise and our first-hand experience of such a preliminary survey. This account will help future researchers with similar scientific objectives who would want to survey the deep ocean archives of this region for evidence of extreme events like major earthquakes.

Odd-even effect in the elastic modulii of alpha,omega-Alkanedicarboxylic acids

Nanoindentation studies on alpha,omega-alkanedicarboxylic acids reveal that the elastic modulus, E, shows an odd-even alternation in exactly the same manner as the melting temperature, T-m. The results are consistent with the hypothesis that the strained molecular conformations in the odd diacids are the reasons for these alternations in T-m. The same packing features that lower T-m in the odd acids lead to easy accommodation of the deformation during nanoindentation and hence their low E.

Warped alternatives to Froggatt-Nielsen models

We consider the Randall-Sundrum (RS) setup to be a theory of flavor, as an alternative to Froggatt-Nielsen models instead of as a solution to the hierarchy problem. The RS framework is modified by taking the low-energy brane to be at the grand unified theory (GUT) scale. This also alleviates constraints from flavor physics. Fermion masses and mixing angles are fit at the GUT scale. The ranges of the bulk mass parameters are determined using a chi(2) fit taking into consideration the variation in O(1) parameters. In the hadronic sector, the heavy top quark requires large bulk mass parameters localizing the right-handed top quark close to the IR brane. Two cases of neutrino masses are considered: (a) Planck scale lepton number violation and (b) Dirac neutrino masses. Contrary to the case of weak scale RS models, both these cases give reasonable fits to the data, with the Planck scale lepton number violation fitting slightly better compared to the Dirac case. In the supersymmetric version, the fits are not significantly different except for the variation in tan beta. If the Higgs superfields and the supersymmetry breaking spurion are localized on the same brane, then the structure of the sfermion masses are determined by the profiles of the zero modes of the hypermultiplets in the bulk. Trilinear terms have the same structure as the Yukawa matrices. The resultant squark spectrum is around similar to 2-3 TeV required by the light Higgs mass to be around 125 GeV and to satisfy the flavor violating constraints.

Numerical simulations of CO2 migration during charnockite genesis

Charnockite is considered to be generated either through the dehydration of granitic magma by CO2 purging or by solid-state dehydration through CO2 metasomatism during granulite facies metamorphism. To understand the extent of dehydration, CO2 migration is quantitatively modeled in silicate melt and metasomatic fluid as a function of temperature, H2O wt%, pressure, basal CO2 flux and dynamic viscosity. Numerical simulations show that CO2 advection through porous and permeable high-grade metamorphic rocks can generate dehydrated patches close to the CO2 flow path, as illustrated by the occurrences of ``incipient charnockites.'' CO2 reaction-front velocity constrained by field observations is 0.69 km/m.y., a reasonable value, which matches well with other studies. On the other hand, temperature, rate of cooling, and basal CO2 flux are the critical parameters affecting CO2 diffusion through a silicate melt. CO2 diffusion through silicate melt can only occur at temperature greater than 840 degrees C and during slow cooling (<= 3.7 x 10(-5) degrees C/yr), features that are typical of magma emplacement in the lower crust. Stalling of CO2 fluxing at similar to 840 degrees C explains why some deep-level plutons contain both hydrous and anhydrous (charnockitic) mineral assemblages. CO2 diffusion through silicate melt is virtually insensitive to pressure. Addition of CO2 basal flux facilitates episodic dehydrated melt migration by generating fracture pathways.

Scalable flow-sensitive pointer analysis for fava with strong updates

The ability to perform strong updates is the main contributor to the precision of flow-sensitive pointer analysis algorithms. Traditional flow-sensitive pointer analyses cannot strongly update pointers residing in the heap. This is a severe restriction for Java programs. In this paper, we propose a new flow-sensitive pointer analysis algorithm for Java that can perform strong updates on heap-based pointers effectively. Instead of points-to graphs, we represent our points-to information as maps from access paths to sets of abstract objects. We have implemented our analysis and run it on several large Java benchmarks. The results show considerable improvement in precision over the points-to graph based flow-insensitive and flow-sensitive analyses, with reasonable running time.

Predictability of nonstationary time series using wavelet and EMD based ARMA models

Research has been undertaken to ascertain the predictability of non-stationary time series using wavelet and Empirical Mode Decomposition (EMD) based time series models. Methods have been developed in the past to decompose a time series into components. Forecasting of these components combined with random component could yield predictions. Using this ideology, wavelet and EMD analyses have been incorporated separately which decomposes a time series into independent orthogonal components with both time and frequency localizations. The component series are fit with specific auto-regressive models to obtain forecasts which are later combined to obtain the actual predictions. Four non-stationary streamflow sites (USGS data resources) of monthly total volumes and two non-stationary gridded rainfall sites (IMD) of monthly total rainfall are considered for the study. The predictability is checked for six and twelve months ahead forecasts across both the methodologies. Based on performance measures, it is observed that wavelet based method has better prediction capabilities over EMD based method despite some of the limitations of time series methods and the manner in which decomposition takes place. Finally, the study concludes that the wavelet based time series algorithm can be used to model events such as droughts with reasonable accuracy. Also, some modifications that can be made in the model have been discussed that could extend the scope of applicability to other areas in the field of hydrology. (C) 2013 Elesvier B.V. All rights reserved.

Effect of SST Variation on ITCZ in APE Simulations

The effect of meridional variation of sea surface temperature (SST) on tropical atmospheric circulation is analyzed using Aqua-planet Experiment (APE) simulations. The meridional SST gradient around the narrow SST peak in CONTROL simulation favours a strong and single equatorial Intertropical Convergence Zone (ITCZ, defined by the maximum of zonally averaged total precipitation) in all APE models. In contrast, flat equatorial SST peak (FLAT simulation) favours split/double ITCZs flanking the SST maximum, in the majority of the APE models. Although there is reasonable agreement for SST sensitivity of ITCZ among the APE models in CONTROL, there exists disparity among them in FLAT case. Similarly, while the total and convective precipitation responses are consistent among the models, the large-scale precipitation response shows considerable inter-model variations in FLAT case. The APE intercomparison indicates that the occurrence and positioning of the ITCZ are primarily related to boundary layer moisture convergence as a response to the meridional variation of SST. Furthermore, the meridional gradient of tropospheric temperature is found to be an important factor that can influence the positioning of ITCZ. FLAT SST distribution is found to be similar to the observed distribution over the Indian region during summer season. Models that yield double ITCZs in this case simulate an easterly jet over the equatorial region (similar to 15 degrees equatorward of the ITCZ). This is analogous to the Tropical Easterly Jet (TEJ), which is a unique feature observed over the Indian region during summer monsoon season, with its core at 12 degrees N, equatorward of the seasonal convergence zone centered along 25 degrees N. In these models, positive meridional temperature gradient and the associated easterly shear in the atmosphere strengthened by moisture convergence penetrate up to the upper troposphere, with which TEJ is in thermal wind balance.

Joint Scheduling & Jamming for Data Secrecy in Wireless Networks

The broadcast nature of the wireless medium jeopardizes secure transmissions. Cryptographic measures fail to ensure security when eavesdroppers have superior computational capability; however, it can be assured from information theoretic security approaches. We use physical layer security to guarantee non-zero secrecy rate in single source, single destination multi-hop networks with eavesdroppers for two cases: when eavesdropper locations and channel gains are known and when their positions are unknown. We propose a two-phase solution which consists of finding activation sets and then obtaining transmit powers subject to SINR constraints for the case when eavesdropper locations are known. We introduce methods to find activation sets and compare their performance. Necessary but reasonable approximations are made in power minimization formulations for tractability reasons. For scenarios with no eavesdropper location information, we suggest vulnerability region (the area having zero secrecy rate) minimization over the network. Our results show that in the absence of location information average number of eavesdroppers who have access to data is reduced.

Physics and chemistry of CdTe/CdS thin film heterojunction photovoltaic devices: fundamental and critical aspects

Among the armoury of photovoltaic materials, thin film heterojunction photovoltaics continue to be a promising candidate for solar energy conversion delivering a vast scope in terms of device design and fabrication. Their production does not require expensive semiconductor substrates and high temperature device processing, which allows reduced cost per unit area while maintaining reasonable efficiency. In this regard, superstrate CdTe/CdS solar cells are extensively investigated because of their suitable bandgap alignments, cost effective methods of production at large scales and stability against proton/electron irradiation. The conversion efficiencies in the range of 6-20% are achieved by structuring the device by varying the absorber/window layer thickness, junction activation/annealing steps, with more suitable front/back contacts, preparation techniques, doping with foreign ions, etc. This review focuses on fundamental and critical aspects like: (a) choice of CdS window layer and CdTe absorber layer; (b) drawbacks associated with the device including environmental problems, optical absorption losses and back contact barriers; (c) structural dynamics at CdS-CdTe interface; (d) influence of junction activation process by CdCl2 or HCF2Cl treatment; (e) interface and grain boundary passivation effects; (f) device degradation due to impurity diffusion and stress; (g) fabrication with suitable front and back contacts; (h) chemical processes occurring at various interfaces; (i) strategies and modifications developed to improve their efficiency. The complexity involved in understanding the multiple aspects of tuning the solar cell efficiency is reviewed in detail by considering the individual contribution from each component of the device. It is expected that this review article will enrich the materials aspects of CdTe/CdS devices for solar energy conversion and stimulate further innovative research interest on this intriguing topic.

Mechanical properties of glasses and TiO2 nanocrystal glass composites in BaO-TiO2-B2O3 system

Glasses and glass-nanocrystal (anatase TiO2) composites in BaO-TiO2-B2O3 system were fabricated by conventional melt-quenching technique and controlled heat treatment respectively. Poisson's ratio and Young's moduli were predicted through Makishima-Mackenzie theoretical equation for the as-quenched glasses by taking the four and three coordinated borons into account. Mechanical properties of the glasses and glass-nanocrystal composites were investigated in detail through nanoindentation and microindentation studies. Predicted Young's moduli of glasses were found to be in reasonable agreement with nanoindentation Measurements. Hardness and Young's modulus were enhanced with increasing volume fraction of nanocrystallites of TiO2 in glass matrix whereas fracture toughness was found susceptible to the surface features. The results were correlated to the structural units and nanocrystals present in the glasses. (C) 2013 Elsevier B.V. All rights reserved.

Experimental and simulation studies on the performance of standing wave thermoacoustic prime mover for pulse tube refrigerator

The thermoacoustic prime mover (TAPM) has gained considerable attention as a pressure wave generator to drive pulse tube refrigerator (PTR) due to no moving parts, reasonable efficiency, use of environmental friendly working fluids etc. To drive PTCs, lower frequencies (f) with larger pressure amplitudes (Delta P) are essential, which are affected by geometric and operating parameters of TAPM as well as working fluids. For driving PTRs, a twin standing wave TAPM is built and studied by using different working fluids such as helium, argon, nitrogen and their binary mixtures. Simulation results of DeltaEc are compared with experimental data wherever possible. DeltaEc predicts slightly increased resonance frequencies, but gives larger Delta P and lower temperature difference Delta T across stack. High mass number working fluid leads to lower frequency with larger Delta P, but higher Delta T. Studies indicate that the binary gas mixture of right composition with lower Delta T can be arrived at to drive TAPM of given geometry. (C) 2013 Elsevier Ltd and IIR. All rights reserved.

A New Voltage Stability Index based on the Tangent Vector of the Power Flow Jacobian

This paper presents a new voltage stability index based on the tangent vector of the power flow jacobian. This index is capable of providing the relative vulnerability information of the system buses from the point of view of voltage collapse. In an effort to compare this index with a similar index, the popular voltage stability index L is studied and it is shown through system studies that the L index is not a very consistent indicator of the voltage collapse point of the system but is only a reasonable indicator of the vulnerability of the system buses to voltage collapse. We also show that the new index can be used in the voltage stability analysis of radial systems which is not possible with the L index. This is a significant result of this investigation since there is a lot of contemporary interest in distributed generation and microgrids which are by and large radial in nature. Simulation results considering several test systems are provided to validate the results and the computational needs of the proposed scheme is assessed in comparison with other schemes

Minimizing total weighted tardiness on heterogeneous batch processors with incompatible job families

In this paper, we address a scheduling problem for minimizing total weighted tardiness. The background for the paper is derived from the automobile gear manufacturing process. We consider the bottleneck operation of heat treatment stage of gear manufacturing. Real-life scenarios like unequal release times, incompatible job families, nonidentical job sizes, heterogeneous batch processors, and allowance for job splitting have been considered. We have developed a mathematical model which takes into account dynamic starting conditions. The problem considered in this study is NP-hard in nature, and hence heuristic algorithms have been proposed to address it. For real-life large-size problems, the performance of the proposed heuristic algorithms is evaluated using the method of estimated optimal solution available in literature. Extensive computational analyses reveal that the proposed heuristic algorithms are capable of consistently obtaining near-optimal statistically estimated solutions in very reasonable computational time.