aspcore.fouriertransform.rdft_weighting

aspcore.fouriertransform.rdft_weighting(num_real_freqs, dft_len, freqs_to_remove_low=0)

The weighting required for the real DFT to be the same as the complex DFT.

np.diag(c_diag) corresponds to C in the paper [brunnstromTimedomain2025]

Parameters:

  • num_freqs (int) – The number of frequencies in the resulting real-valued DFT. Technically redundant information as this is the same as dft_len // 2 + 1 - freqs_to_remove, so will not be necessary later.

  • dft_len (int) – The length of the DFT.

  • freqs_to_remove_low (int) – The number of frequencies that were remove from the lower end of the DFT.

Returns:

c_diag

Return type:

np.ndarray of shape (num_freqs,)