Comparison of the computational cost of a Monte Carlo and deterministic algorithm for computing bilinear forms of matrix powers

In this paper we consider bilinear forms of matrix polynomials and show that these polynomials can be used to construct solutions for the problems of solving systems of linear algebraic equations, matrix inversion and finding extremal eigenvalues. An almost Optimal Monte Carlo (MAO) algorithm for computing bilinear forms of matrix polynomials is presented. Results for the computational costs of a balanced algorithm for computing the bilinear form of a matrix power is presented, i.e., an algorithm for which probability and systematic errors are of the same order, and this is compared with the computational cost for a corresponding deterministic method.

Uniformization of affine dependence programs for parallel embedded system design

This paper is concerned with the uniformization of a system of afine recurrence equations. This transformation is used in the design (or compilation) of highly parallel embedded systems (VLSI systolic arrays, signal processing filters, etc.). In this paper, we present and implement an automatic system to achieve uniformization of systems of afine recurrence equations. We unify the results from many earlier papers, develop some theoretical extensions, and then propose effective uniformization algorithms. Our results can be used in any high level synthesis tool based on polyhedral representation of nested loop computations.

The relative order and inverses of recurrent networks

Differential geometry is used to investigate the structure of neural-network-based control systems. The key aspect is relative order—an invariant property of dynamic systems. Finite relative order allows the specification of a minimal architecture for a recurrent network. Any system with finite relative order has a left inverse. It is shown that a recurrent network with finite relative order has a local inverse that is also a recurrent network with the same weights. The results have implications for the use of recurrent networks in the inverse-model-based control of nonlinear systems.

Analysis of disputed single-parent/child and sibling relationships using 16 STR loci

This study describes the validation of short tandem repeat (STR) systems for the resolution of cases of disputed parentage where only a single parent is available for testing or where the claimed relationship of both parents is in doubt and also cases where sibship must be tested. Three separate multiplex systems the Second Generation Multiplex, Powerplex 1.2 and FFFL have been employed, giving a total of 16 STR loci. Both empirical and theoretical approaches to the validation have been adopted. Appropriate equations have been derived to calculate likelihood ratios for different relationships, incorporating a correction for subpopulation effects. An F(ST) point estimate of 1% has been applied throughout. Empirically, 101 cases of alleged father, alleged mother and child where analysed using six SLP systems and also using the three multiplex STR systems. Of the 202 relationships tested, 197 were independently resolved by both systems, providing either clear evidence of non-parentage or strong support for the relationship.

Dynamic data reconciliation for sequential modular simulators: application to a mixing process

This paper describes a method for dynamic data reconciliation of nonlinear systems that are simulated using the sequential modular approach, and where individual modules are represented by a class of differential algebraic equations. The estimation technique consists of a bank of extended Kalman filters that are integrated with the modules. The paper reports a study based on experimental data obtained from a pilot scale mixing process.

Finite departure from convective quasi-equilibrium: periodic cycle and discharge-recharge mechanism

A simple self–contained theory is proposed for describing life cycles of convective systems as a discharge–recharge process. A closed description is derived for the dynamics of an ensemble of convective plumes based on an energy cycle. The system consists of prognostic equations for the cloud work function and the convective kinetic energy. The system can be closed by intro ducing a functional relationship between the convective kinetic energy and the cloud–base mass flux. The behaviour of this system is considered under a bulk simplification. Previous cloud–resolving mo delling as well as bulk statistical theories for ensemble convective systems suggest that a plausible relationship would be to assume that the convective kinetic energy is linearly proportional to the cloud–base mass flux. As a result, the system reduces to a nonlinear dynamical system with two dependent variables, the cloud–base mass flux and the cloud work function. The fully nonlinear solution of this system always represents a periodic cycle regardless of the initial condition under constant large–scale forcing. Importantly, the inclusion of energy dissipation in this model does not in itself lead the system to an equilibrium.

The Mathematics of Control Theory

This text contains papers presented at the Institute of Mathematics and its Applications Conference on Control Theory, held at the University of Strathclyde in Glasgow. The contributions cover a wide range of topics of current interest to theoreticians and practitioners including algebraic systems theory, nonlinear control systems, adaptive control, robustness issues, infinite dimensional systems, applications studies and connections to mathematical aspects of information theory and data-fusion.

Numerical Methods for Stiff Two-Point Boundary Value Problems

We consider the two-point boundary value problem for stiff systems of ordinary differential equations. For systems that can be transformed to essentially diagonally dominant form with appropriate smoothness conditions, a priori estimates are obtained. Problems with turning points can be treated with this theory, and we discuss this in detail. We give robust difference approximations and present error estimates for these schemes. In particular we give a detailed description of how to transform a general system to essentially diagonally dominant form and then stretch the independent variable so that the system will satisfy the correct smoothness conditions. Numerical examples are presented for both linear and nonlinear problems.

Natal nest locations of the Asian house rat (Rattus tanezumi) in lowland rice–coconut cropping systems: a coconut penthouse or rice bunds with water frontage?

CONTEXT. Rattus tanezumi is a serious crop pest within the island of Luzon, Philippines. In intensive flood-irrigated rice field ecosystems of Luzon, female R. tanezumi are known to primarily nest within the tillers of ripening rice fields and along the banks of irrigation canals. The nesting habits of R. tanezumi in complex rice–coconut cropping systems are unknown. AIMS. To identify the natal nest locations of R. tanezumi females in rice–coconut systems of the Sierra Madre Biodiversity Corridor (SMBC), Luzon, during the main breeding season to develop a management strategy that specifically targets their nesting habitat. METHODS. When rice was at the booting to ripening stage, cage-traps were placed in rice fields adjacent to coconut habitat. Thirty breeding adult R. tanezumi females were fitted with radio-collars and successfully tracked to their nest sites. KEY RESULTS. Most R. tanezumi nests (66.7%) were located in coconut groves, five nests (16.7%) were located in rice fields and five nests (16.7%) were located on the rice field edge. All nests were located above ground level and seven nests were located in coconut tree crowns. The median distance of nest sites to the nearest rice field was 22.5m. Most nest site locations had good cover of ground vegetation and understorey vegetation, but low canopy cover. Only one nest location had an understorey vegetation height of less than 20 cm. CONCLUSIONS. In the coastal lowland rice–coconut cropping systems of the SMBC, female R. tanezumi showed a preference for nesting in adjacent coconut groves. This is contrary to previous studies in intensive flood-irrigated rice ecosystems of Luzon, where the species nests mainly in the banks of irrigation canals. It is important to understand rodent breeding ecology in a specific ecosystem before implementing appropriate management strategies. IMPLICATIONS. In lowland rice–coconut cropping systems, coconut groves adjacent to rice fields should be targeted for the 20 management of R. tanezumi nest sites during the main breeding season as part of an integrated ecologically based approach to rodent pest management.

Correlations between two sets of angular relation equations

It is shown here that the angular relation equations between direct and reciprocal vectors are very similar to the angular relation equations in Euler's theorem. These two sets of equations are usually treated separately as unrelated equations in different fields. In this careful study, the connection between the two sets of angular equations is revealed by considering the cosine rule for the spherical triangle. It is found that understanding of the correlation is hindered by the facts that the same variables are defined differently and different symbols are used to represent them in the two fields. Understanding the connection between different concepts is not only stimulating and beneficial, but also a fundamental tool in innovation and research, and has historical significance. The background of the work presented here contains elements of many scientific disciplines. This work illustrates the common ground of two theories usually considered separately and is therefore of benefit not only for its own sake but also to illustrate a general principle that a theory relevant to one discipline can often be used in another. The paper works with chemistry related concepts using mathematical methodologies unfamiliar to the usual audience of mainstream experimental and theoretical chemists.

Dynamical complexity in a mean-field model of human EEG

A recently proposed mean-field theory of mammalian cortex rhythmogenesis describes the salient features of electrical activity in the cerebral macrocolumn, with the use of inhibitory and excitatory neuronal populations (Liley et al 2002). This model is capable of producing a range of important human EEG (electroencephalogram) features such as the alpha rhythm, the 40 Hz activity thought to be associated with conscious awareness (Bojak & Liley 2007) and the changes in EEG spectral power associated with general anesthetic effect (Bojak & Liley 2005). From the point of view of nonlinear dynamics, the model entails a vast parameter space within which multistability, pseudoperiodic regimes, various routes to chaos, fat fractals and rich bifurcation scenarios occur for physiologically relevant parameter values (van Veen & Liley 2006). The origin and the character of this complex behaviour, and its relevance for EEG activity will be illustrated. The existence of short-lived unstable brain states will also be discussed in terms of the available theoretical and experimental results. A perspective on future analysis will conclude the presentation.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat “patchy” connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a “lattice-directed” traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

A generalized collectively compact operator theory with an application to integral equations on unbounded domains

In this paper a generalization of collectively compact operator theory in Banach spaces is developed. A feature of the new theory is that the operators involved are no longer required to be compact in the norm topology. Instead it is required that the image of a bounded set under the operator family is sequentially compact in a weaker topology. As an application, the theory developed is used to establish solvability results for a class of systems of second kind integral equations on unbounded domains, this class including in particular systems of Wiener-Hopf integral equations with L1 convolutions kernels

Scattering of electromagnetic waves by rough interfaces and inhomogeneous layers

We consider a two-dimensional problem of scattering of a time-harmonic electromagnetic plane wave by an infinite inhomogeneous conducting or dielectric layer at the interface between semi-infinite homogeneous dielectric half-spaces. The magnetic permeability is assumed to be a fixed positive constant. The material properties of the media are characterized completely by an index of refraction, which is a bounded measurable function in the layer and takes positive constant values above and below the layer, corresponding to the homogeneous dielectric media. In this paper, we examine only the transverse magnetic (TM) polarization case. A radiation condition appropriate for scattering by infinite rough surfaces is introduced, a generalization of the Rayleigh expansion condition for diffraction gratings. With the help of the radiation condition the problem is reformulated as an equivalent mixed system of boundary and domain integral equations, consisting of second-kind integral equations over the layer and interfaces within the layer. Assumptions on the variation of the index of refraction in the layer are then imposed which prove to be sufficient, together with the radiation condition, to prove uniqueness of solution and nonexistence of guided wave modes. Recent, general results on the solvability of systems of second kind integral equations on unbounded domains establish existence of solution and continuous dependence in a weighted norm of the solution on the given data. The results obtained apply to the case of scattering by a rough interface between two dielectric media and to many other practical configurations.

Asymptotic behavior at infinity of solutions of multidimensional second kind integral equations

We consider second kind integral equations of the form x(s) - (abbreviated x - K x = y ), in which Ω is some unbounded subset of Rn. Let Xp denote the weighted space of functions x continuous on Ω and satisfying x (s) = O(|s|-p ),s → ∞We show that if the kernel k(s,t) decays like |s — t|-q as |s — t| → ∞ for some sufficiently large q (and some other mild conditions on k are satisfied), then K ∈ B(XP) (the set of bounded linear operators on Xp), for 0 ≤ p ≤ q. If also (I - K)-1 ∈ B(X0) then (I - K)-1 ∈ B(XP) for 0 < p < q, and (I- K)-1∈ B(Xq) if further conditions on k hold. Thus, if k(s, t) = O(|s — t|-q). |s — t| → ∞, and y(s)=O(|s|-p), s → ∞, the asymptotic behaviour of the solution x may be estimated as x (s) = O(|s|-r), |s| → ∞, r := min(p, q). The case when k(s,t) = к(s — t), so that the equation is of Wiener-Hopf type, receives especial attention. Conditions, in terms of the symbol of I — K, for I — K to be invertible or Fredholm on Xp are established for certain cases (Ω a half-space or cone). A boundary integral equation, which models three-dimensional acoustic propaga-tion above flat ground, absorbing apart from an infinite rigid strip, illustrates the practical application and sharpness of the above results. This integral equation mod-els, in particular, road traffic noise propagation along an infinite road surface sur-rounded by absorbing ground. We prove that the sound propagating along the rigid road surface eventually decays with distance at the same rate as sound propagating above the absorbing ground.