| 911 | } |
| 912 | |
| 913 | dft_fields fields::add_dft_fields(component *components, int num_components, const volume where, |
| 914 | const double *freq, size_t Nfreq, bool use_centered_grid, |
| 915 | int decimation_factor, bool persist) { |
| 916 | bool include_dV_and_interp_weights = false; |
| 917 | bool sqrt_dV_and_interp_weights = false; // default option from meep.hpp (expose to user?) |
| 918 | std::complex<double> extra_weight = 1.0; // default option from meep.hpp (expose to user?) |
| 919 | complex<double> stored_weight = 1.0; |
| 920 | dft_chunk *chunks = NULL; |
| 921 | for (int nc = 0; nc < num_components; nc++) |
| 922 | chunks = add_dft(components[nc], where, freq, Nfreq, include_dV_and_interp_weights, |
| 923 | stored_weight, chunks, sqrt_dV_and_interp_weights, extra_weight, |
| 924 | use_centered_grid, 0, decimation_factor, persist); |
| 925 | |
| 926 | return dft_fields(chunks, freq, Nfreq, where); |
| 927 | } |
| 928 | |
| 929 | /***************************************************************/ |
| 930 | /* chunk-level processing for fields::process_dft_component. */ |