Riemannian metric and geometric mean for positive semidefinite matrices of fixed rank

This paper introduces a new metric and mean on the set of positive semidefinite matrices of fixed-rank. The proposed metric is derived from a well-chosen Riemannian quotient geometry that generalizes the reductive geometry of the positive cone and the associated natural metric. The resulting Riemannian space has strong geometrical properties: it is geodesically complete, and the metric is invariant with respect to all transformations that preserve angles (orthogonal transformations, scalings, and pseudoinversion). A meaningful approximation of the associated Riemannian distance is proposed, that can be efficiently numerically computed via a simple algorithm based on SVD. The induced mean preserves the rank, possesses the most desirable characteristics of a geometric mean, and is easy to compute. © 2009 Society for Industrial and Applied Mathematics.

Weighted geometric mean of large-scale matrices: numerical analysis and algorithms

Computing the weighted geometric mean of large sparse matrices is an operation that tends to become rapidly intractable, when the size of the matrices involved grows. However, if we are not interested in the computation of the matrix function itself, but just in that of its product times a vector, the problem turns simpler and there is a chance to solve it even when the matrix mean would actually be impossible to compute. Our interest is motivated by the fact that this calculation has some practical applications, related to the preconditioning of some operators arising in domain decomposition of elliptic problems. In this thesis, we explore how such a computation can be efficiently performed. First, we exploit the properties of the weighted geometric mean and find several equivalent ways to express it through real powers of a matrix. Hence, we focus our attention on matrix powers and examine how well-known techniques can be adapted to the solution of the problem at hand. In particular, we consider two broad families of approaches for the computation of f(A) v, namely quadrature formulae and Krylov subspace methods, and generalize them to the pencil case f(A\B) v. Finally, we provide an extensive experimental evaluation of the proposed algorithms and also try to assess how convergence speed and execution time are influenced by some characteristics of the input matrices. Our results suggest that a few elements have some bearing on the performance and that, although there is no best choice in general, knowing the conditioning and the sparsity of the arguments beforehand can considerably help in choosing the best strategy to tackle the problem.

Quadratic Mean Radius of a Polynomial in C(Z)

* Dedicated to the memory of Prof. N. Obreshkoff

A Model for Bundled Conductors considering a Non-Homogeneous Distribution of the Current through Subconductors

This paper describes an alternative procedure to obtain an equivalent conductor from a bundled conductor, taking into account the distribution of the current in subcondutors of the bundle. Firstly, it is introduced a brief background about the concept of Geometric Mean Radius (GMR) and how this methodology is applied to define an equivalent conductor and its electric parameters. Emphasizing that the classical procedure, using GMR, is limited to premise which the current is equally distributed through subconductors. Afterwards, it is described the development of proposed method and applications for an equivalent conductor obtained from a conventional transmission line bundled conductor and from an equivalent conductor based on a bundle with compressed SF(6) insulation system, where the current is unequally distributed through subconductors.

Características dos parâmetros do condutor equivalente a um feixe de subcondutores de linhas de transmissão: análise inicial

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Um modelo para condutores múltiplos considerando a distribuição da corrente nos subcondutores

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Rank-preserving geometric means of positive semi-definite matrices

The generalization of the geometric mean of positive scalars to positive definite matrices has attracted considerable attention since the seminal work of Ando. The paper generalizes this framework of matrix means by proposing the definition of a rank-preserving mean for two or an arbitrary number of positive semi-definite matrices of fixed rank. The proposed mean is shown to be geometric in that it satisfies all the expected properties of a rank-preserving geometric mean. The work is motivated by operations on low-rank approximations of positive definite matrices in high-dimensional spaces.© 2012 Elsevier Inc. All rights reserved.

On the development of summer thunderstorms in the city of Sao Paulo: Mean meteorological characteristics and pollution effect

This study investigates how the summer thunderstorms developed over the city of Sao Paulo and if the pollution might affect its development or characteristics during the austral summer (December-January-February-March, DJFM months). A total of 605 days from December 1999 to March 2004 was separated as 241 thunderstorms days (TDs) and 364 non-thunderstorm days (NTDs). The analyses are performed by using hourly measurements of air temperature (T), web-bulb temperature (Tw), surface atmospheric pressure (P), wind velocity and direction, rainfall and thunder and lightning observations collected at the Meteorological Station of the University of Sao Paulo in conjunction with aerosol measurements obtained by AERONET (Aerosol Robotic Network), and the NCEP-DOE (National Centers for Environmental Prediction Department of Energy) reanalysis and radiosondes. The wind diurnal cycle shows that for TDs the morning flow is from the northwest rotating to the southeast after 16: 00 local time (LT) and it remains from the east until the night. For the NTDs, the wind is well characterized by the sea-breeze circulation that in the morning has the wind blowing from the northeast and in the afternoon from the southeast. The TDs show that the air temperature diurnal cycle presents higher amplitude and the maximum temperature of the day is 3.2 degrees C higher than in NTDs. Another important factor found is the difference between moisture that is higher during TDs. In terms of precipitation, the TDs represent 40% of total of days analyzed and those days are responsible for more than 60% of the total rain accumulation during the summer, for instance 50% of the TDs had more than 15.5mm day(-1) while the NTDs had 4 mm day(-1). Moreover, the rainfall distribution shows that TDs have higher rainfall rate intensities and an afternoon precipitation maximum; while in the NTDs there isn`t a defined precipitation diurnal cycle. The wind and temperature fields from NCEP reanalysis concur with the local weather station and radiosonde observations. The NCEP composites show that TDs are controlled by synoptic circulation characterized by a pre-frontal situation, with a baroclinic zone situated at southern part of Sao Paulo. In terms of pollution, this study employed the AERONET data to obtain the main aerosol characteristics in the atmospheric column for both TDs and NTDs. The particle size distribution and particle volume size distribution have similar concentrations for both TDs and NTDs and present a similar fine and coarse mode mean radius. In respect to the atmospheric loading, the aerosol optical depth (AOD) at different frequencies presented closed mean values for both TDs and NTDs that were statistically significant at 95% level. The spectral dependency of those values in conjunction with the Angstrom parameter reveal the higher concentration of the fine mode particles that are more likely to be hygroscopic and from urban areas. In summary, no significant aerosol effect could be found on the development of summer thunderstorms, suggesting the strong synoptic control by the baroclinic forcing for deep convective development. (C) 2010 Published by Elsevier B. V.

Selectivity in separation using π electron-rich stationary phases for the comprehensive two-dimensional analysis of cafe espresso

The multidimensional high-performance liquid chromatography separations of the complex sample matrix found in café espresso coffee were completed on the propyl phenyl and butyl phenyl columns that contain 3 and 4 carbon atoms in the spacer chain, respectively. Phenyl type stationary phases are able to undergo unique π–π interactions with aromatic compounds. Previous works have found that there are differences in retention characteristics between these chain lengths and this was explored further here. It was found that when analysing the separations by quadrants, using a geometric approach to factor analysis and by measuring the normalised mean radius, subtle differences in the separations were observed and the butyl phenyl phase was more selective for the high to medium polarity species. However, there was very little difference in separation behaviour for the hydrophobic components within the coffee sample. Overall, the analysis of the entire separation showed very little difference.

Characterization of expiration air jets and droplet size distributions immediately at the mouth opening

Size distributions of expiratory droplets expelled during coughing and speaking and the velocities of the expiration air jets of healthy volunteers were measured. Droplet size was measured using the Interferometric Mie imaging (IMI) technique while the Particle Image Velocimetry (PIV) technique was used for measuring air velocity. These techniques allowed measurements in close proximity to the mouth and avoided air sampling losses. The average expiration air velocity was 11.7 m/s for coughing and 3.9 m/s for speaking. Under the experimental setting, evaporation and condensation effects had negligible impact on the measured droplet size. The geometric mean diameter of droplets from coughing was 13.5m and it was 16.0m for speaking (counting 1 to 100). The estimated total number of droplets expelled ranged from 947 – 2085 per cough and 112 – 6720 for speaking. The estimated droplet concentrations for coughing ranged from 2.4 - 5.2cm-3 per cough and 0.004 – 0.223 cm-3 for speaking.

Radiological aspects of petroleum exploration and production in the sultanate of Oman

This thesis is a study of naturally occurring radioactive materials (NORM) activity concentration, gamma dose rate and radon (222Rn) exhalation from the waste streams of large-scale onshore petroleum operations. Types of activities covered included; sludge recovery from separation tanks, sludge farming, NORM storage, scaling in oil tubulars, scaling in gas production and sedimentation in produced water evaporation ponds. Field work was conducted in the arid desert terrain of an operational oil exploration and production region in the Sultanate of Oman. The main radionuclides found were 226Ra and 210Pb (238U - series), 228Ra and 228Th (232Th - series), and 227Ac (235U - series), along with 40K. All activity concentrations were higher than the ambient soil level and varied over several orders of magnitude. The range of gamma dose rates at a 1 m height above ground for the farm treated sludge had a range of 0.06 0.43 µSv h 1, and an average close to the ambient soil mean of 0.086 ± 0.014 µSv h 1, whereas the untreated sludge gamma dose rates had a range of 0.07 1.78 µSv h 1, and a mean of 0.456 ± 0.303 µSv h 1. The geometric mean of ambient soil 222Rn exhalation rate for area surrounding the sludge was mBq m 2 s 1. Radon exhalation rates reported in oil waste products were all higher than the ambient soil value and varied over three orders of magnitude. This study resulted in some unique findings including: (i) detection of radiotoxic 227Ac in the oil scales and sludge, (ii) need of a new empirical relation between petroleum sludge activity concentrations and gamma dose rates, and (iii) assessment of exhalation of 222Rn from oil sludge. Additionally the study investigated a method to determine oil scale and sludge age by the use of inherent behaviour of radionuclides as 228Ra:226Ra and 228Th:228Ra activity ratios.

The impact of delayed blood centrifuging, choice of collection tube, and type of assay on 25-hydroxyvitamin D concentrations

Studies have examined the associations between cancers and circulating 25-hydroxyvitamin D [25(OH)D], but little is known about the impact of different laboratory practices on 25(OH)D concentrations. We examined the potential impact of delayed blood centrifuging, choice of collection tube, and type of assay on 25(OH)D concentrations. Blood samples from 20 healthy volunteers underwent alternative laboratory procedures: four centrifuging times (2, 24, 72, and 96 h after blood draw); three types of collection tubes (red top serum tube, two different plasma anticoagulant tubes containing heparin or EDTA); and two types of assays (DiaSorin radioimmunoassay [RIA] and chemiluminescence immunoassay [CLIA/LIAISON®]). Log-transformed 25(OH)D concentrations were analyzed using the generalized estimating equations (GEE) linear regression models. We found no difference in 25(OH)D concentrations by centrifuging times or type of assay. There was some indication of a difference in 25(OH)D concentrations by tube type in CLIA/LIAISON®-assayed samples, with concentrations in heparinized plasma (geometric mean, 16.1 ng ml−1) higher than those in serum (geometric mean, 15.3 ng ml−1) (p = 0.01), but the difference was significant only after substantial centrifuging delays (96 h). Our study suggests no necessity for requiring immediate processing of blood samples after collection or for the choice of a tube type or assay.