sarkas.tools.observables.calc_Sk#

sarkas.tools.observables.calc_Sk(nkt, k_list, k_counts, species_np, no_dumps)[source]#

Calculate \(S_{ij}(k)\) at each saved timestep.

Parameters
  • nkt (numpy.ndarray, complex) – Density fluctuations of all species. Shape = ( no_species, no_dumps, no_ka_values)

  • k_list – List of \(k\) indices in each direction with corresponding magnitude and index of k_counts. Shape=(no_ka_values, 5)

  • k_counts (numpy.ndarray) – Number of times each \(k\) magnitude appears.

  • species_np (numpy.ndarray) – Array with number of particles of each species.

  • no_dumps (int) – Number of dumps.

Returns

Sk_all (numpy.ndarray) – Array containing \(S_{ij}(k)\). Shape=(no_Sk, no_ka_values, no_dumps)