From 9ded6e5063e1ab3c012a3a53844b96d06a894223 Mon Sep 17 00:00:00 2001 From: Lachlan Grose Date: Wed, 2 Oct 2024 11:20:18 +1000 Subject: [PATCH] use interpolation propagation --- .../interpolators/supports/_multiscale_structured_grid.py | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/LoopStructural/interpolators/supports/_multiscale_structured_grid.py b/LoopStructural/interpolators/supports/_multiscale_structured_grid.py index 38f7e152..c3511ab5 100644 --- a/LoopStructural/interpolators/supports/_multiscale_structured_grid.py +++ b/LoopStructural/interpolators/supports/_multiscale_structured_grid.py @@ -74,11 +74,11 @@ def propagate_property(self, property, to_level): yi = np.arange(0, self.levels[to_level]['nsteps'][1], 2 ** (to_level - from_level)) zi = np.arange(0, self.levels[to_level]['nsteps'][2], 2 ** (to_level - from_level)) ii, jj, kk = np.meshgrid(xi, yi, zi, indexing='ij') - new_nodes_idx = np.vstack( - [ii.flatten(order='f'), jj.flatten(order='f'), kk.flatten(order='f')] - ).T + # new_nodes_idx = np.vstack( + # [ii.flatten(order='f'), jj.flatten(order='f'), kk.flatten(order='f')] + # ).T self.nsteps = self.levels[to_level]['nsteps'] - gi = self.global_node_indices(new_nodes_idx).astype(int) + # gi = self.global_node_indices(new_nodes_idx).astype(int) self.set_level(from_level) # reset to the original level # get property value on all nodes print(property.shape, self.n_nodes)