Dynamic load balancing of distributed memory parallel computational mechanics using unstructured meshes for multi-physical modelling

As the complexity of parallel applications increase, the performance limitations resulting from computational load imbalance become dominant. Mapping the problem space to the processors in a parallel machine in a manner that balances the workload of each processors will typically reduce the run-time. In many cases the computation time required for a given calculation cannot be predetermined even at run-time and so static partition of the problem returns poor performance. For problems in which the computational load across the discretisation is dynamic and inhomogeneous, for example multi-physics problems involving fluid and solid mechanics with phase changes, the workload for a static subdomain will change over the course of a computation and cannot be estimated beforehand. For such applications the mapping of loads to process is required to change dynamically, at run-time in order to maintain reasonable efficiency. The issue of dynamic load balancing are examined in the context of PHYSICA, a three dimensional unstructured mesh multi-physics continuum mechanics computational modelling code.

Mesh partitioning and load balancing for distributed memory parallel systems

A method is outlined for optimising graph partitions which arise in mapping unstructured mesh calculations to parallel computers. The method employs a relative gain iterative technique to both evenly balance the workload and minimise the number and volume of interprocessor communications. A parallel graph reduction technique is also briefly described and can be used to give a global perspective to the optimisation. The algorithms work efficiently in parallel as well as sequentially and when combined with a fast direct partitioning technique (such as the Greedy algorithm) to give an initial partition, the resulting two-stage process proves itself to be both a powerful and flexible solution to the static graph-partitioning problem. Experiments indicate that the resulting parallel code can provide high quality partitions, independent of the initial partition, within a few seconds. The algorithms can also be used for dynamic load-balancing, reusing existing partitions and in this case the procedures are much faster than static techniques, provide partitions of similar or higher quality and, in comparison, involve the migration of a fraction of the data.

Parallel dynamic graph-partitioning for unstructured meshes

A parallel method for the dynamic partitioning of unstructured meshes is described. The method introduces a new iterative optimisation technique known as relative gain optimisation which both balances the workload and attempts to minimise the interprocessor communications overhead. Experiments on a series of adaptively refined meshes indicate that the algorithm provides partitions of an equivalent or higher quality to static partitioners (which do not reuse the existing partition) and much more rapidly. Perhaps more importantly, the algorithm results in only a small fraction of the amount of data migration compared to the static partitioners.

Dynamic multi-partitioning for parallel finite element applications

The central product of the DRAMA (Dynamic Re-Allocation of Meshes for parallel Finite Element Applications) project is a library comprising a variety of tools for dynamic re-partitioning of unstructured Finite Element (FE) applications. The input to the DRAMA library is the computational mesh, and corresponding costs, partitioned into sub-domains. The core library functions then perform a parallel computation of a mesh re-allocation that will re-balance the costs based on the DRAMA cost model. We discuss the basic features of this cost model, which allows a general approach to load identification, modelling and imbalance minimisation. Results from crash simulations are presented which show the necessity for multi-phase/multi-constraint partitioning components.

A parallel numerical solver using hierarchically tiled arrays

Solving linear systems is an important problem for scientific computing. Exploiting parallelism is essential for solving complex systems, and this traditionally involves writing parallel algorithms on top of a library such as MPI. The SPIKE family of algorithms is one well-known example of a parallel solver for linear systems. The Hierarchically Tiled Array data type extends traditional data-parallel array operations with explicit tiling and allows programmers to directly manipulate tiles. The tiles of the HTA data type map naturally to the block nature of many numeric computations, including the SPIKE family of algorithms. The higher level of abstraction of the HTA enables the same program to be portable across different platforms. Current implementations target both shared-memory and distributed-memory models. In this thesis we present a proof-of-concept for portable linear solvers. We implement two algorithms from the SPIKE family using the HTA library. We show that our implementations of SPIKE exploit the abstractions provided by the HTA to produce a compact, clean code that can run on both shared-memory and distributed-memory models without modification. We discuss how we map the algorithms to HTA programs as well as examine their performance. We compare the performance of our HTA codes to comparable codes written in MPI as well as current state-of-the-art linear algebra routines.

A unified model for the long and high jump

A simple model based on, the maximum energy that an athlete can produce in a small time interval is used to describe the high and long jump. Conservation of angular momentum is used to explain why an athlete should, run horizontally to perform a vertical jump. Our results agree with world records. (c) 2005 American Association of Physics Teachers.

Implementación eficiente de redes neuronales autoorganizadas en CUDA C/C++.

En este documento se expondrá una implementación del problema del viajante de comercio usando una implementación personalizada de un mapa auto-organizado basándose en soluciones anteriores y adaptándolas a la arquitectura CUDA, haciendo a la vez una comparativa de la implementación eficiente en CUDA C/C++ con la implementación de las funciones de GPU incluidas en el Parallel Computing Toolbox de Matlab. La solución que se da reduce en casi un cuarto las iteraciones necesarias para llegar a una solución buena del problema mencionado, además de la mejora inminente del uso de las arquitecturas paralelas. En esta solución se estudia la mejora en tiempo que se consigue con el uso específico de la memoria compartida, siendo esta una de las herramientas más potentes para mejorar el rendimiento. En lo referente a los tiempos de ejecución, se llega a concluir que la mejor solución es el lanzamiento de un kernel de CUDA desde Matlab a través de la funcionalidad incluida en el Parallel Computing Toolbox.

Sistema integrado de generación de mapas de profundidad para videoconferencia 3D en tiempo real

This paper presents the implementation of a high quality real-time 3D video system intended for 3D videoconferencing -- Basically, the system is able to extract depth information from a pair of images coming from a short-baseline camera setup -- The system is based on the use of a variant of the adaptive support-weight algorithm to be applied on GPU-based architectures -- The reason to do it is to get real-time results without compromising accuracy and also to reduce costs by using commodity hardware -- The complete system runs over the GStreamer multimedia software platform to make it even more flexible -- Moreover, an autoestereoscopic display has been used as the end-up terminal for 3D content visualization

Drying of rye (Lolium miltiflorum l.) in fixed bed with parallel air flow.

This work aimed to study the drying of ryegrass seeds (Lolium multiflorum L.) in fixed bed dryer with parallel air flow.

Geometrico-static modelling of continuum parallel robots

In this thesis, we explore three methods for the geometrico-static modelling of continuum parallel robots. Inspired by biological trunks, tentacles and snakes, continuum robot designs can reach confined spaces, manipulate objects in complex environments and conform to curvilinear paths in space. In addition, parallel continuum manipulators have the potential to inherit some of the compactness and compliance of continuum robots while retaining some of the precision, stability and strength of rigid-links parallel robots. Subsequently, the foundation of our work is performed on slender beam by applying the Cosserat rod theory, appropriate to model continuum robots. After that, three different approaches are developed on a case study of a planar parallel continuum robot constituted of two connected flexible links. We solve the forward and inverse geometrico-static problem namely by using (a) shooting methods to obtain a numerical solution, (b) an elliptic method to find a quasi-analytical solution, and (c) the Corde model to perform further model analysis. The performances of each of the studied methods are evaluated and their limits are highlighted. This thesis is divided as follows. Chapter one gives the introduction on the field of the continuum robotics and introduce the parallel continuum robots that is studied in this work. Chapter two describe the geometrico-static problem and gives the mathematical description of this problem. Chapter three explains the numerical approach with the shooting method and chapter four introduce the quasi-analytical solution. Then, Chapter five introduce the analytic method inspired by the Corde model and chapter six gives the conclusions of this work.

Islet Neogenesis-associated Protein (ingap): The Role Of Its Endogenous Production As A Positive Modulator Of Insulin Secretion.

Islet neogenesis-associated protein (INGAP) is a peptide found in pancreatic exocrine-, duct- and islet- non-β-cells from normal hamsters. Its increase induced by either its exogenous administration or by the overexpression of its gene enhances β-cell secretory function and increases β-cell mass by a combination of stimulation of cell replication and islet neogenesis and reduction of β-cell apoptosis. We studied the potential modulatory role of endogenous INGAP in insulin secretion using two different experimental approaches. Hamster islets transfected with INGAP-small interfering RNA (INGAP-siRNA) were used to study glucose-stimulated insulin secretion (GSIS). In parallel, freshly isolated islets were incubated with high glucose and the same concentration of either a specific anti-INGAP rabbit serum or normal rabbit serum. INGAP-siRNA transfected islets reduced their INGAP mRNA and protein content by 35.1% and 47.2%, respectively whereas GSIS decreased by 25.8%. GSIS by transfected islets attained levels comparable to those recorded in control islets when INGAP pentadecapeptide (INGAP-PP) was added to the culture medium. INGAP antibody in the medium decreased significantly GSIS in a dose-dependent manner. These results indicate that endogenous INGAP plays a physiological positive modulatory role in insulin secretion, supporting its possible use in the treatment of prediabetes and Type 2 diabetes.

Hdl Levels And Oxidizability During Myocardial Infarction Are Associated With Reduced Endothelial-mediated Vasodilation And Nitric Oxide Bioavailability.

Acute phase response modifies high-density lipoprotein (HDL) into a dysfunctional particle that may favor oxidative/inflammatory stress and eNOS dysfunction. The present study investigated the impact of this phenomenon on patients presenting ST-elevation myocardial infarction (STEMI). Plasma was obtained from 180 consecutive patients within the first 24-h of onset of STEMI symptoms (D1) and after 5 days (D5). Nitrate/nitrite (NOx) and lipoproteins were isolated by gradient ultracentrifugation. The oxidizability of low-density lipoprotein incubated with HDL (HDLaoxLDL) and the HDL self-oxidizability (HDLautox) were measured after CuSO4 co-incubation. Anti-inflammatory activity of HDL was estimated by VCAM-1 secretion by human umbilical vein endothelial cells after incubation with TNF-α. Flow-mediated dilation (FMD) was assessed at the 30(th) day (D30) after STEMI. Among patients in the first tertile of admission HDL-Cholesterol (<33 mg/dL), the increment of NOx from D1 to D5 [6.7(2; 13) vs. 3.2(-3; 10) vs. 3.5(-3; 12); p = 0.001] and the FMD adjusted for multiple covariates [8.4(5; 11) vs 6.1(3; 10) vs. 5.2(3; 10); p = 0.001] were higher than in those in the second (33-42 mg/dL) or third (>42 mg/dL) tertiles, respectively. From D1 to D5, there was a decrease in HDL size (-6.3 ± 0.3%; p < 0.001) and particle number (-22.0 ± 0.6%; p < 0.001) as well as an increase in both HDLaoxLDL (33%(23); p < 0.001) and HDLautox (65%(25); p < 0.001). VCAM-1 secretion after TNF-a stimulation was reduced after co-incubation with HDL from healthy volunteers (-24%(33); p = 0.009), from MI patients at D1 (-23%(37); p = 0.015) and at D30 (-22%(24); p = 0.042) but not at D5 (p = 0.28). During STEMI, high HDL-cholesterol is associated with a greater decline in endothelial function. In parallel, structural and functional changes in HDL occur reducing its anti-inflammatory and anti-oxidant properties.

Frequency Of The Allelic Variant C.1150t > C In Exon 10 Of The Fibroblast Growth Factor Receptor 3 (fgfr3) Gene Is Not Increased In Patients With Pathogenic Mutations And Related Chondrodysplasia Phenotypes.

Mutations in the FGFR3 gene cause the phenotypic spectrum of FGFR3 chondrodysplasias ranging from lethal forms to the milder phenotype seen in hypochondroplasia (Hch). The p.N540K mutation in the FGFR3 gene occurs in ∼70% of individuals with Hch, and nearly 30% of individuals with the Hch phenotype have no mutations in the FGFR3, which suggests genetic heterogeneity. The identification of a severe case of Hch associated with the typical mutation c.1620C > A and the occurrence of a c.1150T > C change that resulted in a p.F384L in exon 10, together with the suspicion that this second change could be a modulator of the phenotype, prompted us to investigate this hypothesis in a cohort of patients. An analysis of 48 patients with FGFR3 chondrodysplasia phenotypes and 330 healthy (control) individuals revealed no significant difference in the frequency of the C allele at the c.1150 position (p = 0.34). One patient carrying the combination `pathogenic mutation plus the allelic variant c.1150T > C' had a typical achondroplasia (Ach) phenotype. In addition, three other patients with atypical phenotypes showed no association with the allelic variant. Together, these results do not support the hypothesis of a modulatory role for the c.1150T > C change in the FGFR3 gene.

Proinflammatory Activity Of Primarily Infected Endodontic Content Against Macrophages After Different Phases Of The Root Canal Therapy.

This study investigated the presence of target bacterial species and the levels of endotoxins in teeth with apical periodontitis. Levels of inflammatory mediators (interleukin [IL]-1β and tumor necrosis factor [TNF]-α) were determined after macrophage stimulation with endodontic content after different phases of endodontic therapy using different irrigants. Thirty primarily infected root canals were randomly assigned into 3 groups according to the irrigant used for root canal preparation (n = 10 per group): GI: 2.5% sodium hypochlorite, GII: 2% chlorhexidine gel, and GIII (control group): saline solution. Root canal samples were taken by using paper points before (s1) and after root canal instrumentation (s2), subsequently to 17% EDTA (s3), after 30 days of intracanal medication (Ca[OH]2 + saline solution) (s4), and before root canal obturation (s5). Polymerase chain reaction (16S recombinant DNA) and limulus amebocyte lysate assay were used for bacterial and endotoxin detection, respectively. Macrophages were stimulated with the root canal contents for IL-1β/TNF-α measurement using enzyme-linked immunosorbent assay. Porphyromonas gingivalis (17/30), Porphyromonas endodontalis (15/30), and Prevotella nigrescens (11/30) were the most prevalent bacterial species. At s1, endotoxins were detected in 100% of the root canals (median = 32.43 EU/mL). In parallel, substantial amounts of IL-1β and TNF-α were produced by endodontic content-stimulated macrophages. At s2, a significant reduction in endotoxin levels was observed in all groups, with GI presenting the greatest reduction (P < .05). After a root canal rinse with EDTA (s3), intracanal medication (s4), and before root canal obturation (s5), endotoxin levels reduced without differences between groups (P < .05). IL-1β and TNF-α release decreased proportionally to the levels of residual endotoxin (P < .05). Regardless of the use of sodium hypochlorite or CHX, the greatest endotoxin reduction occurs after chemomechanical preparation. Increasing steps of root canal therapy associated with intracanal medication enhances endotoxin reduction, leading to a progressively lower activation of proinflammatory cells such as macrophages.

The C-terminal Region Of The Human P23 Chaperone Modulates Its Structure And Function.

The p23 protein is a chaperone widely involved in protein homeostasis, well known as an Hsp90 co-chaperone since it also controls the Hsp90 chaperone cycle. Human p23 includes a β-sheet domain, responsible for interacting with Hsp90; and a charged C-terminal region whose function is not clear, but seems to be natively unfolded. p23 can undergo caspase-dependent proteolytic cleavage to form p19 (p231-142), which is involved in apoptosis, while p23 has anti-apoptotic activity. To better elucidate the function of the human p23 C-terminal region, we studied comparatively the full-length human p23 and three C-terminal truncation mutants: p23₁₋₁₁₇; p23₁₋₁₃₁ and p23₁₋₁₄₂. Our data indicate that p23 and p19 have distinct characteristics, whereas the other two truncations behave similarly, with some differences to p23 and p19. We found that part of the C-terminal region can fold in an α-helix conformation and slightly contributes to p23 thermal-stability, suggesting that the C-terminal interacts with the β-sheet domain. As a whole, our results suggest that the C-terminal region of p23 is critical for its structure-function relationship. A mechanism where the human p23 C-terminal region behaves as an activation/inhibition module for different p23 activities is proposed.