Function dft_ldos

python/simulation.py:6007–6044

`dft_ldos(fcen=None, df=None, nfreq=None, freq=None, ldos=None)` ##sig Compute the power spectrum of the sources (usually a single point dipole source), normalized to correspond to the LDOS, in either a frequency bandwidth `df` centered at `fcen` and `nfreq` equally spaced freque

(*args, **kwargs)

def dft_ldos(*args, **kwargs):
    """
    `dft_ldos(fcen=None, df=None, nfreq=None, freq=None, ldos=None)`   ##sig

    Compute the power spectrum of the sources (usually a single point dipole source),
    normalized to correspond to the LDOS, in either a frequency bandwidth `df` centered at
    `fcen` and `nfreq` equally spaced frequency points or an array/list `freq` for
    arbitrarily spaced frequencies. One can also pass in an `Ldos` object as
    `dft_ldos(ldos=my_Ldos)`.

    The resulting spectrum is outputted as comma-delimited text, prefixed by `ldos:,`, and
    is also stored in the `ldos_data` variable of the `Simulation` object after the `run`
    is complete. The Fourier-transformed electric field and current source are stored in
    the `ldos_Fdata` and `ldos_Jdata` of the `Simulation` object, respectively.
    """
    ldos = kwargs.get("ldos", None)
    if ldos is None:
        args = fix_dft_args(args, 0)
        freq = args[0]
        if isinstance(freq, (np.ndarray, list)):
            ldos = mp._dft_ldos(freq)
        else:
            raise TypeError(
                "dft_ldos only accepts freq_min,freq_max,nfreq (3 numbers) or freq (array/list) or ldos (keyword argument)"
            )

    def _ldos(sim, todo):
        if todo == "step":
            ldos.update(sim.fields)
        else:
            sim.ldos_data = mp._dft_ldos_ldos(ldos)
            sim.ldos_Fdata = mp._dft_ldos_F(ldos)
            sim.ldos_Jdata = mp._dft_ldos_J(ldos)
            sim.ldos_scale = ldos.overall_scale()
            if verbosity.meep > 0:
                display_csv(sim, "ldos", zip(mp.get_ldos_freqs(ldos), sim.ldos_data))

    return _ldos


def scale_flux_fields(s, flux):