Dear all,
I am currently trying to extract a time averaged and spatially resolved map
of the local dielectric energies from a simulation in which I am looking at
a systems response to a Gaussian wave package. To my understanding meep does
not support internal calculations (ie an integration) of field data.
therefore, exporting the data after each time step (or sufficiently small
amout of time) seems to be the only solution.
As i do not necessarily need the maps to be of the same resolution as the
field, i thought of minimizing the amount of output data by reducing the
resolution. for this I tried the approach below, which however was
unsuccessful.
My question: Does anyone know of a better solution for this. I thought of a
matrix of energy_boxes as well but i am not very sure if this would be
practical, as i would need about 10^4 of these boxes.
Any ideas would be great!
Best, Moritz
...
double resolution = 1.0; // pixels per distance
volume v = vol2d(sx,sy, resolution); // 5x10 2d cell
volume v_out = vol2d(sx,sy, resolution/10.0);
structure s(v, eps, pml(20.0,Y));
fields f(&s);
f.output_hdf5(Dielectric, v_out.surroundings(),0,false,true,0);
...
f0.output_hdf5(D_EnergyDensity, v_out.surroundings(),0,false,true,0);
...
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
