A numerical study of large sparse matrix exponentials arising in Markov chains

Krylov subspace techniques have been shown to yield robust methods for the numerical computation of large sparse matrix exponentials and especially the transient solutions of Markov Chains. The attractiveness of these methods results from the fact that they allow us to compute the action of a matrix exponential operator on an operand vector without having to compute, explicitly, the matrix exponential in isolation. In this paper we compare a Krylov-based method with some of the current approaches used for computing transient solutions of Markov chains. After a brief synthesis of the features of the methods used, wide-ranging numerical comparisons are performed on a power challenge array supercomputer on three different models. (C) 1999 Elsevier Science B.V. All rights reserved.AMS Classification: 65F99; 65L05; 65U05.

Hypoxia stimulates cerebral blood flow in the estuarine crocodile (Crocodylus porosus)

The effect of N-2 respiration on cerebral blood flow (CBF) velocity on the dorsal surface of cerebellum was examined in the estuarine crocodile, Crocodylus porosus, using epi-illumination microscopy. Twelve minutes of N-2 respiration resulted in a 126% increase in CBF velocity. N-2 respiration had no effect on blood pressure, indicating an underlying cerebral vasodilation. In addition, heart rate increased significantly. Systemic injections of aminophylline and the NO synthase (NOS) inhibitor nitro-L-arginine (L-NA) did not affect the hypoxia induced increase in CBF. We conclude that C. porosus responds to hypoxia with adenosine and nitric oxide (NO) independent cerebral vasodilation, and that this is likely to be a mechanism protecting the brain from energy deficiency during prolonged dives. (C) 1999 Elsevier Science Ireland Ltd. All rights reserved.

Determining forest structural attributes using an inverted geometric-optical model in mixed eucalypt forests, Southeast Queensland, Australia

The Montreal Process indicators are intended to provide a common framework for assessing and reviewing progress toward sustainable forest management. The potential of a combined geometrical-optical/spectral mixture analysis model was assessed for mapping the Montreal Process age class and successional age indicators at a regional scale using Landsat Thematic data. The project location is an area of eucalyptus forest in Emu Creek State Forest, Southeast Queensland, Australia. A quantitative model relating the spectral reflectance of a forest to the illumination geometry, slope, and aspect of the terrain surface and the size, shape, and density, and canopy size. Inversion of this model necessitated the use of spectral mixture analysis to recover subpixel information on the fractional extent of ground scene elements (such as sunlit canopy, shaded canopy, sunlit background, and shaded background). Results obtained fron a sensitivity analysis allowed improved allocation of resources to maximize the predictive accuracy of the model. It was found that modeled estimates of crown cover projection, canopy size, and tree densities had significant agreement with field and air photo-interpreted estimates. However, the accuracy of the successional stage classification was limited. The results obtained highlight the potential for future integration of high and moderate spatial resolution-imaging sensors for monitoring forest structure and condition. (C) Elsevier Science Inc., 2000.

Response programming in dementia of the Alzheimer type: A kinematic analysis

Although planning is important for the functioning of patients with dementia of the Alzheimer Type (DAT), little is known about response programming in DAT. This study used a cueing paradigm coupled with quantitative kinematic analysis to document the preparation and execution of movements made by a group of 12 DAT patients and their age and sex matched controls. Participants connected a series of targets placed upon a WACOM SD420 graphics tablet, in response to the pattern of illumination of a set of light emitting diodes (LEDs). In one condition, participants could programme the upcoming movement, whilst in another they were forced to reprogramme this movement on-line (i.e. they were not provided with advance information about the location of the upcoming target). DAT patients were found to have programming deficits, taking longer to initiate movements; particularly in the absence of cues. While problems spontaneously programming a movement might cause a greater reliance upon on-line guidance, when both groups were required to guide the movement on-line, DAT patients continued to show slower and less efficient movements implying declining sensori-motor function; these differences were not simply due to strategy or medication status. (C) 1997 Elsevier Science Ltd.

A low-temperature insulating phase at v=1.5 for 2D holes in high-mobility SiSi1-xGex heterostructures with Landau level degeneracy

Magneto-transport measurements of the 2D hole system (2DHS) in p-type Si-Si1-xGex heterostructures identify the integer quantum Hall effect (IQHE) at dominantly odd-integer filling factors v and two low-temperature insulating phases (IPs) at v = 1.5 and v less than or similar to 0.5, with re-entrance to the quantum Hall effect at v = 1. The temperature dependence, current-voltage characteristics, and tilted field and illumination responses of the IP at v = 1.5 indicate that the important physics is associated with an energy degeneracy of adjacent Landau levels of opposite spin, which provides a basis for consideration of an intrinsic, many-body origin.

Stability of linear difference and differential inclusions

Any given n X n matrix A is shown to be a restriction, to the A-invariant subspace, of a nonnegative N x N matrix B of spectral radius p(B) arbitrarily close to p(A). A difference inclusion x(k+1) is an element of Ax(k), where A is a compact set of matrices, is asymptotically stable if and only if A can be extended to a set B of nonnegative matrices B with \ \B \ \ (1) < 1 or \ \B \ \ (infinity) < 1. Similar results are derived for differential inclusions.

The visual ecology of avian photoreceptors

The spectral sensitivities of avian retinal photoreceptors are examined with respect to microspectrophotometric measurements of single cells, spectrophotometric measurements of extracted or in vitro regenerated visual pigments, and molecular genetic analyses of visual pigment opsin protein sequences. Bird species from diverse orders are compared in relation to their evolution, their habitats and the multiplicity of visual tasks they must perform. Birds have five different types of visual pigment and seven different types of photo receptor-rods, double (uneven twin) cones and four types of single cone. The spectral locations of the wavelengths of maximum absorbance (lambda (max)) of the different visual pigments, and the spectral transmittance characteristics of the intraocular spectral filters (cone oil droplets) that also determine photoreceptor spectral sensitivity, vary according to both habitat and phylogenetic relatedness. The primary influence on avian retinal design appears to be the range of wavelengths available for vision, regardless of whether that range is determined by the spectral distribution of the natural illumination or the spectral transmittance of the ocular media (cornea, aqueous humour, lens, vitreous humour). Nevertheless, other variations in spectral sensitivity exist that reflect the variability and complexity of avian visual ecology. (C) 2001 Elsevier Science Ltd. All rights reserved.

Colourful objects through animal eyes

To understand how bees, birds, and fish may use colour vision for food selection and mate choice, we reconstructed views of biologically important objects taking into account the receptor spectral sensitivities. Reflectance spectra a of flowers, bird plumage, and fish skin were used to calculate receptor quantum catches. The quantum catches were then coded by red, green, and blue of a computer monitor; and powers, birds, and fish were visualized in animal colours. Calculations were performed for different illumination conditions. To simulate colour constancy, we used a von Kries algorithm, i.e., the receptor quantum catches were scaled so that the colour of illumination remained invariant. We show that on land this algorithm compensates reasonably well for changes of object appearance caused by natural changes of illumination, while in water failures of von Kries colour constancy are prominent. (C) 2000 John Wiley & Sons, Inc.

An ESR study on gamma-irradiated poly(vinyl alcohol)

The formation of radicals in poly(vinyl alcohol), PVA, powder irradiated at 77 K by gamma -rays and the transformations of these radicals during photolysis with visible wavelengths and on thermal annealing have been studied. After irradiation a four-line ESR spectrum was observed. It was assigned to a triplet of the C-alpha-radical (38%), with a splitting of 3.27 mT, superimposed on a doublet (62%) with a splitting of 2.7 mT. The doublet appears to be composed of two radicals, one of which is photo-bleachable (58%) and the other which is not photo-bleachable (42%). This suggests that the latter radical is a neutral radical. The photo-bleachable component of the doublet has been assigned to a carbonyl anion radical. but the second doublet due to a neutral radical is unassigned. The total G-value for formation of radicals at 77 K was found to be 2.41 +/- 0.03. Upon illumination with visible light, the anion radicals were removed and the doublet components or the spectrum diminished in intensity, while the three-line spectrum of the C-alpha-radical became more clearly visible. This transition was due to the photo-detachment of electrons from traps which were proposed to be located on carbonyl groups in the polymer resulting from incomplete hydrolysis of the vinyl acetate. The photo-decay of the anion radicals could be satisfactorily described by a two-stage process. The first stage comprised the decay of approximately 80% of the anion radicals present, while the second stage was associated with the decay of the remaining 20%. Subsequent thermal annealing of a photolysed sample to 290 K led to a change in the shape of the spectrum to form a more clearly defined triplet, As the doublet of the neutral radical decays on thermal annealing between 150 and 250K, the C-alpha-radical is formed. (C) 2001 Elsevier Science Ltd. All rights reserved.

A systematic approach to error isolation in computerized wastewater simulation models

Activated sludge models are used extensively in the study of wastewater treatment processes. While various commercial implementations of these models are available, there are many people who need to code models themselves using the simulation packages available to them, Quality assurance of such models is difficult. While benchmarking problems have been developed and are available, the comparison of simulation data with that of commercial models leads only to the detection, not the isolation of errors. To identify the errors in the code is time-consuming. In this paper, we address the problem by developing a systematic and largely automated approach to the isolation of coding errors. There are three steps: firstly, possible errors are classified according to their place in the model structure and a feature matrix is established for each class of errors. Secondly, an observer is designed to generate residuals, such that each class of errors imposes a subspace, spanned by its feature matrix, on the residuals. Finally. localising the residuals in a subspace isolates coding errors. The algorithm proved capable of rapidly and reliably isolating a variety of single and simultaneous errors in a case study using the ASM 1 activated sludge model. In this paper a newly coded model was verified against a known implementation. The method is also applicable to simultaneous verification of any two independent implementations, hence is useful in commercial model development.

Manual asymmetries in the preparation and control of goal-directed movements

The primary purpose of this experiment was to determine if left hand reaction time advantages in manual aiming result from a right hemisphere attentional advantage or an early right hemisphere role in movement preparation. Right-handed participants were required to either make rapid goal-directed movements to small targets or simply lift their hand upon target illumination. The amount of advance information about the target for a particular trial was manipulated by precuing a subset of potential targets prior to the reaction time interval. When participants were required to make aiming movements to targets in left space, the left hand enjoyed a reaction advantage that was not present for aiming in right space: or simple finger lifts. This advantage was independent of the amount or type of advance information provided by the precue. This finding supports the movement planning hypothesis. With respect to movement execution, participants completed their aiming movements more quickly when aiming with their right hand, particularly in right space. This right hand advantage in right space was due to the time required to decelerate the movement and to make feedback-based adjustments late in the movement trajectory. (C) 2001 Academic Press.

A comparison of low-storage strategies for spectral analysis in dissipative systems: filter diagonalisation in the Lanczos representation and harmonic inversion of the Chebychev-order-domain autocorrelation function

We compare the performance of two different low-storage filter diagonalisation (LSFD) strategies in the calculation of complex resonance energies of the HO2, radical. The first is carried out within a complex-symmetric Lanczos subspace representation [H. Zhang, S.C. Smith, Phys. Chem. Chem. Phys. 3 (2001) 2281]. The second involves harmonic inversion of a real autocorrelation function obtained via a damped Chebychev recursion [V.A. Mandelshtam, H.S. Taylor, J. Chem. Phys. 107 (1997) 6756]. We find that while the Chebychev approach has the advantage of utilizing real algebra in the time-consuming process of generating the vector recursion, the Lanczos, method (using complex vectors) requires fewer iterations, especially for low-energy part of the spectrum. The overall efficiency in calculating resonances for these two methods is comparable for this challenging system. (C) 2001 Elsevier Science B.V. All rights reserved.

Transmit array of transistor amplifiers illuminated by a patch array in the reactive near-field region

In this paper, a small transmit array of transistor amplifiers illuminated by a passive array of microstrip patches in the reactive near-field region is investigated as a power-combining structure. The two cases considered are when the transmit array radiates in a free space and when a passive array similar to the one used for illumination collects the radiated power. A comparison of the performance of the proposed structure against the alternative one, which uses a conventional horn antenna as a power-launching/receiving device, is also presented.

Investigations into a power-combining structure using a reflect array of dual-feed aperture-coupled microstrip patch antennas

This paper details an investigation of a power combiner that uses a reflect array of dual-feed aperture-coupled microstrip patch antennas and a corporate-fed dual-polarized array as a signal distributing/combining device. In this configuration, elements of the reflect array receive a linearly polarized wave and retransmit it with an orthogonal polarization using variable-length sections of microstrip lines connecting receive and transmit ports. By applying appropriate lengths of these delay lines, the array focuses the transmitted wave onto the feed array. The operation of the combiner is investigated for a small-size circular reflect array for the cases of -3 dB, -6 dB and -10 dB edge illumination by the 2 x 2-element dual-polarized array.

Spectroscopic study of the penetration depth of grafted polystyrene onto poly (tetrafluoroethylene-co-perfluoropropylvinylether) substrates. 1. Effect of grafting conditions

This study concerns the radiation grafting of styrene onto poly(tetrafluoroethylene-co-perfluoropropylvinylether) (PFA) substrates and the penetration depth of the graft. Grafting was obtained by the simultaneous irradiation method, and the spectroscopic analysis was made with the micro-Raman technique. Effects of grafting conditions such as the type of solvent, dose rate, and irradiation dose on the grafting yield were investigated. Of the different solvents used, the most efficient in terms of increasing grafting yield were dichloromethane, benzene, and methanol, respectively. A mixture of methanol and dichloromethane used as a solvent for styrene achieved a higher degree of grafting and concentration of grafted polystyrene onto the surface of PFA substrates than solutions of the monomer in the separate solvents. The degree of grafting increased with increasing radiation dose up to 500 kGy, stabilizing above this dose. However, the grafting yield decreased with an increase in the dose rate. The increase in the overall grafting yield was accompanied by a proportional increase in the penetration depth of the grafts into the substrate. (C) 2002 Wiley Periodicals, Inc.