devito icon indicating copy to clipboard operation
devito copied to clipboard

DSL and compiler framework for automated finite-differences and stencil computation

Results 157 devito issues
Sort by recently updated
recently updated
newest added

Updated this notebook with a better free surface implementation in the vein of Mathias' implementation in the acoustic examples.

examples

Hi guys, This PR contains the implementations of elastic adjoint equations. The implementation is based on Feng and Schuster (2017): Elastic least-squares reverse time migration (https://doi.org/10.1190/geo2016-0254.1). I hope that you...

examples

MPI `data_gather` broken for `TimeFunction`. MFE below ``` # env DEVITO_LOGGING=DEBUG DEVITO_MPI=1 mpirun -n 2 python mfe_gather2.py import numpy as np from devito import Grid, TimeFunction, Eq, Operator, MPI grid...

MPI

If you try to use the absolute value of an index (for example when implementing a free-surface as in the acoustic examples) with a function that is staggered, the staggering...

Hello everybody. This PR tries to suggest reconstructing the source wave field from the one saved in the border. This should be an option for the gradient computation, in which...

examples

I have been using Devito for a while and found myself using code set up in examples/ for multiple reasons. I propose this issue so as to aid people in...

I made a notebook with a simple implementation of a toy variable-z curvilinear acoustic model with a cosine hill. Still WIP, needs tidy-up and explanations, but lmk your thoughts

examples

This should improve runtime performance of TTI staggered. The basic infrastructure is already in the codebase (used for example by the CIRE algorithm), but it currently doesn't support the case...

Opening for discussion. This helps to work around a bug when injecting to more complex fields such as f(t, u(x,y,z)). Helps to implement Listing 3 of https://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=9460483

compiler

I propose this issue as a venue to report error messages that weren't helpful in the past, you or that they didn't fit your usecase. This will most likely spawn,...