Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDevito (v3.1.0): an embedded domain-specific language for finite differences and geophysical exploration
153
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We introduce Devito, a new domain-specific language for implementing high-performance finite difference partial differential equation solvers. The motivating application is exploration seismology where methods such as Full-Waveform Inversion and Reverse-Time Migration are used to invert terabytes of seismic data to create images of the earths subsurface. Even using modern supercomputers, it can take weeks to process a single seismic survey and create a useful subsurface image. The computational cost is dominated by the numerical solution of wave equations and their corresponding adjoints. Therefore, a great deal of effort is invested in aggressively optimizing the performance of these wave-equation propagators for different computer architectures. Additionally, the actual set of partial differential equations being solved and their numerical discretization is under constant innovation as increasingly realistic representations of the physics are developed, further ratcheting up the cost of practical solvers. By embedding a domain-specific language within Python and making heavy use of SymPy, a symbolic mathematics library, we make it possible to develop finite difference simulators quickly using a syntax that strongly resembles the mathematics. The Devito compiler reads this code and applies a wide range of analysis to generate highly optimized and parallel code. This approach can reduce the development time of a verified and optimized solver from months to days.
rate research
Read More
The SPaCIFY project, which aims at bringing advances in MDE to the satellite flight software industry, advocates a top-down approach built on a domain-specific modeling language named Synoptic. In line with previous approaches to real-time modeling such as Statecharts and Simulink, Synoptic features hierarchical decomposition of application and control modules in synchronous block diagrams and state machines. Its semantics is described in the polychronous model of computation, which is that of the synchronous language Signal.
Computer simulations are used widely across the engineering and science disciplines, including in the research and development of magnetic devices using computational micromagnetics. In this work, we identify and review different approaches to configuring simulation runs: (i) the re-compilation of source code, (ii) the use of configuration files, (iii) the graphical user interface, and (iv) embedding the simulation specification in an existing programming language to express the computational problem. We identify the advantages and disadvantages of different approaches and discuss their implications on effectiveness and reproducibility of computational studies and results. Following on from this, we design and describe a domain specific language for micromagnetics that is embedded in the Python language, and allows users to define the micromagnetic simulations they want to carry out in a flexible way. We have implemented this micromagnetic simulation description language together with a computational backend that executes the simulation task using the Object Oriented MicroMagnetic Framework (OOMMF). We illustrate the use of this Python interface for OOMMF by solving the micromagnetic standard problem 4. All the code is publicly available and is open source.
Orbital geophysical investigations of Enceladus are critical to understanding its energy balance. We identified key science questions for the geophysical exploration of Enceladus, answering which would support future assessment of Enceladus astrobiological potential. Using a Bayesian framework, we explored how science requirements map to measurement requirements. We performed mission simulations to study the sensitivity of a single spacecraft and dual spacecraft configurations to static gravity and tidal Love numbers of Enceladus. We find that mapping Enceladus gravity field, improving the accuracy of the physical libration amplitude, and measuring Enceladus tidal response would provide critical constraints on the internal structure, and establish a framework for assessing Enceladus long-term habitability. This kind of investigation could be carried out as part of a life search mission at little additional resource requirements.
The ANTAREX project relies on a Domain Specific Language (DSL) based on Aspect Oriented Programming (AOP) concepts to allow applications to enforce extra functional properties such as energy-efficiency and performance and to optimize Quality of Service (QoS) in an adaptive way. The DSL approach allows the definition of energy-efficiency, performance, and adaptivity strategies as well as their enforcement at runtime through application autotuning and resource and power management. In this paper, we present an overview of the key outcome of the project, the ANTAREX DSL, and some of its capabilities through a number of examples, including how the DSL is applied in the context of the project use cases.
After reviewing the problematic behavior of some previously suggested finite interval spatial operators of the symmetric Riesz type, we create a wish list leading toward a new spatial operator suitable to use in the space-time fractional differential equation of anomalous diffusion when the transport of material is strictly restricted to a bounded domain. Based on recent studies of wall effects, we introduce a new definition of the spatial operator and illustrate its favorable characteristics. We provide two numerical methods to solve the modified space-time fractional differential equation and show particular results illustrating compliance to our established list of requirements, most important to the conservation principle and the second law of thermodynamics.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
