Skip to content

Instantly share code, notes, and snippets.

@stetie
Created December 7, 2017 15:31
Show Gist options
  • Save stetie/03a5a436e259c00e3bae16145fc51c36 to your computer and use it in GitHub Desktop.
Save stetie/03a5a436e259c00e3bae16145fc51c36 to your computer and use it in GitHub Desktop.
postpic/postpic/_field_calc.py in kspace_Ex(self, **kwargs)
196 fields['By'] = self.By()
197
--> 198 return self._kspace('Ex', fields, **kwargs)
199
200 def kspace_Ey(self, **kwargs):
postpic/postpic/_field_calc.py in _kspace(self, component, fields, alignment, solver, **kwargs)
183
184 if alignment == 'epoch':
--> 185 return helper.kspace_epoch_like(component, fields, dt, align_to='B', **kwargs)
186
187 if alignment == 'epoch-final':
postpic/postpic/helper.py in kspace_epoch_like(component, fields, dt, extent, omega_func, align_to)
536 # print(k, v.extent, [(a[0], a[-1]) for a in v._conjugate_grid().values()])
537
--> 538 return kspace(component, fields, interpolation='fourier')
539
540
postpic/postpic/helper.py in kspace(component, fields, extent, interpolation, omega_func)
721 if interpolation == 'fourier':
722 # print('apply linear phase')
--> 723 field = field._apply_linear_phase(dict(enumerate(grid_shift)))
724
725 # add the field to the result with the appropriate prefactor
postpic/postpic/datahandling.py in _apply_linear_phase(self, dx)
1647 "have transformed_axes_origins not None")
1648
-> 1649 exp_ikdx = helper.linear_phase(self, dx)
1650
1651 ret = self * exp_ikdx
postpic/postpic/helper.py in linear_phase(field, dx)
744 if transform_state is True:
745 # center the axes around 0 to eliminate global phase
--> 746 axes[i] -= axes[i][gridlen//2]
747 else:
748 # start the axes at 0 to eliminate global phase
IndexError: index 0 is out of bounds for axis 0 with size 0
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment