bug fix in rgrid.get_grid_traveltimes

groupeLIAMG · Dec 26, 2024 · 92e15f0 · 92e15f0
1 parent 9792254
commit 92e15f0
Show file tree

Hide file tree

Showing 2 changed files with 80 additions and 28 deletions.
diff --git a/examples/example_Grid3d.ipynb b/examples/example_Grid3d.ipynb
diff --git a/src/ttcrpy/rgrid.pyx b/src/ttcrpy/rgrid.pyx
@@ -407,11 +407,11 @@ cdef class Grid3d:
         cdef vector[double] tmp
         cdef int n
         self.grid.getTT(tmp, thread_no)
-        tt = np.empty((tmp.size(),), order='F')
+        tt = np.empty((tmp.size(),))
         for n in range(tmp.size()):
             tt[n] = tmp[n]
         shape = (self._x.size(), self._y.size(), self._z.size())
-        return tt.reshape(shape)
+        return tt.reshape(shape, order='F')
 
     def ind(self, int i, int j, int k):
         """