Skip to content

switch the spectrum code over to using AMReX's FFT #4

@zingale

Description

@zingale

see AMReX-Codes/amrex#3713

we can do:

int level = 0;
MultiFab mf = plotfile.get(level);
FFT::R2C<Real,FFT::Direction::forward> r2c(plotfile.probDomain(level));
auto const& [cba, cdm] = r2c.getSpectralDataLayout();
cMultiFab spectral_data(cba, cdm, 1, 0);
int icomp = 3; // component 3
MultiFab real_data(mf, amrex::make_alias, icomp, 1);
r2c.forward(read_data, spectral_data);
// spectral_data is now the result of real-to-complex forward FFT of the level 0 component 3 data.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions