Skip to content

Commit 5fccc9f

Browse files
Reduced the number of redistributions in spectral helper (#494)
1 parent 85dc966 commit 5fccc9f

File tree

1 file changed

+4
-6
lines changed

1 file changed

+4
-6
lines changed

pySDC/helpers/spectral_helper.py

Lines changed: 4 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -1560,14 +1560,13 @@ def transform_single_component(self, u, axes=None, padding=None):
15601560
if self.comm is not None:
15611561
_out *= np.prod([self.axes[i].N for i in _axes])
15621562

1563-
axes_next_base = axes_collapsed[(trf + 1) % len(axes_collapsed)]
1563+
axes_next_base = (axes_collapsed + [(-1,)])[trf + 1]
15641564
alignment = alignment if len(axes_next_base) == 0 else self.ndim + axes_next_base[-1]
15651565
result = self.get_aligned(
15661566
_out, axis_in=self.ndim + _axes[0], axis_out=alignment, fft=fft, forward=True, shape=shape
15671567
)
15681568

1569-
fft = self.get_fft(axes=axes, padding=padding)
1570-
return self.get_aligned(result, axis_in=alignment, axis_out=self.ndim - 1, fft=fft, forward=True, shape=shape)
1569+
return result
15711570

15721571
def transform(self, u, axes=None, padding=None):
15731572
"""
@@ -1720,14 +1719,13 @@ def itransform_single_component(self, u, axes=None, padding=None):
17201719
else:
17211720
shape[_ax] = _out.shape[_ax]
17221721

1723-
axes_next_base = axes_collapsed[(trf + 1) % len(axes_collapsed)]
1722+
axes_next_base = (axes_collapsed + [(-1,)])[trf + 1]
17241723
alignment = alignment if len(axes_next_base) == 0 else self.ndim + axes_next_base[0]
17251724
result = self.get_aligned(
17261725
_out, axis_in=self.ndim + _axes[-1], axis_out=alignment, fft=fft, forward=False, shape=shape
17271726
)
17281727

1729-
fft = self.get_fft(axes=axes, padding=padding)
1730-
return self.get_aligned(result, axis_in=alignment, axis_out=self.ndim - 1, fft=fft, shape=shape)
1728+
return result
17311729

17321730
def get_aligned(self, u, axis_in, axis_out, fft=None, forward=False, **kwargs):
17331731
"""

0 commit comments

Comments
 (0)