Multicore Scheduling of Real-Time Irregular Parallel Algorithms in Linux

Face à estagnação da tecnologia uniprocessador registada na passada década, aos principais fabricantes de microprocessadores encontraram na tecnologia multi-core a resposta `as crescentes necessidades de processamento do mercado. Durante anos, os desenvolvedores de software viram as suas aplicações acompanhar os ganhos de performance conferidos por cada nova geração de processadores sequenciais, mas `a medida que a capacidade de processamento escala em função do número de processadores, a computação sequencial tem de ser decomposta em várias partes concorrentes que possam executar em paralelo, para que possam utilizar as unidades de processamento adicionais e completar mais rapidamente. A programação paralela implica um paradigma completamente distinto da programação sequencial. Ao contrário dos computadores sequenciais tipificados no modelo de Von Neumann, a heterogeneidade de arquiteturas paralelas requer modelos de programação paralela que abstraiam os programadores dos detalhes da arquitectura e simplifiquem o desenvolvimento de aplicações concorrentes. Os modelos de programação paralela mais populares incitam os programadores a identificar instruções concorrentes na sua lógica de programação, e a especificá-las sob a forma de tarefas que possam ser atribuídas a processadores distintos para executarem em simultâneo. Estas tarefas são tipicamente lançadas durante a execução, e atribuídas aos processadores pelo motor de execução subjacente. Como os requisitos de processamento costumam ser variáveis, e não são conhecidos a priori, o mapeamento de tarefas para processadores tem de ser determinado dinamicamente, em resposta a alterações imprevisíveis dos requisitos de execução. `A medida que o volume da computação cresce, torna-se cada vez menos viável garantir as suas restrições temporais em plataformas uniprocessador. Enquanto os sistemas de tempo real se começam a adaptar ao paradigma de computação paralela, há uma crescente aposta em integrar execuções de tempo real com aplicações interativas no mesmo hardware, num mundo em que a tecnologia se torna cada vez mais pequena, leve, ubíqua, e portável. Esta integração requer soluções de escalonamento que simultaneamente garantam os requisitos temporais das tarefas de tempo real e mantenham um nível aceitável de QoS para as restantes execuções. Para tal, torna-se imperativo que as aplicações de tempo real paralelizem, de forma a minimizar os seus tempos de resposta e maximizar a utilização dos recursos de processamento. Isto introduz uma nova dimensão ao problema do escalonamento, que tem de responder de forma correcta a novos requisitos de execução imprevisíveis e rapidamente conjeturar o mapeamento de tarefas que melhor beneficie os critérios de performance do sistema. A técnica de escalonamento baseado em servidores permite reservar uma fração da capacidade de processamento para a execução de tarefas de tempo real, e assegurar que os efeitos de latência na sua execução não afectam as reservas estipuladas para outras execuções. No caso de tarefas escalonadas pelo tempo de execução máximo, ou tarefas com tempos de execução variáveis, torna-se provável que a largura de banda estipulada não seja consumida por completo. Para melhorar a utilização do sistema, os algoritmos de partilha de largura de banda (capacity-sharing) doam a capacidade não utilizada para a execução de outras tarefas, mantendo as garantias de isolamento entre servidores. Com eficiência comprovada em termos de espaço, tempo, e comunicação, o mecanismo de work-stealing tem vindo a ganhar popularidade como metodologia para o escalonamento de tarefas com paralelismo dinâmico e irregular. O algoritmo p-CSWS combina escalonamento baseado em servidores com capacity-sharing e work-stealing para cobrir as necessidades de escalonamento dos sistemas abertos de tempo real. Enquanto o escalonamento em servidores permite partilhar os recursos de processamento sem interferências a nível dos atrasos, uma nova política de work-stealing que opera sobre o mecanismo de capacity-sharing aplica uma exploração de paralelismo que melhora os tempos de resposta das aplicações e melhora a utilização do sistema. Esta tese propõe uma implementação do algoritmo p-CSWS para o Linux. Em concordância com a estrutura modular do escalonador do Linux, ´e definida uma nova classe de escalonamento que visa avaliar a aplicabilidade da heurística p-CSWS em circunstâncias reais. Ultrapassados os obstáculos intrínsecos `a programação da kernel do Linux, os extensos testes experimentais provam que o p-CSWS ´e mais do que um conceito teórico atrativo, e que a exploração heurística de paralelismo proposta pelo algoritmo beneficia os tempos de resposta das aplicações de tempo real, bem como a performance e eficiência da plataforma multiprocessador.

Going with the wind: The time for time-of-use tariffs

A Work Project, presented as part of the requirements for the Award of a Masters Degree in Management from the NOVA – School of Business and Economics

Avaliação epidemiológica dos internamentos através de um serviço de urgência de segunda linha (SII)

RESUMO: Introdução: A Comissão Nacional para a reestruturação dos Serviços de Saúde Mental em 2004, fez uma proposta de âmbito regional, ao nível da região de Saúde do Norte, levando a uma alteração da rede de referenciação hospitalar dos internamentos em psiquiatria. Métodos: Realizou-se um estudo epidemiológico longitudinal para avaliar factores preditivos de internamento através de um serviço de Urgência de Segunda Linha (SII), que de algum modo reflectisse o funcionamento dos serviços de internamento na região de Saúde do Norte, ao longo de 12 anos, nomeadamente relacionando com factores organizacionais contemplados no Plano Nacional de Saúde Mental e na reorganização da rede de referenciação hospitalar. Resultados: Durante os 12 anos do estudo, verificou-se um aumento estatisticamente significativo do número e duração de internamentos através do SII, com ponto de partida no ano 2008-2009, e de novo a partir do ano 2010-2011 (nº de internamentos), para o qual contribuiu a alteração da rede hospitalar na região de saúde do Norte, nomeadamente pelo facto do HMLemos, assumir a responsabilidade de novo, dos internamentos das áreas de Famalicão, Gondomar e Santa Maria da Feira. Em relação ao número de internamentos, e na análise exploratória, encontramos nas áreas hospitalares fora da área de influência do HMLemos, uma contribuição positiva significativa para o aumento do nº de internamentos ao longo dos anos com os Dx (290, 296, 297, 291, 309). Em relação à área do HMLemos restrita (PVC, STT, Matosinhos, Porto Ocidental), de referir a contribuição positiva significativa dos Dx 309 e 301, para o aumento do número de internamentos ao longo do tempo, sendo que a prevalência maior se mantém relacionada às Psicoses (Dx 295, 296 e 297). Não se concluiu por uma contribuição estatisticamente significativa ( positiva ou negativa), das variáveis independentes idade, sexo ou natureza do internamento em relação à variável dependente ( duração de internamentos/ano). Em relação á variável dependente (nº de internamentos/ano), relativamente aos doentes fora de área de influência do HMLemos, concluiu-se uma contribuição positiva estatisticamente significativa da variável independente idade. Conclusões: Através da análise exploratória foi possível perceber o esforço realizado pelos hospitais no sentido de melhorar a equidade e acessibilidade dos doentes à Saúde Mental, a par da reorganização da rede hospitalar. De destacar a necessidade de encontrar alternativas às situações de internamento, com menos critérios de gravidade diagnóstica, nomeadamente reforçar a importância da criação de consultas de crises nos respectivos Hospitais de Dia dos DPSM.----------------ABSTRACT:Introduction : The National Commission for the restructuring of mental health services in 2004 , has proposed at a regional level ( North Health Region), a change in the network of hospital referrals of admissions in psychiatry. Methods: We conducted a longitudinal epidemiological study to assess predictors of hospitalization through a Second Line Emergency Service ( SII) , that somehow reflect the operation of inpatient services in North Health Region, over 12 years, particularly relating to organizational factors included in the National Mental Health Plan and reorganization of the hospital referral network. Results: During the period of the study, there was a statistically significant increase in the number and duration of hospitalizations through the SII, with starting point in the year 2008-2009 and again from 2010-2011 (number of admissions) , for which counted the change of the hospital network referral in Northern health region , in particular because Hospital Magalhães Lemos (HMLemos) , took the new responsibility of admissions from areas of Famalicão, Gondomar and Santa Maria da Feira . Regarding the number of hospitalizations, in the exploratory analysis , we found in hospital areas outside the area of influence of HMLemos , a significant positive contribution to the increase in number of admissions over the years with Diagnosis of 290, 296, 297, 291 , 309 in the ICD-9. With respect to the restricted area of HMLemos (PVC, STT , Porto Ocidental and Matosinhos) , we found a significant positive contribution of Diagnosis 309 and 301, to increase the number of hospitalizations over time, with higher prevalence rates remaining the psychoses ( Dx 295, 296 and 297 ) . Did not conclude for any statistically significant contribution (positive or negative) of the independent variables age, sex and nature of admission to the dependent variable (duration of hospitalization / year). In relation to the dependent variable (number of admissions / year) relative to patients outside the area of influence of HMLemos, it was found a statistically significant positive contribution of the independent variable age . Conclusions: Through the exploratory analysis, it was possible to see the efforts made by hospitals to improve the accessibility of patients to Mental Health, throughout the hospital network reorganization. Its important to highlight the need to find alternatives to inpatient admissions in those with less gravity diagnostic criteria, reinforcing the importance of creating specific crisis consultations in Day Hospital regime.

Wood specific gravity of trees and forest types in the southern Peruvian Amazon

Estimates of terrestrial biomass depend critically on reliable information about the specific gravity of the wood of forest trees. The study reported on here was carried out in the southern Peruvian Amazon and involved collection of wood samples from trees (70 spp.) in intact forest stands. Results demonstrate the high degree of variability in specific gravity (ovendry weight/green volume) in trees at single locations. Three forest types (swamp, high terrace forest with alluvial soil, and sandy-soil forest) had values close to the average reported for tropical forest woods (.69). Two early successional forest types, which make up as much as 12% of the total vegetated area in this part of the Amazon, had values significantly lower (.40). An increase in specific gravity with increasing age of the tree, which has been reported in some spe cies of tropical-forest woods, is seen in a positive relationship between specific gravity and di ameter for a species prevalent in one plot. Increases in specific gravity with tree and forest age may be significant in estimating changes in carbon stores over time.

Forecasting time series combining Holt-Winters and bootstrap approaches

Publicado em "AIP Conference Proceedings", Vol. 1648

Time-variations of equivalent water heights'from Grace Mission and in-situ river stages in the Amazon basin

Gravity Recovery and Climate Experiment (GRACE) mission is dedicated to measuring temporal variations of the Earth's gravity field. In this study, the Stokes coefficients made available by Groupe de Recherche en Géodésie Spatiale (GRGS) at a 10-day interval were converted into equivalent water height (EWH) for a ~4-year period in the Amazon basin (from July-2002 to May-2006). The seasonal amplitudes of EWH signal are the largest on the surface of Earth and reach ~ 1250mm at that basin's center. Error budget represents ~130 mm of EWH, including formal errors on Stokes coefficient, leakage errors (12 ~ 21 mm) and spectrum truncation (10 ~ 15 mm). Comparison between in situ river level time series measured at 233 ground-based hydrometric stations (HS) in the Amazon basin and vertically-integrated EWH derived from GRACE is carried out in this paper. Although EWH and HS measure different water bodies, in most of the cases a high correlation (up to ~80%) is detected between the HS series and EWH series at the same site. This correlation allows adjusting linear relationships between in situ and GRACE-based series for the major tributaries of the Amazon river. The regression coefficients decrease from up to down stream along the rivers reaching the theoretical value 1 at the Amazon's mouth in the Atlantic Ocean. The variation of the regression coefficients versus the distance from estuary is analysed for the largest rivers in the basin. In a second step, a classification of the proportionality between in situ and GRACE time-series is proposed.

Specification and testing of multiplicative time-varying GARCH models with applications

In this article, we develop a specification technique for building multiplicative time-varying GARCH models of Amado and Teräsvirta (2008, 2013). The variance is decomposed into an unconditional and a conditional component such that the unconditional variance component is allowed to evolve smoothly over time. This nonstationary component is defined as a linear combination of logistic transition functions with time as the transition variable. The appropriate number of transition functions is determined by a sequence of specification tests. For that purpose, a coherent modelling strategy based on statistical inference is presented. It is heavily dependent on Lagrange multiplier type misspecification tests. The tests are easily implemented as they are entirely based on auxiliary regressions. Finite-sample properties of the strategy and tests are examined by simulation. The modelling strategy is illustrated in practice with two real examples: an empirical application to daily exchange rate returns and another one to daily coffee futures returns.

Predicting SO2 pollution incidents by means of additive models with optimum variable selection

The aim of this paper is to predict time series of SO2 concentrations emitted by coal-fired power stations in order to estimate in advance emission episodes and analyze the influence of some meteorological variables in the prediction. An emission episode is said to occur when the series of bi-hourly means of SO2 is greater than a specific level. For coal-fired power stations it is essential to predict emission epi- sodes sufficiently in advance so appropriate preventive measures can be taken. We proposed a meth- odology to predict SO2 emission episodes based on using an additive model and an algorithm for variable selection. The methodology was applied to the estimation of SO2 emissions registered in sampling lo- cations near a coal-fired power station located in Northern Spain. The results obtained indicate a good performance of the model considering only two terms of the time series and that the inclusion of the meteorological variables in the model is not significant.

Blow-up and finite time extinction for p(x, t)-curl systems arising in electromagnetism

"Available online 22 March 2016"

Evaluation of the determinism of the actuator sensor interface (ASI) on programmable logic controllers (PLC), personal computers (PC) and real-time operating systems (RTOS)

The purpose of this study was to evaluate the determinism of the AS-lnterface network and the 3 main families of control systems, which may use it, namely PLC, PC and RTOS. During the course of this study the PROFIBUS and Ethernet field level networks were also considered in order to ensure that they would not introduce unacceptable latencies into the overall control system. This research demonstrated that an incorrectly configured Ethernet network introduces unacceptable variable duration latencies into the control system, thus care must be exercised if the determinism of a control system is not to be compromised. This study introduces a new concept of using statistics and process capability metrics in the form of CPk values, to specify how suitable a control system is for a given control task. The PLC systems, which were tested, demonstrated extremely deterministic responses, but when a large number of iterations were introduced in the user program, the mean control system latency was much too great for an AS-I network. Thus the PLC was found to be unsuitable for an AS-I network if a large, complex user program Is required. The PC systems, which were tested were non-deterministic and had latencies of variable duration. These latencies became extremely exaggerated when a graphing ActiveX was included in the control application. These PC systems also exhibited a non-normal frequency distribution of control system latencies, and as such are unsuitable for implementation with an AS-I network. The RTOS system, which was tested, overcame the problems identified with the PLC systems and produced an extremely deterministic response, even when a large number of iterations were introduced in the user program. The RTOS system, which was tested, is capable of providing a suitable deterministic control system response, even when an extremely large, complex user program is required.

Polynomial-time complexity for instances of the endomorphism problem in free groups

We say the endomorphism problem is solvable for an element W in a free group F if it can be decided effectively whether, given U in F, there is an endomorphism Φ of F sending W to U. This work analyzes an approach due to C. Edmunds and improved by C. Sims. Here we prove that the approach provides an efficient algorithm for solving the endomorphism problem when W is a two- generator word. We show that when W is a two-generator word this algorithm solves the problem in time polynomial in the length of U. This result gives a polynomial-time algorithm for solving, in free groups, two-variable equations in which all the variables occur on one side of the equality and all the constants on the other side.

Real time GPS precise positioning of LEO satellites

Report for the scientific sojourn at the German Aerospace Center (DLR) , Germany, during June and July 2006. The main objective of the two months stay has been to apply the techniques of LEO (Low Earth Orbiters) satellites GPS navigation which DLR currently uses in real time navigation. These techniques comprise the use of a dynamical model which takes into account the precise earth gravity field and models to account for the effects which perturb the LEO’s motion (such as drag forces due to earth’s atmosphere, solar pressure, due to the solar radiation impacting on the spacecraft, luni-solar gravity, due to the perturbation of the gravity field for the sun and moon attraction, and tidal forces, due to the ocean and solid tides). A high parameterized software was produced in the first part of work, which has been used to asses which accuracy could be reached exploring different models and complexities. The objective was to study the accuracy vs complexity, taking into account that LEOs at different heights have different behaviors. In this frame, several LEOs have been selected in a wide range of altitudes, and several approaches with different complexity have been chosen. Complexity is a very important issue, because processors onboard spacecrafts have very limited computing and memory resources, so it is mandatory to keep the algorithms simple enough to let the satellite process it by itself.

Variable viral clearance despite adequate ganciclovir plasma levels during valganciclovir treatment for cytomegalovirus disease in D+/R- transplant recipients.

ABSTRACT: BACKGROUND: Valganciclovir, the oral prodrug of ganciclovir, has been demonstrated equivalent to iv ganciclovir for CMV disease treatment in solid organ transplant recipients. Variability in ganciclovir exposure achieved with valganciclovir could be implicated as a contributing factor for explaining variations in the therapeutic response. This prospective observational study aimed to correlate clinical and cytomegalovirus (CMV) viral load response (DNAemia) with ganciclovir plasma concentrations in patients treated with valganciclovir for CMV infection/disease. METHODS: Seven CMV D+/R- transplant recipients (4 kidney, 2 liver and 1 heart) were treated with valganciclovir (initial dose was 900-1800 mg/day for 3-6.5 weeks, followed by 450-900 mg/day for 2-9 weeks). DNAemia was monitored by real time quantitative PCR and ganciclovir plasma concentration was measured at trough (Ctrough) and 3 h after drug administration (C3h) by HPLC. RESULTS: Four patients presented with CMV syndrome, two had CMV tissue-invasive disease after prophylaxis discontinuation, and one liver recipient was treated pre-emptively for asymptomatic rising CMV viral load 5 weeks post-transplantation in the absence of prophylaxis. CMV DNAemia decreased during the first week of treatment in all recipients except in one patient (median decrease: -1.2 log copies/mL, range: -1.8 to 0) despite satisfactory ganciclovir exposure (AUC0-12 = 48 mg.h/L, range for the 7 patients: 40-118 mg.h/L). Viral clearance was obtained in five patients after a median of time of 34 days (range: 28-82 days). Two patients had recurrent CMV disease despite adequate ganciclovir exposure (65 mg.h/L, range: 44-118 mg.h/L). CONCLUSIONS: Valganciclovir treatment for CMV infection/disease in D+/R- transplant recipients can thus result in variable viral clearance despite adequate ganciclovir plasma concentrations, probably correlating inversely with anti-CMV immune responses after primary infection.

Antibody-mediated status epilepticus: a retrospective multicenter survey.

Background: In recent years, an increasing number of auto-antibodies (AB) have been detected in the CSF and serum of patients with new onset epilepsy. Some of these patients develop convulsive or nonconvulsive status epilepticus (AB-SE), necessitating intensive medical care and administration of multiple antiepileptic and immunomodulatory treatments of uncertain effectiveness. Objectives: In this retrospective multicenter survey we aimed to determine the spectrum of gravity, the duration and the prognosis of the disorder. In addition, we sought to identify the antibodies associated with this condition, as well as determine whether there is a most effective treatment regime. Methods: 12 European Neurology University Clinics, with extensive experience in the treatment of SE patients, were sent a detailed questionnaire regarding symptoms and treatment of AB-SE patients. Seven centers responded positively, providing a total of 13 patients above the age of 16. Results: AB-SE affects mainly women (12/13, 92%) with a variable age at onset (17-69 years, median: 25 years). The duration of the disease is also variable (10 days to 12 years, median: 2 months). Only the 3 oldest patients died (55-69 years). Most patients were diagnosed with anti NMDAR encephalitis (8/13) and had oligoclonal bands in the CSF (9/13). No specific treatment regimen (antiepileptic, immunomodulatory) was found to be clearly superior. Most of the surviving 10 patients (77%) recovered completely or nearly so within 2 years of index poststatus. Conclusion: AB-SE is a severe but potentially reversible condition. Long duration does not seem to imply fatal outcome; however, age older than 50 years at time of onset appears to be a risk factor for death. There was no evidence for an optimal antiepileptic or immunomodulatory treatment. A prospective multicenter study is warranted in order to stratify the optimal treatment algorithm, determine clear risk factors of unfavorable outcome and long-term prognosis.

Reliable acetabular cup orientation with a new gravity-assisted guidance system.

Acetabular cup orientation is a key factor determining hip stability, and standard mechanical guides have shown little help in improving alignment. An in vitro study was carried out to compare the accuracy and precision of a new gravity-assisted guidance system with a standard mechanical guide. Three hundred ten cups were impacted by 5 surgeons, and the final cup orientation was measured. With the new guide, the average error in anteversion was 0.4 degrees , compared with 10.4 degrees with the standard guide and 0.3 degrees and -4.7 degrees , respectively, for abduction angles. The average time required for orienting the cups was similar for both guides. The accuracy and reproducibility obtained with the new guide were better (P < .0001). These good results would require a clinical validation.