A parallel multistep predictor-corrector algorithm for solving ordinary differential equations

In this paper we give a generalized predictor-corrector algorithm for solving ordinary differential equations with specified initial values. The method uses multiple correction steps which can be carried out in parallel with a prediction step. The proposed method gives a larger stability interval compared to the existing parallel predictor-corrector methods. A method has been suggested to implement the algorithm in multiple processor systems with efficient utilization of all the processors.

TCP: thread contention predictor for parallel programs

With proliferation of chip multicores (CMPs) on desktops and embedded platforms, multi-threaded programs have become ubiquitous. Existence of multiple threads may cause resource contention, such as, in on-chip shared cache and interconnects, depending upon how they access resources. Hence, we propose a tool - Thread Contention Predictor (TCP) to help quantify the number of threads sharing data and their sharing pattern. We demonstrate its use to predict a more profitable shared, last level on-chip cache (LLC) access policy on CMPs. Our cache configuration predictor is 2.2 times faster compared to the cycle-accurate simulations. We also demonstrate its use for identifying hot data structures in a program which may cause performance degradation due to false data sharing. We fix layout of such data structures and show up-to 10% and 18% improvement in execution time and energy-delay product (EDP), respectively.

STAT-1 expression is regulated by IGFBP-3 in malignant glioma cells and is a strong predictor of poor survival in patients with glioblastoma Laboratory investigation

Object. Insulin-like growth factor binding proteins (IGEBPs) have been implicated in the pathogenesis of glioma. In a previous study the authors demonstrated that IGFBP-3 is a novel glioblastoma biomarker associated with poor survival. Since signal transducer and activator of transcription 1 (STAT-1) has been shown to be regulated by IGFBP-3 during chondrogenesis and is a prosurvival and radioresistant molecule in different tumors, the aim in the present study was to explore the functional significance of IGFBP-3 in malignant glioma cells, to determine if STAT-1 is indeed regulated by IGFBP-3, and to study the potential of STAT-1 as a biomarker in glioblastoma. Methods. The functional significance of IGFBP-3 was investigated using the short hairpin (sh)RNA gene knockdown approach on U251MG cells. STAT-1 regulation by IGFBP-3 was tested on U251MG and U87MG cells by shRNA gene knockdown and exogenous treatment with recombinant IGFBP-3 protein. Subsequently, the expression of STAT-1 was analyzed with real-time reverse transcription polymerase chain reaction (RT-PCR) and immunohistochemistry (IHC) in glioblastoma and control brain tissues. Survival analyses were done on a uniformly treated prospective cohort of adults with newly diagnosed glioblastoma (136 patients) using Kaplan-Meier and Cox regression models. Results. IGFBP-3 knockdown significantly impaired proliferation, motility, migration, and invasive capacity of U251MG cells in vitro (p < 0.005). Exogenous overexpression of IGFBP-3 in U251MG and U87MG cells demonstrated STAT-1 regulation. The mean transcript levels (by real-time RT-PCR) and the mean labeling index of STAT-1 (by IHC) were significantly higher in glioblastoma than in control brain tissues (p = 0.0239 and p < 0.001, respectively). Multivariate survival analysis revealed that STAT-1 protein expression (HR 1.015, p = 0.033, 95% CI 1.001-1.029) along with patient age (HR 1.025, p = 0.005, 95% CI 1.008-1.042) were significant predictors of shorter survival in patients with glioblastoma. Conclusions. IGFBP-3 influences tumor cell proliferation, migration, and invasion and regulates STAT-1 expression in malignant glioma cells. STAT-1 is overexpressed in human glioblastoma tissues and emerges as a novel prognostic biomarker.

Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions

We present a generalization of the finite volume evolution Galerkin scheme [M. Lukacova-Medvid'ova,J. Saibertov'a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533-562; M. Luacova-Medvid'ova, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1-30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor-corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme.

Role of predictors in downscaling surface temperature to river basin in India for IPCC SRES scenarios using support vector machine

In this paper, downscaling models are developed using a support vector machine (SVM) for obtaining projections of monthly mean maximum and minimum temperatures (T-max and T-min) to river-basin scale. The effectiveness of the model is demonstrated through application to downscale the predictands for the catchment of the Malaprabha reservoir in India, which is considered to be a climatically sensitive region. The probable predictor variables are extracted from (1) the National Centers for Environmental Prediction (NCEP) reanalysis dataset for the period 1978-2000, and (2) the simulations from the third-generation Canadian Coupled Global Climate Model (CGCM3) for emission scenarios A1B, A2, B1 and COMMIT for the period 1978-2100. The predictor variables are classified into three groups, namely A, B and C. Large-scale atmospheric variables Such as air temperature, zonal and meridional wind velocities at 925 nib which are often used for downscaling temperature are considered as predictors in Group A. Surface flux variables such as latent heat (LH), sensible heat, shortwave radiation and longwave radiation fluxes, which control temperature of the Earth's surface are tried as plausible predictors in Group B. Group C comprises of all the predictor variables in both the Groups A and B. The scatter plots and cross-correlations are used for verifying the reliability of the simulation of the predictor variables by the CGCM3 and to Study the predictor-predictand relationships. The impact of trend in predictor variables on downscaled temperature was studied. The predictor, air temperature at 925 mb showed an increasing trend, while the rest of the predictors showed no trend. The performance of the SVM models that are developed, one for each combination of predictor group, predictand, calibration period and location-based stratification (land, land and ocean) of climate variables, was evaluated. In general, the models which use predictor variables pertaining to land surface improved the performance of SVM models for downscaling T-max and T-min

Comparison of Prediction Based LSF Quantization Methods using Split VQ

Further improvement in performance, to achieve near transparent quality LSF quantization, is shown to be possible by using a higher order two dimensional (2-D) prediction in the coefficient domain. The prediction is performed in a closed-loop manner so that the LSF reconstruction error is the same as the quantization error of the prediction residual. We show that an optimum 2-D predictor, exploiting both inter-frame and intra-frame correlations, performs better than existing predictive methods. Computationally efficient split vector quantization technique is used to implement the proposed 2-D prediction based method. We show further improvement in performance by using weighted Euclidean distance.

Middle aged wasps mate through most of the year, without regard to body size, ovarian development and nestmateship: a laboratory study of the primitively eusocial wasp Ropalidia marginata

We studied the mating behaviour of the primi-tively eusocial wasp Ropalidia marginata and the factors that may influence sperm transfer. By introducing a male and a female R. marginata into ventilated transparent plastic boxes, we were able to observe mating behaviour, and it involved mounting and short or long conjugation of the wasps. Dissection of female wasps after the observation indicated that long conjugation is a good behavioural predictor of sperm transfer. This finding makes it possible to obtain mated females without dissecting them every time. We tested the effect of age, season, relatedness, body size and female's ovarian status on mating. Under laboratory conditions, mating success declined rapidly below and above the ages 5-20 days. Within this age range mating success was significantly low in December compared to other months tested. There was no nestmate discrimination, and there was no influence of male and female body size or of the ovarian state of the female on the probability of sperm transfer.

Direct numerical simulation of autoignition in a non-premixed, turbulent medium

In this paper, direct numerical simulation of autoignition in an initially non-premixed medium under isotropic, homogeneous, and decaying turbulence is presented. The pressure-based method developed herein is a spectral implementation of the sequential steps followed in the predictor-corrector type of algorithms; it includes the effects of density fluctuations caused by spatial inhomogeneities ill temperature and species. The velocity and pressure field are solved in the spectral space while the scalars and density field are solved in the physical space. The presented results reveal that the autoignition spots originate and evolve at locations where (1) the composition corresponds to a small range around a specific mixture fraction, and (2) the conditional scaler dissipation rate is low. A careful examination of the data obtained indicates that the autoignition spots originate in the vortex cores, and the hot gases travel outward as combustion progresses. Hence, the applicability of the transient laminar flamelet model for this problem is questioned. The dependence of autoignition characteristics on parameters such as (1) die initial eddy-turnover time and (2) the initial ratio of length scale of scalars to that of velocities are investigated. Certain implications of new results on the conditional moment closure modeling are discussed.

A statistical analysis of audience sound absorption (A)

The absorption produced by the audience in concert halls is considered a random variable. Beranek's proposal [L. L. Beranek, Music, Acoustics and Architecture (Wiley, New York, 1962), p. 543] that audience absorption is proportional to the area they occupy and not to their number is subjected to a statistical hypothesis test. A two variable linear regression model of the absorption with audience area and residual area as regressor variables is postulated for concert halls without added absorptive materials. Since Beranek's contention amounts to the statement that audience absorption is independent of the seating density, the test of the hypothesis lies in categorizing halls by seating density and examining for significant differences among slopes of regression planes of the different categories. Such a test shows that Beranek's hypothesis can be accepted. It is also shown that the audience area is a better predictor of the absorption than the audience number. The absorption coefficients and their 95% confidence limits are given for the audience and residual areas. A critique of the regression model is presented.

A Ten-microRNA Expression Signature Predicts Survival in Glioblastoma

Glioblastoma (GBM) is the most common and aggressive primary brain tumor with very poor patient median survival. To identify a microRNA (miRNA) expression signature that can predict GBM patient survival, we analyzed the miRNA expression data of GBM patients (n = 222) derived from The Cancer Genome Atlas (TCGA) dataset. We divided the patients randomly into training and testing sets with equal number in each group. We identified 10 significant miRNAs using Cox regression analysis on the training set and formulated a risk score based on the expression signature of these miRNAs that segregated the patients into high and low risk groups with significantly different survival times (hazard ratio HR] = 2.4; 95% CI = 1.4-3.8; p < 0.0001). Of these 10 miRNAs, 7 were found to be risky miRNAs and 3 were found to be protective. This signature was independently validated in the testing set (HR = 1.7; 95% CI = 1.1-2.8; p = 0.002). GBM patients with high risk scores had overall poor survival compared to the patients with low risk scores. Overall survival among the entire patient set was 35.0% at 2 years, 21.5% at 3 years, 18.5% at 4 years and 11.8% at 5 years in the low risk group, versus 11.0%, 5.5%, 0.0 and 0.0% respectively in the high risk group (HR = 2.0; 95% CI = 1.4-2.8; p < 0.0001). Cox multivariate analysis with patient age as a covariate on the entire patient set identified risk score based on the 10 miRNA expression signature to be an independent predictor of patient survival (HR = 1.120; 95% CI = 1.04-1.20; p = 0.003). Thus we have identified a miRNA expression signature that can predict GBM patient survival. These findings may have implications in the understanding of gliomagenesis, development of targeted therapy and selection of high risk cancer patients for adjuvant therapy.

Endophytic fungal communities in woody perennials of three tropical forest types of the Western Ghats, southern India

Fungal endophytes of tropical trees are expected to be exceptionally species rich as a consequence of high tree diversity in the tropics and the purported host restriction among the endophytes. Based on this premise, endophytes have been regarded as a focal group for estimating fungal numbers because their possible hyperdiverse nature would reflect significantly global fungal diversity. We present our consolidated ten-year work on 75 dicotyledonous tree hosts belonging to 33 families and growing in three different types of tropical forests of the NBR in the Western Ghats, southern India. We conclude that endophyte diversity in these forests is limited due to loose host affiliations among endophytes. Some endophytes have a wide host range and colonize taxonomically disparate hosts suggesting adaptations in them to counter a variety of defense chemicals in their hosts. Furthermore, such polyphagous endophytes dominate the endophyte assemblages of different tree hosts. Individual leaves may be densely colonized but only by a few endophyte species. It appears that the environment (the type of forest in this case) has a larger role in determining the endophyte assemblage of a plant host than the taxonomy of the host plant. Thus, different tropical plant communities have to be studied for their endophyte diversity to test the generalization that endophytes are hyperdiverse in the tropics, estimate their true species richness, and use them as a predictor group for more accurate assessment of global fungal diversity.

Post-seizure EEG fractal dimension of first ECT predicts antidepressant response at two weeks

Background: Duration of seizure by itself is an insufficient criterion for a therapeutically adequate seizure in ECT. Therefore, measures of seizure EEG other than its duration need to be explored as indices of seizure adequacy and predictors of treatment response. We measured the EEG seizure using a geometrical method-fractal dimension (FD) and examined if this measure predicted remission. Methods: Data from an efficacy study on melancholic depressives (n = 40) is used for the present exploration. They received thrice or once weekly ECTs, each schedule at two energy levels - high or low energy level. FD was computed for early-, mid- and post-seizure phases of the ictal EEG. Average of the two channels was used for analysis. Results: Two-thirds of the patients (n = 25) were remitted at the end of 2 weeks. As expected, a significantly higher proportion of patients receiving thrice weekly ECT remitted than in patients receiving once weekly ECT. Smaller post-seizure FD at first ECT is the only variable which predicted remission status after six ECTs. within the once weekly ECT group too, smaller post-seizure FD was associated with remission status. Conclusions: Post-seizure FD is proposed as a novel measure of seizure adequacy and predictor of treatment response. Clinical implications: Seizure measures at first ECT may guide selection of ECT schedule to optimize ECT. Limitations: The study examined short term antidepressant effects only. The results may not be generalized to medication-resistant depressives. (C) 1999 Elsevier Science B.V. All rights reserved.

A new parallel overlapped domain decomposition method for nonlinear dynamic finite element analysis

In this paper a new parallel algorithm for nonlinear transient dynamic analysis of large structures has been presented. An unconditionally stable Newmark-beta method (constant average acceleration technique) has been employed for time integration. The proposed parallel algorithm has been devised within the broad framework of domain decomposition techniques. However, unlike most of the existing parallel algorithms (devised for structural dynamic applications) which are basically derived using nonoverlapped domains, the proposed algorithm uses overlapped domains. The parallel overlapped domain decomposition algorithm proposed in this paper has been formulated by splitting the mass, damping and stiffness matrices arises out of finite element discretisation of a given structure. A predictor-corrector scheme has been formulated for iteratively improving the solution in each step. A computer program based on the proposed algorithm has been developed and implemented with message passing interface as software development environment. PARAM-10000 MIMD parallel computer has been used to evaluate the performances. Numerical experiments have been conducted to validate as well as to evaluate the performance of the proposed parallel algorithm. Comparisons have been made with the conventional nonoverlapped domain decomposition algorithms. Numerical studies indicate that the proposed algorithm is superior in performance to the conventional domain decomposition algorithms. (C) 2003 Elsevier Ltd. All rights reserved.

Correction for Bias in Downscaling GCM Simulations for Hydrologic Impact Assessment

In a statistical downscaling model, it is important to remove the bias of General Circulations Model (GCM) outputs resulting from various assumptions about the geophysical processes. One conventional method for correcting such bias is standardisation, which is used prior to statistical downscaling to reduce systematic bias in the mean and variances of GCM predictors relative to the observations or National Centre for Environmental Prediction/ National Centre for Atmospheric Research (NCEP/NCAR) reanalysis data. A major drawback of standardisation is that it may reduce the bias in the mean and variance of the predictor variable but it is much harder to accommodate the bias in large-scale patterns of atmospheric circulation in GCMs (e.g. shifts in the dominant storm track relative to observed data) or unrealistic inter-variable relationships. While predicting hydrologic scenarios, such uncorrected bias should be taken care of; otherwise it will propagate in the computations for subsequent years. A statistical method based on equi-probability transformation is applied in this study after downscaling, to remove the bias from the predicted hydrologic variable relative to the observed hydrologic variable for a baseline period. The model is applied in prediction of monsoon stream flow of Mahanadi River in India, from GCM generated large scale climatological data.

3D Image Reconstruction From Truncated Helical Cone Beam Projection Data - A Linear Prediction Approach

With the introduction of 2D flat-panel X-ray detectors, 3D image reconstruction using helical cone-beam tomography is fast replacing the conventional 2D reconstruction techniques. In 3D image reconstruction, the source orbit or scanning geometry should satisfy the data sufficiency or completeness condition for exact reconstruction. The helical scan geometry satisfies this condition and hence can give exact reconstruction. The theoretically exact helical cone-beam reconstruction algorithm proposed by Katsevich is a breakthrough and has attracted interest in the 3D reconstruction using helical cone-beam Computed Tomography.In many practical situations, the available projection data is incomplete. One such case is where the detector plane does not completely cover the full extent of the object being imaged in lateral direction resulting in truncated projections. This result in artifacts that mask small features near to the periphery of the ROI when reconstructed using the convolution back projection (CBP) method assuming that the projection data is complete. A number of techniques exist which deal with completion of missing data followed by the CBP reconstruction. In 2D, linear prediction (LP)extrapolation has been shown to be efficient for data completion, involving minimal assumptions on the nature of the data, producing smooth extensions of the missing projection data.In this paper, we propose to extend the LP approach for extrapolating helical cone beam truncated data. The projection on the multi row flat panel detectors has missing columns towards either ends in the lateral direction in truncated data situation. The available data from each detector row is modeled using a linear predictor. The available data is extrapolated and this completed projection data is backprojected using the Katsevich algorithm. Simulation results show the efficacy of the proposed method.