devito
devito copied to clipboard
devitocodes
DSL and compiler framework for automated finite-differences and stencil computation
When using a conditional dimension that contains `Ge` / `Le` types of conditions and a saved buffer defined over that conditional dimension, there is a conflict in the calculation of...
Symbols generated within Devito currently do not have guaranteed unique names, potentially producing inadvertent clashes which result in cryptic errors being raised. Devito currently makes use of counters (among other...
hi FabioLuporini, are there any examples of forward modelling using a 3D velocity model ?
I'm trying to solve the wave equation in $\mathbb{R}^2$ and use the abc methods which add a damping term to the original equation. The amplitude of the damping term is...
MFE ``` import devito shape = (12,13,14) nfreq = 5 nrec = 6 space_order = 8 grid = devito.Grid(shape=shape) freq_dim = devito.DefaultDimension(name="freq_dim", default_value=nfreq) p0re = devito.Function(name="p0re", grid=grid, dimensions=(*grid.dimensions, freq_dim), shape=(*shape,...
Cross-rank assignments fail, while expected to work ```python from devito import Grid, TimeFunction, Eq, Operator grid = Grid(shape=(4, 4)) u = TimeFunction(name="u", grid=grid, space_order=2) u.data[0, 1:-1, 1:-1] = 1. u.data[0,...
Reproducer ``` import numpy as np from devito import Eq, Grid, TimeFunction, Operator grid = Grid(shape=(4, 4)) u = TimeFunction(name='u', grid=grid) usave = TimeFunction(name='usave', grid=grid, save=5) eqns = [Eq(u.forward, u...
Reproducer: ``` import numpy as np from devito import * from devito.types import Symbol configuration['language'] = 'openmp' grid = Grid(shape=(10,)) x, = grid.dimensions f = Function(name='f', grid=grid, space_order=0) s =...
Metadata
Owner
Metadata
DSL and compiler framework for automated finite-differences and stencil computation