diff --git a/src/mpb.cpp b/src/mpb.cpp index 324922e51..a4aab11c4 100644 --- a/src/mpb.cpp +++ b/src/mpb.cpp @@ -432,6 +432,8 @@ void *fields::get_eigenmode(double omega_src, /*--------------------------------------------------------------*/ complex *cdata = (complex *) mdata->fft_data; + broadcast(0, (double*) H.data, 2 * H.n * H.p); + maxwell_compute_h_from_H(mdata, H, (scalar_complex*)cdata, band_num - 1, 1); /* choose deterministic phase, maximizing power in real part; see fix_field_phase routine in MPB.*/