A fast cross-entropy method for estimating buffer overflows in queuing networks

In this paper, we propose a fast adaptive importance sampling method for the efficient simulation of buffer overflow probabilities in queueing networks. The method comprises three stages. First, we estimate the minimum cross-entropy tilting parameter for a small buffer level; next, we use this as a starting value for the estimation of the optimal tilting parameter for the actual (large) buffer level. Finally, the tilting parameter just found is used to estimate the overflow probability of interest. We study various properties of the method in more detail for the M/M/1 queue and conjecture that similar properties also hold for quite general queueing networks. Numerical results support this conjecture and demonstrate the high efficiency of the proposed algorithm.

Automated regional myocardial displacement for facilitating the interpretation of dobutamine echocardiography

Quantification of stress echocardiography may overcome the training requirements and subjective nature of visual wall motion score (WMS) assessment, but quantitative approaches may be difficult to apply and require significant time for image processing. The integral of long-axis myocardial velocity is displacement, which may be represented as a color map over the left ventricular myocardium. This study was designed to explore the feasibility and accuracy of measuring long-axis myocardial displacement, derived from tissue Doppler, for the detection of coronary artery disease (CAD) during dobutamine stress echocardiography (DBE). One hundred thirty patients underwent standard DBE, including 30 patients at low risk of CAD, 30 patients with normal coronary angiography (both groups studied to define normal ranges of displacement), and 70 patients who underwent coronary angiography in whom the accuracy of normal ranges was tested. Regional myocardial displacement was obtained by analysis of color tissue Doppler apical images acquired at peak stress. Displacement was compared with WMS, and with the presence of CAD by angiography. The analysis time was 3.2 +/- 1.5 minutes per patient. Segmental displacement was correlated with wall motion (normal 7.4 +/- 3.2 mm, ischemia 5.8 +/- 4.2 mm, viability 4.6 +/- 3.0 mm, scar 4.5 +/- 3.5 mm, p <0.001). Reversal of normal base-apex displacement was an insensitive (19%) but specific (90%) marker of CAD. The sum of displacements within each vascular territory had a sensitivity and specificity of 89% and 79%, respectively, for prediction of significant CAD, compared with 86% and 78%, respectively, for WMS (p = NS). The displacements in the basal segments had a sensitivity and specificity of 83% and 78%, respectively (p = NS). Regional myocardial displacement during DBE is feasible and offers a fast and accurate method for the diagnosis of CAD. (C),2002 by Excerpta Medica, Inc.

Fast Structure-Based Assignment of 15N HSQC Spectra of Selectively 15N-Labeled Paramagnetic Proteins

A novel strategy for fast NMR resonance assignment of N-15 HSQC spectra of proteins is presented. It requires the structure coordinates of the protein, a paramagnetic center, and one or more residue-selectively N-15-labeled samples. Comparison of sensitive undecoupled N-15 HSQC spectra recorded of paramagnetic and diamagnetic samples yields data for every cross-peak on pseudocontact shift, paramagnetic relaxation enhancement, cross-correlation between Curie-spin and dipole-dipole relaxation, and residual dipolar coupling. Comparison of these four different paramagnetic quantities with predictions from the three-dimensional structure simultaneously yields the resonance assignment and the anisotropy of the susceptibility tensor of the paramagnetic center. The method is demonstrated with the 30 kDa complex between the N-terminal domain of the epsilon subunit and the theta subunit of Escherichia Coll DNA polymerase III. The program PLATYPUS was developed to perform the assignment, provide a measure of reliability of the assignment, and determine the susceptibility tensor anisotropy.

Operational regimes for a simplified one-step artificial chaperone refolding method

The artificial chaperone method for protein refolding developed by Rozema et al. (Rozema, D.; Gellman, S. H. J. Am. Chem. Soc. 1995, 117 (8), 2373-2374) involves the sequential dilution of denatured protein into a buffer containing detergent (cetyltrimethylammonium bromide, CTAB) and then into a refolding buffer containing cyclodextrin WD). In this paper a simplified one-step artificial chaperone method is reported, whereby CTAB is added directly to the denatured solution, which is then diluted directly into a refolding buffer containing P-cyclodextrin (P-CD). This new method can be applied at high protein concentrations, resulting in smaller processing volumes and a more concentrated protein solution following refolding. The increase in achievable protein concentration results from the enhanced solubility of CTAB at elevated temperatures in concentrated denaturant. The refolding yields obtained for the new method were significantly higher than for control experiments lacking additives and were comparable to the yields obtained with the classical two-step approach. A study of the effect of beta-CD and CTAB concentrations on refolding yield suggested two operational regimes: slow stripping ( beta-CDXTABsimilar to1), most suited for higher protein concentrations, and fast stripping (beta-CD/CTABsimilar to2.7), best suited for lower protein concentrations. An increased chaotrope concentration resulted in higher refolding yields and an enlarged operational regime.

Fast-convergent iterative scheme for filtering velocity signals and finding Kolmogorov scales

The present fundamental knowledge of fluid turbulence has been established primarily from hot- and cold-wire measurements. Unfortunately, however, these measurements necessarily suffer from contamination by noise since no certain method has previously been available to optimally filter noise from the measured signals. This limitation has impeded our progress of understanding turbulence profoundly. We address this limitation by presenting a simple, fast-convergent iterative scheme to digitally filter signals optimally and find Kolmogorov scales definitely. The great efficacy of the scheme is demonstrated by its application to the instantaneous velocity measured in a turbulent jet.

On the fast approximation of Green's functions in MPIE formulations for planar layered media

The numerical implementation of the complex image approach for the Green's function of a mixed-potential integralequation formulation is examined and is found to be limited to low values of k(0) rho (in this context k(0) rho = 2 pirho/ lambda(0), where rho is the distance between the source and the field points of the Green's function and lambda(0) is the free space wavelength). This is a clear limitation for problems of large dimension or high frequency where this limit is easily exceeded. This paper examines the various strategies and proposes a hybrid method whereby most of the above problems can be avoided. An efficient integral method that is valid for large k(0) rho is combined with the complex image method in order to take advantage of the relative merits of both schemes. It is found that a wide overlapping region exists between the two techniques allowing a very efficient and consistent approach for accurately calculating the Green's functions. In this paper, the method developed for the computation of the Green's function is used for planar structures containing both lossless and lossy media.

A design method for superconducting MRI magnets with ferromagnetic material

The emphasis of this work is on the optimal design of MRI magnets with both superconducting coils and ferromagnetic rings. The work is directed to the automated design of MRI magnet systems containing superconducting wire and both `cold' and `warm' iron. Details of the optimization procedure are given and the results show combined superconducting and iron material MRI magnets with excellent field characteristics. Strong, homogeneous central magnetic fields are produced with little stray or external field leakage. The field calculations are performed using a semi-analytical method for both current coil and iron material sources. Design examples for symmetric, open and asymmetric clinical MRI magnets containing both superconducting coils and ferromagnetic material are presented.

Database technologies for l-system simulations in virtual plant applications on bioinformatics

One of the most important advantages of database systems is that the underlying mathematics is rich enough to specify very complex operations with a small number of statements in the database language. This research covers an aspect of biological informatics that is the marriage of information technology and biology, involving the study of real-world phenomena using virtual plants derived from L-systems simulation. L-systems were introduced by Aristid Lindenmayer as a mathematical model of multicellular organisms. Not much consideration has been given to the problem of persistent storage for these simulations. Current procedures for querying data generated by L-systems for scientific experiments, simulations and measurements are also inadequate. To address these problems the research in this paper presents a generic process for data-modeling tools (L-DBM) between L-systems and database systems. This paper shows how L-system productions can be generically and automatically represented in database schemas and how a database can be populated from the L-system strings. This paper further describes the idea of pre-computing recursive structures in the data into derived attributes using compiler generation. A method to allow a correspondence between biologists' terms and compiler-generated terms in a biologist computing environment is supplied. Once the L-DBM gets any specific L-systems productions and its declarations, it can generate the specific schema for both simple correspondence terminology and also complex recursive structure data attributes and relationships.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Enhancing experiential learning via virtual environments

Research indicates Virtual Reality (VR) is delivering on it's promised potential to provide enhanced training and education outcomes. A significant research project, at the University of Queensland, has constructed a number of virtual contexts in which the phenomena experienced by patients who have psychosis are reproduced for use in psychiatry education. Symptoms of psychosis reproduced include delusions, hallucinations and thought disorder. The new software enables psychiatry students to experience the inner world of a patient with psychosis. Lecturers in psychiatry report VR has the potential to enhance student's abilities to actually 'feel' the types of emotions and physiological reactions a hallucination precipitates in a patient. The current work of the project and stages of software development will be demonstrated. The virtual environments provide a new method of delivering experiential learning opportunities to higher education classrooms.

The Equilibrium Flux Method for the calculation of flows with non-equilbrium chemical reactions

The Equilibrium Flux Method [1] is a kinetic theory based finite volume method for calculating the flow of a compressible ideal gas. It is shown here that, in effect, the method solves the Euler equations with added pseudo-dissipative terms and that it is a natural upwinding scheme. The method can be easily modified so that the flow of a chemically reacting gas mixture can be calculated. Results from the method for a one-dimensional non-equilibrium reacting flow are shown to agree well with a conventional continuum solution. Results are also presented for the calculation of a plane two-dimensional flow, at hypersonic speed, of a dissociating gas around a blunt-nosed body.

A wavelet-based adaptive technique for adsorption problems involving steep gradients

A general, fast wavelet-based adaptive collocation method is formulated for heat and mass transfer problems involving a steep moving profile of the dependent variable. The technique of grid adaptation is based on sparse point representation (SPR). The method is applied and tested for the case of a gas–solid non-catalytic reaction in a porous solid at high Thiele modulus. Accurate and convergent steep profiles are obtained for Thiele modulus as large as 100 for the case of slab and found to match the analytical solution.

Towards realistic simulations of lava dome growth using the level set method

The level set method has been implemented in a computational volcanology context. New techniques are presented to solve the advection equation and the reinitialisation equation. These techniques are based upon an algorithm developed in the finite difference context, but are modified to take advantage of the robustness of the finite element method. The resulting algorithm is tested on a well documented Rayleigh–Taylor instability benchmark [19], and on an axisymmetric problem where the analytical solution is known. Finally, the algorithm is applied to a basic study of lava dome growth.

A rapid single-step multiplex method for discriminating between Trichogramma (Hymenoptera: Trichogrammatidae) species in Australia

Inaccurate species identification confounds insect ecological studies. Examining aspects of Trichogramma ecology pertinent to the novel insect resistance management strategy for future transgenic cotton, Gossypium hirsutum L., production in the Ord River Irrigation Area (ORIA) of Western Australia required accurate differentiation between morphologically similar Trichogramma species. Established molecular diagnostic methods for Trichogramma identification use species-specific sequence difference in the internal transcribed spacer (ITS)-2 chromosomal region; yet, difficulties arise discerning polymerase chain reaction (PCR) fragments of similar base pair length by gel electrophoresis. This necessitates the restriction enzyme digestion of PCR-amplified ITS-2 fragments to readily differentiate Trichogramma australicum Girault and Trichogramma pretiosum Riley. To overcome the time and expense associated with a two-step diagnostic procedure, we developed a “one-step” multiplex PCR technique using species-specific primers designed to the ITS-2 region. This approach allowed for a high-throughput analysis of samples as part of ongoing ecological studies examining Trichogramma biological control potential in the ORIA where these two species occur in sympatry.