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.

Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20\% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

A Comparative Study of Two Semantic Systems: Body Russian and Russian Phraseology

Grigorij Kreidlin (Russia). A Comparative Study of Two Semantic Systems: Body Russian and Russian Phraseology. Mr. Kreidlin teaches in the Department of Theoretical and Applied Linguistics of the State University of Humanities in Moscow and worked on this project from August 1996 to July 1998. The classical approach to non-verbal and verbal oral communication is based on a traditional separation of body and mind. Linguists studied words and phrasemes, the products of mind activities, while gestures, facial expressions, postures and other forms of body language were left to anthropologists, psychologists, physiologists, and indeed to anyone but linguists. Only recently have linguists begun to turn their attention to gestures and semiotic and cognitive paradigms are now appearing that raise the question of designing an integral model for the unified description of non-verbal and verbal communicative behaviour. This project attempted to elaborate lexical and semantic fragments of such a model, producing a co-ordinated semantic description of the main Russian gestures (including gestures proper, postures and facial expressions) and their natural language analogues. The concept of emblematic gestures and gestural phrasemes and of their semantic links permitted an appropriate description of the transformation of a body as a purely physical substance into a body as a carrier of essential attributes of Russian culture - the semiotic process called the culturalisation of the human body. Here the human body embodies a system of cultural values and displays them in a text within the area of phraseology and some other important language domains. The goal of this research was to develop a theory that would account for the fundamental peculiarities of the process. The model proposed is based on the unified lexicographic representation of verbal and non-verbal units in the Dictionary of Russian Gestures, which the Mr. Kreidlin had earlier complied in collaboration with a group of his students. The Dictionary was originally oriented only towards reflecting how the lexical competence of Russian body language is represented in the Russian mind. Now a special type of phraseological zone has been designed to reflect explicitly semantic relationships between the gestures in the entries and phrasemes and to provide the necessary information for a detailed description of these. All the definitions, rules of usage and the established correlations are written in a semantic meta-language. Several classes of Russian gestural phrasemes were identified, including those phrasemes and idioms with semantic definitions close to those of the corresponding gestures, those phraseological units that have lost touch with the related gestures (although etymologically they are derived from gestures that have gone out of use), and phrasemes and idioms which have semantic traces or reflexes inherited from the meaning of the related gestures. The basic assumptions and practical considerations underlying the work were as follows. (1) To compare meanings one has to be able to state them. To state the meaning of a gesture or a phraseological expression, one needs a formal semantic meta-language of propositional character that represents the cognitive and mental aspects of the codes. (2) The semantic contrastive analysis of any semiotic codes used in person-to-person communication also requires a single semantic meta-language, i.e. a formal semantic language of description,. This language must be as linguistically and culturally independent as possible and yet must be open to interpretation through any culture and code. Another possible method of conducting comparative verbal-non-verbal semantic research is to work with different semantic meta-languages and semantic nets and to learn how to combine them, translate from one to another, etc. in order to reach a common basis for the subsequent comparison of units. (3) The practical work in defining phraseological units and organising the phraseological zone in the Dictionary of Russian Gestures unexpectedly showed that semantic links between gestures and gestural phrasemes are reflected not only in common semantic elements and syntactic structure of semantic propositions, but also in general and partial cognitive operations that are made over semantic definitions. (4) In comparative semantic analysis one should take into account different values and roles of inner form and image components in the semantic representation of non-verbal and verbal units. (5) For the most part, gestural phrasemes are direct semantic derivatives of gestures. The cognitive and formal techniques can be regarded as typological features for the future functional-semantic classification of gestural phrasemes: two phrasemes whose meaning can be obtained by the same cognitive or purely syntactic operations (or types of operations) over the meanings of the corresponding gestures, belong by definition to one and the same class. The nature of many cognitive operations has not been studied well so far, but the first steps towards its comprehension and description have been taken. The research identified 25 logically possible classes of relationships between a gesture and a gestural phraseme. The calculation is based on theoretically possible formal (set-theory) correlations between signifiers and signified of the non-verbal and verbal units. However, in order to examine which of them are realised in practice a complete semantic and lexicographic description of all (not only central) everyday emblems and gestural phrasemes is required and this unfortunately does not yet exist. Mr. Kreidlin suggests that the results of the comparative analysis of verbal and non-verbal units could also be used in other research areas such as the lexicography of emotions.

Black hair follicular dysplasia in Large Münsterländer dogs: clinical, histological and ultrastructural features

Four Large Münsterländer cross-bred dogs affected with black hair follicular dysplasia (BHFD) and one unaffected control littermate were observed, and skin was sampled weekly over the first 19 weeks of life. Affected dogs were born with silvery grey hair, a consequence of melanin clumping in the hair shafts. Hair bulb melanocytes were densely pigmented, and contained abundant stage IV melanosomes but adjacent matrix keratinocytes lacked melanosomes. Melanin clumping was not prominent in epidermal melanocytes in the haired skin but occurred in the foot pads. Follicular changes progressed from bulbar clumping, clumping in the isthmus/infundibulum and finally to dysplastic hair shafts. Alopecia developed progressively in pigmented areas. Silver-grey hair, melanin clumping, accumulation of stage IV melanosomes within melanocytes and insufficient melanin transfer to adjacent keratinocytes are also classic features of human Griscelli syndrome. The underlying cause in Griscelli syndrome is a defect of melanocytic intracellular transport proteins leading to inadequate and disorganized melanosome transfer to keratinocytes with resultant melanin clumping. In view of the correlation in the phenotype, histology and ultrastructure between both disorders, a defect in intracellular melanosome transport is postulated as the pathogenic mechanism in BHFD.

Differences across health care systems in outcome and cost-utility of surgical and conservative treatment of chronic low back pain: a study protocol

BACKGROUND: There is little evidence on differences across health care systems in choice and outcome of the treatment of chronic low back pain (CLBP) with spinal surgery and conservative treatment as the main options. At least six randomised controlled trials comparing these two options have been performed; they show conflicting results without clear-cut evidence for superior effectiveness of any of the evaluated interventions and could not address whether treatment effect varied across patient subgroups. Cost-utility analyses display inconsistent results when comparing surgical and conservative treatment of CLBP. Due to its higher feasibility, we chose to conduct a prospective observational cohort study. METHODS: This study aims to examine if1. Differences across health care systems result in different treatment outcomes of surgical and conservative treatment of CLBP2. Patient characteristics (work-related, psychological factors, etc.) and co-interventions (physiotherapy, cognitive behavioural therapy, return-to-work programs, etc.) modify the outcome of treatment for CLBP3. Cost-utility in terms of quality-adjusted life years differs between surgical and conservative treatment of CLBP.This study will recruit 1000 patients from orthopaedic spine units, rehabilitation centres, and pain clinics in Switzerland and New Zealand. Effectiveness will be measured by the Oswestry Disability Index (ODI) at baseline and after six months. The change in ODI will be the primary endpoint of this study.Multiple linear regression models will be used, with the change in ODI from baseline to six months as the dependent variable and the type of health care system, type of treatment, patient characteristics, and co-interventions as independent variables. Interactions will be incorporated between type of treatment and different co-interventions and patient characteristics. Cost-utility will be measured with an index based on EQol-5D in combination with cost data. CONCLUSION: This study will provide evidence if differences across health care systems in the outcome of treatment of CLBP exist. It will classify patients with CLBP into different clinical subgroups and help to identify specific target groups who might benefit from specific surgical or conservative interventions. Furthermore, cost-utility differences will be identified for different groups of patients with CLBP. Main results of this study should be replicated in future studies on CLBP.

Electronic medical record systems, data quality and loss to follow-up: survey of antiretroviral therapy programmes in resource-limited settings

OBJECTIVE: To describe the electronic medical databases used in antiretroviral therapy (ART) programmes in lower-income countries and assess the measures such programmes employ to maintain and improve data quality and reduce the loss of patients to follow-up. METHODS: In 15 countries of Africa, South America and Asia, a survey was conducted from December 2006 to February 2007 on the use of electronic medical record systems in ART programmes. Patients enrolled in the sites at the time of the survey but not seen during the previous 12 months were considered lost to follow-up. The quality of the data was assessed by computing the percentage of missing key variables (age, sex, clinical stage of HIV infection, CD4+ lymphocyte count and year of ART initiation). Associations between site characteristics (such as number of staff members dedicated to data management), measures to reduce loss to follow-up (such as the presence of staff dedicated to tracing patients) and data quality and loss to follow-up were analysed using multivariate logit models. FINDINGS: Twenty-one sites that together provided ART to 50 060 patients were included (median number of patients per site: 1000; interquartile range, IQR: 72-19 320). Eighteen sites (86%) used an electronic database for medical record-keeping; 15 (83%) such sites relied on software intended for personal or small business use. The median percentage of missing data for key variables per site was 10.9% (IQR: 2.0-18.9%) and declined with training in data management (odds ratio, OR: 0.58; 95% confidence interval, CI: 0.37-0.90) and weekly hours spent by a clerk on the database per 100 patients on ART (OR: 0.95; 95% CI: 0.90-0.99). About 10 weekly hours per 100 patients on ART were required to reduce missing data for key variables to below 10%. The median percentage of patients lost to follow-up 1 year after starting ART was 8.5% (IQR: 4.2-19.7%). Strategies to reduce loss to follow-up included outreach teams, community-based organizations and checking death registry data. Implementation of all three strategies substantially reduced losses to follow-up (OR: 0.17; 95% CI: 0.15-0.20). CONCLUSION: The quality of the data collected and the retention of patients in ART treatment programmes are unsatisfactory for many sites involved in the scale-up of ART in resource-limited settings, mainly because of insufficient staff trained to manage data and trace patients lost to follow-up.

Computational simulations of latent heat thermal energy storage systems - with innovative and first-principles based simulation for the underlying unsteady melting (and solidification) processes

This thesis develops an effective modeling and simulation procedure for a specific thermal energy storage system commonly used and recommended for various applications (such as an auxiliary energy storage system for solar heating based Rankine cycle power plant). This thermal energy storage system transfers heat from a hot fluid (termed as heat transfer fluid - HTF) flowing in a tube to the surrounding phase change material (PCM). Through unsteady melting or freezing process, the PCM absorbs or releases thermal energy in the form of latent heat. Both scientific and engineering information is obtained by the proposed first-principle based modeling and simulation procedure. On the scientific side, the approach accurately tracks the moving melt-front (modeled as a sharp liquid-solid interface) and provides all necessary information about the time-varying heat-flow rates, temperature profiles, stored thermal energy, etc. On the engineering side, the proposed approach is unique in its ability to accurately solve – both individually and collectively – all the conjugate unsteady heat transfer problems for each of the components of the thermal storage system. This yields critical system level information on the various time-varying effectiveness and efficiency parameters for the thermal storage system.

The Influence of Enterprise Systems on Business and Information Technology

Business strategy is important to all organizations. Nearly all Fortune 500 firms are implementing Enterprise Resource Planning (ERP) systems to improve the execution of their business strategy and to improve integration with its information technology (IT) strategy. Successful implementation of these multi-million dollar software systems are requiring new emphasis on change management and on Business and IT strategic alignment. This paper examines business and IT strategic alignment and seeks to explore whether an ERP implementation can drive business process reengineering and business and IT strategic alignment. An overview of business strategy and strategic alignment are followed by an analysis of ERP. The “As-Is/To-Be” process model is then presented and explained as a simple, but vital tool for improving business strategy, strategic alignment, and ERP implementation success.

Over two decades of blended and online physics courses at Michigan State University

In Fall 1992, our first physics course offered online homework. Over two decades later, we have seven physics courses online, spanning the whole range of introductory course offerings, with a total of over 1600 students in 2014. We found that several of the the purely online courses had better learning success than traditional lecture courses, as measured by exam scores. Particularly successful were online materials with embedded assessment. This result can be interpreted in different ways, but may serve as an indicator that during in-class lectures, we are oftentimes not taking advantage of the fact that we have the students on-site.

Evaluation of Year-round Forage Management Systems for Spring- and Fall-Calving Beef Cows (A Progress Report)

A year-round grazing system for spring- and fall-calving cows was developed to compare animal production and performance, hay production and feeding, winter forage composition changes, and summer pasture yield and nutrient composition to that from a conventional, or minimal land system. Systems compared forage from smooth bromegrass-orchardgrass-birdsfoot trefoil pastures for both systems in the summer and corn crop residues and stockpiled grass-legume pastures for the year-round system to drylot hay feeding during winter for the minimal land system. The year-round grazing system utilized 1.67 acres of smooth bromegrassorchardgrass- birdsfoot trefoil (SB-O-T) pasture per cow in the summer, compared with 3.33 acres of (SB-O-T) pasture per cow in the control (minimal land) system. In addition to SB-O-T pastures, the year-round grazing system utilized 2.5 acres of tall fescue-red clover (TFRC) and 2.5 acres of smooth bromegrass-red clover (SBRC) per cow for grazing in both mid-summer and winter for fall- and spring-calving cows, respectively. First-cutting hay was harvested from the TF-RC and SB-RC pastures, and regrowth was grazed for approximately 45 days in the summer. These pastures were then fertilized with 40 lbs N/acre and stockpiled for winter grazing. Also utilized during the winter for spring-calving cows in the year-round grazing system were corn crop residue (CCR) pastures at an allowance of 2.5 acres per cow. In the minimal land system, hay was harvested from three-fourths of the area in SB-O-T pastures and stored for feeding in a drylot through the winter. Summer grazing was managed with rotational stocking for both systems, and winter grazing of stockpiled forages and corn crop residues by year-round system cows was managed by strip-stocking. Hay was fed to maintain a body condition score of 5 on a 9 point scale for spring-calving cows in both systems. Hay was supplemented as needed to maintain a body condition score of 3 for fall-calving cows nursing calves through the winter. Although initial condition scores for cows in both systems were different at the initiation of grazing for both winter and summer, there were no significant differences (P > .05) in overall condition score changes throughout both grazing seasons. In year 1, fall-calving cows in the year-round grazing system lost more (P < .05) body weight during winter than spring-calving cows in either system. In year 2, there were no differences seen in weight changes over winter for any group of cows. Average daily gains of fall calves in the yearround system were 1.9 lbs/day compared with weight gains of 2.5 lbs/day for spring calves from both systems. Yearly growing animal production from pastures for both years did not differ between systems when weight gains of stockers that grazed summer pastures in the year-round grazing system were added to weight gains of suckling calves. Carcass characteristics for all calves finished in the feedlot for both systems were similar. There were no significant differences in hay production between systems for year 1; however, amounts of hay needed to maintain cows were 923, 1373, 4732 lbs dry matter/cow for year-round fall-calving, year-round spring-calving, and minimal land spring-calving cows, respectively. In year 2, hay production per acre in the minimal land system was greater (P < .05) than for the year-round system, but the amounts of hay required per cow were 0, 0, and 4720 lbs dry matter/cow for yearround fall-calving, year-round spring-calving, and minimal land spring-calving cows, respectively.

Water–montmorillonite systems: Neutron scattering and tracer throughdiffusion studies

Designs for deep geological respositories of nuclear waste include bentonite as a hydraulic and chemisorption buffer material to protect the biosphere from leakage of radionuclides. Bentonite is chosen because it is a cheap, naturally occurring material with the required properties. It consists essentially of montmorillonite, a swelling clay mineral. Upon contact with groundwater such clays can seal the repository by incorporating water in the interlayers of their crystalline structure. The intercalated water exhibits significantly different properties to bulk water in the surrounding interparticle pores, such as lower diffusion coefficients (González Sánchez et. al. 2008). This doctoral thesis presents water distribution and diffusion behavior on various time and space scales in montmorillonite. Experimental results are presented for Na- and Cs-montmorillonite samples with a range of bulk dry densities (0.8 to 1.7 g/cm3). The experimental methods employed were neutron scattering (backscattering, diffraction, time-of-flight), adsorption measurements (water, nitrogen) and tracer-through diffusion. For the tracer experiments the samples were fully saturated via the liquid phase under volume-constrained conditions. In contrast, for the neutron scattering experiments, the samples were hydrated via the vapor phase and subsequently compacted, leaving a significant fraction of interparticle pores unfilled with water. Owing to these differences in saturation, the water contents of the samples for neutron scattering were characterized by gravimetry whereas those for the tracer experiments were obtained from the bulk dry density. The amount of surface water in interlayer pores could be successfully discriminated from the amount of bulk-like water in interparticle pores in Na- and Csmontmorillonite using neutron spectroscopy. For the first time in the literature, the distribution of water between these two pore environments was deciphered as a function of gravimetric water content. The amount was compared to a geometrical estimation of the amount of interlayer and interparticle water determined by neutron diffraction and adsorption measurements. The relative abundances of the 1 to 4 molecular water layers in the interlayer were determined from the area ratios of the (001)-diffraction peaks. Depending on the characterization method, different fractions of surface water and interlayer water were obtained. Only surface and interlayer water exists in amontmorillonite with water contents up to 0.18 g/g according to spectroscopic measurements and up to 0.32 g/g according to geometrical estimations, respectively. At higher water contents, bulk-like and interparticle water also exists. The amounts increase monotonically, but not linearly, from zero to 0.33 g/g for bulk-like water and to 0.43 g/g for interparticle water. It was found that water most likely redistributes between the surface and interlayer sites during the spectroscopic measurements and therefore the reported fraction is relevant only below about -10 ºC (Anderson, 1967). The redistribution effect can explain the discrepancy in fractions between the methods. In a novel approach the fractions of water in different pore environments were treated as a fixed parameter to derive local diffusion coefficients for water from quasielastic neutron scattering data, in particular for samples with high water contents. Local diffusion coefficients were obtained for the 1 to 4 molecular water layers in the interlayer of 0.5·10–9, 0.9·10–9, 1.5·10–9 and 1.4·10–9 m²/s, respectively, taking account of the different water fractions (molecular water layer, bulk-like water). The diffusive transport of 22Na and HTO through Na-montmorillonite was measured on the laboratory experimental scale (i.e. cm, days) by tracer through-diffusion experiments. We confirmed that diffusion of HTO is independent of the ionic strength of the external solution in contact with the clay sample but dependent on the bulk dry density. In contrast, the diffusion of 22Na was found to depend on both the ionic strength of the pore solution and on the bulk dry density. The ratio of the pore and surface diffusion could be experimentally determined for 22Na from the dependence of the diffusion coefficient on the ionic strength. Activation energies were derived from the temperaturedependent diffusion coefficients via the Arrhenius relation. In samples with high bulk dry density the activation energies are slightly higher than those of bulk water whereas in low density samples they are lower. The activation energies as a function of ionic strengths of the pore solutions are similar for 22Na and HTO. The facts that (i) the slope of the logarithmic effective diffusion coefficients as a function of the logarithmic ionic strength is less than unity for low bulk dry densities and (ii) two water populations can be observed for high gravimetric water contents (low bulk dry densities) support the interlayer and interparticle porosity model proposed by Glaus et al. (2007), Bourg et al. (2006, 2007) and Gimmi and Kosakowski (2011).