rhodrin
rhodrin
The problem appears to be within the `_setitem_` method when setting part/all of one Devito `Function`'s MPI distributed `Data` to another `Function`'s MPI distributed `Data`. Certain slicing short hands do...
Fix is under way here: https://github.com/devitocodes/devito/tree/mpi_slicing_fix In this branch, the above MFE is fixed but below is a new failing example with `-np 4`: ``` import numpy as np from...
Just a follow up to this owing to the recent chat in #qestions. So below is code for a lid driven cavity problem: ``` import numpy as np from devito...
Note: For the above script to now work with `master` the following changes must be made: ---------------------------------- ``` eq_u = Eq(u.dt + u*u.dx + v*u.dy - nu*u.laplace, -1./rho*p.dx) eq_v =...
Note: There's a also a native solution to the [lid-driven cavity](https://github.com/devitocodes/devito/blob/break_devito/examples/cfd/lid_driven_cavity.py) problem on branch [break_devito](https://github.com/devitocodes/devito/tree/break_devito) c-code produced is: ``` #define _POSIX_C_SOURCE 200809L #include "stdlib.h" #include "math.h" #include "sys/time.h" #include "xmmintrin.h"...
`Data` objects are not designed to be 'stripped' away from their `Function` in this manner. Is there any reason why: ``` from devito import Grid, Function import pickle g =...
Yeah, I see what you mean. So creating a 'new' `Data` object clearly works, but I'll have a think re. the best way to get 'that' attribute behaviour in `data.py`...
For my reference: seems like a solution could be in here https://stackoverflow.com/questions/26598109/preserve-custom-attributes-when-pickling-subclass-of-numpy-array
So the changes required to properly pickle our `Data` type seem non-trivial - for a start `mpi4py.MPI.Comm` objects can't be pickled. Therefore I'd suggest casting to `numpy` arrays prior to...
After a quick discussion, we'll keep this open but it's not clear to us (aside from for user convenience) in what situations this functionality is particularly important? (Pickling `Data` may...