devito
devito copied to clipboard
DSL and compiler framework for automated finite-differences and stencil computation
The following MFE fails with `DEVITO_MPI='full'`, even on a single rank. MFE: ```python from devito import Grid, SubDomain, Function, Eq, Operator class SD0(SubDomain): name = 'sd0' def define(self, dimensions): return...
Issues with the operator construction time for multi-stage methods. The ideia of the multi-stage methods implemented in Devito is to write a PDE system as first order in time: where...
See https://github.com/devitocodes/devito/blob/v4.8.11/tests/test_symbolics.py#L257-L260 for example of use. Spitting out the generated code we can see that in the case where the constant value is an `np.float64` sympy eagerly turns this into...
Undeclared var, possibly due to wrong scheduling (Initial thread from @felipeguedessilva, here: https://devitocodes.slack.com/archives/C7JMLMSG0/p1736908013546609) mfe.py ```python from devito import Grid, TimeFunction, Function, Eq, Operator # Setup Configuration grid = Grid(shape=(30, 30))...
MFE: ```python from devito.types.basic import DataSymbol s0 = DataSymbol(name='s') s1 = DataSymbol(name='s') s2 = s0._rebuild(name='s2') s3 = s1._rebuild(name='s3') # Causes s1 to be incorrectly substituted # mapper = {s0: s2,...
`Operator._prepare_arguments` currently contains the following lines: ```python # Process Thicknesses for p in self.parameters: if isinstance(p, Thickness): args.update(p._arg_values(grid=grid, **kwargs)) ``` which could be removed by making `Thickness` have `is_Input =...
Updates the requirements on [dask-cuda](https://github.com/rapidsai/dask-cuda) to permit the latest version. Changelog Sourced from dask-cuda's changelog. dask-cuda 25.12.00 (10 Dec 2025) 🐛 Bug Fixes Print test_all_to_all debug info when assertion fails...