Please explain Meep simulation output

[SilvieLuisa]

Hi,
I am running a set of very similar simulations with parallel Meep, where I use 8 cores.
Majority of the simulations run perfectly, but some of them provide this output:

Initializing structure...
Splitting into 8 chunks by voxels
time for choose_chunkdivision = 0.000213185 s
Working in 3D dimensions.
Computational cell is 2 x 2 x 9.06 with resolution 50
block, center = (0,0,0.9)
size (2,2,1.8)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon diagonal = (1,1,1)
block, center = (0,0,0.15)
size (2,2,0.3)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon diagonal = (1.1,1.1,1.1)
block, center = (0,0,0.05)
size (2,2,0.1)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon diagonal = (1.1,1.1,1.1)
block, center = (0,0,1.75)
size (2,2,3.5)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon diagonal = (1.1,1.1,1.1)
block, center = (0,0,0.03)
size (2,2,0.06)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon diagonal = (1.1,1.1,1.1)
block, center = (0,0,0.65)
size (2,2,1.3)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon diagonal = (1.1,1.1,1.1)
block, center = (-0.9,0.9,0.05)
size (0.2,0.2,0.1)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon diagonal = (1.1,1.1,1.1)
block, center = (-0.9,0.7,0.05)
size (0.2,0.2,0.1)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon diagonal = (1.1,1.1,1.1)
block, center = (-0.9,0.3,0.05)
size (0.2,0.2,0.1)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon diagonal = (1.1,1.1,1.1)
block, center = (-0.9,0.1,0.05)
size (0.2,0.2,0.1)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon diagonal = (1.1,1.1,1.1)
...(+ 25 objects not shown)...
time for set_epsilon = 2.60473 s
lorentzian susceptibility: frequency=4.84853, gamma=0
lorentzian susceptibility: frequency=5.67964, gamma=0
lorentzian susceptibility: frequency=7.35789, gamma=0.0017962
lorentzian susceptibility: frequency=3.1146, gamma=0
lorentzian susceptibility: frequency=2.73203, gamma=0.1614
lorentzian susceptibility: frequency=5.01686, gamma=0
lorentzian susceptibility: frequency=3.06792, gamma=3.6912e-51
lorentzian susceptibility: frequency=5.88666, gamma=0
lorentzian susceptibility: frequency=5.75657, gamma=3.26362e-08
lorentzian susceptibility: frequency=9.22293, gamma=1.24813e-07
lorentzian susceptibility: frequency=9.67528, gamma=0
lorentzian susceptibility: frequency=10.9301, gamma=0

meep: Total dft_chunks size does not agree with size allocated for output array.
meep: Total dft_chunks size does not agree with size allocated for output array.
meep: Total dft_chunks size does not agree with size allocated for output array.
meep: Total dft_chunks size does not agree with size allocated for output array.
meep: Total dft_chunks size does not agree with size allocated for output array.
meep: Total dft_chunks size does not agree with size allocated for output array.
meep: Total dft_chunks size does not agree with size allocated for output array.
meep: Total dft_chunks size does not agree with size allocated for output array.
Abort(1) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
Abort(1) on node 1 (rank 1 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 1
Abort(1) on node 3 (rank 3 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 3
Abort(1) on node 4 (rank 4 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 4
Abort(1) on node 5 (rank 5 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 5
Abort(1) on node 6 (rank 6 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 6
Abort(1) on node 7 (rank 7 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 7
Abort(1) on node 2 (rank 2 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 2

The differences between simulations running fine and those with errors are just simulation volume size, and some difference in the embedded structure.

Does anyone have experience with this 'total dft_chunks size does not agree with size allocated for output array' message? In what cases it can happen? I cannot put all my code here, it is very complex. So I am not expecting somebody could tell me where exactly the issue is. But maybe, if you could just help me identifying potential weak points, of typical situations when it could happen, it would be easier to debug it...

thanks a lot!

[SilvieLuisa]

The attempt of this simulation was to compute transmission spectra, as in the tutorial. And this particular error happened due to the fact that the "normalization simulation" was run with different resolution than the "true" simulation. Hopefully this will help somebody in the future, the error message is not very clear.