SEARCH · Search NASA

Results for “Parallel in Time”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,621 records · Page 90

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Dynamic Phasor Modeling of Multi-Converter Systems

Generalized form of Dynamic Phasor (DP)-based modeling of multi-converter systems containing high order harmonics has not been proposed in the literature yet due to the complexity and large number of variables that should be included in the model. In this work a generalized form of a system-level model comprising high number of single-phase converters is presented. The proposed method can be applied to any number of single-phase voltage source converters (VSIs). Test cases of this modeling method containing 5 to 100 parallel connected single-phase voltage source inverters are modeled and simulation times are recorded. The results of the proposed method are compared and validated with conventional average models as well as detailed switching models. A systematic approach for comparing the accuracy and timestep between dynamic phasor modeling method and detailed switching model is illustrated. Advantage of DP models over conventional average models for stability assessment are also discussed at the end of the paper.

Xue, Yaosuo↗

DOE OSTI2012

Development of flow network analysis code for block type VHTR core by linear theory method

VHTR (Very High Temperature Reactor) is high-efficiency nuclear reactor which is capable of generating hydrogen with high temperature of coolant. PMR (Prismatic Modular Reactor) type reactor consists of hexagonal prismatic fuel blocks and reflector blocks. The flow paths in the prismatic VHTR core consist of coolant holes, bypass gaps and cross gaps. Complicated flow paths are formed in the core since the coolant holes and bypass gap are connected by the cross gap. Distributed coolant was mixed in the core through the cross gap so that the flow characteristics could not be modeled as a simple parallel pipe system. It requires lot of effort and takes very long time to analyze the core flow with CFD analysis. Hence, it is important to develop the code for VHTR core flow which can predict the core flow distribution fast and accurate. In this study, steady state flow network analysis code is developed using flow network algorithm. Developed flow network analysis code was named as FLASH code and it was validated with the experimental data and CFD simulation results. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2017

Savannah River Site HB-Line Plutonium Oxide Production Process Simulation Model - 17458

The AFS II/Plutonium Oxide Production process is in its initial operational phase at the Savannah River Site's (SRS) HB-Line. The process is new for the facility where it is being performed, but the facility itself has been in use on and off for over 50 years. Beginning this process inside the facility has resulted in multiple process delays, due to both equipment issues and process restrictions. In order to better understand the process and to identify ways to improve facility performance, Savannah River Nuclear Solutions (SRNS), in conjuncture with Newport News Shipbuilding (NNS), developed the HB-Line Process Simulation Model. This tool completely replicates the operations that comprise the HB-Line process and simulates all of the constraints associated with the process. The model realistically accounts for the number of personnel required for a process step, the time required for a process step to occur, rules related to the parallel operation of process steps, and realistic equipment and operational delay data. The result is a tool that allows SRNS to test different scenarios for improving the process throughput, while realistically accounting for potential process delays and restrictions. The tool is currently in use and has allowed SRNS to present to their customer, the National Nuclear Security Administration (NNSA), estimated production capabilities over wide-ranging timetables and assumptions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2018

Distinguishing Research Reactor Fuels Using Neutron Multiplicity - 25615

Neutron measurements are one of the reliable means of interrogation of fissile material. Neutron measurements are used for total counting (pure metal fissile samples) or counting coincidence/multiplicity (for impure fissile samples). Various flavors of neutron measurement are used routinely by the IAEA to determine the fissile mass of an unknown sample. The objective of this work is to design an experiment to distinguish the two different fuel types using neutron multiplicity method. The design of experiment is performed using Monte Carlo N- Particle (MCNP) transport code. General Electric (GE) designed 1-MWt pool-type research reactor at UMass Lowell (UMLRR) uses Materials Testing Reactor (MTR) plate-type fuel consisting of 19 full and 2 partial fuel assemblies, controlled by four boron control blades and one regulating control rod. Two different fuel types are available at the UMLRR facility. In 2008, UMass Lowell obtained slightly burnt fuel assemblies from the decommissioned WPI research reactor. The simulation results discussed in this paper show that within the limits of statistical error, the single and double moments for the WPI and UMLRR fuels can be measured. The two fuels have very similar composition and enrichments resulting in very similar neutron multiplicity moments. This implies that neutron multiplicity measurements alone would not be sufficient to distinguish between the two fuel types at low burnups. The change in the observed single-double ratio between each fuel type is less than 6 % which falls below the uncertainty in the efficiency of the JCC71 slab detectors. The neutron signature from low burnup and low enriched uranium fuel are not strong, however, these can be measured under low background and long measurement times. This would be an important consideration in designing the experiments and parallel gamma measurements will also be planned. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Enabling real-time adaptation of machine learning models at x-ray Free Electron Laser facilities with high-speed training optimized computational hardware

The emergence of novel computational hardware is enabling a new paradigm for rapid machine learning model training. For the Department of Energy’s major research facilities, this developing technology will enable a highly adaptive approach to experimental sciences. In this manuscript we present the per-epoch and end-to-end training times for an example of a streaming diagnostic that is planned for the upcoming high-repetition rate x-ray Free Electron Laser, the Linac Coherent Light Source-II. We explore the parameter space of batch size and data parallel training across multiple Graphics Processing Units and Reconfigurable Dataflow Units. We show the landscape of training times with a goal of full model retraining in under 15 min. Although a full from scratch retraining of a model may not be required in all cases, we nevertheless present an example of the application of emerging computational hardware for adapting machine learning models to changing environments in real-time, during streaming data acquisition, at the rates expected for the data fire hoses of accelerator-based user facilities.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1993

Parametric instabilities of a large-amplitude circularly polarized Alfven wave - Linear growth in two-dimensional geometries

The growth of parametric instabilities, which may lead to the development of a turbulent cascade, is studied using an MHD code that permits nonlinear couplings in the parallel direction to the ambient magnetic field and one perpendicular direction. Compressibility is included in the form of a polytropic equation of state. Parametric instabilities associated with a parallel-propagating decay instability are found to dominate the low-beta case. An obliquely propagating filamentationlike instability dominates the high-beta case. The nonlinear growth of the nth harmonic of a daughter wave growing as a factor of n times the fundamental's growth rate is found in both cases. Nonlinear saturation is caused by the parallel decay instability in the low-beta case and by the oblique filamentationlike instability in the high-beta case.

Ghosh, S.↗

DOE OSTI2019

Software for Sparse Tensor Decomposition on Emerging Computing Architectures

In this paper, we develop software for decomposing sparse tensors that is portable to and performant on a variety of multicore, manycore, and GPU computing architectures. The result is a single code whose performance matches optimized architecture-specific implementations. The key to a portable approach is to determine multiple levels of parallelism that can be mapped in different ways to different architectures, and we explain how to do this for the matricized tensor times Khatri--Rao product (MTTKRP), which is the key kernel in canonical polyadic tensor decomposition. Our implementation leverages the Kokkos framework, which enables a single code to achieve high performance across multiple architectures that differ in how they approach fine-grained parallelism. We also introduce a new construct for portable thread-local arrays, which we call compile-time polymorphic arrays. Not only are the specifics of our approaches and implementation interesting for tuning tensor computations, but they also provide a roadmap for developing other portable high-performance codes. As a last step in optimizing performance, we modify the MTTKRP algorithm itself to do a permuted traversal of tensor nonzeros to reduce atomic-write contention. Lastly, we test the performance of our implementation on 16- and 68-core Intel CPUs and the K80 and P100 NVIDIA GPUs, showing that we are competitive with state-of-the-art architecture-specific codes while having the advantage of being able to run on a variety of architectures.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

Real-time data system: Incorporating new technology in mission critical environments

If the Space Station Freedom is to remain viable over its 30-year life span, it must be able to incorporate new information systems technologies. These technologies are necessary to enhance mission effectiveness and to enable new NASA missions, such as supporting the Lunar-Mars Initiative. Hi-definition television (HDTV), neural nets, model-based reasoning, advanced languages, CPU designs, and computer networking standards are areas which have been forecasted to make major strides in the next 30 years. A major challenge to NASA is to bring these technologies online without compromising mission safety. In past programs, NASA managers have been understandably reluctant to rely on new technologies for mission critical activities until they are proven in noncritical areas. NASA must develop strategies to allow inflight confidence building and migration of technologies into the trusted tool base. NASA has successfully met this challenge and developed a winning strategy in the Space Shuttle Mission Control Center. This facility, which is clearly among NASA's most critical, is based on 1970's mainframe architecture. Changes to the mainframe are very expensive due to the extensive testing required to prove that changes do not have unanticipated impact on critical processes. Systematic improvement efforts in this facility have been delayed due to this 'risk to change.' In the real-time data system (RTDS) we have introduced a network of engineering computer workstations which run in parallel to the mainframe system. These workstations are located next to flight controller operating positions in mission control and, in some cases, the display units are mounted in the traditional mainframe consoles. This system incorporates several major improvements over the mainframe consoles including automated fault detection by real-time expert systems and color graphic animated schematics of subsystems driven by real-time telemetry. The workstations have the capability of recording telemetry data and providing 'instant replay' for flight controllers. RTDS also provides unique graphics animated by real-time telemetry such as workstation emulation of the shuttle's flight instruments and displays of the remote manipulator system (RMS) position. These systems have been used successfully as prime operational tools since STS-26 and have supported seven shuttle missions.

Muratore, John F.↗

DOE OSTI2021

Transient cycling of nitrogen, organic carbon and oxygen within the free-flowing Columbia River corridor: Linking exposure time dependent biogeochemical reactions to river stage fluctuations (Final project report)

The objective of this project was to develop next-generation techniques for representing the transformations of complex reaction networks like those impacting transient river-corridors, then apply them to develop upscaling tools. The project contributed significantly to the development of Lagrangian “mass transfer particle tracking” (MTPT) tools that are the first numerical framework that explicitly separates mixing and spreading processes; a crucial distinction for accurate representations of reaction rates across scales. The MTPT approach uses an operator splitting scheme where physical transport processes are simulated using classical random walk methods and each particle is treated as a “container” that carries with it any number of chemical species. Mixing between particles is simulated using a colocation-probability based mass transfer kernel and reactions are evaluated on each particle after mixing. The method is stable, accurate, and also allows the explicit modeling of residence time distributions of the individual containers of mass. The MTPT scheme allowed us to explore several research questions related to reactive transport in river corridors and complex reaction networks. We found that it is possible to quantify the exposure time distributions (ETDs) of reactants and that this does suggest some pathways to upscaling. However, comparable approximations can be obtained using simplified (reduced dimensionality) MTPT simulations in less time with similar uncertainty, so development of extensive ETD-based methods was not productive. The major findings are that i) the MTPT schemes are robust and highly accurate across a range of arbitrarily complex reactions, ii) the methods can be efficiently parallelized and the parallel performance characteristics of MTPT are predictable, and iii) the MTPT tools allow simultaneous tracking of residence time. The development of these research tools into efficient software packages continues. The methods are already available to the community because open-source, working examples have been included with all publications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Multipinhole collimator with 20 apertures for a brain SPECT application

Purpose: Several new technologies for single photon emission computed tomography (SPECT) instrumentation with parallel-hole collimation have been proposed to improve detector sensitivity and signal collection efficiency. Benefits from improved signal efficiency include shorter acquisition times and lower dose requirements. In this paper, the authors show a possibility of over an order of magnitude enhancement in photon detection efficiency (from 7.6 × 10{sup −5} to 1.6 × 10{sup −3}) for dopamine transporter (DaT) imaging of the striatum over the conventional SPECT parallel-hole collimators by use of custom-designed 20 multipinhole (20-MPH) collimators with apertures of 0.75 cm diameter. Methods: Quantifying specific binding ratio (SBR) of {sup 123}I-ioflupane or {sup 123}I-iometopane’s signal at the striatal region is a common brain imaging method to confirm the diagnosis of the Parkinson’s disease. The authors performed imaging of a striatal phantom filled with aqueous solution of I-123 and compared camera recovery ratios of SBR acquired between low-energy high-resolution (LEHR) parallel-hole collimators and 20-MPH collimators. Results: With only two-thirds of total acquisition time (20 min against 30 min), a comparable camera recovery ratio of SBR was achieved using 20-MPH collimators in comparison to that from the LEHR collimator study. Conclusions: Their systematic analyses showed that the 20-MPH collimator could be a promising alternative for the DaT SPECT imaging for brain over the traditional LEHR collimator, which could give both shorter scan time and improved diagnostic accuracy.

60 APPLIED LIFE SCIENCES↗

NASA NTRS1983

Parallel solution of finite element equations

The paper examines several parallel processing solution algorithms for finite element equations arising in linear equilibrium problems. Two basic groups of algorithms, direct and iterative, are investigated with respect to a number of parallel computer architectures and associated selection criteria. The direct algorithms include: LR-Gauss, Crout, Cholesky, Cyclic Reduction and WZ-factorization. The iterative methods examined are: Accelerated Gauss-Seidel, Surrogate Stiffness, Jacobi, Series Expansion, and Energy Monte Carlo. For real-time applications, where the object is to minimize the execution time, Cyclic Reduction appears to be best suited. This assumes a computer with an unlimited number of parallel processors. However, for computers with a limited number of parallel processors that must be used efficiently, both Gauss factorization and Jacobi-like iterative methods rank favorably.

Salama, M.↗

NASA NTRS1997

HARP: A Dynamic Inertial Spectral Partitioner

Partitioning unstructured graphs is central to the parallel solution of computational science and engineering problems. Spectral partitioners, such recursive spectral bisection (RSB), have proven effecfive in generating high-quality partitions of realistically-sized meshes. The major problem which hindered their wide-spread use was their long execution times. This paper presents a new inertial spectral partitioner, called HARP. The main objective of the proposed approach is to quickly partition the meshes at runtime in a manner that works efficiently for real applications in the context of distributed-memory machines. The underlying principle of HARP is to find the eigenvectors of the unpartitioned vertices and then project them onto the eigerivectors of the original mesh. Results for various meshes ranging in size from 1000 to 100,000 vertices indicate that HARP can indeed partition meshes rapidly at runtime. Experimental results show that our largest mesh can be partitioned sequentially in only a few seconds on an SP2 which is several times faster than other spectral partitioners while maintaining the solution quality of the proven RSB method. A parallel WI version of HARP has also been implemented on IBM SP2 and Cray T3E. Parallel HARP, running on 64 processors SP2 and T3E, can partition a mesh containing more than 100,000 vertices into 64 subgrids in about half a second. These results indicate that graph partitioning can now be truly embedded in dynamically-changing real-world applications.

Simon, Horst D.↗

NASA NTRS2006

Implementing Access to Data Distributed on Many Processors

A reference architecture is defined for an object-oriented implementation of domains, arrays, and distributions written in the programming language Chapel. This technology primarily addresses domains that contain arrays that have regular index sets with the low-level implementation details being beyond the scope of this discussion. What is defined is a complete set of object-oriented operators that allows one to perform data distributions for domain arrays involving regular arithmetic index sets. What is unique is that these operators allow for the arbitrary regions of the arrays to be fragmented and distributed across multiple processors with a single point of access giving the programmer the illusion that all the elements are collocated on a single processor. Today's massively parallel High Productivity Computing Systems (HPCS) are characterized by a modular structure, with a large number of processing and memory units connected by a high-speed network. Locality of access as well as load balancing are primary concerns in these systems that are typically used for high-performance scientific computation. Data distributions address these issues by providing a range of methods for spreading large data sets across the components of a system. Over the past two decades, many languages, systems, tools, and libraries have been developed for the support of distributions. Since the performance of data parallel applications is directly influenced by the distribution strategy, users often resort to low-level programming models that allow fine-tuning of the distribution aspects affecting performance, but, at the same time, are tedious and error-prone. This technology presents a reusable design of a data-distribution framework for data parallel high-performance applications. Distributions are a means to express locality in systems composed of large numbers of processor and memory components connected by a network. Since distributions have a great effect on the performance of applications, it is important that the distribution strategy is flexible, so its behavior can change depending on the needs of the application. At the same time, high productivity concerns require that the user be shielded from error-prone, tedious details such as communication and synchronization.

James, Mark↗

NASA NTRS1993

Operations automation using the Link Monitor and Control Operator Assistant

The Link Monitor and Control Operator Assistant (LMC OA) is a knowledge-based prototype system which uses AI techniques to provide semiautomated monitor and control functions to support operations of the Deep Space Network (DSN) 70-m antenna at the Goldstone Deep Space Communications Complex (DSCC). The manual and time-consuming process of configuring the 70-m antenna and its associated communications and processing equipment, known as precalibration, is an overhead activity; the time spent in precalibration is time which cannot be spent supporting actual mission operations. Therefore, the major goal of the LMC OA task is to demonstrate techniques that reduce precalibration time, decrease operations overhead, and increase the availability of this valuable and oversubscribed NASA resource. The LMC OA prototype was tested in a parallel, experimental mode at the Goldstone DSCC performing semiautomated precalibration using the actual operational equipment. This test demonstrated that a reduction of 40 percent in precalibration time can be achieved with the LMC OA prototype.

Lee, Lorrine F.↗

NASA NTRS2002

Maximizing TDRS Command Load Lifetime

The GNC software onboard ISS utilizes TORS command loads, and a simplistic model of TORS orbital motion to generate onboard TORS state vectors. Each TORS command load contains five "invariant" orbital elements which serve as inputs to the onboard propagation algorithm. These elements include semi-major axis, inclination, time of last ascending node crossing, right ascension of ascending node, and mean motion. Running parallel to the onboard software is the TORS Command Builder Tool application, located in the JSC Mission Control Center. The TORS Command Builder Tool is responsible for building the TORS command loads using a ground TORS state vector, mirroring the onboard propagation algorithm, and assessing the fidelity of current TORS command loads onboard ISS. The tool works by extracting a ground state vector at a given time from a current TORS ephemeris, and then calculating the corresponding "onboard" TORS state vector at the same time using the current onboard TORS command load. The tool then performs a comparison between these two vectors and displays the relative differences in the command builder tool GUI. If the RSS position difference between these two vectors exceeds the tolerable lim its, a new command load is built using the ground state vector and uplinked to ISS. A command load's lifetime is therefore defined as the time from when a command load is built to the time the RSS position difference exceeds the tolerable limit. From the outset of TORS command load operations (STS-98), command load lifetime was limited to approximately one week due to the simplicity of both the onboard propagation algorithm, and the algorithm used by the command builder tool to generate the invariant orbital elements. It was soon desired to extend command load lifetime in order to minimize potential risk due to frequent ISS commanding. Initial studies indicated that command load lifetime was most sensitive to changes in mean motion. Finding a suitable value for mean motion was therefore the key to achieving this goal. This goal was eventually realized through development of an Excel spreadsheet tool called EMMIE (Excel Mean Motion Interactive Estimation). EMMIE utilizes ground ephemeris nodal data to perform a least-squares fit to inferred mean anomaly as a function of time, thus generating an initial estimate for mean motion. This mean motion in turn drives a plot of estimated downtrack position difference versus time. The user can then manually iterate the mean motion, and determine an optimal value that will maximize command load lifetime. Once this optimal value is determined, the mean motion initially calculated by the command builder tool is overwritten with the new optimal value, and the command load is built for uplink to ISS. EMMIE also provides the capability for command load lifetime to be tracked through multiple TORS ephemeris updates. Using EMMIE, TORS command load lifetimes of approximately 30 days have been achieved.

Brown, Aaron J.↗

DOE OSTI2015

Measurement of Scintillation and Ionization Yield and Scintillation Pulse Shape from Nuclear Recoils in Liquid Argon

We have measured the scintillation and ionization yield of recoiling nuclei in liquid argon as a function of applied electric field by exposing a dual-phase liquid argon time projection chamber (LAr-TPC) to a low energy pulsed narrow band neutron beam produced at the Notre Dame Institute for Structure and Nuclear Astrophysics. Liquid scintillation counters were arranged to detect and identify neutrons scattered in the TPC and to select the energy of the recoiling nuclei. We also report measurements of the scintillation yields for nuclear recoils with energies from 10.3 to 57.3 keV and for median applied electric fields from 0 to 970 V/cm. For the ionization yields, we report measurements from 16.9 to 57.3 keV and for electric fields from 96.4 to 486 V/cm. Furthermore, we report the observation of an anticorrelation between scintillation and ionization from nuclear recoils, which is similar to the anticorrelation between scintillation and ionization from electron recoils. Assuming that the energy loss partitions into excitons and ion pairs from 83m Kr internal conversion electrons is comparable to that from 207 Bi conversion electrons, we obtained the numbers of excitons (N ex ) and ion pairs (N i ) and their ratio (N ex /N i ) produced by nuclear recoils from 16.9 to 57.3 keV. Motivated by arguments suggesting direction sensitivity in LAr-TPC signals due to columnar recombination, a comparison of the light and charge yield of recoils parallel and perpendicular to the applied electric field is presented for the first time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2016

THE IBEX RIBBON AND THE PICKUP ION RING STABILITY IN THE OUTER HELIOSHEATH. II. MONTE-CARLO AND PARTICLE-IN-CELL MODEL RESULTS

The nearly circular ribbon of energetic neutral atom (ENA) emission discovered by NASA’s Interplanetary Boundary EXplorer satellite ( IBEX ), is most commonly attributed to the effect of charge exchange of secondary pickup ions (PUIs) gyrating about the magnetic field in the outer heliosheath (OHS) and the interstellar space beyond. The first paper in the series (Paper I) presented a theoretical analysis of the pickup process in the OHS and hybrid-kinetic simulations, revealing that the kinetic properties of freshly injected proton rings depend sensitively on the details of their velocity distribution. It was demonstrated that only rings that are not too narrow (parallel thermal spread above a few km s{sup −1}) and not too wide (parallel temperature smaller than the core plasma temperature) could remain stable for a period of time long enough to generate ribbon ENAs. This paper investigates the role of electron dynamics and the extra spatial degree of freedom in the ring ion scattering process with the help of two-dimensional full particle-in-cell (PIC) kinetic simulations. A good agreement is observed between ring evolution under unstable conditions in hybrid and PIC models, and the dominant modes are found to propagate parallel to the magnetic field. We also present more realistic ribbon PUI distributions generated using Monte Carlo simulations of atomic hydrogen in the global heliosphere and examine the effect of both the cold ring-like and the hot “halo” PUIs produced from heliosheath ENAs on the ring stability. It is shown that the second PUI population enhances the fluctuation growth rate, leading to faster isotropization of the solar-wind-derived ring ions.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2017

Numerical Simulation of Illumination and Thermal Conditions at the Lunar Poles Using LOLA DTMs

We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.

Glaser, P.↗