Dirac structures and boundary control systems associated with skew-symmetric differential operators

Associated with a skew-symmetric linear operator on the spatial domain $[a,b]$ we define a Dirac structure which includes the port variables on the boundary of this spatial domain. This Dirac structure is a subspace of a Hilbert space. Naturally, associated to this Dirac structure is infinite dimensional system. We parameterize the boundary port variables for which the $C_{0}$-semigroup associated to this system is contractive or unitary. Furthermore, this parameterization is used to split the boundary port variables into inputs and outputs. Similarly, we define a linear port controlled Hamiltonian system associated with the previously defined Dirac structure and a symmetric positive operator defining the energy of the system. We illustrate this theory on the example of the Timoshenko Beam. \u

Rice Intensification in a Changing Environment: Impact on Water Availability in Inland Valley Landscapes in Benin

This study assesses the impact of climate change on hydrological processes under rice intensification in three headwater inland valley watersheds characterized by different land conditions. The Soil and Water Assessment Tool was used to simulate the combined impacts of two land use scenarios defined as converting 25% and 75% of lowland savannah into rice cultivation, and two climate scenarios (A1B and B1) of the Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios. The simulations were performed based on the traditional and the rainfed-bunded rice cultivation systems and analyzed up to the year 2049 with a special focus on the period of 2030–2049. Compared to land use, climate change impact on hydrological processes was overwhelming at all watersheds. The watersheds with a high portion of cultivated areas are more sensitive to changes in climate resulting in a decrease of water yield of up to 50% (145 mm). Bunded fields cause a rise in surface runoff projected to be up to 28% (18 mm) in their lowlands, while processes were insignificantly affected at the vegetation dominated-watershed. Analyzing three watersheds instead of one as is usually done provides further insight into the natural variability and therefore gives more evidence of possible future processes and management strategie

A framework for evaluating statistical dependencies and rank correlations in power law graphs

We analyze dependencies in power law graph data (Web sample, Wikipedia sample and a preferential attachment graph) using statistical inference for multivariate regular variation. To the best of our knowledge, this is the first attempt to apply the well developed theory of regular variation to graph data. The new insights this yields are striking: the three above-mentioned data sets are shown to have a totally different dependence structure between different graph parameters, such as in-degree and PageRank. Based on the proposed methodology, we suggest a new measure for rank correlations. Unlike most known methods, this measure is especially sensitive to rank permutations for topranked nodes. Using this method, we demonstrate that the PageRank ranking is not sensitive to moderate changes in the damping factor

On the accuracy of phase-type approximations of heavy-tailed risk models

Numerical evaluation of ruin probabilities in the classical risk model is an important problem. If claim sizes are heavy-tailed, then such evaluations are challenging. To overcome this, an attractive way is to approximate the claim sizes with a phase-type distribution. What is not clear though is how many phases are enough in order to achieve a specific accuracy in the approximation of the ruin probability. The goals of this paper are to investigate the number of phases required so that we can achieve a pre-specified accuracy for the ruin probability and to provide error bounds. Also, in the special case of a completely monotone claim size distribution we develop an algorithm to estimate the ruin probability by approximating the excess claim size distribution with a hyperexponential one. Finally, we compare our approximation with the heavy traffic and heavy tail approximations.Comment: 24 pages, 13 figures, 8 tables, 38 reference

Parallel Algorithm for Solving Kepler's Equation on Graphics Processing Units: Application to Analysis of Doppler Exoplanet Searches

[Abridged] We present the results of a highly parallel Kepler equation solver using the Graphics Processing Unit (GPU) on a commercial nVidia GeForce 280GTX and the "Compute Unified Device Architecture" programming environment. We apply this to evaluate a goodness-of-fit statistic (e.g., chi^2) for Doppler observations of stars potentially harboring multiple planetary companions (assuming negligible planet-planet interactions). We tested multiple implementations using single precision, double precision, pairs of single precision, and mixed precision arithmetic. We find that the vast majority of computations can be performed using single precision arithmetic, with selective use of compensated summation for increased precision. However, standard single precision is not adequate for calculating the mean anomaly from the time of observation and orbital period when evaluating the goodness-of-fit for real planetary systems and observational data sets. Using all double precision, our GPU code outperforms a similar code using a modern CPU by a factor of over 60. Using mixed-precision, our GPU code provides a speed-up factor of over 600, when evaluating N_sys > 1024 models planetary systems each containing N_pl = 4 planets and assuming N_obs = 256 observations of each system. We conclude that modern GPUs also offer a powerful tool for repeatedly evaluating Kepler's equation and a goodness-of-fit statistic for orbital models when presented with a large parameter space.Comment: 19 pages, to appear in New Astronom

Core Formation by a Population of Massive Remnants

Core radii of globular clusters in the Large and Small Magellanic Clouds show an increasing trend with age. We propose that this trend is a dynamical effect resulting from the accumulation of massive stars and stellar-mass black holes at the cluster centers. The black holes are remnants of stars with initial masses exceeding 20-25 solar masses; as their orbits decay by dynamical friction, they heat the stellar background and create a core. Using analytical estimates and N-body experiments, we show that the sizes of the cores so produced and their growth rates are consistent with what is observed. We propose that this mechanism is responsible for the formation of cores in all globular clusters and possibly in other systems as well.Comment: 5 page

Corrected phase-type approximations of heavy-tailed queueing models in a Markovian environment

Significant correlations between arrivals of load-generating events make the numerical evaluation of the workload of a system a challenging problem. In this paper, we construct highly accurate approximations of the workload distribution of the MAP/G/1 queue that capture the tail behavior of the exact workload distribution and provide a bounded relative error. Motivated by statistical analysis, we consider the service times as a mixture of a phase-type and a heavy-tailed distribution. With the aid of perturbation analysis, we derive our approximations as a sum of the workload distribution of the MAP/PH/1 queue and a heavy-tailed component that depends on the perturbation parameter. We refer to our approximations as corrected phase-type approximations, and we exhibit their performance with a numerical study.Comment: Received the Marcel Neuts Student Paper Award at the 8th International Conference on Matrix Analytic Methods in Stochastic Models 201