A Massively-Parallel Mixed-Mode Computer Designed to Support

Total Page:16

File Type:pdf, Size:1020Kb

A Massively-Parallel Mixed-Mode Computer Designed to Support This paper appeared in th International Parallel Processing Symposium Proc of nd Work shop on Heterogeneous Processing pages NewportBeach CA April Triton A MassivelyParallel MixedMo de Computer Designed to Supp ort High Level Languages Christian G Herter Thomas M Warschko Walter F Tichy and Michael Philippsen University of Karlsruhe Dept of Informatics Postfach D Karlsruhe Germany Mo dula Abstract Mo dula pronounced Mo dulastar is a small ex We present the architectureofTriton a scalable tension of Mo dula for massively parallel program mixedmode SIMDMIMD paral lel computer The ming The programming mo del of Mo dula incor novel features of Triton are p orates b oth data and control parallelism and allows hronous and asynchronous execution mixed sync Support for highlevel machineindependent pro Mo dula is problemorientedinthesensethatthe gramming languages programmer can cho ose the degree of parallelism and mix the control mo de SIMD or MIMDlike as need Fast SIMDMIMD mode switching ed bytheintended algorithm Parallelism maybe nested to arbitrary depth Pro cedures may b e called Special hardware for barrier synchronization of from sequential or parallel contexts and can them multiple process groups selves generate parallel activity without any restric tions Most Mo dula programs can b e translated into ecient co de for b oth SIMD and MIMD archi A selfrouting deadlockfreeperfect shue inter tectures connect with latency hiding Overview of language extensions The architecture is the outcomeofanintegrated de Mo dula extends Mo dula with the following two sign in which a machineindependent programming language constructs language optimizing compiler and paral lel computer were designed handinhand Parallelism can b e created in Mo dula pro grams only by means of the FORALL statement There are twoversions of this statement a syn Intro duction chronous and an asynchronous one The goal of adequate programmability of parallel ma chines can b est b e achieved by tightly coupling the de The distribution of array data may b e optional sign of machineindep endent programming languages ly sp ecied by socalled al locators in a machine compilers and parallel hardware In the past the de indep endentway Allo cators do not haveany se velopment of parallel computers has b een mainly driv mantic meaning they are hints for the compiler en by hardware considerations without regard to the law of the weakest link Because of the compactness and simplicity of the ex The overall p erformance of a parallel computer is de tensions they could easily b e incorp orated in other termined by the p erformance of its slowest part with imp erative programming languages suchasFortran resp ect to the requirements of the software C or Ada Ignoring software requirements has resulted in unsat A synchronous FORALL statement creates a set of isfactory p erformance of these machines on machine synchronously running pro cesses The number of the indep endent parallel programs Toavoid these short pro cesses is determined bytheFORALL statement and comings and to show that highlevel parallel program is not limited to the numberofthePEsofthema ming do es not necessarily lead to p o or p erformance chine As long as there is no branch in the control we sp ecically analyzed the requirements of program ow of the statements inside a FORALL statementthe ming languages and compilers b efore designing the semantics of the execution is equivalent to a SIMD ma hardware chine executing those statements If there are branch The following section outlines the parallel program es in the control ow that dene alternatives likeIF ming language Mo dula and derives general require THEN ELSE or CASE the pro cesses are partitioned ments for parallel computers In section we describ e into several groups Each of those groups is executing the architecture of Triton We emphasize those fea one branch in the control ow The pro cesses b elong tures which arose from software requirements ing to one group execute synchronously in SIMD style but groups are allowed to execute concurrently with are signicant and therefore require that the compil resp ect to each other There is no assumption ab out er detect the faster cases However it is often imp ossi the relative sp eeds of two pro cesses in dierent groups ble to know statically for which case to optimize For instance we found that in many cases it was imp ossi In contrast to the synchronous FORALL statement ble to determine in the compiler whether a pro cedure an asynchronous FORALL statement creates a set of in would access lo cal or nonlo cal memory The generat dep endent pro cesses running at an unsp ecied relative ed co de thus has to check all three cases at runtime sp eed Common to b oth variants is that the termi Such simple frequent and dynamic analyses could b e nation of a FORALL statement is determined bythe done more eciently in hardware termination of the last pro cess created by the FORALL Barrier synchronization is extremely frequent Basi statement The end of a FORALL statement always callyevery communication requires a barrier synchro denes a synchronization barrier For further details nization The reasoning is as follows Communication ab out the language see for detailed discussion of sends or receives data Unless communication is re compilation techniques and optimization see dundant there must b e a write b etween two succes Software requirements for parallel sive communication calls to the same cell It follows computers that a synchronization op eration must b e placed some On distributed memory machines the distribution of where b etween one of the communication calls and the array data over the available pro cessors is a central write in order to avoid race conditions If many pro problem Two conicting goals data lo calityand cesses communicate at once as in massively parallel maximum degree of parallelism must b e recon machines this typ e of synchronization amounts to a ciled Data lo cality means that data elements should barrier No pro cess may pro ceed until all pro cesses b e stored lo cal to the pro cessors that need them to have completed either their write or communication minimize communication costs Perfect data lo cali op eration Because of its frequencywe need ty could b e achieved by employing a single pro ces fast barrier synchronization sor Parallelism which reduces runtime may unfor In principle the communication network could b e used tunately reduce lo cality and increase communication for barrier synchronization However communication costs Additional goals for data distribution are networks usually have high latency whichmake them exploiting sp ecial comm unication patterns supp orted to o slow for fast barrier synchronization These nets by hardware and generating simple address calcu are optimized for transp orting data while barrier syn lations to prevent addressing from b ecoming a domi chronization requires the transp ort of only one or two nant cost bits but must also implement a reduction or scan op Even with optimal layout of the data there will eration on these bits still b e communication in the general case In fact in An additional complication is that there are usu most massively parallel applications that are not triv ally several groups of pro cesses that need to com ially parallel communication is almost as frequentas municate among themselves necessitating multiple computation Therefore nonoverlapping barriers Consider for instance a the network MIPS measure must approach the CPU pip elined architecture in which one set of pro cesses MIPS measure passes data to another set Eachsetmayhaveto A second p erformanceoriented recommendation is an synchronize internally indep endent of the other Sim ilarly an IFstatement within a FORALL divides a indep endently op erating network with asynchronous set of pro cesses into two subsets whichmayhaveto message delivery synchronize indep endentlyThus weneed since it allows the delivery of packets concurrently barrier synchronization for multiple indep endent with computation That would enable the compiler to sets of pro cesses interleave computation and communication and thus to hide some of the communication latency On the other hand we recommend a Triton shared address space The p o or programmabilityoftodays parallel ma chines is a consequence of the fact that the design A shared address space do es not imply shared memo of these machines has b een driven mostly by hard ry it only means that every pro cessor can generate ad ware considerations Programmability seems to have dresses for the entire memory in the system System b een a secondary issue resulting in languages designed wide addresses are esp ecially imp ortant for p ointers sp ecically for a particular machine mo del Suchlan b ecause otherwise they would havetobesimulated guages do not satisfy the needs of programmers who quite ineciently in software Even the memory of need to write machineindep endent applications the control pro cessor eg the frontend of a SIMD Triton matches most of the recommendations of machine should b e part of the shared address space the previous section Lo oking from a general p ointof Furthermore and similar to the ab ove wecallfora view Triton is determined by the following state uniform memory access instructions ments Many parallel machines to dayprovide a set of instruc General Architecture Triton is a SAMD tions for accessing lo cal memory a second one for ac synchronousasynchronous instruction streams mut cessing memory in neighb ors and a third
Recommended publications
  • 2.5 Classification of Parallel Computers
    52 // Architectures 2.5 Classification of Parallel Computers 2.5 Classification of Parallel Computers 2.5.1 Granularity In parallel computing, granularity means the amount of computation in relation to communication or synchronisation Periods of computation are typically separated from periods of communication by synchronization events. • fine level (same operations with different data) ◦ vector processors ◦ instruction level parallelism ◦ fine-grain parallelism: – Relatively small amounts of computational work are done between communication events – Low computation to communication ratio – Facilitates load balancing 53 // Architectures 2.5 Classification of Parallel Computers – Implies high communication overhead and less opportunity for per- formance enhancement – If granularity is too fine it is possible that the overhead required for communications and synchronization between tasks takes longer than the computation. • operation level (different operations simultaneously) • problem level (independent subtasks) ◦ coarse-grain parallelism: – Relatively large amounts of computational work are done between communication/synchronization events – High computation to communication ratio – Implies more opportunity for performance increase – Harder to load balance efficiently 54 // Architectures 2.5 Classification of Parallel Computers 2.5.2 Hardware: Pipelining (was used in supercomputers, e.g. Cray-1) In N elements in pipeline and for 8 element L clock cycles =) for calculation it would take L + N cycles; without pipeline L ∗ N cycles Example of good code for pipelineing: §doi =1 ,k ¤ z ( i ) =x ( i ) +y ( i ) end do ¦ 55 // Architectures 2.5 Classification of Parallel Computers Vector processors, fast vector operations (operations on arrays). Previous example good also for vector processor (vector addition) , but, e.g. recursion – hard to optimise for vector processors Example: IntelMMX – simple vector processor.
    [Show full text]
  • Massively Parallel Computing with CUDA
    Massively Parallel Computing with CUDA Antonino Tumeo Politecnico di Milano 1 GPUs have evolved to the point where many real world applications are easily implemented on them and run significantly faster than on multi-core systems. Future computing architectures will be hybrid systems with parallel-core GPUs working in tandem with multi-core CPUs. Jack Dongarra Professor, University of Tennessee; Author of “Linpack” Why Use the GPU? • The GPU has evolved into a very flexible and powerful processor: • It’s programmable using high-level languages • It supports 32-bit and 64-bit floating point IEEE-754 precision • It offers lots of GFLOPS: • GPU in every PC and workstation What is behind such an Evolution? • The GPU is specialized for compute-intensive, highly parallel computation (exactly what graphics rendering is about) • So, more transistors can be devoted to data processing rather than data caching and flow control ALU ALU Control ALU ALU Cache DRAM DRAM CPU GPU • The fast-growing video game industry exerts strong economic pressure that forces constant innovation GPUs • Each NVIDIA GPU has 240 parallel cores NVIDIA GPU • Within each core 1.4 Billion Transistors • Floating point unit • Logic unit (add, sub, mul, madd) • Move, compare unit • Branch unit • Cores managed by thread manager • Thread manager can spawn and manage 12,000+ threads per core 1 Teraflop of processing power • Zero overhead thread switching Heterogeneous Computing Domains Graphics Massive Data GPU Parallelism (Parallel Computing) Instruction CPU Level (Sequential
    [Show full text]
  • CS 677: Parallel Programming for Many-Core Processors Lecture 1
    1 CS 677: Parallel Programming for Many-core Processors Lecture 1 Instructor: Philippos Mordohai Webpage: mordohai.github.io E-mail: [email protected] Objectives • Learn how to program massively parallel processors and achieve – High performance – Functionality and maintainability – Scalability across future generations • Acquire technical knowledge required to achieve the above goals – Principles and patterns of parallel programming – Processor architecture features and constraints – Programming API, tools and techniques 2 Important Points • This is an elective course. You chose to be here. • Expect to work and to be challenged. • If your programming background is weak, you will probably suffer. • This course will evolve to follow the rapid pace of progress in GPU programming. It is bound to always be a little behind… 3 Important Points II • At any point ask me WHY? • You can ask me anything about the course in class, during a break, in my office, by email. – If you think a homework is taking too long or is wrong. – If you can’t decide on a project. 4 Logistics • Class webpage: http://mordohai.github.io/classes/cs677_s20.html • Office hours: Tuesdays 5-6pm and by email • Evaluation: – Homework assignments (40%) – Quizzes (10%) – Midterm (15%) – Final project (35%) 5 Project • Pick topic BEFORE middle of the semester • I will suggest ideas and datasets, if you can’t decide • Deliverables: – Project proposal – Presentation in class – Poster in CS department event – Final report (around 8 pages) 6 Project Examples • k-means • Perceptron • Boosting – General – Face detector (group of 2) • Mean Shift • Normal estimation for 3D point clouds 7 More Ideas • Look for parallelizable problems in: – Image processing – Cryptanalysis – Graphics • GPU Gems – Nearest neighbor search 8 Even More… • Particle simulations • Financial analysis • MCMC • Games/puzzles 9 Resources • Textbook – Kirk & Hwu.
    [Show full text]
  • Core Processors
    UNIVERSITY OF CALIFORNIA Los Angeles Parallel Algorithms for Medical Informatics on Data-Parallel Many-Core Processors A dissertation submitted in partial satisfaction of the requirements for the degree Doctor of Philosophy in Computer Science by Maryam Moazeni 2013 © Copyright by Maryam Moazeni 2013 ABSTRACT OF THE DISSERTATION Parallel Algorithms for Medical Informatics on Data-Parallel Many-Core Processors by Maryam Moazeni Doctor of Philosophy in Computer Science University of California, Los Angeles, 2013 Professor Majid Sarrafzadeh, Chair The extensive use of medical monitoring devices has resulted in the generation of tremendous amounts of data. Storage, retrieval, and analysis of such data require platforms that can scale with data growth and adapt to the various behavior of the analysis and processing algorithms. In recent years, many-core processors and more specifically many-core Graphical Processing Units (GPUs) have become one of the most promising platforms for high performance processing of data, due to the massive parallel processing power they offer. However, many of the algorithms and data structures used in medical and bioinformatics systems do not follow a data-parallel programming paradigm, and hence cannot fully benefit from the parallel processing power of ii data-parallel many-core architectures. In this dissertation, we present three techniques to adapt several non-data parallel applications in different dwarfs to modern many-core GPUs. First, we present a load balancing technique to maximize parallelism in non-serial polyadic Dynamic Programming (DP), which is a family of dynamic programming algorithms with more non-uniform data access pattern. We show that a bottom-up approach to solving the DP problem exploits more parallelism and therefore yields higher performance.
    [Show full text]
  • Massively Parallel Computers: Why Not Prirallel Computers for the Masses?
    Maslsively Parallel Computers: Why Not Prwallel Computers for the Masses? Gordon Bell Abstract In 1989 I described the situation in high performance computers including several parallel architectures that During the 1980s the computers engineering and could deliver teraflop power by 1995, but with no price science research community generally ignored parallel constraint. I felt SIMDs and multicomputers could processing. With a focus on high performance computing achieve this goal. A shared memory multiprocessor embodied in the massive 1990s High Performance looked infeasible then. Traditional, multiple vector Computing and Communications (HPCC) program that processor supercomputers such as Crays would simply not has the short-term, teraflop peak performanc~:goal using a evolve to a teraflop until 2000. Here's what happened. network of thousands of computers, everyone with even passing interest in parallelism is involted with the 1. During the first half of 1992, NEC's four processor massive parallelism "gold rush". Funding-wise the SX3 is the fastest computer, delivering 90% of its situation is bright; applications-wise massit e parallelism peak 22 Glops for the Linpeak benchmark, and Cray's is microscopic. While there are several programming 16 processor YMP C90 has the greatest throughput. models, the mainline is data parallel Fortran. However, new algorithms are required, negating th~: decades of 2. The SIMD hardware approach of Thinking Machines progress in algorithms. Thus, utility will no doubt be the was abandoned because it was only suitable for a few, Achilles Heal of massive parallelism. very large scale problems, barely multiprogrammed, and uneconomical for workloads. It's unclear whether The Teraflop: 1992 large SIMDs are "generation" scalable, and they are clearly not "size" scalable.
    [Show full text]
  • Introduction to Parallel Computing
    INTRODUCTION TO PARALLEL COMPUTING Plamen Krastev Office: 38 Oxford, Room 117 Email: [email protected] FAS Research Computing Harvard University OBJECTIVES: To introduce you to the basic concepts and ideas in parallel computing To familiarize you with the major programming models in parallel computing To provide you with with guidance for designing efficient parallel programs 2 OUTLINE: Introduction to Parallel Computing / High Performance Computing (HPC) Concepts and terminology Parallel programming models Parallelizing your programs Parallel examples 3 What is High Performance Computing? Pravetz 82 and 8M, Bulgarian Apple clones Image credit: flickr 4 What is High Performance Computing? Pravetz 82 and 8M, Bulgarian Apple clones Image credit: flickr 4 What is High Performance Computing? Odyssey supercomputer is the major computational resource of FAS RC: • 2,140 nodes / 60,000 cores • 14 petabytes of storage 5 What is High Performance Computing? Odyssey supercomputer is the major computational resource of FAS RC: • 2,140 nodes / 60,000 cores • 14 petabytes of storage Using the world’s fastest and largest computers to solve large and complex problems. 5 Serial Computation: Traditionally software has been written for serial computations: To be run on a single computer having a single Central Processing Unit (CPU) A problem is broken into a discrete set of instructions Instructions are executed one after another Only one instruction can be executed at any moment in time 6 Parallel Computing: In the simplest sense, parallel
    [Show full text]
  • Massively Parallel Processor Architectures for Resource-Aware Computing
    Massively Parallel Processor Architectures for Resource-aware Computing Vahid Lari, Alexandru Tanase, Frank Hannig, and J¨urgen Teich Hardware/Software Co-Design, Department of Computer Science Friedrich-Alexander University Erlangen-N¨urnberg (FAU), Germany {vahid.lari, alexandru-petru.tanase, hannig, teich}@cs.fau.de Abstract—We present a class of massively parallel processor severe, when considering massively parallel architectures, with architectures called invasive tightly coupled processor arrays hundreds to thousands of resources, all working with real-time (TCPAs). The presented processor class is a highly parame- requirements. terizable template, which can be tailored before runtime to fulfill costumers’ requirements such as performance, area cost, As a remedy, we present a domain-specific class of mas- and energy efficiency. These programmable accelerators are well sively parallel processor architectures called invasive tightly suited for domain-specific computing from the areas of signal, coupled processor arrays (TCPA), which offer built-in and image, and video processing as well as other streaming processing scalable resource management. The term “invasive” stems from applications. To overcome future scaling issues (e. g., power con- sumption, reliability, resource management, as well as application a novel paradigm called invasive computing [3], for designing parallelization and mapping), TCPAs are inherently designed and programming future massively parallel computing systems in a way to support self-adaptivity and resource awareness at (e.g., heterogeneous MPSoCs). The main idea and novelty of hardware level. Here, we follow a recently introduced resource- invasive computing is to introduce resource-aware program- aware parallel computing paradigm called invasive computing ming support in the sense that a given application gets the where an application can dynamically claim, execute, and release ability to explore and dynamically spread its computations resources.
    [Show full text]
  • Real-Time Network Traffic Simulation Methodology with a Massively Parallel Computing Architecture
    TRANSPORTATION RESEA RCH RECORD 1358 13 Advanced Traffic Management System: Real-Time Network Traffic Simulation Methodology with a Massively Parallel Computing Architecture THANAVAT }UNCHAYA, GANG-LEN CHANG, AND ALBERTO SANTIAGO The advent of parallel computing architectures presents an op­ of these ATMS developments would ensure optimal net­ portunity fo r lra nspon a tion professionals 1 imulate a large- ca le workwise performance. tra ffi c network with sufficienily fa t response time for real-time The anticipated benefits of these control methods depend opcrarion. Howeve r, ii neccssira tc ·. a fundament al change in the on the complex interactions among principal traffic system modeling algorithm LO tuke full ad va ntage of parallel computing. • uch a methodology t imulare tra ffic n twork with the Con­ components. These systems include driver behavior, level of nection Machine, a massively parallel computer, is described . The congestion, dynamic nature of traffic patterns, and the net­ basic parallel computing architectures are introdu ed, along with work's geometric configuration. It is crucial to the design of a list of commercially available parall el comput ers. This is fol ­ these strategies that a comprehensive understanding of the lowed by an in-depth presentation of the proposed simulation complex interrelations between these key system components meth odology with a massively parallel computer. The propo ed traffic simulation model ha. an inherent path-proces ing capa­ be established. Because it is often difficult for theoretical bilit y to represent drivers ' roure choice behavior at the individual­ formulations to take all such complexities into account, traffic vehicle level.
    [Show full text]
  • Massively Parallel Message Passing: Time to Start from Scratch?
    Massively Parallel Message Passing: Time to Start from Scratch? Anthony Skjellum, PhD University of Alabama at Birmingham Department of Computer and Information Sciences College of Arts and Sciences Birmingham, AL, USA [email protected] Abstract primitive abstractions do not do enough to abstract how operations are performed compared to what is to be performed. Extensions Yes: New requirements have invalidated the use of MPI and to and attempts to standardize with distributed shared memory (cf, other similar library-type message passing systems at full scale on MPI-2 one-sided primitives [10]) have been included, but are systems of 100's of millions of cores. Uses on substantial subsets currently somewhat incompatible with PGAS languages and sys- of these systems should remain viable, subject to limitations. tems [15], and furthermore MPI isn't easy to mutually convert No: If the right set of implementation and application deci- with BSP-style notations or semantics either [16]. Connection sions are made, and runtime strategies put in place, existing MPI with publish-subscribe models [17] used in military / space-based notations can used within exascale systems at least to a degree. fault-tolerant message passing systems are absent (see also [18]). Maybe: Implications of instantiating “MPI worlds” as islands The hypothesis of this paper is that new requirements arising inside exascale systems cannot damage overall scalability too in exascale computation have at least partly invalidated the use of much, or they will inhibit exascalability, and thus fail. the parallel model expressed and implied by MPI and other simi- lar library-type message passing systems at full scale on systems Categories and Subject Descriptors F.1.2.
    [Show full text]
  • An Event-Driven Massively Parallel Fine-Grained Processor Array
    An Event-Driven Massively Parallel Fine-Grained Processor Array Declan Walsh and Piotr Dudek School of Electrical and Electronic Engineering The University of Manchester, UK [email protected]; [email protected] Abstract—A multi-core event-driven parallel processor array connected to its four nearest direct neighbours and only local design is presented. Using relatively simple 8-bit processing cores communications are performed. This network topology and a 2D mesh network topology, the architecture focuses on reduces the long wire interconnects which also reduces power reducing the area occupation of a single processor core. A large and enables operations at higher frequencies. number of these processor cores can be implemented on a single The SIMD implementations are useful for performing tasks integrated chip to create a MIMD architecture capable of providing a powerful processing performance. Each processor which require data-level parallelism. However, for tasks which core is an event-driven processor which can enter an idle mode require more complex processing than simple pixel-parallel when no data is changing locally. An 8 × 8 prototype processor operators, e.g., event-driven networks [12], message-passing array is implemented in a 65 nm CMOS process in 1,875 µm × inference networks [13], or even some low-level image 1,875 µm. This processor array is capable of performing 5.12 processing operations [14, 15], SIMD processing may not be GOPS operating at 80 MHz with an average power consumption optimal and a more flexible MIMD (Multiple Instruction of 75.4 mW. Multiple Data) solution is more suitable.
    [Show full text]
  • A Massively Parallel Digital Learning Processor
    A Massively Parallel Digital Learning Processor Hans Peter Graf Srihari Cadambi Igor Durdanovic [email protected] [email protected] [email protected] Venkata Jakkula Murugan Sankardadass Eric Cosatto Srimat Chakradhar [email protected] [email protected] [email protected] [email protected] NEC Laboratories, America 4 Independence Way, Suite 200; Princeton, NJ 07738, USA Abstract We present a new, massively parallel architecture for accelerating machine learning algorithms, based on arrays of vector processing elements (VPEs) with variable-resolution arithmetic. Groups of VPEs operate in SIMD (single instruction multiple data) mode, and each group is connected to an independent memory bank. The memory bandwidth thus scales with the number of VPEs, while the main data flows are local, keeping power dissipation low. With 256 VPEs, implemented on two FPGAs (field programmable gate array) chips, we obtain a sustained speed of 19 GMACS (billion multiply- accumulate per sec.) for SVM training, and 86 GMACS for SVM classification. This performance is more than an order of magnitude higher than that of any FPGA implementation reported so far. The speed on one FPGA is similar to the fastest speeds published on a Graphics Processor for the MNIST problem, despite a clock rate that is an order of magnitude lower. Tests with Convolutional Neural Networks show similar compute performances. This massively parallel architecture is particularly attractive for embedded applications, where low power dissipation is critical. 1 Introduction Machine learning demands higher and higher compute-performance, but serial processors are not improving that much anymore - at least not as quickly as they used to.
    [Show full text]
  • Implementation of a Message Passing Interface Into a Cloud-Resolving Model for Massively Parallel Computing
    Implementation of a Message Passing Interface into a Cloud-Resolving Model for Massively Parallel Computing Hann-Ming Henry Juang’, Wei-Kuo Tao’, Xiping Zen$’3, Chung-Lin Shie23, Joanne Simpson2and Steve hg4, ’EnvironmentalModeling Center NCEP, NOR4 Washington, DC 2Laboratolyfor Atmospheres NASNGoddard Space Flight Center Greenbelt, MLI 20771, USA 3GoddardEarth Sciences and Technology Center University of Maryland, Baltimore County Baltimore, MD 4Science Systems and Applications Inc. Lanham, MD 20706 June 2,2004 Mon. Wea. Rev. 1 Corresponding author address: Dr. Ham-Ming Henry Juang, Environmental Modeling Center, NCEP, NOAA, W/NP5 Room 204,5200 Auth Road, Camp Springs, MD 20746 Abstract The capability for massively parallel programming (MPP) using a message passing interface (MPI) has been implemented into a three-dimensional version of the Goddard Cumulus Ensemble (GCE) model. The design for the MPP with MPI uses the concept of maintaining similar code structure between the whole domain as well as the portions after decomposition. Hence the model follows the same integration for single and multiple tasks (CPUs). Also, it provides for minimal changes to the on@ code, so it is easily modified andor managed by the model developers and users who have little knowledge of MPP. The entire model domain could be sliced into one- or two-dimensional decomposition with a halo regime, which is overlaid on partial domains. The halo regime requires that no data be fetched across tasks during the computational stage, but it must be updated before the next computational stage through data exchange via MPI. For reproducible purposes, transposing data among tasks is required for spectral transform (Fast Fourier Transform, FFT), which is used in the anelastic version of the model for solving the pressure equation.
    [Show full text]