We turn to a similar but slightly different example for which there exists an analytic solution via the Fresnel equations: computing the broadband reflectance spectrum of a planar air-dielectric interface for an incident planewave over a range of angles. Similar to the previous example, we will need to run two simulations: (1) an empty cell with air/vacuum ($n$=1) everywhere to obtain the incident flux, and (2) with the dielectric ($n$=3.5) interface to obtain the reflected flux. For each angle of the incident planewave, a separate simulation is necessary.
A 1d cell must be used since a higher-dimensional cell will introduce artificial modes due to band folding. We will use a Gaussian source spanning visible wavelengths of 0.4 to 0.8 μm. Unlike a continuous-wave (CW) source, a pulsed source turns off. This enables a termination condition of when there are no fields remaining in the cell (due to absorption by the PMLs) via the run function stop_when_fields_decayed
, similar to the previous example.
Creating an oblique planewave source typically requires specifying two parameters: (1) for periodic structures, the Bloch-periodic wavevector $\vec{k}$ via k_point
, and (2) the source amplitude function amp_func
for setting the $e^{i\vec{k} \cdot \vec{r}}$ spatial dependence ($\vec{r}$ is the position vector). Since we have a 1d cell and the source is at a single point, it is not necessary to specify the source amplitude (see this 2d example for how this is done). The magnitude of the Bloch-periodic wavevector is specified according to the dispersion relation formula for a planewave in homogeneous media with index $n$: $ω=c|\vec{k}|/n$. As the source in this example is incident from air, $|\vec{k}|$ is simply equal to the frequency $ω$ (the minimum frequency of the pulse which excludes the 2π factor). Note that a fixed wavevector only applies to a single frequency. Any broadband source is therefore incident at a specified angle for only a single frequency. This is described in more detail in Section 4.5 ("Efficient Frequency-Angle Coverage") in Chapter 4 ("Electromagnetic Wave Source Conditions") of the book Advances in FDTD Computational Electrodynamics: Photonics and Nanotechnology.
In this example, the plane of incidence which contains $\vec{k}$ and the surface normal vector is $xz$. The source angle θ is defined in degrees in the counterclockwise (CCW) direction around the $y$ axis with 0 degrees along the +$z$ axis. In Meep, a 1d cell is defined along the $z$ direction. When $\vec{k}$ is not set, only the $E_x$ and $H_y$ field components are permitted. A non-zero $\vec{k}$ results in a 3d simulation where all field components are allowed and are complex (the fields are real, by default). A current source with $E_x$ polarization lies in the plane of incidence and corresponds to the convention of $\mathcal{P}$-polarization. In order to model the $\mathcal{S}$-polarization, we must use an $E_y$ source. This example involves just the $\mathcal{P}$-polarization.
import meep as mp
import math
import numpy as np
import numpy.matlib
import matplotlib.pyplot as plt
resolution = 50 # pixels/um
dpml = 1.0 # PML thickness
sz = 10 + 2 * dpml
cell_size = mp.Vector3(z=sz)
pml_layers = [mp.PML(dpml)]
wvl_min = 0.4 # min wavelength
wvl_max = 0.8 # max wavelength
fmin = 1 / wvl_max # min frequency
fmax = 1 / wvl_min # max frequency
fcen = 0.5 * (fmin + fmax) # center frequency
df = fmax - fmin # frequency width
nfreq = 50 # number of frequency bins
def planar_reflectance(theta):
# rotation angle (in degrees) of source: CCW around Y axis, 0 degrees along +Z axis
theta_r = math.radians(theta)
# plane of incidence is XZ; rotate counter clockwise (CCW) about y-axis
k = mp.Vector3(z=fmin).rotate(mp.Vector3(y=1), theta_r)
# if normal incidence, force number of dimensions to be 1
if theta_r == 0:
dimensions = 1
else:
dimensions = 3
sources = [
mp.Source(
mp.GaussianSource(fcen, fwidth=df),
component=mp.Ex,
center=mp.Vector3(z=-0.5 * sz + dpml),
)
]
sim = mp.Simulation(
cell_size=cell_size,
boundary_layers=pml_layers,
sources=sources,
k_point=k,
dimensions=dimensions,
resolution=resolution,
)
refl_fr = mp.FluxRegion(center=mp.Vector3(z=-0.25 * sz))
refl = sim.add_flux(fcen, df, nfreq, refl_fr)
sim.run(
until_after_sources=mp.stop_when_fields_decayed(
50, mp.Ex, mp.Vector3(z=-0.5 * sz + dpml), 1e-9
)
)
empty_flux = mp.get_fluxes(refl)
empty_data = sim.get_flux_data(refl)
sim.reset_meep()
# add a block with n=3.5 for the air-dielectric interface
geometry = [
mp.Block(
mp.Vector3(mp.inf, mp.inf, 0.5 * sz),
center=mp.Vector3(z=0.25 * sz),
material=mp.Medium(index=3.5),
)
]
sim = mp.Simulation(
cell_size=cell_size,
geometry=geometry,
boundary_layers=pml_layers,
sources=sources,
k_point=k,
dimensions=dimensions,
resolution=resolution,
)
refl = sim.add_flux(fcen, df, nfreq, refl_fr)
sim.load_minus_flux_data(refl, empty_data)
sim.run(
until_after_sources=mp.stop_when_fields_decayed(
50, mp.Ex, mp.Vector3(z=-0.5 * sz + dpml), 1e-9
)
)
refl_flux = mp.get_fluxes(refl)
freqs = mp.get_flux_freqs(refl)
wvls = np.empty(nfreq)
theta_out = np.empty(nfreq)
R = np.empty(nfreq)
for i in range(nfreq):
wvls[i] = 1 / freqs[i]
theta_out[i] = math.degrees(math.asin(k.x / freqs[i]))
R[i] = -refl_flux[i] / empty_flux[i]
print("refl:, {}, {}, {}, {}".format(k.x, wvls[i], theta_out[i], R[i]))
return k.x * np.ones(nfreq), wvls, theta_out, R
The function planar_reflectance
computes the reflectance at each frequency as well as the wavevector component $k_x$ and the corresponding angle for the ($k_x$, ω) pair. For those frequencies not equal to the minimum frequency of the source, this is not the same as the specified angle of the incident planewave, but rather sin-1(kx/ω).
The reflectance spectrum is generated over the angular range of 0$^\circ$ to 80$^\circ$ in increments of 5$^\circ$.
theta_in = np.arange(0, 85, 5)
wvl = np.empty(nfreq)
kxs = np.empty((nfreq, theta_in.size))
thetas = np.empty((nfreq, theta_in.size))
Rmeep = np.empty((nfreq, theta_in.size))
for j in range(theta_in.size):
kxs[:, j], wvl, thetas[:, j], Rmeep[:, j] = planar_reflectance(theta_in[j])
# create a 2d matrix for the wavelength by repeating the column vector for each angle
wvls = np.transpose(np.matlib.repmat(wvl, theta_in.size, 1))
----------- Initializing structure... time for choose_chunkdivision = 2.40803e-05 s Working in 1D dimensions. Computational cell is 0 x 0 x 12 with resolution 50 time for set_epsilon = 7.39098e-05 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.25332329653323415 / 0.25332329653323415 = 1.0 field decay(t = 100.01): 6.806395978139866e-16 / 0.25332329653323415 = 2.686841704370019e-15 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 2.71797e-05 s Working in 1D dimensions. Computational cell is 0 x 0 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.000118971 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.25332329652480207 / 0.25332329652480207 = 1.0 field decay(t = 100.01): 1.9736380723733672e-11 / 0.25332329652480207 = 7.790985272371642e-11 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.0, 0.8, 0.0, 0.2946488869770689 refl:, 0.0, 0.784, 0.0, 0.29411278345905817 refl:, 0.0, 0.7686274509803922, 0.0, 0.2935428929532476 refl:, 0.0, 0.7538461538461539, 0.0, 0.29293473702048123 refl:, 0.0, 0.739622641509434, 0.0, 0.2923073383425624 refl:, 0.0, 0.7259259259259259, 0.0, 0.29164890565286355 refl:, 0.0, 0.7127272727272727, 0.0, 0.29098366907954704 refl:, 0.0, 0.7, 0.0, 0.29030209388401956 refl:, 0.0, 0.6877192982456141, 0.0, 0.28959812074303476 refl:, 0.0, 0.6758620689655173, 0.0, 0.28888122149723194 refl:, 0.0, 0.664406779661017, 0.0, 0.28815027904146656 refl:, 0.0, 0.6533333333333333, 0.0, 0.2874049842378321 refl:, 0.0, 0.6426229508196721, 0.0, 0.2866449790098622 refl:, 0.0, 0.632258064516129, 0.0, 0.28586718403268657 refl:, 0.0, 0.6222222222222222, 0.0, 0.2850720751223629 refl:, 0.0, 0.6124999999999999, 0.0, 0.28425994421351475 refl:, 0.0, 0.6030769230769231, 0.0, 0.2834305476307166 refl:, 0.0, 0.593939393939394, 0.0, 0.2825862116808428 refl:, 0.0, 0.5850746268656717, 0.0, 0.28172703444497454 refl:, 0.0, 0.5764705882352942, 0.0, 0.28085005219605824 refl:, 0.0, 0.5681159420289855, 0.0, 0.2799551048425921 refl:, 0.0, 0.56, 0.0, 0.2790440893576681 refl:, 0.0, 0.552112676056338, 0.0, 0.27811614989969974 refl:, 0.0, 0.5444444444444444, 0.0, 0.27716906652214396 refl:, 0.0, 0.536986301369863, 0.0, 0.2762029568966099 refl:, 0.0, 0.5297297297297298, 0.0, 0.27521916171170663 refl:, 0.0, 0.5226666666666667, 0.0, 0.2742172838003258 refl:, 0.0, 0.5157894736842105, 0.0, 0.2731957001654285 refl:, 0.0, 0.509090909090909, 0.0, 0.27215409786728667 refl:, 0.0, 0.5025641025641026, 0.0, 0.2710932430411507 refl:, 0.0, 0.4962025316455696, 0.0, 0.27001246262335565 refl:, 0.0, 0.49, 0.0, 0.268910054299966 refl:, 0.0, 0.4839506172839506, 0.0, 0.2677861833693871 refl:, 0.0, 0.47804878048780486, 0.0, 0.2666421048931029 refl:, 0.0, 0.47228915662650606, 0.0, 0.26547679611705527 refl:, 0.0, 0.4666666666666667, 0.0, 0.2642883298672654 refl:, 0.0, 0.4611764705882353, 0.0, 0.2630777642724503 refl:, 0.0, 0.4558139534883721, 0.0, 0.261846739625148 refl:, 0.0, 0.4505747126436782, 0.0, 0.2605925093678602 refl:, 0.0, 0.4454545454545454, 0.0, 0.2593113496920778 refl:, 0.0, 0.44044943820224725, 0.0, 0.2580054281160953 refl:, 0.0, 0.43555555555555553, 0.0, 0.25667922159327894 refl:, 0.0, 0.4307692307692308, 0.0, 0.2553302969233571 refl:, 0.0, 0.4260869565217391, 0.0, 0.25395270771108786 refl:, 0.0, 0.421505376344086, 0.0, 0.25254987709773263 refl:, 0.0, 0.41702127659574467, 0.0, 0.25113151895670055 refl:, 0.0, 0.4126315789473684, 0.0, 0.24969310724781182 refl:, 0.0, 0.4083333333333333, 0.0, 0.24821384020578532 refl:, 0.0, 0.4041237113402062, 0.0, 0.24668825253510232 refl:, 0.0, 0.4, 0.0, 0.24514471927281825 ----------- Initializing structure... time for choose_chunkdivision = 3.91006e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00191998 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.25242167342001054 / 0.25242167342001054 = 1.0 field decay(t = 100.01): 1.8867425273501205e-14 / 0.25242167342001054 = 7.474566275498553e-14 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 3.88622e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00448012 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.2524216734361254 / 0.2524216734361254 = 1.0 field decay(t = 100.01): 2.0310851045086525e-11 / 0.2524216734361254 = 8.046397430379974e-11 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.1089446784345727, 0.8, 5.0, 0.2933049593316432 refl:, 0.1089446784345727, 0.784, 4.899752997934953, 0.2928302248191875 refl:, 0.1089446784345727, 0.7686274509803922, 4.803451415694315, 0.2923055731561396 refl:, 0.1089446784345727, 0.7538461538461539, 4.710866569098618, 0.2917466563012179 refl:, 0.1089446784345727, 0.739622641509434, 4.621787131270349, 0.2911629722039186 refl:, 0.1089446784345727, 0.7259259259259259, 4.536017514803643, 0.2905470818015815 refl:, 0.1089446784345727, 0.7127272727272727, 4.45337643175598, 0.28992287302292685 refl:, 0.1089446784345727, 0.7, 4.373695609047488, 0.28927810074884625 refl:, 0.1089446784345727, 0.6877192982456141, 4.296818640028281, 0.28860968806846743 refl:, 0.1089446784345727, 0.6758620689655173, 4.222599955651541, 0.28792678919854764 refl:, 0.1089446784345727, 0.664406779661017, 4.150903900954991, 0.287228177914529 refl:, 0.1089446784345727, 0.6533333333333333, 4.081603904475513, 0.28651426279788256 refl:, 0.1089446784345727, 0.6426229508196721, 4.014581729857545, 0.2857836503745603 refl:, 0.1089446784345727, 0.632258064516129, 3.9497268003123533, 0.2850332352892976 refl:, 0.1089446784345727, 0.6222222222222222, 3.886935587780479, 0.2842645360517432 refl:, 0.1089446784345727, 0.6124999999999999, 3.826111059675504, 0.2834779258184347 refl:, 0.1089446784345727, 0.6030769230769231, 3.7671621769699377, 0.2826728763462015 refl:, 0.1089446784345727, 0.593939393939394, 3.7100034381453693, 0.2818517562980362 refl:, 0.1089446784345727, 0.5850746268656717, 3.6545544641874574, 0.2810148164875731 refl:, 0.1089446784345727, 0.5764705882352942, 3.6007396203768605, 0.28015923094614514 refl:, 0.1089446784345727, 0.5681159420289855, 3.5484876711229183, 0.2792846976539734 refl:, 0.1089446784345727, 0.56, 3.4977314645183832, 0.2783930023602693 refl:, 0.1089446784345727, 0.552112676056338, 3.4484076436699387, 0.2774836544079885 refl:, 0.1089446784345727, 0.5444444444444444, 3.4004563821883713, 0.2765546012121586 refl:, 0.1089446784345727, 0.536986301369863, 3.353821141510494, 0.27560563009266964 refl:, 0.1089446784345727, 0.5297297297297298, 3.3084484479779683, 0.2746381229823153 refl:, 0.1089446784345727, 0.5226666666666667, 3.2642876878206, 0.27365213822849105 refl:, 0.1089446784345727, 0.5157894736842105, 3.221290918387684, 0.2726459819587968 refl:, 0.1089446784345727, 0.509090909090909, 3.1794126941439105, 0.2716189126246391 refl:, 0.1089446784345727, 0.5025641025641026, 3.1386099060992034, 0.27057188214742395 refl:, 0.1089446784345727, 0.4962025316455696, 3.098841633477275, 0.2695046844668159 refl:, 0.1089446784345727, 0.49, 3.060069006547709, 0.26841546175441744 refl:, 0.1089446784345727, 0.4839506172839506, 3.0222550796530636, 0.26730405423702647 refl:, 0.1089446784345727, 0.47804878048780486, 2.985364713557375, 0.2661719545046955 refl:, 0.1089446784345727, 0.47228915662650606, 2.9493644663270095, 0.26501846138724444 refl:, 0.1089446784345727, 0.4666666666666667, 2.914222492030252, 0.2638414562416327 refl:, 0.1089446784345727, 0.4611764705882353, 2.879908446609464, 0.2626417282749142 refl:, 0.1089446784345727, 0.4558139534883721, 2.84639340033999, 0.2614210758850605 refl:, 0.1089446784345727, 0.4505747126436782, 2.8136497563440477, 0.26017691547552313 refl:, 0.1089446784345727, 0.4454545454545454, 2.7816511746763815, 0.25890544613245936 refl:, 0.1089446784345727, 0.44044943820224725, 2.750372501542037, 0.2576089367857666 refl:, 0.1089446784345727, 0.43555555555555553, 2.719789703245832, 0.25629198685456567 refl:, 0.1089446784345727, 0.4307692307692308, 2.689879804508449, 0.2549519690830257 refl:, 0.1089446784345727, 0.4260869565217391, 2.660620830815881, 0.2535828658853722 refl:, 0.1089446784345727, 0.421505376344086, 2.631991754497782, 0.2521885223227569 refl:, 0.1089446784345727, 0.41702127659574467, 2.603972444256211, 0.25077854364690916 refl:, 0.1089446784345727, 0.4126315789473684, 2.576543617889841, 0.24934708585148915 refl:, 0.1089446784345727, 0.4083333333333333, 2.549686797979975, 0.24787364791621788 refl:, 0.1089446784345727, 0.4041237113402062, 2.5233842703240543, 0.2463554261780657 refl:, 0.1089446784345727, 0.4, 2.4976190449198983, 0.24482106271892765 ----------- Initializing structure... time for choose_chunkdivision = 3.91006e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00183392 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.24974540035413884 / 0.24974540035413884 = 1.0 field decay(t = 100.01): 6.006906599019063e-14 / 0.24974540035413884 = 2.405212104207434e-13 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.19617e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00453806 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.24974540044608917 / 0.24974540044608917 = 1.0 field decay(t = 100.01): 2.1519308385459792e-11 / 0.24974540044608917 = 8.616498380759976e-11 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.2170602220836629, 0.8, 10.0, 0.2892656432100493 refl:, 0.2170602220836629, 0.784, 9.798006528153513, 0.288965065351793 refl:, 0.2170602220836629, 0.7686274509803922, 9.604050171837292, 0.2885818690170257 refl:, 0.2170602220836629, 0.7538461538461539, 9.417658416993296, 0.28817212668801895 refl:, 0.2170602220836629, 0.739622641509434, 9.238395240840497, 0.28771910655279653 refl:, 0.2170602220836629, 0.7259259259259259, 9.06585764090149, 0.2872339286740473 refl:, 0.2170602220836629, 0.7127272727272727, 8.899672554443574, 0.2867321699029375 refl:, 0.2170602220836629, 0.7, 8.739494117841588, 0.2861987897084524 refl:, 0.2170602220836629, 0.6877192982456141, 8.585001222725978, 0.28564000163622033 refl:, 0.2170602220836629, 0.6758620689655173, 8.435895331947279, 0.2850609513737311 refl:, 0.2170602220836629, 0.664406779661017, 8.291898523577625, 0.2844608473128053 refl:, 0.2170602220836629, 0.6533333333333333, 8.152751735551202, 0.2838419805391297 refl:, 0.2170602220836629, 0.6426229508196721, 8.018213187256704, 0.28320005732998216 refl:, 0.2170602220836629, 0.632258064516129, 7.88805695754783, 0.28253271331372304 refl:, 0.2170602220836629, 0.6222222222222222, 7.762071701325296, 0.28184393698477367 refl:, 0.2170602220836629, 0.6124999999999999, 7.640059489140416, 0.2811344534194306 refl:, 0.2170602220836629, 0.6030769230769231, 7.521834756238996, 0.28040341723697604 refl:, 0.2170602220836629, 0.593939393939394, 7.407223349155971, 0.27965263406145413 refl:, 0.2170602220836629, 0.5850746268656717, 7.296061659428915, 0.2788825964767311 refl:, 0.2170602220836629, 0.5764705882352942, 7.188195835257705, 0.27809153810572734 refl:, 0.2170602220836629, 0.5681159420289855, 7.083481063027822, 0.27727879066309463 refl:, 0.2170602220836629, 0.56, 6.981780911561047, 0.2764453605102568 refl:, 0.2170602220836629, 0.552112676056338, 6.882966732780441, 0.2755918264938716 refl:, 0.2170602220836629, 0.5444444444444444, 6.786917113194022, 0.274717044631823 refl:, 0.2170602220836629, 0.536986301369863, 6.693517371228444, 0.2738198627903547 refl:, 0.2170602220836629, 0.5297297297297298, 6.602659095992853, 0.2729014096931157 refl:, 0.2170602220836629, 0.5226666666666667, 6.514239723534344, 0.27196303613184486 refl:, 0.2170602220836629, 0.5157894736842105, 6.428162147069652, 0.2710031700214432 refl:, 0.2170602220836629, 0.509090909090909, 6.3443343580501015, 0.2700199119299729 refl:, 0.2170602220836629, 0.5025641025641026, 6.262669115245526, 0.26901445255671513 refl:, 0.2170602220836629, 0.4962025316455696, 6.1830836393232005, 0.2679878933990742 refl:, 0.2170602220836629, 0.49, 6.105499330654852, 0.2669382978356776 refl:, 0.2170602220836629, 0.4839506172839506, 6.029841508312739, 0.2658645437475429 refl:, 0.2170602220836629, 0.47804878048780486, 5.956039168418179, 0.26476841577725835 refl:, 0.2170602220836629, 0.47228915662650606, 5.884024760185935, 0.26365019399207895 refl:, 0.2170602220836629, 0.4666666666666667, 5.813733978168244, 0.2625075913337705 refl:, 0.2170602220836629, 0.4611764705882353, 5.745105569345401, 0.26134055474931606 refl:, 0.2170602220836629, 0.4558139534883721, 5.678081153837623, 0.2601509047307578 refl:, 0.2170602220836629, 0.4505747126436782, 5.612605058127373, 0.2589367681695608 refl:, 0.2170602220836629, 0.4454545454545454, 5.5486241597837695, 0.2576946836654639 refl:, 0.2170602220836629, 0.44044943820224725, 5.486087742772728, 0.256426818319966 refl:, 0.2170602220836629, 0.43555555555555553, 5.424947362519055, 0.2551375207965109 refl:, 0.2170602220836629, 0.4307692307692308, 5.365156719961124, 0.2538238582454513 refl:, 0.2170602220836629, 0.4260869565217391, 5.306671543905553, 0.25248026000140456 refl:, 0.2170602220836629, 0.421505376344086, 5.249449481049872, 0.25111138429336477 refl:, 0.2170602220836629, 0.41702127659574467, 5.19344999309547, 0.24972542744836562 refl:, 0.2170602220836629, 0.4126315789473684, 5.138634260422533, 0.24831408973125035 refl:, 0.2170602220836629, 0.4083333333333333, 5.08496509184317, 0.24685936094756866 refl:, 0.2170602220836629, 0.4041237113402062, 5.032406839989342, 0.24536500901018427 refl:, 0.2170602220836629, 0.4, 4.980925321928872, 0.24385884403865477 ----------- Initializing structure... time for choose_chunkdivision = 3.69549e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00182509 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.24537918139687429 / 0.24537918139687429 = 1.0 field decay(t = 100.01): 1.2755953350203158e-13 / 0.24537918139687429 = 5.198466013941004e-13 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.3869e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00469923 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.2453791816184804 / 0.2453791816184804 = 1.0 field decay(t = 100.01): 2.2796495624696677e-11 / 0.2453791816184804 = 9.290313658369374e-11 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.3235238063781509, 0.8, 14.999999999999998, 0.28250329291621734 refl:, 0.3235238063781509, 0.784, 14.693171512000124, 0.2824643918340546 refl:, 0.3235238063781509, 0.7686274509803922, 14.398780921441814, 0.2823434344182011 refl:, 0.3235238063781509, 0.7538461538461539, 14.116078899389818, 0.2821736091292865 refl:, 0.3235238063781509, 0.739622641509434, 13.844375746673084, 0.2819428557616974 refl:, 0.3235238063781509, 0.7259259259259259, 13.583035518835887, 0.28168405270433 refl:, 0.3235238063781509, 0.7127272727272727, 13.331470838933798, 0.28138555830619527 refl:, 0.3235238063781509, 0.7, 13.089138305160036, 0.28104620398023994 refl:, 0.3235238063781509, 0.6877192982456141, 12.855534414514684, 0.2806783763575039 refl:, 0.3235238063781509, 0.6758620689655173, 12.630191935533823, 0.28027812041674527 refl:, 0.3235238063781509, 0.664406779661017, 12.412676672931594, 0.27984738142823545 refl:, 0.3235238063781509, 0.6533333333333333, 12.202584575236058, 0.27938821203729947 refl:, 0.3235238063781509, 0.6426229508196721, 11.999539143408532, 0.27889571979963745 refl:, 0.3235238063781509, 0.632258064516129, 11.803189104258525, 0.2783707318747749 refl:, 0.3235238063781509, 0.6222222222222222, 11.613206317390016, 0.27781778136422264 refl:, 0.3235238063781509, 0.6124999999999999, 11.429283888592414, 0.277238647976418 refl:, 0.3235238063781509, 0.6030769230769231, 11.25113446614539, 0.2766335972563 refl:, 0.3235238063781509, 0.593939393939394, 11.078488699542484, 0.2760021980086433 refl:, 0.3235238063781509, 0.5850746268656717, 10.91109384273803, 0.2753445983679576 refl:, 0.3235238063781509, 0.5764705882352942, 10.748712486253877, 0.2746618737814714 refl:, 0.3235238063781509, 0.5681159420289855, 10.591121404404543, 0.27395346744418553 refl:, 0.3235238063781509, 0.56, 10.438110505558328, 0.273218482798473 refl:, 0.3235238063781509, 0.552112676056338, 10.289481874787974, 0.272458491647809 refl:, 0.3235238063781509, 0.5444444444444444, 10.145048899510067, 0.27167442241968504 refl:, 0.3235238063781509, 0.536986301369863, 10.004635469795673, 0.27086454108766905 refl:, 0.3235238063781509, 0.5297297297297298, 9.868075245978739, 0.2700288346581343 refl:, 0.3235238063781509, 0.5226666666666667, 9.735210987013529, 0.2691696679864803 refl:, 0.3235238063781509, 0.5157894736842105, 9.60589393375409, 0.26828660843247504 refl:, 0.3235238063781509, 0.509090909090909, 9.479983241961918, 0.26737703033509175 refl:, 0.3235238063781509, 0.5025641025641026, 9.3573454604044, 0.26644151143474176 refl:, 0.3235238063781509, 0.4962025316455696, 9.237854049896542, 0.26548236142478776 refl:, 0.3235238063781509, 0.49, 9.121388939570695, 0.26449868322216635 refl:, 0.3235238063781509, 0.4839506172839506, 9.00783611704105, 0.2634884411439653 refl:, 0.3235238063781509, 0.47804878048780486, 8.897087249467761, 0.26245254376131855 refl:, 0.3235238063781509, 0.47228915662650606, 8.789039332825531, 0.2613924415313307 refl:, 0.3235238063781509, 0.4666666666666667, 8.683594366947904, 0.2603070110071445 refl:, 0.3235238063781509, 0.4611764705882353, 8.5806590541555, 0.25919500218989733 refl:, 0.3235238063781509, 0.4558139534883721, 8.480144519487721, 0.2580569212415155 refl:, 0.3235238063781509, 0.4505747126436782, 8.381966050745921, 0.2568924101657985 refl:, 0.3235238063781509, 0.4454545454545454, 8.286042856724524, 0.2557000732428357 refl:, 0.3235238063781509, 0.44044943820224725, 8.192297842157469, 0.2544810421270595 refl:, 0.3235238063781509, 0.43555555555555553, 8.100657398042404, 0.25323728280622937 refl:, 0.3235238063781509, 0.4307692307692308, 8.011051206126568, 0.25196641359587585 refl:, 0.3235238063781509, 0.4260869565217391, 7.923412056447124, 0.2506655440913078 refl:, 0.3235238063781509, 0.421505376344086, 7.837675676917115, 0.2493388062231876 refl:, 0.3235238063781509, 0.41702127659574467, 7.753780574036344, 0.24799024185005186 refl:, 0.3235238063781509, 0.4126315789473684, 7.671667883886533, 0.24661141415700796 refl:, 0.3235238063781509, 0.4083333333333333, 7.591281232641979, 0.24519197231834583 refl:, 0.3235238063781509, 0.4041237113402062, 7.512566605892175, 0.2437413870063278 refl:, 0.3235238063781509, 0.4, 7.435472226131853, 0.24228102145492145 ----------- Initializing structure... time for choose_chunkdivision = 3.69549e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00182009 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.2394607243108125 / 0.2394607243108125 = 1.0 field decay(t = 100.01): 2.1433673054356211e-13 / 0.2394607243108125 = 8.950809413963009e-13 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.19617e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00456595 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.23946072470582733 / 0.23946072470582733 = 1.0 field decay(t = 100.01): 2.292441064702724e-11 / 0.23946072470582733 = 9.57334889685539e-11 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.4275251791570859, 0.8, 20.0, 0.27291140903435845 refl:, 0.4275251791570859, 0.784, 19.583468236198428, 0.2732447943229538 refl:, 0.4275251791570859, 0.7686274509803922, 19.184283570762837, 0.2735018552296093 refl:, 0.4275251791570859, 0.7538461538461539, 18.80136282780684, 0.27366446189146665 refl:, 0.4275251791570859, 0.739622641509434, 18.433712921276232, 0.273781793501152 refl:, 0.4275251791570859, 0.7259259259259259, 18.080421505473552, 0.27384858957084335 refl:, 0.4275251791570859, 0.7127272727272727, 17.74064878668832, 0.27385095593564335 refl:, 0.4275251791570859, 0.7, 17.413620328364388, 0.27380597771345955 refl:, 0.4275251791570859, 0.6877192982456141, 17.098620709730763, 0.2737136095477295 refl:, 0.4275251791570859, 0.6758620689655173, 16.794987920274295, 0.273573243892546 refl:, 0.4275251791570859, 0.664406779661017, 16.50210839086207, 0.2733876760611016 refl:, 0.4275251791570859, 0.6533333333333333, 16.21941257752279, 0.27315442100604564 refl:, 0.4275251791570859, 0.6426229508196721, 15.946371026493122, 0.2728767468514616 refl:, 0.4275251791570859, 0.632258064516129, 15.682490859619808, 0.27255955817684324 refl:, 0.4275251791570859, 0.6222222222222222, 15.427312627971409, 0.27220355920057415 refl:, 0.4275251791570859, 0.6124999999999999, 15.18040748886694, 0.2718114250177165 refl:, 0.4275251791570859, 0.6030769230769231, 14.941374667722823, 0.2713856767538751 refl:, 0.4275251791570859, 0.593939393939394, 14.709839171355561, 0.2709245485567433 refl:, 0.4275251791570859, 0.5850746268656717, 14.485449723819883, 0.27042763367941924 refl:, 0.4275251791570859, 0.5764705882352942, 14.267876899642369, 0.2698980098701011 refl:, 0.4275251791570859, 0.5681159420289855, 14.056811432538773, 0.2693367446910556 refl:, 0.4275251791570859, 0.56, 13.851962680467771, 0.2687424115671181 refl:, 0.4275251791570859, 0.552112676056338, 13.653057230248155, 0.2681156337689324 refl:, 0.4275251791570859, 0.5444444444444444, 13.459837627011849, 0.26745868029520964 refl:, 0.4275251791570859, 0.536986301369863, 13.272061215531336, 0.2667718443134419 refl:, 0.4275251791570859, 0.5297297297297298, 13.089499081989581, 0.2660542100897866 refl:, 0.4275251791570859, 0.5226666666666667, 12.911935086088024, 0.26530649401202394 refl:, 0.4275251791570859, 0.5157894736842105, 12.73916497454351, 0.2645300141858129 refl:, 0.4275251791570859, 0.509090909090909, 12.570995568032437, 0.2637243046129653 refl:, 0.4275251791570859, 0.5025641025641026, 12.407244014521105, 0.26288850346304676 refl:, 0.4275251791570859, 0.4962025316455696, 12.247737102692609, 0.2620237658492669 refl:, 0.4275251791570859, 0.49, 12.092310629857906, 0.2611315394355365 refl:, 0.4275251791570859, 0.4839506172839506, 11.94080881933422, 0.26021080661473406 refl:, 0.4275251791570859, 0.47804878048780486, 11.793083782798943, 0.259260103032999 refl:, 0.4275251791570859, 0.47228915662650606, 11.648995023590732, 0.25828076720993753 refl:, 0.4275251791570859, 0.4666666666666667, 11.508408977339421, 0.2572746105587765 refl:, 0.4275251791570859, 0.4611764705882353, 11.37119858666983, 0.25624008547337146 refl:, 0.4275251791570859, 0.4558139534883721, 11.23724290704701, 0.2551747627993682 refl:, 0.4275251791570859, 0.4505747126436782, 11.106426741117268, 0.254079676475692 refl:, 0.4275251791570859, 0.4454545454545454, 10.978640299154753, 0.2529572394711134 refl:, 0.4275251791570859, 0.44044943820224725, 10.853778883451254, 0.25180680229382696 refl:, 0.4275251791570859, 0.43555555555555553, 10.731742594690267, 0.25062573390797804 refl:, 0.4275251791570859, 0.4307692307692308, 10.61243605852877, 0.2494136021716301 refl:, 0.4275251791570859, 0.4260869565217391, 10.495768170772996, 0.24817270591418827 refl:, 0.4275251791570859, 0.421505376344086, 10.381651859681224, 0.24690457932715454 refl:, 0.4275251791570859, 0.41702127659574467, 10.270003864057921, 0.24560713839504952 refl:, 0.4275251791570859, 0.4126315789473684, 10.160744525922071, 0.24427619924007887 refl:, 0.4275251791570859, 0.4083333333333333, 10.053797596639106, 0.24291142234897387 refl:, 0.4275251791570859, 0.4041237113402062, 9.949090055502005, 0.24152100364291515 refl:, 0.4275251791570859, 0.4, 9.846551939834079, 0.2401173303544212 ----------- Initializing structure... time for choose_chunkdivision = 3.60012e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00180292 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.2321756751554577 / 0.2321756751554577 = 1.0 field decay(t = 100.01): 3.189584689866576e-13 / 0.2321756751554577 = 1.373780732081829e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.1008e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00448799 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.23217567572618025 / 0.23217567572618025 = 1.0 field decay(t = 100.01): 2.5113865846466273e-11 / 0.23217567572618025 = 1.0816751482650867e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.5282728271758743, 0.8, 25.0, 0.26019389637751095 refl:, 0.5282728271758743, 0.784, 24.46679999189225, 0.26109407463193457 refl:, 0.5282728271758743, 0.7686274509803922, 23.95662859236144, 0.2618692444514501 refl:, 0.5282728271758743, 0.7538461538461539, 23.46797596156662, 0.26257268253476557 refl:, 0.5282728271758743, 0.739622641509434, 22.99946566139384, 0.26318360986138806 refl:, 0.5282728271758743, 0.7259259259259259, 22.54983979533518, 0.2636884917611983 refl:, 0.5282728271758743, 0.7127272727272727, 22.117946146436154, 0.2641153624871549 refl:, 0.5282728271758743, 0.7, 21.70272699899909, 0.26446400131053266 refl:, 0.5282728271758743, 0.6877192982456141, 21.303209386217617, 0.2647361502445863 refl:, 0.5282728271758743, 0.6758620689655173, 20.91849655102253, 0.26494127762665987 refl:, 0.5282728271758743, 0.664406779661017, 20.54776044367494, 0.26507870572225684 refl:, 0.5282728271758743, 0.6533333333333333, 20.19023510896709, 0.2651506958652906 refl:, 0.5282728271758743, 0.6426229508196721, 19.84521083974841, 0.26516452813944175 refl:, 0.5282728271758743, 0.632258064516129, 19.51202899301187, 0.2651243633954837 refl:, 0.5282728271758743, 0.6222222222222222, 19.190077380825585, 0.2650321194109797 refl:, 0.5282728271758743, 0.6124999999999999, 18.878786161658894, 0.26488903211488374 refl:, 0.5282728271758743, 0.6030769230769231, 18.577624168665142, 0.2646966666638105 refl:, 0.5282728271758743, 0.593939393939394, 18.28609562066815, 0.2644576983261531 refl:, 0.5282728271758743, 0.5850746268656717, 18.003737169291743, 0.26417368977600525 refl:, 0.5282728271758743, 0.5764705882352942, 17.730115242139682, 0.26384517442033967 refl:, 0.5282728271758743, 0.5681159420289855, 17.46482364739322, 0.2634743009887046 refl:, 0.5282728271758743, 0.56, 17.207481409818527, 0.2630632173645374 refl:, 0.5282728271758743, 0.552112676056338, 16.957730812108316, 0.2626116463801122 refl:, 0.5282728271758743, 0.5444444444444444, 16.715235618835717, 0.26212001269997814 refl:, 0.5282728271758743, 0.536986301369863, 16.479679463167955, 0.2615912946581222 refl:, 0.5282728271758743, 0.5297297297297298, 16.250764378950514, 0.26102663469509957 refl:, 0.5282728271758743, 0.5226666666666667, 16.028209462892256, 0.26042439141314616 refl:, 0.5282728271758743, 0.5157894736842105, 15.811749653412038, 0.2597853758930952 refl:, 0.5282728271758743, 0.509090909090909, 15.60113461429131, 0.25911262730709367 refl:, 0.5282728271758743, 0.5025641025641026, 15.396127712651408, 0.2584060888970194 refl:, 0.5282728271758743, 0.4962025316455696, 15.196505081970068, 0.2576642746546754 refl:, 0.5282728271758743, 0.49, 15.002054761894165, 0.2568888948983748 refl:, 0.5282728271758743, 0.4839506172839506, 14.81257590751694, 0.25608173235353965 refl:, 0.5282728271758743, 0.47804878048780486, 14.627878061586326, 0.25524132391655996 refl:, 0.5282728271758743, 0.47228915662650606, 14.44778048381156, 0.254367293135576 refl:, 0.5282728271758743, 0.4666666666666667, 14.272111532051735, 0.2534621115579794 refl:, 0.5282728271758743, 0.4611764705882353, 14.100708090713411, 0.25252576242439695 refl:, 0.5282728271758743, 0.4558139534883721, 13.933415042164041, 0.25155529986985325 refl:, 0.5282728271758743, 0.4505747126436782, 13.770084777392734, 0.2505511615005892 refl:, 0.5282728271758743, 0.4454545454545454, 13.610576742526133, 0.2495168318010228 refl:, 0.5282728271758743, 0.44044943820224725, 13.454757018141454, 0.24845156231536186 refl:, 0.5282728271758743, 0.43555555555555553, 13.30249792861589, 0.2473511491305419 refl:, 0.5282728271758743, 0.4307692307692308, 13.153677679016674, 0.24621635842933665 refl:, 0.5282728271758743, 0.4260869565217391, 13.008180017272103, 0.2450522325472059 refl:, 0.5282728271758743, 0.421505376344086, 12.865893919575512, 0.24385875474669194 refl:, 0.5282728271758743, 0.41702127659574467, 12.726713297162874, 0.24263075122017086 refl:, 0.5282728271758743, 0.4126315789473684, 12.590536722774571, 0.24136640786665386 refl:, 0.5282728271758743, 0.4083333333333333, 12.457267175263851, 0.24006853393330466 refl:, 0.5282728271758743, 0.4041237113402062, 12.326811800951434, 0.23874368705732368 refl:, 0.5282728271758743, 0.4, 12.199081690448809, 0.23740409018120856 ----------- Initializing structure... time for choose_chunkdivision = 3.8147e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00182295 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.22375100718842017 / 0.22375100718842017 = 1.0 field decay(t = 100.01): 4.405483227742841e-13 / 0.22375100718842017 = 1.9689221885973432e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 3.91006e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00455689 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.22375100783150073 / 0.22375100783150073 = 1.0 field decay(t = 100.01): 2.6045908524047717e-11 / 0.22375100783150073 = 1.1640577075595567e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.6249999999999999, 0.8, 29.999999999999993, 0.2441304698440292 refl:, 0.6249999999999999, 0.784, 29.34058157502373, 0.2458848744990323 refl:, 0.6249999999999999, 0.7686274509803922, 28.711017527148794, 0.2474738617848535 refl:, 0.6249999999999999, 0.7538461538461539, 28.10922128260952, 0.24886145900455134 refl:, 0.6249999999999999, 0.739622641509434, 27.53330580109674, 0.2500790166042368 refl:, 0.6249999999999999, 0.7259259259259259, 26.98155921981659, 0.2511706573365747 refl:, 0.6249999999999999, 0.7127272727272727, 26.452424118557172, 0.25213606232996266 refl:, 0.6249999999999999, 0.7, 25.944479772370002, 0.25298936027517865 refl:, 0.6249999999999999, 0.6877192982456141, 25.45642688540588, 0.2537413324767768 refl:, 0.6249999999999999, 0.6758620689655173, 24.98707439783817, 0.2543902864127213 refl:, 0.6249999999999999, 0.664406779661017, 24.53532803474182, 0.25494503656419537 refl:, 0.6249999999999999, 0.6533333333333333, 24.10018032643974, 0.25541541144304014 refl:, 0.6249999999999999, 0.6426229508196721, 23.680701877992792, 0.2558038470011666 refl:, 0.6249999999999999, 0.632258064516129, 23.276033704033082, 0.25611422492454605 refl:, 0.6249999999999999, 0.6222222222222222, 22.885380476158563, 0.25635226359269697 refl:, 0.6249999999999999, 0.6124999999999999, 22.508004555237044, 0.2565212690513844 refl:, 0.6249999999999999, 0.6030769230769231, 22.14322070144805, 0.25662363175125674 refl:, 0.6249999999999999, 0.593939393939394, 21.79039137167391, 0.25666316405820044 refl:, 0.6249999999999999, 0.5850746268656717, 21.448922527676533, 0.256643853974939 refl:, 0.6249999999999999, 0.5764705882352942, 21.118259889941786, 0.25656744769664874 refl:, 0.6249999999999999, 0.5681159420289855, 20.797885581592283, 0.2564348059692535 refl:, 0.6249999999999999, 0.56, 20.487315114722662, 0.2562488813390074 refl:, 0.6249999999999999, 0.552112676056338, 20.186094678183196, 0.2560129738604307 refl:, 0.6249999999999999, 0.5444444444444444, 19.89379869145801, 0.25572798412991815 refl:, 0.6249999999999999, 0.536986301369863, 19.610027594036424, 0.2553946757545203 refl:, 0.6249999999999999, 0.5297297297297298, 19.33440584370919, 0.2550156190782023 refl:, 0.6249999999999999, 0.5226666666666667, 19.066580100655866, 0.2545925551142987 refl:, 0.6249999999999999, 0.5157894736842105, 18.806217577124535, 0.25412540391688965 refl:, 0.6249999999999999, 0.509090909090909, 18.553004535020655, 0.25361528880754325 refl:, 0.6249999999999999, 0.5025641025641026, 18.306644915884704, 0.25306481668131153 refl:, 0.6249999999999999, 0.4962025316455696, 18.066859089603533, 0.25247461257562936 refl:, 0.6249999999999999, 0.49, 17.833382709813034, 0.2518436268795026 refl:, 0.6249999999999999, 0.4839506172839506, 17.605965665348197, 0.2511731587087493 refl:, 0.6249999999999999, 0.47804878048780486, 17.38437111831227, 0.25046601975628274 refl:, 0.6249999999999999, 0.47228915662650606, 17.168374620396104, 0.24972209992178607 refl:, 0.6249999999999999, 0.4666666666666667, 16.957763300004142, 0.24893969262180818 refl:, 0.6249999999999999, 0.4611764705882353, 16.752335113553887, 0.24811981145545228 refl:, 0.6249999999999999, 0.4558139534883721, 16.551898155026578, 0.24726440526511775 refl:, 0.6249999999999999, 0.4505747126436782, 16.35627001847215, 0.24637290574011259 refl:, 0.6249999999999999, 0.4454545454545454, 16.165277208722518, 0.24544450395766448 refl:, 0.6249999999999999, 0.44044943820224725, 15.978754596053776, 0.24448027110760395 refl:, 0.6249999999999999, 0.43555555555555553, 15.796544910968182, 0.2434805663315552 refl:, 0.6249999999999999, 0.4307692307692308, 15.618498275648548, 0.2424444686589383 refl:, 0.6249999999999999, 0.4260869565217391, 15.44447176897603, 0.24137271247893222 refl:, 0.6249999999999999, 0.421505376344086, 15.274329022304034, 0.24026658833730827 refl:, 0.6249999999999999, 0.41702127659574467, 15.107939843449019, 0.23912295725148944 refl:, 0.6249999999999999, 0.4126315789473684, 14.94517986659885, 0.2379359418651507 refl:, 0.6249999999999999, 0.4083333333333333, 14.78593022605342, 0.23670966886718617 refl:, 0.6249999999999999, 0.4041237113402062, 14.630077251904048, 0.23546524192379745 refl:, 0.6249999999999999, 0.4, 14.477512185929921, 0.23421631294354242 ----------- Initializing structure... time for choose_chunkdivision = 3.69549e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00170493 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.21444706390029775 / 0.21444706390029775 = 1.0 field decay(t = 100.01): 5.749510873093161e-13 / 0.21444706390029775 = 2.681086310310252e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 3.98159e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00455904 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.21444706445458084 / 0.21444706445458084 = 1.0 field decay(t = 100.01): 2.737725033432578e-11 / 0.21444706445458084 = 1.2766437444110684e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.7169705454388076, 0.8, 35.0, 0.22495734929489236 refl:, 0.7169705454388076, 0.784, 34.201491482224874, 0.2276898082643635 refl:, 0.7169705454388076, 0.7686274509803922, 33.4413597291606, 0.23011693205976114 refl:, 0.7169705454388076, 0.7538461538461539, 32.71669424141867, 0.23234222429680992 refl:, 0.7169705454388076, 0.739622641509434, 32.02489215108529, 0.23437969108695692 refl:, 0.7169705454388076, 0.7259259259259259, 31.363616172201723, 0.23622265510793064 refl:, 0.7169705454388076, 0.7127272727272727, 30.73075960358516, 0.23789643333674684 refl:, 0.7169705454388076, 0.7, 30.12441698853809, 0.23940633852434773 refl:, 0.7169705454388076, 0.6877192982456141, 29.542859352190995, 0.24075746810602938 refl:, 0.7169705454388076, 0.6758620689655173, 28.984513173223203, 0.2419699286741594 refl:, 0.7169705454388076, 0.664406779661017, 28.447942424882076, 0.24305307471225007 refl:, 0.7169705454388076, 0.6533333333333333, 27.931833156170086, 0.2440117580982575 refl:, 0.7169705454388076, 0.6426229508196721, 27.43498018882287, 0.24485799959527518 refl:, 0.7169705454388076, 0.632258064516129, 26.95627558715663, 0.24559973015569014 refl:, 0.7169705454388076, 0.6222222222222222, 26.49469862174186, 0.246240896723184 refl:, 0.7169705454388076, 0.6124999999999999, 26.049306998353728, 0.2467892899221832 refl:, 0.7169705454388076, 0.6030769230769231, 25.619229163859046, 0.24725182835118295 refl:, 0.7169705454388076, 0.593939393939394, 25.20365753294407, 0.24763129782088683 refl:, 0.7169705454388076, 0.5850746268656717, 24.80184250561243, 0.24793098018363321 refl:, 0.7169705454388076, 0.5764705882352942, 24.41308716651858, 0.24815599140237504 refl:, 0.7169705454388076, 0.5681159420289855, 24.036742574466203, 0.24831042769942835 refl:, 0.7169705454388076, 0.56, 23.672203564580393, 0.24839661296784196 refl:, 0.7169705454388076, 0.552112676056338, 23.318904997368072, 0.24841763613734327 refl:, 0.7169705454388076, 0.5444444444444444, 22.976318398592397, 0.2483776953869765 refl:, 0.7169705454388076, 0.536986301369863, 22.643948941980693, 0.2482787140153355 refl:, 0.7169705454388076, 0.5297297297297298, 22.321332733561132, 0.24812129287197915 refl:, 0.7169705454388076, 0.5226666666666667, 22.008034362119382, 0.24790881905698 refl:, 0.7169705454388076, 0.5157894736842105, 21.703644685073712, 0.2476447465584906 refl:, 0.7169705454388076, 0.509090909090909, 21.407778823139967, 0.2473290950751843 refl:, 0.7169705454388076, 0.5025641025641026, 21.120074340620857, 0.24696260849110413 refl:, 0.7169705454388076, 0.4962025316455696, 20.840189591108683, 0.24654850371687256 refl:, 0.7169705454388076, 0.49, 20.56780221092012, 0.24608799987122265 refl:, 0.7169705454388076, 0.4839506172839506, 20.302607744753697, 0.24558075371549817 refl:, 0.7169705454388076, 0.47804878048780486, 20.044318389931625, 0.24502869978373035 refl:, 0.7169705454388076, 0.47228915662650606, 19.79266184720379, 0.24443385311897003 refl:, 0.7169705454388076, 0.4666666666666667, 19.54738026749198, 0.24379544974149034 refl:, 0.7169705454388076, 0.4611764705882353, 19.308229285168313, 0.24311337479355408 refl:, 0.7169705454388076, 0.4558139534883721, 19.07497712952083, 0.24239050667955775 refl:, 0.7169705454388076, 0.4505747126436782, 18.847403806983923, 0.24162858318113342 refl:, 0.7169705454388076, 0.4454545454545454, 18.62530034751981, 0.24082594853100417 refl:, 0.7169705454388076, 0.44044943820224725, 18.408468109247096, 0.23998192730485585 refl:, 0.7169705454388076, 0.43555555555555553, 18.196718136035763, 0.23909849717321205 refl:, 0.7169705454388076, 0.4307692307692308, 17.98987056333767, 0.23817577312539148 refl:, 0.7169705454388076, 0.4260869565217391, 17.787754068005896, 0.23721144993308976 refl:, 0.7169705454388076, 0.421505376344086, 17.590205358285754, 0.2362040994387872 refl:, 0.7169705454388076, 0.41702127659574467, 17.39706870053973, 0.23515140286483463 refl:, 0.7169705454388076, 0.4126315789473684, 17.20819547960629, 0.23405349223820357 refl:, 0.7169705454388076, 0.4083333333333333, 17.02344378999232, 0.23292711235344504 refl:, 0.7169705454388076, 0.4041237113402062, 16.842678055366356, 0.23179557860673924 refl:, 0.7169705454388076, 0.4, 16.665768674058118, 0.2306452929106223 ----------- Initializing structure... time for choose_chunkdivision = 3.88622e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.001827 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.20454876789491913 / 0.20454876789491913 = 1.0 field decay(t = 100.01): 7.178635176478305e-13 / 0.20454876789491913 = 3.5094981262200107e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.1008e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00457597 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.20454876813558653 / 0.20454876813558653 = 1.0 field decay(t = 100.01): 2.817457244887082e-11 / 0.20454876813558653 = 1.3774012283562185e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.8034845121081741, 0.8, 39.99999999999999, 0.20174574732399359 refl:, 0.8034845121081741, 0.784, 39.04509528455022, 0.20589297541852558 refl:, 0.8034845121081741, 0.7686274509803922, 38.139646206191365, 0.20966853236552122 refl:, 0.8034845121081741, 0.7538461538461539, 37.27949726994566, 0.21305998243366164 refl:, 0.8034845121081741, 0.739622641509434, 36.46098976062103, 0.2161315876607331 refl:, 0.8034845121081741, 0.7259259259259259, 35.680884053707004, 0.21892504829192913 refl:, 0.8034845121081741, 0.7127272727272727, 34.93629688867566, 0.22145999607107203 refl:, 0.8034845121081741, 0.7, 34.22465020933728, 0.22376578875003042 refl:, 0.8034845121081741, 0.6877192982456141, 33.54362905484555, 0.22586507657858915 refl:, 0.8034845121081741, 0.6758620689655173, 32.89114661043291, 0.22776751087226974 refl:, 0.8034845121081741, 0.664406779661017, 32.265314978971055, 0.22948642344304926 refl:, 0.8034845121081741, 0.6533333333333333, 31.664420565727855, 0.23104021748754944 refl:, 0.8034845121081741, 0.6426229508196721, 31.086903214629274, 0.23243958627884664 refl:, 0.8034845121081741, 0.632258064516129, 30.531338419096144, 0.2336946954040348 refl:, 0.8034845121081741, 0.6222222222222222, 29.996422070856152, 0.23482007111973577 refl:, 0.8034845121081741, 0.6124999999999999, 29.480957317802616, 0.23582329675340516 refl:, 0.8034845121081741, 0.6030769230769231, 28.983843185365757, 0.2367091706706393 refl:, 0.8034845121081741, 0.593939393939394, 28.504064681021994, 0.2374872020062935 refl:, 0.8034845121081741, 0.5850746268656717, 28.0406841528977, 0.238164188713723 refl:, 0.8034845121081741, 0.5764705882352942, 27.592833714170844, 0.23874383693611143 refl:, 0.8034845121081741, 0.5681159420289855, 27.159708577552564, 0.2392325677866777 refl:, 0.8034845121081741, 0.56, 26.740561170354134, 0.23963653341947785 refl:, 0.8034845121081741, 0.552112676056338, 26.33469592188663, 0.23995882233761215 refl:, 0.8034845121081741, 0.5444444444444444, 25.94146463225064, 0.2402024441340643 refl:, 0.8034845121081741, 0.536986301369863, 25.560262345758233, 0.24037225972376797 refl:, 0.8034845121081741, 0.5297297297297298, 25.190523663916437, 0.24047264456383183 refl:, 0.8034845121081741, 0.5226666666666667, 24.831719442577725, 0.24050529500914392 refl:, 0.8034845121081741, 0.5157894736842105, 24.483353825914165, 0.24047249027879555 refl:, 0.8034845121081741, 0.509090909090909, 24.144961576600384, 0.2403786028523382 refl:, 0.8034845121081741, 0.5025641025641026, 23.816105667237974, 0.24022523323130834 refl:, 0.8034845121081741, 0.4962025316455696, 23.496375102814216, 0.24001232236295097 refl:, 0.8034845121081741, 0.49, 23.185382948015043, 0.23974362785116177 refl:, 0.8034845121081741, 0.4839506172839506, 22.882764536632674, 0.23942247402605424 refl:, 0.8034845121081741, 0.47804878048780486, 22.58817584322306, 0.23904770280781232 refl:, 0.8034845121081741, 0.47228915662650606, 22.301291999661412, 0.23862003627448614 refl:, 0.8034845121081741, 0.4666666666666667, 22.021805941382443, 0.23814371157526962 refl:, 0.8034845121081741, 0.4611764705882353, 21.749427169932794, 0.23762004855970983 refl:, 0.8034845121081741, 0.4558139534883721, 21.483880620051718, 0.2370480859309985 refl:, 0.8034845121081741, 0.4505747126436782, 21.224905620871482, 0.2364296182998171 refl:, 0.8034845121081741, 0.4454545454545454, 20.972254942022712, 0.23576666887769465 refl:, 0.8034845121081741, 0.44044943820224725, 20.725693916468796, 0.23505766085138496 refl:, 0.8034845121081741, 0.43555555555555553, 20.48499963280006, 0.2343015839900727 refl:, 0.8034845121081741, 0.4307692307692308, 20.249960190511292, 0.2335015406593809 refl:, 0.8034845121081741, 0.4260869565217391, 20.020374012481035, 0.23265820727265243 refl:, 0.8034845121081741, 0.421505376344086, 19.79604920948226, 0.23176538827635734 refl:, 0.8034845121081741, 0.41702127659574467, 19.576802992091316, 0.2308222435106821 refl:, 0.8034845121081741, 0.4126315789473684, 19.362461125837058, 0.22984486603425214 refl:, 0.8034845121081741, 0.4083333333333333, 19.15285742585141, 0.22884869903534621 refl:, 0.8034845121081741, 0.4041237113402062, 18.9478332876545, 0.22781887957732183 refl:, 0.8034845121081741, 0.4, 18.747237251037504, 0.22672147973970627 ----------- Initializing structure... time for choose_chunkdivision = 3.69549e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00182915 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.1943564114319852 / 0.1943564114319852 = 1.0 field decay(t = 100.01): 8.613369502151169e-13 / 0.1943564114319852 = 4.431739317828168e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.1008e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.0045619 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.19435641121678224 / 0.19435641121678224 = 1.0 field decay(t = 100.01): 2.929219089204079e-11 / 0.19435641121678224 = 1.507137876680009e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.8838834764831843, 0.8, 44.99999999999999, 0.17502993854504276 refl:, 0.8838834764831843, 0.784, 43.865246854678944, 0.1808719938626426 refl:, 0.8838834764831843, 0.7686274509803922, 42.79498686880298, 0.1861559520009496 refl:, 0.8838834764831843, 0.7538461538461539, 41.78306957363623, 0.19097169679998788 refl:, 0.8838834764831843, 0.739622641509434, 40.82419653278361, 0.19534107187624397 refl:, 0.8838834764831843, 0.7259259259259259, 39.913765803587594, 0.1993129093290473 refl:, 0.8838834764831843, 0.7127272727272727, 39.047751315763335, 0.2029343832170293 refl:, 0.8838834764831843, 0.7, 38.2226079814105, 0.20622282929854657 refl:, 0.8838834764831843, 0.6877192982456141, 37.435196089128915, 0.20921161798306928 refl:, 0.8838834764831843, 0.6758620689655173, 36.682720369951866, 0.21193294983781294 refl:, 0.8838834764831843, 0.664406779661017, 35.962680378531154, 0.2144100930400042 refl:, 0.8838834764831843, 0.6533333333333333, 35.272829708778914, 0.21666480283186992 refl:, 0.8838834764831843, 0.6426229508196721, 34.61114218453038, 0.21871574748298628 refl:, 0.8838834764831843, 0.632258064516129, 33.975783613579615, 0.2205809372467119 refl:, 0.8838834764831843, 0.6222222222222222, 33.36508802077569, 0.22227155758719339 refl:, 0.8838834764831843, 0.6124999999999999, 32.77753751829758, 0.22379898710754023 refl:, 0.8838834764831843, 0.6030769230769231, 32.21174515294742, 0.2251782492582617 refl:, 0.8838834764831843, 0.593939393939394, 31.666440208036935, 0.22641733524917368 refl:, 0.8838834764831843, 0.5850746268656717, 31.14045554291693, 0.22752444186490398 refl:, 0.8838834764831843, 0.5764705882352942, 30.6327166347442, 0.22851133567283208 refl:, 0.8838834764831843, 0.5681159420289855, 30.142232050688204, 0.22938387252712697 refl:, 0.8838834764831843, 0.56, 29.66808512880701, 0.2301469724038025 refl:, 0.8838834764831843, 0.552112676056338, 29.20942668547734, 0.23080897055556673 refl:, 0.8838834764831843, 0.5444444444444444, 28.765468598924116, 0.23137563774572306 refl:, 0.8838834764831843, 0.536986301369863, 28.33547814384704, 0.23185070033374858 refl:, 0.8838834764831843, 0.5297297297297298, 27.91877297273401, 0.2322392190088606 refl:, 0.8838834764831843, 0.5226666666666667, 27.514716656212954, 0.23254636390929412 refl:, 0.8838834764831843, 0.5157894736842105, 27.12271470851516, 0.23277522069046336 refl:, 0.8838834764831843, 0.509090909090909, 26.742211035417114, 0.2329277669161597 refl:, 0.8838834764831843, 0.5025641025641026, 26.372684751370738, 0.233008303204484 refl:, 0.8838834764831843, 0.4962025316455696, 26.013647320299032, 0.2330211687709118 refl:, 0.8838834764831843, 0.49, 25.66463998102006, 0.23296698717426295 refl:, 0.8838834764831843, 0.4839506172839506, 25.32523142370283, 0.23284755035583113 refl:, 0.8838834764831843, 0.47804878048780486, 24.99501568834094, 0.2326674384773264 refl:, 0.8838834764831843, 0.47228915662650606, 24.673610260104642, 0.23242802131165444 refl:, 0.8838834764831843, 0.4666666666666667, 24.360654339720863, 0.2321295593079615 refl:, 0.8838834764831843, 0.4611764705882353, 24.055807269832574, 0.23177582439583932 refl:, 0.8838834764831843, 0.4558139534883721, 23.75874710068368, 0.2313690613476082 refl:, 0.8838834764831843, 0.4505747126436782, 23.46916928052956, 0.23090791937820976 refl:, 0.8838834764831843, 0.4454545454545454, 23.18678545794031, 0.23039322072655252 refl:, 0.8838834764831843, 0.44044943820224725, 22.911322384688706, 0.2298284080169662 refl:, 0.8838834764831843, 0.43555555555555553, 22.64252090923418, 0.22921510204235154 refl:, 0.8838834764831843, 0.4307692307692308, 22.380135051959574, 0.2285521481969639 refl:, 0.8838834764831843, 0.4260869565217391, 22.123931154313652, 0.22783726183812486 refl:, 0.8838834764831843, 0.421505376344086, 21.873687094881706, 0.22707071257620573 refl:, 0.8838834764831843, 0.41702127659574467, 21.629191566166806, 0.22625950710942383 refl:, 0.8838834764831843, 0.4126315789473684, 21.390243406530644, 0.22541142020660002 refl:, 0.8838834764831843, 0.4083333333333333, 21.156650982328358, 0.22451918835662907 refl:, 0.8838834764831843, 0.4041237113402062, 20.928231615787418, 0.22356313569076627 refl:, 0.8838834764831843, 0.4, 20.704811054635428, 0.2225460542244332 ----------- Initializing structure... time for choose_chunkdivision = 3.8147e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00170398 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.18417632037527296 / 0.18417632037527296 = 1.0 field decay(t = 100.01): 1.0178427596032929e-12 / 0.18417632037527296 = 5.526458328244166e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.19617e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00469804 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.1841763197552185 / 0.1841763197552185 = 1.0 field decay(t = 100.01): 2.871301261786922e-11 / 0.1841763197552185 = 1.5589958935019742e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.9575555538987225, 0.8, 50.0, 0.14423305713624324 refl:, 0.9575555538987225, 0.784, 48.6530933185261, 0.15238194243131842 refl:, 0.9575555538987225, 0.7686274509803922, 47.39207738436693, 0.15968064236197976 refl:, 0.9575555538987225, 0.7538461538461539, 46.207396358038395, 0.166276599415515 refl:, 0.9575555538987225, 0.739622641509434, 45.091066324150106, 0.172235759713748 refl:, 0.9575555538987225, 0.7259259259259259, 44.036334278275795, 0.17762992926260002 refl:, 0.9575555538987225, 0.7127272727272727, 43.037427456538545, 0.18252108862299557 refl:, 0.9575555538987225, 0.7, 42.089365210987516, 0.18696921001018088 refl:, 0.9575555538987225, 0.6877192982456141, 41.18781524255363, 0.1910237249088231 refl:, 0.9575555538987225, 0.6758620689655173, 40.32898196453324, 0.1947139212974769 refl:, 0.9575555538987225, 0.664406779661017, 39.50951857909891, 0.1980824834747826 refl:, 0.9575555538987225, 0.6533333333333333, 38.726456948394606, 0.20115657663131353 refl:, 0.9575555538987225, 0.6426229508196721, 37.97715101972185, 0.2039606150416146 refl:, 0.9575555538987225, 0.632258064516129, 37.25923071465613, 0.2065188140121264 refl:, 0.9575555538987225, 0.6222222222222222, 36.57056399547784, 0.2088504663121866 refl:, 0.9575555538987225, 0.6124999999999999, 35.909225393217994, 0.21097639595491705 refl:, 0.9575555538987225, 0.6030769230769231, 35.273469693564564, 0.2129096520737003 refl:, 0.9575555538987225, 0.593939393939394, 34.6617097783444, 0.21466551482576457 refl:, 0.9575555538987225, 0.5850746268656717, 34.07249784378941, 0.21625947963710546 refl:, 0.9575555538987225, 0.5764705882352942, 33.504509384470694, 0.2176988158642594 refl:, 0.9575555538987225, 0.5681159420289855, 32.95652945897567, 0.21899518230958334 refl:, 0.9575555538987225, 0.56, 32.42744085087325, 0.2201600193505694 refl:, 0.9575555538987225, 0.552112676056338, 31.91621381391884, 0.2211986038256294 refl:, 0.9575555538987225, 0.5444444444444444, 31.42189714930561, 0.2221192201942536 refl:, 0.9575555538987225, 0.536986301369863, 30.943610409083956, 0.22293005892732626 refl:, 0.9575555538987225, 0.5297297297297298, 30.480537056602333, 0.22363546953064717 refl:, 0.9575555538987225, 0.5226666666666667, 30.031918444163345, 0.22424057729239344 refl:, 0.9575555538987225, 0.5157894736842105, 29.597048491686717, 0.22475131722696567 refl:, 0.9575555538987225, 0.509090909090909, 29.17526896927061, 0.22517270735852093 refl:, 0.9575555538987225, 0.5025641025641026, 28.76596530209651, 0.22550829130360708 refl:, 0.9575555538987225, 0.4962025316455696, 28.36856282886028, 0.22576168560272664 refl:, 0.9575555538987225, 0.49, 27.982523455399697, 0.22593769488637538 refl:, 0.9575555538987225, 0.4839506172839506, 27.60734265386835, 0.22603879926079012 refl:, 0.9575555538987225, 0.47804878048780486, 27.24254676502359, 0.22606671684807417 refl:, 0.9575555538987225, 0.47228915662650606, 26.88769056722596, 0.22602615267005247 refl:, 0.9575555538987225, 0.4666666666666667, 26.54235508080712, 0.2259196867892714 refl:, 0.9575555538987225, 0.4611764705882353, 26.206145580726076, 0.22574707111781514 refl:, 0.9575555538987225, 0.4558139534883721, 25.878689794039474, 0.22551123719165922 refl:, 0.9575555538987225, 0.4505747126436782, 25.55963626177365, 0.22521534349380676 refl:, 0.9575555538987225, 0.4454545454545454, 25.248652847395146, 0.2248599992539298 refl:, 0.9575555538987225, 0.44044943820224725, 24.945425376307522, 0.22444833597837902 refl:, 0.9575555538987225, 0.43555555555555553, 24.64965639271632, 0.2239834956983558 refl:, 0.9575555538987225, 0.4307692307692308, 24.361064021851533, 0.2234627996648224 refl:, 0.9575555538987225, 0.4260869565217391, 24.079380926958596, 0.22288475662006782 refl:, 0.9575555538987225, 0.421505376344086, 23.804353351700417, 0.22225599445518354 refl:, 0.9575555538987225, 0.41702127659574467, 23.535740239681235, 0.221581310984325 refl:, 0.9575555538987225, 0.4126315789473684, 23.27331242373393, 0.22085233490014997 refl:, 0.9575555538987225, 0.4083333333333333, 23.016851878423925, 0.22005799710068308 refl:, 0.9575555538987225, 0.4041237113402062, 22.76615102993319, 0.2192050777671104 refl:, 0.9575555538987225, 0.4, 22.521012118111, 0.21831293144861783 ----------- Initializing structure... time for choose_chunkdivision = 3.69549e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00181794 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.17431196343224817 / 0.17431196343224817 = 1.0 field decay(t = 100.01): 1.2265996190910864e-12 / 0.17431196343224817 = 7.036806854440848e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.1008e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00457215 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.1743119626418479 / 0.1743119626418479 = 1.0 field decay(t = 100.01): 2.9177363352545384e-11 / 0.1743119626418479 = 1.6738589199695374e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.0239400553612397, 0.8, 54.99999999999999, 0.11017776540639429 refl:, 1.0239400553612397, 0.784, 53.39534223243568, 0.12101496561527146 refl:, 1.0239400553612397, 0.7686274509803922, 51.90867481651075, 0.1306771391774596 refl:, 1.0239400553612397, 0.7538461538461539, 50.524208813748004, 0.13937263755149704 refl:, 1.0239400553612397, 0.739622641509434, 49.22931269112476, 0.14718623029207292 refl:, 1.0239400553612397, 0.7259259259259259, 48.013685667151286, 0.15426411470561957 refl:, 1.0239400553612397, 0.7127272727272727, 46.86879216202182, 0.16065288211741083 refl:, 1.0239400553612397, 0.7, 45.787462696319906, 0.16645430349254092 refl:, 1.0239400553612397, 0.6877192982456141, 44.7636047370949, 0.17172280994961617 refl:, 1.0239400553612397, 0.6758620689655173, 43.79198840588911, 0.17652234609536085 refl:, 1.0239400553612397, 0.664406779661017, 42.868084512177944, 0.18089733371336672 refl:, 1.0239400553612397, 0.6533333333333333, 41.98794000564178, 0.18488560841000873 refl:, 1.0239400553612397, 0.6426229508196721, 41.148080731197126, 0.18853239113664866 refl:, 1.0239400553612397, 0.632258064516129, 40.345434464399176, 0.1918629583011057 refl:, 1.0239400553612397, 0.6222222222222222, 39.57726925315509, 0.19490950570711735 refl:, 1.0239400553612397, 0.6124999999999999, 38.841143478370526, 0.197693493807791 refl:, 1.0239400553612397, 0.6030769230769231, 38.13486500383326, 0.2002377008340049 refl:, 1.0239400553612397, 0.593939393939394, 37.45645745905701, 0.20256263625319026 refl:, 1.0239400553612397, 0.5850746268656717, 36.80413218014341, 0.20468094633969736 refl:, 1.0239400553612397, 0.5764705882352942, 36.176264682924305, 0.2066127525275387 refl:, 1.0239400553612397, 0.5681159420289855, 35.57137479946729, 0.20837010022232308 refl:, 1.0239400553612397, 0.56, 34.98810980028866, 0.2099625922722152 refl:, 1.0239400553612397, 0.552112676056338, 34.42522996870841, 0.21140498333943816 refl:, 1.0239400553612397, 0.5444444444444444, 33.881596203503065, 0.21270480827886393 refl:, 1.0239400553612397, 0.536986301369863, 33.35615931039313, 0.2138698898904216 refl:, 1.0239400553612397, 0.5297297297297298, 32.847950708397136, 0.2149102942440705 refl:, 1.0239400553612397, 0.5226666666666667, 32.3560743283589, 0.21583196386549577 refl:, 1.0239400553612397, 0.5157894736842105, 31.87969952142079, 0.21664120771803203 refl:, 1.0239400553612397, 0.509090909090909, 31.41805482739263, 0.2173447451325115 refl:, 1.0239400553612397, 0.5025641025641026, 30.97042247873359, 0.21794797391235327 refl:, 1.0239400553612397, 0.4962025316455696, 30.536133536635912, 0.21845532089955091 refl:, 1.0239400553612397, 0.49, 30.114563572550164, 0.21887060442537454 refl:, 1.0239400553612397, 0.4839506172839506, 29.70512882224086, 0.21919938384338775 refl:, 1.0239400553612397, 0.47804878048780486, 29.307282750744406, 0.21944555930181503 refl:, 1.0239400553612397, 0.47228915662650606, 28.920512975908146, 0.21961061853019698 refl:, 1.0239400553612397, 0.4666666666666667, 28.544338505905493, 0.21969923529515722 refl:, 1.0239400553612397, 0.4611764705882353, 28.178307252549438, 0.2197148009196021 refl:, 1.0239400553612397, 0.4558139534883721, 27.821993787605123, 0.2196570179143748 refl:, 1.0239400553612397, 0.4505747126436782, 27.474997313821564, 0.21953044408992073 refl:, 1.0239400553612397, 0.4454545454545454, 27.13693982621647, 0.21934102718237625 refl:, 1.0239400553612397, 0.44044943820224725, 26.80746444237855, 0.2190876250731135 refl:, 1.0239400553612397, 0.43555555555555553, 26.486233883298294, 0.21876881981339763 refl:, 1.0239400553612397, 0.4307692307692308, 26.172929088582443, 0.2183889418977696 refl:, 1.0239400553612397, 0.4260869565217391, 25.867247951913267, 0.21795300571421844 refl:, 1.0239400553612397, 0.421505376344086, 25.56890416433822, 0.2174619823180691 refl:, 1.0239400553612397, 0.41702127659574467, 25.277626154460066, 0.21691277583436266 refl:, 1.0239400553612397, 0.4126315789473684, 24.993156115881433, 0.21630289003159692 refl:, 1.0239400553612397, 0.4083333333333333, 24.715249113370163, 0.21563472861478358 refl:, 1.0239400553612397, 0.4041237113402062, 24.443672260178424, 0.21491428731801748 refl:, 1.0239400553612397, 0.4, 24.178203959791162, 0.21414224338441656 ----------- Initializing structure... time for choose_chunkdivision = 3.69549e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00193405 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.16505586379118406 / 0.16505586379118406 = 1.0 field decay(t = 100.01): 1.6796657965795737e-12 / 0.16505586379118406 = 1.0176347316594321e-11 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.1008e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00460005 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.16505586311279843 / 0.16505586311279843 = 1.0 field decay(t = 100.01): 3.3223715036954263e-11 / 0.16505586311279843 = 2.0128769987558289e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.0825317547305482, 0.8, 59.99999999999999, 0.0735591806237639 refl:, 1.0825317547305482, 0.784, 58.07108488083593, 0.08751326272230432 refl:, 1.0825317547305482, 0.7686274509803922, 56.311309260669624, 0.09991502299592077 refl:, 1.0825317547305482, 0.7538461538461539, 54.69254540505164, 0.11103055503068811 refl:, 1.0825317547305482, 0.739622641509434, 53.19365179333026, 0.12095912989771326 refl:, 1.0825317547305482, 0.7259259259259259, 51.798237196004614, 0.12988938721094026 refl:, 1.0825317547305482, 0.7127272727272727, 50.49327371888257, 0.13795276336817966 refl:, 1.0825317547305482, 0.7, 49.268194122622916, 0.14525663007066364 refl:, 1.0825317547305482, 0.6877192982456141, 48.11428082873218, 0.15187367313781958 refl:, 1.0825317547305482, 0.6758620689655173, 47.024238635081495, 0.15789709160960827 refl:, 1.0825317547305482, 0.664406779661017, 45.99188753936469, 0.1633722181385885 refl:, 1.0825317547305482, 0.6533333333333333, 45.01193662107726, 0.16837298204602447 refl:, 1.0825317547305482, 0.6426229508196721, 44.07981414374115, 0.17292974882236836 refl:, 1.0825317547305482, 0.632258064516129, 43.191537588393835, 0.17709784614849625 refl:, 1.0825317547305482, 0.6222222222222222, 42.34361264736431, 0.1809078709487836 refl:, 1.0825317547305482, 0.6124999999999999, 41.53295361356457, 0.18439488499615878 refl:, 1.0825317547305482, 0.6030769230769231, 40.756819839298586, 0.1875870463998401 refl:, 1.0825317547305482, 0.593939393939394, 40.01276444425686, 0.1905058338431826 refl:, 1.0825317547305482, 0.5850746268656717, 39.298592486049785, 0.19318024997825003 refl:, 1.0825317547305482, 0.5764705882352942, 38.61232652957572, 0.19562392092427208 refl:, 1.0825317547305482, 0.5681159420289855, 37.9521780657052, 0.19785775786220766 refl:, 1.0825317547305482, 0.56, 37.31652360111115, 0.19989832950774553 refl:, 1.0825317547305482, 0.552112676056338, 36.70388451304266, 0.20175483810150383 refl:, 1.0825317547305482, 0.5444444444444444, 36.112909964600526, 0.20344366885684412 refl:, 1.0825317547305482, 0.536986301369863, 35.54236232752769, 0.20497463489294676 refl:, 1.0825317547305482, 0.5297297297297298, 34.991104674473426, 0.20635674385823924 refl:, 1.0825317547305482, 0.5226666666666667, 34.4580899908132, 0.20760148101525444 refl:, 1.0825317547305482, 0.5157894736842105, 33.94235182430718, 0.2087154840835457 refl:, 1.0825317547305482, 0.509090909090909, 33.442996144127704, 0.2097054505580977 refl:, 1.0825317547305482, 0.5025641025641026, 32.95919422270091, 0.21057842554599443 refl:, 1.0825317547305482, 0.4962025316455696, 32.4901763870572, 0.21134139847096645 refl:, 1.0825317547305482, 0.49, 32.035226512949734, 0.2120001650839906 refl:, 1.0825317547305482, 0.4839506172839506, 31.593677156370052, 0.21255799660930014 refl:, 1.0825317547305482, 0.47804878048780486, 31.164905234389202, 0.2130209470984357 refl:, 1.0825317547305482, 0.47228915662650606, 30.748328181342448, 0.2133940857681709 refl:, 1.0825317547305482, 0.4666666666666667, 30.343400517915825, 0.21367846587930897 refl:, 1.0825317547305482, 0.4611764705882353, 29.949610780196405, 0.21387912459727104 refl:, 1.0825317547305482, 0.4558139534883721, 29.566478763613983, 0.21400172267644174 refl:, 1.0825317547305482, 0.4505747126436782, 29.193553043244272, 0.2140466258675452 refl:, 1.0825317547305482, 0.4454545454545454, 28.830408737409662, 0.21401666370463251 refl:, 1.0825317547305482, 0.44044943820224725, 28.47664548610104, 0.21391675944760216 refl:, 1.0825317547305482, 0.43555555555555553, 28.131885619609587, 0.21374717744561747 refl:, 1.0825317547305482, 0.4307692307692308, 27.79577249602797, 0.21350966862372817 refl:, 1.0825317547305482, 0.4260869565217391, 27.46796898905743, 0.2132082256317414 refl:, 1.0825317547305482, 0.421505376344086, 27.14815610992485, 0.21284457069627338 refl:, 1.0825317547305482, 0.41702127659574467, 26.836031749237822, 0.21242137686646906 refl:, 1.0825317547305482, 0.4126315789473684, 26.531309526343414, 0.2119364849641077 refl:, 1.0825317547305482, 0.4083333333333333, 26.23371773525153, 0.21138678961842441 refl:, 1.0825317547305482, 0.4041237113402062, 25.94299837747481, 0.210778113086434 refl:, 1.0825317547305482, 0.4, 25.65890627325528, 0.21012210961113728 ----------- Initializing structure... time for choose_chunkdivision = 3.79086e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00197196 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.15668265201002066 / 0.15668265201002066 = 1.0 field decay(t = 100.01): 2.244682331122525e-12 / 0.15668265201002066 = 1.432629778936194e-11 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 3.91006e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00452209 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.1566826516265942 / 0.1566826516265942 = 1.0 field decay(t = 100.01): 3.452976999440086e-11 / 0.1566826516265942 = 2.2038030143051282e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.1328847337958123, 0.8, 64.99999999999999, 0.037128559862592 refl:, 1.1328847337958123, 0.784, 62.645633431180784, 0.0536330451598086 refl:, 1.1328847337958123, 0.7686274509803922, 60.547811567584965, 0.06876671939925749 refl:, 1.1328847337958123, 0.7538461538461539, 58.65172469781183, 0.08235712625730884 refl:, 1.1328847337958123, 0.739622641509434, 56.919780417243786, 0.09463354363126596 refl:, 1.1328847337958123, 0.7259259259259259, 55.32480276439999, 0.10558417457801342 refl:, 1.1328847337958123, 0.7127272727272727, 53.84635215217889, 0.11545216623293213 refl:, 1.1328847337958123, 0.7, 52.468573108244414, 0.12431529050720325 refl:, 1.1328847337958123, 0.6877192982456141, 51.17885777022349, 0.13235617584005718 refl:, 1.1328847337958123, 0.6758620689655173, 49.96697533403755, 0.1396274674220686 refl:, 1.1328847337958123, 0.664406779661017, 48.82448244349733, 0.14625610319510743 refl:, 1.1328847337958123, 0.6533333333333333, 47.744310665402864, 0.15227959565841903 refl:, 1.1328847337958123, 0.6426229508196721, 46.720469815981424, 0.15778582117204662 refl:, 1.1328847337958123, 0.632258064516129, 45.7478295128327, 0.16280641813158056 refl:, 1.1328847337958123, 0.6222222222222222, 44.82195500406973, 0.16739947756546966 refl:, 1.1328847337958123, 0.6124999999999999, 43.93898156096055, 0.17159829305125546 refl:, 1.1328847337958123, 0.6030769230769231, 43.09551684594588, 0.17544425185043788 refl:, 1.1328847337958123, 0.593939393939394, 42.288563952640665, 0.17897035579329929 refl:, 1.1328847337958123, 0.5850746268656717, 41.51545997434592, 0.1821986118653532 refl:, 1.1328847337958123, 0.5764705882352942, 40.77382641078957, 0.18516054972374968 refl:, 1.1328847337958123, 0.5681159420289855, 40.06152872079658, 0.18786984191653888 refl:, 1.1328847337958123, 0.56, 39.37664302675111, 0.19035023250949035 refl:, 1.1328847337958123, 0.552112676056338, 38.717428473384835, 0.19261983460766804 refl:, 1.1328847337958123, 0.5444444444444444, 38.08230410219336, 0.19469085220660431 refl:, 1.1328847337958123, 0.536986301369863, 37.4698293655767, 0.19658114960976503 refl:, 1.1328847337958123, 0.5297297297297298, 36.87868759977799, 0.1983001785222384 refl:, 1.1328847337958123, 0.5226666666666667, 36.307671922074825, 0.19985845493543128 refl:, 1.1328847337958123, 0.5157894736842105, 35.75567312877807, 0.2012671524485849 refl:, 1.1328847337958123, 0.509090909090909, 35.22166925577522, 0.2025355459067762 refl:, 1.1328847337958123, 0.5025641025641026, 34.70471652928195, 0.20367251569638592 refl:, 1.1328847337958123, 0.4962025316455696, 34.203941485939055, 0.2046830665024618 refl:, 1.1328847337958123, 0.49, 33.7185340819115, 0.20557518603611066 refl:, 1.1328847337958123, 0.4839506172839506, 33.247741642788455, 0.20635576080782017 refl:, 1.1328847337958123, 0.47804878048780486, 32.790863531764046, 0.20702823770770473 refl:, 1.1328847337958123, 0.47228915662650606, 32.347246434236666, 0.20759985728196045 refl:, 1.1328847337958123, 0.4666666666666667, 31.91628017368787, 0.2080745784918885 refl:, 1.1328847337958123, 0.4611764705882353, 31.497393987322308, 0.20845408884561908 refl:, 1.1328847337958123, 0.4558139534883721, 31.090053201105746, 0.20874542652135622 refl:, 1.1328847337958123, 0.4505747126436782, 30.693756253024553, 0.208954041733091 refl:, 1.1328847337958123, 0.4454545454545454, 30.308032020994055, 0.20908048492966044 refl:, 1.1328847337958123, 0.44044943820224725, 29.9324374181665, 0.2091260332602117 refl:, 1.1328847337958123, 0.43555555555555553, 29.5665552236732, 0.20909503137987662 refl:, 1.1328847337958123, 0.4307692307692308, 29.209992121269647, 0.2089911211638158 refl:, 1.1328847337958123, 0.4260869565217391, 28.86237692208816, 0.2088197236636275 refl:, 1.1328847337958123, 0.421505376344086, 28.523358950864424, 0.20858506970983773 refl:, 1.1328847337958123, 0.41702127659574467, 28.192606577687982, 0.20827627963637158 refl:, 1.1328847337958123, 0.4126315789473684, 27.86980587961545, 0.20789057239525374 refl:, 1.1328847337958123, 0.4083333333333333, 27.554659418441407, 0.20744945601962797 refl:, 1.1328847337958123, 0.4041237113402062, 27.246885122601544, 0.2069758959357581 refl:, 1.1328847337958123, 0.4, 26.946215262627685, 0.2064567398821898 ----------- Initializing structure... time for choose_chunkdivision = 3.91006e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00178099 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.14944315649486128 / 0.14944315649486128 = 1.0 field decay(t = 100.01): 4.987937476668207e-12 / 0.14944315649486128 = 3.337682095091267e-11 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.19617e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00460696 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.14944315644072229 / 0.14944315644072229 = 1.0 field decay(t = 100.01): 3.807851695152499e-11 / 0.14944315644072229 = 2.5480268122300477e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.1746157759823854, 0.8, 70.0, 0.008739943705374597 refl:, 1.1746157759823854, 0.784, 67.05783140972835, 0.024245008080781356 refl:, 1.1746157759823854, 0.7686274509803922, 64.53417775149677, 0.040309353659514456 refl:, 1.1746157759823854, 0.7538461538461539, 62.31059707934406, 0.05557589230945354 refl:, 1.1746157759823854, 0.739622641509434, 60.31629899672389, 0.06970273468511892 refl:, 1.1746157759823854, 0.7259259259259259, 58.50481213939229, 0.08255233289128476 refl:, 1.1746157759823854, 0.7127272727272727, 56.843594941553754, 0.09413192002673795 refl:, 1.1746157759823854, 0.7, 55.30875739675561, 0.10463895552192765 refl:, 1.1746157759823854, 0.6877192982456141, 53.88211612099151, 0.11407537277856838 refl:, 1.1746157759823854, 0.6758620689655173, 52.54943143507829, 0.12264564757179623 refl:, 1.1746157759823854, 0.664406779661017, 51.29929235382265, 0.130398421159986 refl:, 1.1746157759823854, 0.6533333333333333, 50.122379445158494, 0.13744278277017355 refl:, 1.1746157759823854, 0.6426229508196721, 49.010959419193355, 0.14384756806524446 refl:, 1.1746157759823854, 0.632258064516129, 47.958527855818005, 0.1496895084660299 refl:, 1.1746157759823854, 0.6222222222222222, 46.9595500055931, 0.1550183718219218 refl:, 1.1746157759823854, 0.6124999999999999, 46.00926848902442, 0.15990604921366627 refl:, 1.1746157759823854, 0.6030769230769231, 45.103557823882895, 0.16437489242095366 refl:, 1.1746157759823854, 0.593939393939394, 44.23881248012816, 0.16847755182701835 refl:, 1.1746157759823854, 0.5850746268656717, 43.41185942214757, 0.17223577319827998 refl:, 1.1746157759823854, 0.5764705882352942, 42.619888854495166, 0.17567909726641667 refl:, 1.1746157759823854, 0.5681159420289855, 41.860398715356425, 0.17883983394634828 refl:, 1.1746157759823854, 0.56, 41.131149701242414, 0.18173430538711705 refl:, 1.1746157759823854, 0.552112676056338, 40.430128463310105, 0.18439052329945838 refl:, 1.1746157759823854, 0.5444444444444444, 39.75551721886283, 0.18682349494977193 refl:, 1.1746157759823854, 0.536986301369863, 39.10566845306894, 0.18904694018036272 refl:, 1.1746157759823854, 0.5297297297297298, 38.479083699202135, 0.19107722550893155 refl:, 1.1746157759823854, 0.5226666666666667, 37.87439561623927, 0.19292653635113324 refl:, 1.1746157759823854, 0.5157894736842105, 37.29035275442276, 0.19460885764872243 refl:, 1.1746157759823854, 0.509090909090909, 36.72580652885898, 0.19613444229151364 refl:, 1.1746157759823854, 0.5025641025641026, 36.179700019842116, 0.19751218330354978 refl:, 1.1746157759823854, 0.4962025316455696, 35.651058294458764, 0.19875050713230683 refl:, 1.1746157759823854, 0.49, 35.138980002927255, 0.19985575200201458 refl:, 1.1746157759823854, 0.4839506172839506, 34.64263004924646, 0.20083821988564857 refl:, 1.1746157759823854, 0.47804878048780486, 34.161233172132036, 0.20170489252979734 refl:, 1.1746157759823854, 0.47228915662650606, 33.69406830116886, 0.20245705085154594 refl:, 1.1746157759823854, 0.4666666666666667, 33.24046357629389, 0.20310286904519861 refl:, 1.1746157759823854, 0.4611764705882353, 32.79979193741594, 0.20364737629446264 refl:, 1.1746157759823854, 0.4558139534883721, 32.37146720614281, 0.2040929385314762 refl:, 1.1746157759823854, 0.4505747126436782, 31.954940593960732, 0.20445102567587897 refl:, 1.1746157759823854, 0.4454545454545454, 31.549697581364857, 0.20471965623854066 refl:, 1.1746157759823854, 0.44044943820224725, 31.155255120816925, 0.2048962694167634 refl:, 1.1746157759823854, 0.43555555555555553, 30.771159123350074, 0.204990790346435 refl:, 1.1746157759823854, 0.4307692307692308, 30.396982194426784, 0.20500592934593215 refl:, 1.1746157759823854, 0.4260869565217391, 30.032321589495837, 0.20495712156136717 refl:, 1.1746157759823854, 0.421505376344086, 29.67679736376329, 0.2048311249915097 refl:, 1.1746157759823854, 0.41702127659574467, 29.33005069412437, 0.2046077238887242 refl:, 1.1746157759823854, 0.4126315789473684, 28.991742354111988, 0.20432648650919877 refl:, 1.1746157759823854, 0.4083333333333333, 28.66155132518973, 0.20399611274158247 refl:, 1.1746157759823854, 0.4041237113402062, 28.339173529827534, 0.20365744502245528 refl:, 1.1746157759823854, 0.4, 28.024320673604695, 0.20323355142072505 ----------- Initializing structure... time for choose_chunkdivision = 3.8147e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00180006 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.14355855111904067 / 0.14355855111904067 = 1.0 field decay(t = 100.01): 9.667874472848302e-12 / 0.14355855111904067 = 6.73444695386454e-11 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.19617e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00464296 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.14355855132715553 / 0.14355855132715553 = 1.0 field decay(t = 100.01): 5.780725980834341e-11 / 0.14355855132715553 = 4.0267374721974217e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.2074072828613354, 0.8, 75.0, 0.002256637590250018 refl:, 1.2074072828613354, 0.784, 71.19254305121969, 0.004043927155760674 refl:, 1.2074072828613354, 0.7686274509803922, 68.13230108184366, 0.017666840957282458 refl:, 1.2074072828613354, 0.7538461538461539, 65.5329128763873, 0.0333568804620509 refl:, 1.2074072828613354, 0.739622641509434, 63.25596068464814, 0.04865237243003264 refl:, 1.2074072828613354, 0.7259259259259259, 61.22160180379117, 0.06275076172785536 refl:, 1.2074072828613354, 0.7127272727272727, 59.378629248099124, 0.07575594386491534 refl:, 1.2074072828613354, 0.7, 57.691773022139905, 0.08757338578695935 refl:, 1.2074072828613354, 0.6877192982456141, 56.13545775731428, 0.09824417945974033 refl:, 1.2074072828613354, 0.6758620689655173, 54.69040263962552, 0.10794438844280667 refl:, 1.2074072828613354, 0.664406779661017, 53.3416232015308, 0.11670036595288552 refl:, 1.2074072828613354, 0.6533333333333333, 52.0771858356392, 0.12465519843605913 refl:, 1.2074072828613354, 0.6426229508196721, 50.88739414232009, 0.1318734346297872 refl:, 1.2074072828613354, 0.632258064516129, 49.76423655717566, 0.13846190112373255 refl:, 1.2074072828613354, 0.6222222222222222, 48.70099913863975, 0.14446187315076348 refl:, 1.2074072828613354, 0.6124999999999999, 47.69198665804526, 0.14995944852053927 refl:, 1.2074072828613354, 0.6030769230769231, 46.73231696062462, 0.1549820926980573 refl:, 1.2074072828613354, 0.593939393939394, 45.817766249324485, 0.1595835789325249 refl:, 1.2074072828613354, 0.5850746268656717, 44.94465059915575, 0.16380210131720174 refl:, 1.2074072828613354, 0.5764705882352942, 44.109733785683, 0.16766536391558942 refl:, 1.2074072828613354, 0.5681159420289855, 43.31015457776125, 0.17121616499277004 refl:, 1.2074072828613354, 0.56, 42.543368664412704, 0.17447361559923755 refl:, 1.2074072828613354, 0.552112676056338, 41.80710174663229, 0.17746600029237136 refl:, 1.2074072828613354, 0.5444444444444444, 41.099311260706344, 0.1802097333049947 refl:, 1.2074072828613354, 0.536986301369863, 40.418154855029364, 0.18272225722756205 refl:, 1.2074072828613354, 0.5297297297297298, 39.76196420912173, 0.18502166386205582 refl:, 1.2074072828613354, 0.5226666666666667, 39.12922312098109, 0.18712270741981965 refl:, 1.2074072828613354, 0.5157894736842105, 38.51854903626071, 0.18904224669910968 refl:, 1.2074072828613354, 0.509090909090909, 37.928677376424204, 0.19078915859924261 refl:, 1.2074072828613354, 0.5025641025641026, 37.35844816099257, 0.19237299906480013 refl:, 1.2074072828613354, 0.4962025316455696, 36.806794523773725, 0.19380527183932186 refl:, 1.2074072828613354, 0.49, 36.272732803336254, 0.19509371745070436 refl:, 1.2074072828613354, 0.4839506172839506, 35.755353950220794, 0.19624819583455538 refl:, 1.2074072828613354, 0.47804878048780486, 35.253816041990945, 0.1972786745002075 refl:, 1.2074072828613354, 0.47228915662650606, 34.76733773550362, 0.1981865257168301 refl:, 1.2074072828613354, 0.4666666666666667, 34.295192516152845, 0.19897642469530394 refl:, 1.2074072828613354, 0.4611764705882353, 33.83670362811667, 0.19965916344877177 refl:, 1.2074072828613354, 0.4558139534883721, 33.391239589169516, 0.2002382841140244 refl:, 1.2074072828613354, 0.4505747126436782, 32.95821020943854, 0.20072161188836213 refl:, 1.2074072828613354, 0.4454545454545454, 32.537063046367535, 0.2011093588503599 refl:, 1.2074072828613354, 0.44044943820224725, 32.12728023870626, 0.20139529444222973 refl:, 1.2074072828613354, 0.43555555555555553, 31.728375671037586, 0.2015949789880908 refl:, 1.2074072828613354, 0.4307692307692308, 31.33989242755132, 0.20171523571194905 refl:, 1.2074072828613354, 0.4260869565217391, 30.96140049976025, 0.20176487549210922 refl:, 1.2074072828613354, 0.421505376344086, 30.592494717856905, 0.20173258595752414 refl:, 1.2074072828613354, 0.41702127659574467, 30.23279287960783, 0.20158399392540655 refl:, 1.2074072828613354, 0.4126315789473684, 29.88193405422147, 0.20138903983141088 refl:, 1.2074072828613354, 0.4083333333333333, 29.539577041619506, 0.20117661380600185 refl:, 1.2074072828613354, 0.4041237113402062, 29.2053989700849, 0.2009185301554802 refl:, 1.2074072828613354, 0.4, 28.879094017427605, 0.2005570833595033 ----------- Initializing structure... time for choose_chunkdivision = 3.71933e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00173497 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.13921342756108773 / 0.13921342756108773 = 1.0 field decay(t = 100.01): 1.3144200426074514e-11 / 0.13921342756108773 = 9.441761945202273e-11 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 4.31538e-05 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.00474286 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.1392134279350055 / 0.1392134279350055 = 1.0 field decay(t = 100.01): 9.871965143904996e-11 / 0.1392134279350055 = 7.091244925391763e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.23100969126526, 0.8, 79.99999999999994, 0.028668309855084118 refl:, 1.23100969126526, 0.784, 74.82079670263325, 0.002131947030522914 refl:, 1.23100969126526, 0.7686274509803922, 71.11813585396203, 0.004117095029749368 refl:, 1.23100969126526, 0.7538461538461539, 68.12392493650243, 0.01750499625965161 refl:, 1.23100969126526, 0.739622641509434, 65.57213417097068, 0.032817770182446984 refl:, 1.23100969126526, 0.7259259259259259, 63.331955818023665, 0.04791694905496988 refl:, 1.23100969126526, 0.7127272727272727, 61.32721657609733, 0.06192830286716911 refl:, 1.23100969126526, 0.7, 59.508764013939604, 0.07458277694558078 refl:, 1.23100969126526, 0.6877192982456141, 57.84259890532987, 0.08617714917468623 refl:, 1.23100969126526, 0.6758620689655173, 56.30398980710504, 0.0966379797155269 refl:, 1.23100969126526, 0.664406779661017, 54.87424701948212, 0.10612489602225829 refl:, 1.23100969126526, 0.6533333333333333, 53.53881754577338, 0.11477448569944776 refl:, 1.23100969126526, 0.6426229508196721, 52.2860932123638, 0.1226196975060367 refl:, 1.23100969126526, 0.632258064516129, 51.10662968356778, 0.12979595071970462 refl:, 1.23100969126526, 0.6222222222222222, 49.992614952243095, 0.1363423839815458 refl:, 1.23100969126526, 0.6124999999999999, 48.937495982303304, 0.14232506775291698 refl:, 1.23100969126526, 0.6030769230769231, 47.93570930843333, 0.14779669116407043 refl:, 1.23100969126526, 0.593939393939394, 46.98248211140486, 0.15279783449569376 refl:, 1.23100969126526, 0.5850746268656717, 46.07368235996428, 0.1573761848030311 refl:, 1.23100969126526, 0.5764705882352942, 45.205703916101704, 0.16157162995599855 refl:, 1.23100969126526, 0.5681159420289855, 44.37537706820699, 0.16542676815817645 refl:, 1.23100969126526, 0.56, 43.57989789517631, 0.16896454842958672 refl:, 1.23100969126526, 0.552112676056338, 42.81677180338089, 0.1722206134003931 refl:, 1.23100969126526, 0.5444444444444444, 42.083767886770914, 0.17520135452911917 refl:, 1.23100969126526, 0.536986301369863, 41.378881661291835, 0.17793593217612066 refl:, 1.23100969126526, 0.5297297297297298, 40.70030435652892, 0.1804419317536035 refl:, 1.23100969126526, 0.5226666666666667, 40.046397397863686, 0.18273421151574976 refl:, 1.23100969126526, 0.5157894736842105, 39.41567103836258, 0.18483834118872827 refl:, 1.23100969126526, 0.509090909090909, 38.806766338771496, 0.18675576413095382 refl:, 1.23100969126526, 0.5025641025641026, 38.218439871701385, 0.18849963174434395 refl:, 1.23100969126526, 0.4962025316455696, 37.64955065969264, 0.19008108230487103 refl:, 1.23100969126526, 0.49, 37.099048958377544, 0.19150758216597197 refl:, 1.23100969126526, 0.4839506172839506, 36.56596657389412, 0.19279776266416668 refl:, 1.23100969126526, 0.47804878048780486, 36.049408464086156, 0.19395439242202947 refl:, 1.23100969126526, 0.47228915662650606, 35.54854542021205, 0.19498180503114848 refl:, 1.23100969126526, 0.4666666666666667, 35.062607663064846, 0.19588433555808407 refl:, 1.23100969126526, 0.4611764705882353, 34.59087921692149, 0.19666557505940196 refl:, 1.23100969126526, 0.4558139534883721, 34.13269294834013, 0.19734832119524592 refl:, 1.23100969126526, 0.4505747126436782, 33.6874261758209, 0.19793121915708817 refl:, 1.23100969126526, 0.4454545454545454, 33.25449677173539, 0.19840657303015277 refl:, 1.23100969126526, 0.44044943820224725, 32.83335969046958, 0.19877980002307316 refl:, 1.23100969126526, 0.43555555555555553, 32.42350386700291, 0.19904606256565613 refl:, 1.23100969126526, 0.4307692307692308, 32.0244494386133, 0.1992480800037741 refl:, 1.23100969126526, 0.4260869565217391, 31.63574524940934, 0.19939085390804745 refl:, 1.23100969126526, 0.421505376344086, 31.2569666032255, 0.19942098648283285 refl:, 1.23100969126526, 0.41702127659574467, 30.887713235292672, 0.19934811719121265 refl:, 1.23100969126526, 0.4126315789473684, 30.527607477190394, 0.1991752884538809 refl:, 1.23100969126526, 0.4083333333333333, 30.176292593038497, 0.19904014318467667 refl:, 1.23100969126526, 0.4041237113402062, 29.83343126780691, 0.19894516024468498 refl:, 1.23100969126526, 0.4, 29.498704231103652, 0.19854708072325522 Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.007622 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.23946072431081186 / 0.23946072431081186 = 1.0 field decay(t = 100.01): 2.143367304660103e-13 / 0.23946072431081186 = 8.950809410724429e-13 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000111818 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.0173719 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.23946072470582727 / 0.23946072470582727 = 1.0 field decay(t = 100.01): 2.2924410648414684e-11 / 0.23946072470582727 = 9.573348897434795e-11 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.4275251791570859, 0.8, 20.0, 0.27287026974445927 refl:, 0.4275251791570859, 0.784, 19.583468236198428, 0.2732024101218738 refl:, 0.4275251791570859, 0.7686274509803922, 19.184283570762837, 0.2734622444200681 refl:, 0.4275251791570859, 0.7538461538461539, 18.80136282780684, 0.27365329059595683 refl:, 0.4275251791570859, 0.739622641509434, 18.433712921276232, 0.27378076185161915 refl:, 0.4275251791570859, 0.7259259259259259, 18.080421505473552, 0.27384851240423785 refl:, 0.4275251791570859, 0.7127272727272727, 17.74064878668832, 0.2738584846379507 refl:, 0.4275251791570859, 0.7, 17.413620328364388, 0.27381376974644933 refl:, 0.4275251791570859, 0.6877192982456141, 17.098620709730763, 0.27371850538265774 refl:, 0.4275251791570859, 0.6758620689655173, 16.794987920274295, 0.2735750966910372 refl:, 0.4275251791570859, 0.664406779661017, 16.50210839086207, 0.2733847949229058 refl:, 0.4275251791570859, 0.6533333333333333, 16.21941257752279, 0.2731503266875855 refl:, 0.4275251791570859, 0.6426229508196721, 15.946371026493122, 0.2728749017452174 refl:, 0.4275251791570859, 0.632258064516129, 15.682490859619808, 0.2725596258354739 refl:, 0.4275251791570859, 0.6222222222222222, 15.427312627971409, 0.2722049153191827 refl:, 0.4275251791570859, 0.6124999999999999, 15.18040748886694, 0.2718130795196188 refl:, 0.4275251791570859, 0.6030769230769231, 14.941374667722823, 0.2713863599041548 refl:, 0.4275251791570859, 0.593939393939394, 14.709839171355561, 0.2709245158308668 refl:, 0.4275251791570859, 0.5850746268656717, 14.485449723819883, 0.2704275378957415 refl:, 0.4275251791570859, 0.5764705882352942, 14.267876899642369, 0.26989771202120866 refl:, 0.4275251791570859, 0.5681159420289855, 14.056811432538773, 0.26933628796739983 refl:, 0.4275251791570859, 0.56, 13.851962680467771, 0.26874214027712073 refl:, 0.4275251791570859, 0.552112676056338, 13.653057230248155, 0.2681155761222093 refl:, 0.4275251791570859, 0.5444444444444444, 13.459837627011849, 0.2674588673074579 refl:, 0.4275251791570859, 0.536986301369863, 13.272061215531336, 0.2667723151469476 refl:, 0.4275251791570859, 0.5297297297297298, 13.089499081989581, 0.26605454446111854 refl:, 0.4275251791570859, 0.5226666666666667, 12.911935086088024, 0.2653063511402134 refl:, 0.4275251791570859, 0.5157894736842105, 12.73916497454351, 0.2645297343354968 refl:, 0.4275251791570859, 0.509090909090909, 12.570995568032437, 0.2637242645202348 refl:, 0.4275251791570859, 0.5025641025641026, 12.407244014521105, 0.2628885752055156 refl:, 0.4275251791570859, 0.4962025316455696, 12.247737102692609, 0.26202377900308255 refl:, 0.4275251791570859, 0.49, 12.092310629857906, 0.2611315306153532 refl:, 0.4275251791570859, 0.4839506172839506, 11.94080881933422, 0.2602107728685866 refl:, 0.4275251791570859, 0.47804878048780486, 11.793083782798943, 0.25926014291141586 refl:, 0.4275251791570859, 0.47228915662650606, 11.648995023590732, 0.25828106391345235 refl:, 0.4275251791570859, 0.4666666666666667, 11.508408977339421, 0.25727477704208235 refl:, 0.4275251791570859, 0.4611764705882353, 11.37119858666983, 0.2562395153461572 refl:, 0.4275251791570859, 0.4558139534883721, 11.23724290704701, 0.2551740956437644 refl:, 0.4275251791570859, 0.4505747126436782, 11.106426741117268, 0.2540803028372474 refl:, 0.4275251791570859, 0.4454545454545454, 10.978640299154753, 0.25295873317085843 refl:, 0.4275251791570859, 0.44044943820224725, 10.853778883451254, 0.2518070732573948 refl:, 0.4275251791570859, 0.43555555555555553, 10.731742594690267, 0.2506247305646683 refl:, 0.4275251791570859, 0.4307692307692308, 10.61243605852877, 0.249413825082604 refl:, 0.4275251791570859, 0.4260869565217391, 10.495768170772996, 0.2481743775271508 refl:, 0.4275251791570859, 0.421505376344086, 10.381651859681224, 0.24690419443990608 refl:, 0.4275251791570859, 0.41702127659574467, 10.270003864057921, 0.2456037441531416 refl:, 0.4275251791570859, 0.4126315789473684, 10.160744525922071, 0.24427554305470142 refl:, 0.4275251791570859, 0.4083333333333333, 10.053797596639106, 0.2429193465329389 refl:, 0.4275251791570859, 0.4041237113402062, 9.949090055502005, 0.24153322642643202 refl:, 0.4275251791570859, 0.4, 9.846551939834079, 0.2401180039660742 ----------- Initializing structure... time for choose_chunkdivision = 0.000110865 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00934601 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.23217567515545764 / 0.23217567515545764 = 1.0 field decay(t = 100.01): 3.1895846890533654e-13 / 0.23217567515545764 = 1.3737807317315728e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000108957 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.016118 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.23217567572617967 / 0.23217567572617967 = 1.0 field decay(t = 100.01): 2.511386584143695e-11 / 0.23217567572617967 = 1.0816751480484724e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.5282728271758743, 0.8, 25.0, 0.2601481976279196 refl:, 0.5282728271758743, 0.784, 24.46679999189225, 0.2610696372423531 refl:, 0.5282728271758743, 0.7686274509803922, 23.95662859236144, 0.2618770325163737 refl:, 0.5282728271758743, 0.7538461538461539, 23.467975961566616, 0.2625788276586943 refl:, 0.5282728271758743, 0.739622641509434, 22.99946566139384, 0.2631817618626755 refl:, 0.5282728271758743, 0.7259259259259259, 22.54983979533518, 0.2636926798419164 refl:, 0.5282728271758743, 0.7127272727272727, 22.117946146436154, 0.2641189902441729 refl:, 0.5282728271758743, 0.7, 21.70272699899909, 0.2644658127354963 refl:, 0.5282728271758743, 0.6877192982456141, 21.303209386217617, 0.26473707694833154 refl:, 0.5282728271758743, 0.6758620689655173, 20.91849655102253, 0.26493869102120726 refl:, 0.5282728271758743, 0.664406779661017, 20.54776044367494, 0.26507564513451415 refl:, 0.5282728271758743, 0.6533333333333333, 20.19023510896709, 0.2651500052068409 refl:, 0.5282728271758743, 0.6426229508196721, 19.84521083974841, 0.26516501686911953 refl:, 0.5282728271758743, 0.632258064516129, 19.51202899301187, 0.26512536746580095 refl:, 0.5282728271758743, 0.6222222222222222, 19.190077380825585, 0.26503310672375896 refl:, 0.5282728271758743, 0.6124999999999999, 18.878786161658894, 0.26488908132632305 refl:, 0.5282728271758743, 0.6030769230769231, 18.577624168665142, 0.26469630209606193 refl:, 0.5282728271758743, 0.593939393939394, 18.28609562066815, 0.26445767784614754 refl:, 0.5282728271758743, 0.5850746268656717, 18.003737169291743, 0.26417370471931184 refl:, 0.5282728271758743, 0.5764705882352942, 17.730115242139682, 0.26384505077481407 refl:, 0.5282728271758743, 0.5681159420289855, 17.46482364739322, 0.26347427561024556 refl:, 0.5282728271758743, 0.56, 17.207481409818527, 0.26306325708398504 refl:, 0.5282728271758743, 0.552112676056338, 16.957730812108316, 0.262611771262809 refl:, 0.5282728271758743, 0.5444444444444444, 16.715235618835717, 0.2621203557572436 refl:, 0.5282728271758743, 0.536986301369863, 16.479679463167955, 0.2615915488530281 refl:, 0.5282728271758743, 0.5297297297297298, 16.250764378950514, 0.26102646172552363 refl:, 0.5282728271758743, 0.5226666666666667, 16.028209462892256, 0.26042407309811766 refl:, 0.5282728271758743, 0.5157894736842105, 15.811749653412038, 0.25978529806673434 refl:, 0.5282728271758743, 0.509090909090909, 15.60113461429131, 0.2591127130188696 refl:, 0.5282728271758743, 0.5025641025641026, 15.396127712651408, 0.25840617419006323 refl:, 0.5282728271758743, 0.4962025316455696, 15.196505081970068, 0.2576643632602327 refl:, 0.5282728271758743, 0.49, 15.002054761894165, 0.2568889467859571 refl:, 0.5282728271758743, 0.4839506172839506, 14.81257590751694, 0.2560818324629057 refl:, 0.5282728271758743, 0.47804878048780486, 14.627878061586326, 0.2552417285733339 refl:, 0.5282728271758743, 0.47228915662650606, 14.44778048381156, 0.25436772779199046 refl:, 0.5282728271758743, 0.4666666666666667, 14.272111532051735, 0.2534618054882236 refl:, 0.5282728271758743, 0.4611764705882353, 14.100708090713411, 0.2525248442883521 refl:, 0.5282728271758743, 0.4558139534883721, 13.933415042164041, 0.25155500702726324 refl:, 0.5282728271758743, 0.4505747126436782, 13.770084777392734, 0.25055182001671567 refl:, 0.5282728271758743, 0.4454545454545454, 13.610576742526133, 0.24951717022490025 refl:, 0.5282728271758743, 0.44044943820224725, 13.454757018141454, 0.24845117408859577 refl:, 0.5282728271758743, 0.43555555555555553, 13.30249792861589, 0.24735169966979012 refl:, 0.5282728271758743, 0.4307692307692308, 13.153677679016674, 0.2462187296769865 refl:, 0.5282728271758743, 0.4260869565217391, 13.008180017272103, 0.24505436171970374 refl:, 0.5282728271758743, 0.421505376344086, 12.865893919575512, 0.24385825418409707 refl:, 0.5282728271758743, 0.41702127659574467, 12.726713297162874, 0.2426284381212972 refl:, 0.5282728271758743, 0.4126315789473684, 12.590536722774571, 0.24136614460408662 refl:, 0.5282728271758743, 0.4083333333333333, 12.457267175263851, 0.24007405600657153 refl:, 0.5282728271758743, 0.4041237113402062, 12.326811800951434, 0.238751417728497 refl:, 0.5282728271758743, 0.4, 12.199081690448809, 0.2373967326334257 ----------- Initializing structure... time for choose_chunkdivision = 0.000108004 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00868416 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.22375100718842 / 0.22375100718842 = 1.0 field decay(t = 100.01): 4.4054832270861915e-13 / 0.22375100718842 = 1.9689221883038714e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000108957 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.0171471 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.22375100783149995 / 0.22375100783149995 = 1.0 field decay(t = 100.01): 2.604590852088861e-11 / 0.22375100783149995 = 1.1640577074183724e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.6249999999999999, 0.8, 29.999999999999993, 0.2442050337697129 refl:, 0.6249999999999999, 0.784, 29.34058157502373, 0.24592026654412946 refl:, 0.6249999999999999, 0.7686274509803922, 28.711017527148794, 0.2474593095717218 refl:, 0.6249999999999999, 0.7538461538461539, 28.10922128260952, 0.24883736533494075 refl:, 0.6249999999999999, 0.739622641509434, 27.53330580109674, 0.25006920079221073 refl:, 0.6249999999999999, 0.7259259259259259, 26.98155921981659, 0.25116546178265714 refl:, 0.6249999999999999, 0.7127272727272727, 26.452424118557172, 0.2521360799590939 refl:, 0.6249999999999999, 0.7, 25.944479772370002, 0.2529918131418104 refl:, 0.6249999999999999, 0.6877192982456141, 25.45642688540588, 0.2537410511682661 refl:, 0.6249999999999999, 0.6758620689655173, 24.98707439783817, 0.2543901348468849 refl:, 0.6249999999999999, 0.664406779661017, 24.53532803474182, 0.2549460427718285 refl:, 0.6249999999999999, 0.6533333333333333, 24.10018032643974, 0.25541571682136277 refl:, 0.6249999999999999, 0.6426229508196721, 23.680701877992792, 0.25580383669879514 refl:, 0.6249999999999999, 0.632258064516129, 23.276033704033082, 0.2561140687592272 refl:, 0.6249999999999999, 0.6222222222222222, 22.885380476158563, 0.2563515425795029 refl:, 0.6249999999999999, 0.6124999999999999, 22.508004555237044, 0.2565208309064006 refl:, 0.6249999999999999, 0.6030769230769231, 22.14322070144805, 0.2566238743420462 refl:, 0.6249999999999999, 0.593939393939394, 21.79039137167391, 0.25666334314703 refl:, 0.6249999999999999, 0.5850746268656717, 21.448922527676533, 0.256643740595326 refl:, 0.6249999999999999, 0.5764705882352942, 21.118259889941786, 0.25656730578199227 refl:, 0.6249999999999999, 0.5681159420289855, 20.797885581592283, 0.25643465337781024 refl:, 0.6249999999999999, 0.56, 20.487315114722662, 0.256248873599916 refl:, 0.6249999999999999, 0.552112676056338, 20.186094678183196, 0.25601329709669296 refl:, 0.6249999999999999, 0.5444444444444444, 19.89379869145801, 0.25572829145150733 refl:, 0.6249999999999999, 0.536986301369863, 19.610027594036424, 0.255394616552182 refl:, 0.6249999999999999, 0.5297297297297298, 19.33440584370919, 0.2550154350698738 refl:, 0.6249999999999999, 0.5226666666666667, 19.066580100655866, 0.25459252830565576 refl:, 0.6249999999999999, 0.5157894736842105, 18.806217577124535, 0.25412542518907427 refl:, 0.6249999999999999, 0.509090909090909, 18.553004535020655, 0.2536152551480315 refl:, 0.6249999999999999, 0.5025641025641026, 18.306644915884704, 0.2530647648651237 refl:, 0.6249999999999999, 0.4962025316455696, 18.066859089603533, 0.25247451149799194 refl:, 0.6249999999999999, 0.49, 17.833382709813034, 0.251843602956665 refl:, 0.6249999999999999, 0.4839506172839506, 17.605965665348197, 0.25117351956545597 refl:, 0.6249999999999999, 0.47804878048780486, 17.38437111831227, 0.2504665488589744 refl:, 0.6249999999999999, 0.47228915662650606, 17.168374620396104, 0.2497220477831976 refl:, 0.6249999999999999, 0.4666666666666667, 16.957763300004142, 0.24893902021832587 refl:, 0.6249999999999999, 0.4611764705882353, 16.752335113553887, 0.24811948817045582 refl:, 0.6249999999999999, 0.4558139534883721, 16.551898155026578, 0.24726486187036864 refl:, 0.6249999999999999, 0.4505747126436782, 16.35627001847215, 0.24637335126615545 refl:, 0.6249999999999999, 0.4454545454545454, 16.165277208722518, 0.24544445758208924 refl:, 0.6249999999999999, 0.44044943820224725, 15.978754596053776, 0.24448044797494067 refl:, 0.6249999999999999, 0.43555555555555553, 15.796544910968182, 0.24348141406546195 refl:, 0.6249999999999999, 0.4307692307692308, 15.618498275648548, 0.24244508607372425 refl:, 0.6249999999999999, 0.4260869565217391, 15.44447176897603, 0.24137187009594935 refl:, 0.6249999999999999, 0.421505376344086, 15.274329022304034, 0.2402639528246127 refl:, 0.6249999999999999, 0.41702127659574467, 15.107939843449019, 0.23912065251888115 refl:, 0.6249999999999999, 0.4126315789473684, 14.94517986659885, 0.23794038016947253 refl:, 0.6249999999999999, 0.4083333333333333, 14.78593022605342, 0.23672496742401775 refl:, 0.6249999999999999, 0.4041237113402062, 14.630077251904048, 0.23547697979370008 refl:, 0.6249999999999999, 0.4, 14.477512185929921, 0.23419574038931623 ----------- Initializing structure... time for choose_chunkdivision = 0.000112057 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00802183 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.2144470639002981 / 0.2144470639002981 = 1.0 field decay(t = 100.01): 5.749510870152026e-13 / 0.2144470639002981 = 2.6810863089387505e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000110865 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.016397 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.21444706445458073 / 0.21444706445458073 = 1.0 field decay(t = 100.01): 2.737725032991902e-11 / 0.21444706445458073 = 1.276643744205575e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.7169705454388076, 0.8, 35.0, 0.22482605012853205 refl:, 0.7169705454388076, 0.784, 34.201491482224874, 0.2275879297659293 refl:, 0.7169705454388076, 0.7686274509803922, 33.4413597291606, 0.2300861771421993 refl:, 0.7169705454388076, 0.7538461538461539, 32.71669424141867, 0.23234663839124975 refl:, 0.7169705454388076, 0.739622641509434, 32.02489215108529, 0.2343911194485937 refl:, 0.7169705454388076, 0.7259259259259259, 31.363616172201723, 0.23623793864435802 refl:, 0.7169705454388076, 0.7127272727272727, 30.73075960358516, 0.23790584943527182 refl:, 0.7169705454388076, 0.7, 30.12441698853809, 0.2394088473062314 refl:, 0.7169705454388076, 0.6877192982456141, 29.542859352190995, 0.24075882698160853 refl:, 0.7169705454388076, 0.6758620689655173, 28.984513173223203, 0.24196961848162934 refl:, 0.7169705454388076, 0.664406779661017, 28.447942424882076, 0.24305087670480616 refl:, 0.7169705454388076, 0.6533333333333333, 27.931833156170086, 0.24400999332366932 refl:, 0.7169705454388076, 0.6426229508196721, 27.43498018882287, 0.24485687989558083 refl:, 0.7169705454388076, 0.632258064516129, 26.95627558715663, 0.24559913215899493 refl:, 0.7169705454388076, 0.6222222222222222, 26.49469862174186, 0.246241502598612 refl:, 0.7169705454388076, 0.6124999999999999, 26.049306998353728, 0.24679042717466868 refl:, 0.7169705454388076, 0.6030769230769231, 25.619229163859046, 0.2472523465378225 refl:, 0.7169705454388076, 0.593939393939394, 25.20365753294407, 0.24763123330505407 refl:, 0.7169705454388076, 0.5850746268656717, 24.80184250561243, 0.24793067603660995 refl:, 0.7169705454388076, 0.5764705882352942, 24.41308716651858, 0.24815554480891827 refl:, 0.7169705454388076, 0.5681159420289855, 24.036742574466203, 0.24831017863638788 refl:, 0.7169705454388076, 0.56, 23.672203564580393, 0.24839671574872801 refl:, 0.7169705454388076, 0.552112676056338, 23.318904997368072, 0.24841768322180216 refl:, 0.7169705454388076, 0.5444444444444444, 22.976318398592397, 0.24837749388988659 refl:, 0.7169705454388076, 0.536986301369863, 22.643948941980693, 0.24827859346807957 refl:, 0.7169705454388076, 0.5297297297297298, 22.321332733561132, 0.24812142120897318 refl:, 0.7169705454388076, 0.5226666666666667, 22.008034362119382, 0.24790902685773963 refl:, 0.7169705454388076, 0.5157894736842105, 21.703644685073712, 0.24764492885992392 refl:, 0.7169705454388076, 0.509090909090909, 21.407778823139967, 0.24732917319523862 refl:, 0.7169705454388076, 0.5025641025641026, 21.120074340620857, 0.24696243326222161 refl:, 0.7169705454388076, 0.4962025316455696, 20.840189591108683, 0.24654821236606225 refl:, 0.7169705454388076, 0.49, 20.56780221092012, 0.246087935342834 refl:, 0.7169705454388076, 0.4839506172839506, 20.302607744753697, 0.2455807741452668 refl:, 0.7169705454388076, 0.47804878048780486, 20.044318389931625, 0.24502832199568722 refl:, 0.7169705454388076, 0.47228915662650606, 19.79266184720379, 0.24443328421700983 refl:, 0.7169705454388076, 0.4666666666666667, 19.54738026749198, 0.24379556106465222 refl:, 0.7169705454388076, 0.4611764705882353, 19.308229285168313, 0.24311440973651832 refl:, 0.7169705454388076, 0.4558139534883721, 19.07497712952083, 0.2423917032614313 refl:, 0.7169705454388076, 0.4505747126436782, 18.847403806983923, 0.24162921270758458 refl:, 0.7169705454388076, 0.4454545454545454, 18.62530034751981, 0.24082577322829152 refl:, 0.7169705454388076, 0.44044943820224725, 18.408468109247096, 0.23998068018388435 refl:, 0.7169705454388076, 0.43555555555555553, 18.196718136035763, 0.23909605307915388 refl:, 0.7169705454388076, 0.4307692307692308, 17.98987056333767, 0.23817255818097147 refl:, 0.7169705454388076, 0.4260869565217391, 17.787754068005896, 0.23720797538105634 refl:, 0.7169705454388076, 0.421505376344086, 17.590205358285754, 0.236202366418755 refl:, 0.7169705454388076, 0.41702127659574467, 17.39706870053973, 0.23515820136339463 refl:, 0.7169705454388076, 0.4126315789473684, 17.20819547960629, 0.2340750107196232 refl:, 0.7169705454388076, 0.4083333333333333, 17.02344378999232, 0.2329510888033509 refl:, 0.7169705454388076, 0.4041237113402062, 16.842678055366356, 0.23178864386805378 refl:, 0.7169705454388076, 0.4, 16.665768674058118, 0.2305906700422077 ----------- Initializing structure... time for choose_chunkdivision = 0.000108957 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00882411 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.20454876789491885 / 0.20454876789491885 = 1.0 field decay(t = 100.01): 7.178635175771892e-13 / 0.20454876789491885 = 3.509498125874663e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000108004 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.018625 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.20454876813558592 / 0.20454876813558592 = 1.0 field decay(t = 100.01): 2.8174572450318613e-11 / 0.20454876813558592 = 1.3774012284270024e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.8034845121081741, 0.8, 39.99999999999999, 0.20179462162966672 refl:, 0.8034845121081741, 0.784, 39.04509528455022, 0.2059250719323383 refl:, 0.8034845121081741, 0.7686274509803922, 38.139646206191365, 0.20966397204966036 refl:, 0.8034845121081741, 0.7538461538461539, 37.27949726994566, 0.21305427397378518 refl:, 0.8034845121081741, 0.739622641509434, 36.46098976062103, 0.21613269723029083 refl:, 0.8034845121081741, 0.7259259259259259, 35.680884053707004, 0.21892741270212848 refl:, 0.8034845121081741, 0.7127272727272727, 34.93629688867566, 0.22146683922801672 refl:, 0.8034845121081741, 0.7, 34.22465020933728, 0.22377279980470224 refl:, 0.8034845121081741, 0.6877192982456141, 33.54362905484555, 0.22586624586940138 refl:, 0.8034845121081741, 0.6758620689655173, 32.89114661043291, 0.22776517254801126 refl:, 0.8034845121081741, 0.664406779661017, 32.265314978971055, 0.2294835019692795 refl:, 0.8034845121081741, 0.6533333333333333, 31.664420565727855, 0.2310373437439053 refl:, 0.8034845121081741, 0.6426229508196721, 31.086903214629274, 0.23243828166083527 refl:, 0.8034845121081741, 0.632258064516129, 30.531338419096144, 0.23369542692345546 refl:, 0.8034845121081741, 0.6222222222222222, 29.996422070856152, 0.2348212532271984 refl:, 0.8034845121081741, 0.6124999999999999, 29.480957317802616, 0.23582404434823864 refl:, 0.8034845121081741, 0.6030769230769231, 28.983843185365757, 0.23670970324155213 refl:, 0.8034845121081741, 0.593939393939394, 28.504064681021994, 0.2374874618810027 refl:, 0.8034845121081741, 0.5850746268656717, 28.0406841528977, 0.23816414677923553 refl:, 0.8034845121081741, 0.5764705882352942, 27.592833714170844, 0.2387438174110003 refl:, 0.8034845121081741, 0.5681159420289855, 27.159708577552564, 0.239232508084756 refl:, 0.8034845121081741, 0.56, 26.740561170354134, 0.23963612123153596 refl:, 0.8034845121081741, 0.552112676056338, 26.33469592188663, 0.23995828312103898 refl:, 0.8034845121081741, 0.5444444444444444, 25.94146463225064, 0.24020224568501714 refl:, 0.8034845121081741, 0.536986301369863, 25.560262345758233, 0.24037244343095185 refl:, 0.8034845121081741, 0.5297297297297298, 25.190523663916437, 0.2404730149043355 refl:, 0.8034845121081741, 0.5226666666666667, 24.831719442577725, 0.2405057199519815 refl:, 0.8034845121081741, 0.5157894736842105, 24.483353825914165, 0.24047273277761147 refl:, 0.8034845121081741, 0.509090909090909, 24.144961576600384, 0.2403784518025041 refl:, 0.8034845121081741, 0.5025641025641026, 23.816105667237974, 0.24022489861203003 refl:, 0.8034845121081741, 0.4962025316455696, 23.496375102814216, 0.24001212440412875 refl:, 0.8034845121081741, 0.49, 23.185382948015043, 0.2397434129524054 refl:, 0.8034845121081741, 0.4839506172839506, 22.882764536632674, 0.2394219863711472 refl:, 0.8034845121081741, 0.47804878048780486, 22.58817584322306, 0.23904730465890137 refl:, 0.8034845121081741, 0.47228915662650606, 22.301291999661412, 0.23862033792003395 refl:, 0.8034845121081741, 0.4666666666666667, 22.021805941382443, 0.23814466238694304 refl:, 0.8034845121081741, 0.4611764705882353, 21.749427169932794, 0.23762099874005177 refl:, 0.8034845121081741, 0.4558139534883721, 21.483880620051718, 0.23704845905896132 refl:, 0.8034845121081741, 0.4505747126436782, 21.224905620871482, 0.2364290006894749 refl:, 0.8034845121081741, 0.4454545454545454, 20.972254942022712, 0.23576489053070165 refl:, 0.8034845121081741, 0.44044943820224725, 20.725693916468796, 0.23505540516189377 refl:, 0.8034845121081741, 0.43555555555555553, 20.48499963280006, 0.23429986728741037 refl:, 0.8034845121081741, 0.4307692307692308, 20.249960190511292, 0.23350021386178307 refl:, 0.8034845121081741, 0.4260869565217391, 20.020374012481035, 0.23265749904183425 refl:, 0.8034845121081741, 0.421505376344086, 19.79604920948226, 0.23176986721335244 refl:, 0.8034845121081741, 0.41702127659574467, 19.576802992091316, 0.23083693098169072 refl:, 0.8034845121081741, 0.4126315789473684, 19.362461125837058, 0.22986107450928206 refl:, 0.8034845121081741, 0.4083333333333333, 19.15285742585141, 0.2288423958877941 refl:, 0.8034845121081741, 0.4041237113402062, 18.9478332876545, 0.22777913143869027 refl:, 0.8034845121081741, 0.4, 18.747237251037504, 0.22667336344348024 ----------- Initializing structure... time for choose_chunkdivision = 0.000111103 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.006598 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.1943564114319847 / 0.1943564114319847 = 1.0 field decay(t = 100.01): 8.613369505850488e-13 / 0.1943564114319847 = 4.431739319731548e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000108004 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.0163691 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.1943564112167818 / 0.1943564112167818 = 1.0 field decay(t = 100.01): 2.929219089217826e-11 / 0.1943564112167818 = 1.5071378766870859e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.8838834764831843, 0.8, 44.99999999999999, 0.1749506546010766 refl:, 0.8838834764831843, 0.784, 43.865246854678944, 0.18084614262894314 refl:, 0.8838834764831843, 0.7686274509803922, 42.79498686880298, 0.18616490888473136 refl:, 0.8838834764831843, 0.7538461538461539, 41.78306957363623, 0.1909834896778329 refl:, 0.8838834764831843, 0.739622641509434, 40.82419653278361, 0.1953532231998348 refl:, 0.8838834764831843, 0.7259259259259259, 39.913765803587594, 0.19932045676242707 refl:, 0.8838834764831843, 0.7127272727272727, 39.047751315763335, 0.20293092176670602 refl:, 0.8838834764831843, 0.7, 38.2226079814105, 0.20621465898807328 refl:, 0.8838834764831843, 0.6877192982456141, 37.435196089128915, 0.20920558346668153 refl:, 0.8838834764831843, 0.6758620689655173, 36.682720369951866, 0.21192936258982145 refl:, 0.8838834764831843, 0.664406779661017, 35.962680378531154, 0.2144091185242741 refl:, 0.8838834764831843, 0.6533333333333333, 35.272829708778914, 0.21666666414632094 refl:, 0.8838834764831843, 0.6426229508196721, 34.61114218453038, 0.2187182335956619 refl:, 0.8838834764831843, 0.632258064516129, 33.975783613579615, 0.22058228827610923 refl:, 0.8838834764831843, 0.6222222222222222, 33.36508802077569, 0.22227211810191635 refl:, 0.8838834764831843, 0.6124999999999999, 32.77753751829758, 0.22379917822267162 refl:, 0.8838834764831843, 0.6030769230769231, 32.21174515294742, 0.225178014297656 refl:, 0.8838834764831843, 0.593939393939394, 31.666440208036935, 0.22641701701757366 refl:, 0.8838834764831843, 0.5850746268656717, 31.14045554291693, 0.22752429105750807 refl:, 0.8838834764831843, 0.5764705882352942, 30.6327166347442, 0.22851103000945724 refl:, 0.8838834764831843, 0.5681159420289855, 30.142232050688204, 0.2293833309292736 refl:, 0.8838834764831843, 0.56, 29.66808512880701, 0.23014665030198314 refl:, 0.8838834764831843, 0.552112676056338, 29.20942668547734, 0.2308091072410071 refl:, 0.8838834764831843, 0.5444444444444444, 28.765468598924116, 0.23137604508998047 refl:, 0.8838834764831843, 0.536986301369863, 28.33547814384704, 0.23185120425684555 refl:, 0.8838834764831843, 0.5297297297297298, 27.91877297273401, 0.23223967284122596 refl:, 0.8838834764831843, 0.5226666666666667, 27.514716656212954, 0.23254649427327514 refl:, 0.8838834764831843, 0.5157894736842105, 27.12271470851516, 0.23277497605931805 refl:, 0.8838834764831843, 0.509090909090909, 26.742211035417114, 0.2329274877105852 refl:, 0.8838834764831843, 0.5025641025641026, 26.372684751370738, 0.23300817835325607 refl:, 0.8838834764831843, 0.4962025316455696, 26.013647320299032, 0.23302096475520206 refl:, 0.8838834764831843, 0.49, 25.66463998102006, 0.23296661309174865 refl:, 0.8838834764831843, 0.4839506172839506, 25.32523142370283, 0.2328473819749458 refl:, 0.8838834764831843, 0.47804878048780486, 24.99501568834094, 0.23266777547883646 refl:, 0.8838834764831843, 0.47228915662650606, 24.673610260104642, 0.2324286704369094 refl:, 0.8838834764831843, 0.4666666666666667, 24.360654339720863, 0.23213016162250638 refl:, 0.8838834764831843, 0.4611764705882353, 24.055807269832574, 0.23177609164136181 refl:, 0.8838834764831843, 0.4558139534883721, 23.75874710068368, 0.2313686559958002 refl:, 0.8838834764831843, 0.4505747126436782, 23.46916928052956, 0.23090689658521463 refl:, 0.8838834764831843, 0.4454545454545454, 23.18678545794031, 0.23039248392782377 refl:, 0.8838834764831843, 0.44044943820224725, 22.911322384688706, 0.22982847953509086 refl:, 0.8838834764831843, 0.43555555555555553, 22.64252090923418, 0.22921468760037644 refl:, 0.8838834764831843, 0.4307692307692308, 22.380135051959574, 0.2285503145110825 refl:, 0.8838834764831843, 0.4260869565217391, 22.123931154313652, 0.22783728199610298 refl:, 0.8838834764831843, 0.421505376344086, 21.873687094881706, 0.22707701516938034 refl:, 0.8838834764831843, 0.41702127659574467, 21.629191566166806, 0.226267933168201 refl:, 0.8838834764831843, 0.4126315789473684, 21.390243406530644, 0.22540929123705422 refl:, 0.8838834764831843, 0.4083333333333333, 21.156650982328358, 0.22450324859888127 refl:, 0.8838834764831843, 0.4041237113402062, 20.928231615787418, 0.2235503997974889 refl:, 0.8838834764831843, 0.4, 20.704811054635428, 0.22254911459777174 ----------- Initializing structure... time for choose_chunkdivision = 0.000109911 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.006423 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.18417632037527257 / 0.18417632037527257 = 1.0 field decay(t = 100.01): 1.0178427595748204e-12 / 0.18417632037527257 = 5.5264583280895835e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000110865 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.0184629 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.18417631975521828 / 0.18417631975521828 = 1.0 field decay(t = 100.01): 2.8713012611847576e-11 / 0.18417631975521828 = 1.5589958931750262e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 0.9575555538987225, 0.8, 50.0, 0.14429067301668 refl:, 0.9575555538987225, 0.784, 48.6530933185261, 0.15241249489296654 refl:, 0.9575555538987225, 0.7686274509803922, 47.39207738436693, 0.15970237410745958 refl:, 0.9575555538987225, 0.7538461538461539, 46.207396358038395, 0.16628917742963503 refl:, 0.9575555538987225, 0.739622641509434, 45.091066324150106, 0.1722363131532337 refl:, 0.9575555538987225, 0.7259259259259259, 44.036334278275795, 0.1776288135399773 refl:, 0.9575555538987225, 0.7127272727272727, 43.037427456538545, 0.1825236470768417 refl:, 0.9575555538987225, 0.7, 42.089365210987516, 0.1869718164862598 refl:, 0.9575555538987225, 0.6877192982456141, 41.18781524255363, 0.19102552097133876 refl:, 0.9575555538987225, 0.6758620689655173, 40.32898196453324, 0.19471651743120957 refl:, 0.9575555538987225, 0.664406779661017, 39.50951857909891, 0.1980846376246996 refl:, 0.9575555538987225, 0.6533333333333333, 38.726456948394606, 0.20115681830261414 refl:, 0.9575555538987225, 0.6426229508196721, 37.97715101972185, 0.20395988101926588 refl:, 0.9575555538987225, 0.632258064516129, 37.25923071465613, 0.20651807398114963 refl:, 0.9575555538987225, 0.6222222222222222, 36.57056399547784, 0.20884957043569513 refl:, 0.9575555538987225, 0.6124999999999999, 35.909225393217994, 0.21097550520798009 refl:, 0.9575555538987225, 0.6030769230769231, 35.273469693564564, 0.21290924565401623 refl:, 0.9575555538987225, 0.593939393939394, 34.6617097783444, 0.214665416754943 refl:, 0.9575555538987225, 0.5850746268656717, 34.07249784378941, 0.2162593463941121 refl:, 0.9575555538987225, 0.5764705882352942, 33.504509384470694, 0.21769885837479563 refl:, 0.9575555538987225, 0.5681159420289855, 32.95652945897567, 0.2189956192414337 refl:, 0.9575555538987225, 0.56, 32.42744085087325, 0.22016058268759847 refl:, 0.9575555538987225, 0.552112676056338, 31.91621381391884, 0.22119898455404 refl:, 0.9575555538987225, 0.5444444444444444, 31.42189714930561, 0.22211935917204448 refl:, 0.9575555538987225, 0.536986301369863, 30.943610409083956, 0.2229298821236265 refl:, 0.9575555538987225, 0.5297297297297298, 30.480537056602333, 0.22363493657078803 refl:, 0.9575555538987225, 0.5226666666666667, 30.031918444163345, 0.2242400133418431 refl:, 0.9575555538987225, 0.5157894736842105, 29.597048491686717, 0.22475115535497284 refl:, 0.9575555538987225, 0.509090909090909, 29.17526896927061, 0.2251729335780506 refl:, 0.9575555538987225, 0.5025641025641026, 28.76596530209651, 0.22550858491249842 refl:, 0.9575555538987225, 0.4962025316455696, 28.36856282886028, 0.22576196217784258 refl:, 0.9575555538987225, 0.49, 27.982523455399697, 0.22593807040959396 refl:, 0.9575555538987225, 0.4839506172839506, 27.60734265386835, 0.22603917633476883 refl:, 0.9575555538987225, 0.47804878048780486, 27.24254676502359, 0.22606684493527063 refl:, 0.9575555538987225, 0.47228915662650606, 26.88769056722596, 0.22602594884033364 refl:, 0.9575555538987225, 0.4666666666666667, 26.54235508080712, 0.22591913511748662 refl:, 0.9575555538987225, 0.4611764705882353, 26.206145580726076, 0.2257461894573964 refl:, 0.9575555538987225, 0.4558139534883721, 25.878689794039474, 0.2255106337747593 refl:, 0.9575555538987225, 0.4505747126436782, 25.55963626177365, 0.22521611643853037 refl:, 0.9575555538987225, 0.4454545454545454, 25.248652847395146, 0.2248621048576318 refl:, 0.9575555538987225, 0.44044943820224725, 24.945425376307522, 0.2244497571914091 refl:, 0.9575555538987225, 0.43555555555555553, 24.64965639271632, 0.22398274531115286 refl:, 0.9575555538987225, 0.4307692307692308, 24.361064021851533, 0.22346171326862457 refl:, 0.9575555538987225, 0.4260869565217391, 24.079380926958596, 0.22288585592315466 refl:, 0.9575555538987225, 0.421505376344086, 23.804353351700417, 0.22225690920257898 refl:, 0.9575555538987225, 0.41702127659574467, 23.535740239681235, 0.22157652716566317 refl:, 0.9575555538987225, 0.4126315789473684, 23.27331242373393, 0.22084351574148542 refl:, 0.9575555538987225, 0.4083333333333333, 23.016851878423925, 0.220057130518136 refl:, 0.9575555538987225, 0.4041237113402062, 22.76615102993319, 0.21921919579654955 refl:, 0.9575555538987225, 0.4, 22.521012118111, 0.21832995493951468 ----------- Initializing structure... time for choose_chunkdivision = 0.000108004 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00824594 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.17431196343224772 / 0.17431196343224772 = 1.0 field decay(t = 100.01): 1.226599619152291e-12 / 0.17431196343224772 = 7.0368068547919865e-12 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000108957 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.0162461 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.17431196264184795 / 0.17431196264184795 = 1.0 field decay(t = 100.01): 2.9177363356244976e-11 / 0.17431196264184795 = 1.6738589201817764e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.0239400553612397, 0.8, 54.99999999999999, 0.11015750327021501 refl:, 1.0239400553612397, 0.784, 53.39534223243568, 0.12100519786278613 refl:, 1.0239400553612397, 0.7686274509803922, 51.90867481651075, 0.13067356321855755 refl:, 1.0239400553612397, 0.7538461538461539, 50.524208813748004, 0.1393743175816963 refl:, 1.0239400553612397, 0.739622641509434, 49.22931269112476, 0.14718461760941978 refl:, 1.0239400553612397, 0.7259259259259259, 48.013685667151286, 0.15425442212255983 refl:, 1.0239400553612397, 0.7127272727272727, 46.86879216202182, 0.16064116852684127 refl:, 1.0239400553612397, 0.7, 45.787462696319906, 0.1664466525245801 refl:, 1.0239400553612397, 0.6877192982456141, 44.7636047370949, 0.17171830451504444 refl:, 1.0239400553612397, 0.6758620689655173, 43.79198840588911, 0.17651889499651638 refl:, 1.0239400553612397, 0.664406779661017, 42.868084512177944, 0.1808955749197593 refl:, 1.0239400553612397, 0.6533333333333333, 41.98794000564178, 0.18488590864155194 refl:, 1.0239400553612397, 0.6426229508196721, 41.148080731197126, 0.1885334768453189 refl:, 1.0239400553612397, 0.632258064516129, 40.345434464399176, 0.19186399028554912 refl:, 1.0239400553612397, 0.6222222222222222, 39.57726925315509, 0.1949105801621511 refl:, 1.0239400553612397, 0.6124999999999999, 38.841143478370526, 0.19769436000574783 refl:, 1.0239400553612397, 0.6030769230769231, 38.13486500383326, 0.20023789123980396 refl:, 1.0239400553612397, 0.593939393939394, 37.45645745905701, 0.2025624060128163 refl:, 1.0239400553612397, 0.5850746268656717, 36.80413218014341, 0.204680856177325 refl:, 1.0239400553612397, 0.5764705882352942, 36.176264682924305, 0.2066128346051995 refl:, 1.0239400553612397, 0.5681159420289855, 35.57137479946729, 0.20837010840175924 refl:, 1.0239400553612397, 0.56, 34.98810980028866, 0.2099625172331732 refl:, 1.0239400553612397, 0.552112676056338, 34.42522996870841, 0.21140487263605512 refl:, 1.0239400553612397, 0.5444444444444444, 33.881596203503065, 0.2127045522560479 refl:, 1.0239400553612397, 0.536986301369863, 33.35615931039313, 0.21386953098052403 refl:, 1.0239400553612397, 0.5297297297297298, 32.847950708397136, 0.21491013773280657 refl:, 1.0239400553612397, 0.5226666666666667, 32.3560743283589, 0.21583212552616557 refl:, 1.0239400553612397, 0.5157894736842105, 31.87969952142079, 0.2166414381263308 refl:, 1.0239400553612397, 0.509090909090909, 31.41805482739263, 0.21734485598596862 refl:, 1.0239400553612397, 0.5025641025641026, 30.97042247873359, 0.2179480717640861 refl:, 1.0239400553612397, 0.4962025316455696, 30.536133536635912, 0.21845549976557374 refl:, 1.0239400553612397, 0.49, 30.114563572550164, 0.21887075023784372 refl:, 1.0239400553612397, 0.4839506172839506, 29.70512882224086, 0.21919941164026552 refl:, 1.0239400553612397, 0.47804878048780486, 29.307282750744406, 0.21944546389616373 refl:, 1.0239400553612397, 0.47228915662650606, 28.920512975908146, 0.21961022544356806 refl:, 1.0239400553612397, 0.4666666666666667, 28.544338505905493, 0.2196984674330614 refl:, 1.0239400553612397, 0.4611764705882353, 28.178307252549438, 0.21971428386147993 refl:, 1.0239400553612397, 0.4558139534883721, 27.821993787605123, 0.21965764862748693 refl:, 1.0239400553612397, 0.4505747126436782, 27.474997313821564, 0.21953193965864698 refl:, 1.0239400553612397, 0.4454545454545454, 27.13693982621647, 0.21934178652625613 refl:, 1.0239400553612397, 0.44044943820224725, 26.80746444237855, 0.21908700093353123 refl:, 1.0239400553612397, 0.43555555555555553, 26.486233883298294, 0.21876853077786598 refl:, 1.0239400553612397, 0.4307692307692308, 26.172929088582443, 0.2183905958409472 refl:, 1.0239400553612397, 0.4260869565217391, 25.867247951913267, 0.21795488284012238 refl:, 1.0239400553612397, 0.421505376344086, 25.56890416433822, 0.21746063743615796 refl:, 1.0239400553612397, 0.41702127659574467, 25.277626154460066, 0.21690850348145124 refl:, 1.0239400553612397, 0.4126315789473684, 24.993156115881433, 0.21630059246127992 refl:, 1.0239400553612397, 0.4083333333333333, 24.715249113370163, 0.21563737630418256 refl:, 1.0239400553612397, 0.4041237113402062, 24.443672260178424, 0.2149184652144797 refl:, 1.0239400553612397, 0.4, 24.178203959791162, 0.2141434469115452 ----------- Initializing structure... time for choose_chunkdivision = 0.000112057 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.0063448 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.16505586379118362 / 0.16505586379118362 = 1.0 field decay(t = 100.01): 1.6796657955389861e-12 / 0.16505586379118362 = 1.0176347310289891e-11 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000111818 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.020278 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.16505586311279807 / 0.16505586311279807 = 1.0 field decay(t = 100.01): 3.3223715038456985e-11 / 0.16505586311279807 = 2.0128769988468766e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.0825317547305482, 0.8, 59.99999999999999, 0.07360773441443474 refl:, 1.0825317547305482, 0.784, 58.07108488083593, 0.08756183349395384 refl:, 1.0825317547305482, 0.7686274509803922, 56.311309260669624, 0.09993186958234106 refl:, 1.0825317547305482, 0.7538461538461539, 54.69254540505164, 0.11102792920544496 refl:, 1.0825317547305482, 0.739622641509434, 53.19365179333026, 0.12095671757064261 refl:, 1.0825317547305482, 0.7259259259259259, 51.798237196004614, 0.12989111490169047 refl:, 1.0825317547305482, 0.7127272727272727, 50.49327371888257, 0.13795385843502186 refl:, 1.0825317547305482, 0.7, 49.268194122622916, 0.14525656081749277 refl:, 1.0825317547305482, 0.6877192982456141, 48.11428082873218, 0.15187584382721003 refl:, 1.0825317547305482, 0.6758620689655173, 47.024238635081495, 0.15790192109849083 refl:, 1.0825317547305482, 0.664406779661017, 45.99188753936469, 0.16337692994168027 refl:, 1.0825317547305482, 0.6533333333333333, 45.01193662107726, 0.16837576220748254 refl:, 1.0825317547305482, 0.6426229508196721, 44.07981414374115, 0.1729308563822558 refl:, 1.0825317547305482, 0.632258064516129, 43.191537588393835, 0.17709787538099556 refl:, 1.0825317547305482, 0.6222222222222222, 42.34361264736431, 0.18090702518494378 refl:, 1.0825317547305482, 0.6124999999999999, 41.53295361356457, 0.18439374819992094 refl:, 1.0825317547305482, 0.6030769230769231, 40.756819839298586, 0.1875864032458676 refl:, 1.0825317547305482, 0.593939393939394, 40.01276444425686, 0.1905057507798468 refl:, 1.0825317547305482, 0.5850746268656717, 39.298592486049785, 0.19318030905511865 refl:, 1.0825317547305482, 0.5764705882352942, 38.61232652957572, 0.1956239605857035 refl:, 1.0825317547305482, 0.5681159420289855, 37.9521780657052, 0.19785784772625742 refl:, 1.0825317547305482, 0.56, 37.31652360111115, 0.19989839276598212 refl:, 1.0825317547305482, 0.552112676056338, 36.70388451304266, 0.2017547526768604 refl:, 1.0825317547305482, 0.5444444444444444, 36.112909964600526, 0.20344360330889855 refl:, 1.0825317547305482, 0.536986301369863, 35.54236232752769, 0.20497480625666326 refl:, 1.0825317547305482, 0.5297297297297298, 34.991104674473426, 0.20635699484209216 refl:, 1.0825317547305482, 0.5226666666666667, 34.4580899908132, 0.20760151265764157 refl:, 1.0825317547305482, 0.5157894736842105, 33.94235182430718, 0.20871528450622134 refl:, 1.0825317547305482, 0.509090909090909, 33.442996144127704, 0.20970524738496094 refl:, 1.0825317547305482, 0.5025641025641026, 32.95919422270091, 0.2105783514069682 refl:, 1.0825317547305482, 0.4962025316455696, 32.4901763870572, 0.2113414215404359 refl:, 1.0825317547305482, 0.49, 32.035226512949734, 0.21200027362523924 refl:, 1.0825317547305482, 0.4839506172839506, 31.593677156370052, 0.21255807961996584 refl:, 1.0825317547305482, 0.47804878048780486, 31.164905234389202, 0.21302071809138926 refl:, 1.0825317547305482, 0.47228915662650606, 30.748328181342448, 0.21339356617448996 refl:, 1.0825317547305482, 0.4666666666666667, 30.343400517915825, 0.2136782430141135 refl:, 1.0825317547305482, 0.4611764705882353, 29.949610780196405, 0.21387969509149066 refl:, 1.0825317547305482, 0.4558139534883721, 29.566478763613983, 0.2140025648813409 refl:, 1.0825317547305482, 0.4505747126436782, 29.193553043244272, 0.21404664550288702 refl:, 1.0825317547305482, 0.4454545454545454, 28.830408737409662, 0.2140158741099825 refl:, 1.0825317547305482, 0.44044943820224725, 28.47664548610104, 0.2139163515040305 refl:, 1.0825317547305482, 0.43555555555555553, 28.131885619609587, 0.213747755347228 refl:, 1.0825317547305482, 0.4307692307692308, 27.79577249602797, 0.2135103221579689 refl:, 1.0825317547305482, 0.4260869565217391, 27.46796898905743, 0.21320790783224955 refl:, 1.0825317547305482, 0.421505376344086, 27.14815610992485, 0.2128440342182022 refl:, 1.0825317547305482, 0.41702127659574467, 26.836031749237822, 0.21242171734687695 refl:, 1.0825317547305482, 0.4126315789473684, 26.531309526343414, 0.21193780892543815 refl:, 1.0825317547305482, 0.4083333333333333, 26.23371773525153, 0.21138950860192895 refl:, 1.0825317547305482, 0.4041237113402062, 25.94299837747481, 0.21078013578978685 refl:, 1.0825317547305482, 0.4, 25.65890627325528, 0.2101171783272214 ----------- Initializing structure... time for choose_chunkdivision = 0.000109911 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00752401 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.15668265201002152 / 0.15668265201002152 = 1.0 field decay(t = 100.01): 2.2446823295144107e-12 / 0.15668265201002152 = 1.4326297779098349e-11 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000112057 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.0143809 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.1566826516265948 / 0.1566826516265948 = 1.0 field decay(t = 100.01): 3.452976998428892e-11 / 0.1566826516265948 = 2.2038030136597423e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.1328847337958123, 0.8, 64.99999999999999, 0.037170886043112084 refl:, 1.1328847337958123, 0.784, 62.645633431180784, 0.05362783985573467 refl:, 1.1328847337958123, 0.7686274509803922, 60.547811567584965, 0.06875709351287543 refl:, 1.1328847337958123, 0.7538461538461539, 58.65172469781183, 0.082354616478198 refl:, 1.1328847337958123, 0.739622641509434, 56.919780417243786, 0.09462755396270182 refl:, 1.1328847337958123, 0.7259259259259259, 55.32480276439999, 0.10557645638030469 refl:, 1.1328847337958123, 0.7127272727272727, 53.84635215217889, 0.11544950273548941 refl:, 1.1328847337958123, 0.7, 52.468573108244414, 0.1243168343725702 refl:, 1.1328847337958123, 0.6877192982456141, 51.17885777022349, 0.13235703881433808 refl:, 1.1328847337958123, 0.6758620689655173, 49.96697533403755, 0.139626360552005 refl:, 1.1328847337958123, 0.664406779661017, 48.82448244349733, 0.14625471025031767 refl:, 1.1328847337958123, 0.6533333333333333, 47.744310665402864, 0.1522788790416092 refl:, 1.1328847337958123, 0.6426229508196721, 46.720469815981424, 0.1577852862490647 refl:, 1.1328847337958123, 0.632258064516129, 45.7478295128327, 0.16280555175366354 refl:, 1.1328847337958123, 0.6222222222222222, 44.82195500406973, 0.167398621072885 refl:, 1.1328847337958123, 0.6124999999999999, 43.93898156096055, 0.17159777718924688 refl:, 1.1328847337958123, 0.6030769230769231, 43.09551684594588, 0.17544391308281032 refl:, 1.1328847337958123, 0.593939393939394, 42.288563952640665, 0.1789701508099767 refl:, 1.1328847337958123, 0.5850746268656717, 41.51545997434592, 0.18219866439447 refl:, 1.1328847337958123, 0.5764705882352942, 40.77382641078957, 0.185160771243366 refl:, 1.1328847337958123, 0.5681159420289855, 40.06152872079658, 0.18786997681737858 refl:, 1.1328847337958123, 0.56, 39.37664302675111, 0.19035022781311492 refl:, 1.1328847337958123, 0.552112676056338, 38.717428473384835, 0.19261993394590182 refl:, 1.1328847337958123, 0.5444444444444444, 38.08230410219336, 0.19469110743273724 refl:, 1.1328847337958123, 0.536986301369863, 37.4698293655767, 0.1965812755931606 refl:, 1.1328847337958123, 0.5297297297297298, 36.87868759977799, 0.19830004161335185 refl:, 1.1328847337958123, 0.5226666666666667, 36.307671922074825, 0.19985824848269654 refl:, 1.1328847337958123, 0.5157894736842105, 35.75567312877807, 0.20126711563403166 refl:, 1.1328847337958123, 0.509090909090909, 35.22166925577522, 0.20253568251243403 refl:, 1.1328847337958123, 0.5025641025641026, 34.70471652928195, 0.2036726801968963 refl:, 1.1328847337958123, 0.4962025316455696, 34.203941485939055, 0.20468320295652365 refl:, 1.1328847337958123, 0.49, 33.7185340819115, 0.20557517176695547 refl:, 1.1328847337958123, 0.4839506172839506, 33.247741642788455, 0.20635542059630924 refl:, 1.1328847337958123, 0.47804878048780486, 32.790863531764046, 0.2070277254966668 refl:, 1.1328847337958123, 0.47228915662650606, 32.347246434236666, 0.20759964094140226 refl:, 1.1328847337958123, 0.4666666666666667, 31.91628017368787, 0.20807496149136276 refl:, 1.1328847337958123, 0.4611764705882353, 31.497393987322308, 0.20845469225613245 refl:, 1.1328847337958123, 0.4558139534883721, 31.090053201105746, 0.20874562394291427 refl:, 1.1328847337958123, 0.4505747126436782, 30.693756253024553, 0.20895386380475825 refl:, 1.1328847337958123, 0.4454545454545454, 30.308032020994055, 0.20908023176898172 refl:, 1.1328847337958123, 0.44044943820224725, 29.9324374181665, 0.209125738596599 refl:, 1.1328847337958123, 0.43555555555555553, 29.5665552236732, 0.20909456295065154 refl:, 1.1328847337958123, 0.4307692307692308, 29.209992121269647, 0.20899049887482632 refl:, 1.1328847337958123, 0.4260869565217391, 28.86237692208816, 0.2088190433084917 refl:, 1.1328847337958123, 0.421505376344086, 28.523358950864424, 0.2085833129971406 refl:, 1.1328847337958123, 0.41702127659574467, 28.192606577687982, 0.2082742041204916 refl:, 1.1328847337958123, 0.4126315789473684, 27.86980587961545, 0.2078930131964556 refl:, 1.1328847337958123, 0.4083333333333333, 27.554659418441407, 0.20745662776433935 refl:, 1.1328847337958123, 0.4041237113402062, 27.246885122601544, 0.20697561960434535 refl:, 1.1328847337958123, 0.4, 26.946215262627685, 0.20643519600186538 ----------- Initializing structure... time for choose_chunkdivision = 0.000110865 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00959396 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.14944315649486126 / 0.14944315649486126 = 1.0 field decay(t = 100.01): 4.987937476864823e-12 / 0.14944315649486126 = 3.3376820952228335e-11 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000108957 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.0218339 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.14944315644072245 / 0.14944315644072245 = 1.0 field decay(t = 100.01): 3.807851695141726e-11 / 0.14944315644072245 = 2.548026812222836e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.1746157759823854, 0.8, 70.0, 0.008715934286147569 refl:, 1.1746157759823854, 0.784, 67.05783140972835, 0.024210799209089452 refl:, 1.1746157759823854, 0.7686274509803922, 64.53417775149677, 0.04029873111181387 refl:, 1.1746157759823854, 0.7538461538461539, 62.31059707934406, 0.0555678814867085 refl:, 1.1746157759823854, 0.739622641509434, 60.31629899672389, 0.06970125775811391 refl:, 1.1746157759823854, 0.7259259259259259, 58.50481213939229, 0.08255524971343331 refl:, 1.1746157759823854, 0.7127272727272727, 56.843594941553754, 0.09413263635827869 refl:, 1.1746157759823854, 0.7, 55.30875739675561, 0.10463698983404623 refl:, 1.1746157759823854, 0.6877192982456141, 53.88211612099151, 0.11407493400247079 refl:, 1.1746157759823854, 0.6758620689655173, 52.54943143507829, 0.1226472707657511 refl:, 1.1746157759823854, 0.664406779661017, 51.29929235382265, 0.13039930819744874 refl:, 1.1746157759823854, 0.6533333333333333, 50.122379445158494, 0.13744242808130103 refl:, 1.1746157759823854, 0.6426229508196721, 49.010959419193355, 0.1438468426182844 refl:, 1.1746157759823854, 0.632258064516129, 47.958527855818005, 0.14968921112177416 refl:, 1.1746157759823854, 0.6222222222222222, 46.9595500055931, 0.1550186900782557 refl:, 1.1746157759823854, 0.6124999999999999, 46.00926848902442, 0.15990623687364083 refl:, 1.1746157759823854, 0.6030769230769231, 45.103557823882895, 0.1643748532843572 refl:, 1.1746157759823854, 0.593939393939394, 44.23881248012816, 0.16847752780425043 refl:, 1.1746157759823854, 0.5850746268656717, 43.41185942214757, 0.1722356163059057 refl:, 1.1746157759823854, 0.5764705882352942, 42.619888854495166, 0.17567897390885354 refl:, 1.1746157759823854, 0.5681159420289855, 41.860398715356425, 0.17883996497200733 refl:, 1.1746157759823854, 0.56, 41.131149701242414, 0.18173480153232927 refl:, 1.1746157759823854, 0.552112676056338, 40.430128463310105, 0.18439121641848574 refl:, 1.1746157759823854, 0.5444444444444444, 39.75551721886283, 0.18682378756973045 refl:, 1.1746157759823854, 0.536986301369863, 39.10566845306894, 0.18904674104030247 refl:, 1.1746157759823854, 0.5297297297297298, 38.479083699202135, 0.19107686788902303 refl:, 1.1746157759823854, 0.5226666666666667, 37.87439561623927, 0.19292625454975332 refl:, 1.1746157759823854, 0.5157894736842105, 37.29035275442276, 0.19460884663783543 refl:, 1.1746157759823854, 0.509090909090909, 36.72580652885898, 0.1961345571897698 refl:, 1.1746157759823854, 0.5025641025641026, 36.179700019842116, 0.1975122933988443 refl:, 1.1746157759823854, 0.4962025316455696, 35.651058294458764, 0.19875051751344916 refl:, 1.1746157759823854, 0.49, 35.138980002927255, 0.19985541870413245 refl:, 1.1746157759823854, 0.4839506172839506, 34.64263004924646, 0.20083789883449643 refl:, 1.1746157759823854, 0.47804878048780486, 34.161233172132036, 0.2017049645889793 refl:, 1.1746157759823854, 0.47228915662650606, 33.69406830116886, 0.20245751259203246 refl:, 1.1746157759823854, 0.4666666666666667, 33.24046357629389, 0.20310354230391978 refl:, 1.1746157759823854, 0.4611764705882353, 32.79979193741594, 0.2036476720462249 refl:, 1.1746157759823854, 0.4558139534883721, 32.37146720614281, 0.20409307625590972 refl:, 1.1746157759823854, 0.4505747126436782, 31.954940593960732, 0.20445096223605416 refl:, 1.1746157759823854, 0.4454545454545454, 31.549697581364857, 0.20471845410715214 refl:, 1.1746157759823854, 0.44044943820224725, 31.155255120816925, 0.20489444063169535 refl:, 1.1746157759823854, 0.43555555555555553, 30.771159123350074, 0.2049889592193817 refl:, 1.1746157759823854, 0.4307692307692308, 30.396982194426784, 0.20500542736591415 refl:, 1.1746157759823854, 0.4260869565217391, 30.032321589495837, 0.20495768982428397 refl:, 1.1746157759823854, 0.421505376344086, 29.67679736376329, 0.20483056853588974 refl:, 1.1746157759823854, 0.41702127659574467, 29.33005069412437, 0.20461216911883964 refl:, 1.1746157759823854, 0.4126315789473684, 28.991742354111988, 0.2043374679527374 refl:, 1.1746157759823854, 0.4083333333333333, 28.66155132518973, 0.20400220026620766 refl:, 1.1746157759823854, 0.4041237113402062, 28.339173529827534, 0.2036399273981529 refl:, 1.1746157759823854, 0.4, 28.024320673604695, 0.20318419566870832 ----------- Initializing structure... time for choose_chunkdivision = 0.000108957 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00775695 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.14355855111904026 / 0.14355855111904026 = 1.0 field decay(t = 100.01): 9.667874474337293e-12 / 0.14355855111904026 = 6.734446954901761e-11 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000108957 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.017272 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.1435585513271557 / 0.1435585513271557 = 1.0 field decay(t = 100.01): 5.780725980413138e-11 / 0.1435585513271557 = 4.0267374719040156e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.2074072828613354, 0.8, 75.0, 0.002279010064127535 refl:, 1.2074072828613354, 0.784, 71.19254305121969, 0.004052585814133056 refl:, 1.2074072828613354, 0.7686274509803922, 68.13230108184366, 0.017669454881364767 refl:, 1.2074072828613354, 0.7538461538461539, 65.5329128763873, 0.03336169435169541 refl:, 1.2074072828613354, 0.739622641509434, 63.25596068464814, 0.048654509225169895 refl:, 1.2074072828613354, 0.7259259259259259, 61.22160180379117, 0.06275216167717797 refl:, 1.2074072828613354, 0.7127272727272727, 59.378629248099124, 0.07576309243511313 refl:, 1.2074072828613354, 0.7, 57.691773022139905, 0.08758135787027703 refl:, 1.2074072828613354, 0.6877192982456141, 56.13545775731428, 0.098248929360899 refl:, 1.2074072828613354, 0.6758620689655173, 54.69040263962552, 0.10794664183105998 refl:, 1.2074072828613354, 0.664406779661017, 53.3416232015308, 0.11670121483117171 refl:, 1.2074072828613354, 0.6533333333333333, 52.0771858356392, 0.12465585169755833 refl:, 1.2074072828613354, 0.6426229508196721, 50.88739414232009, 0.13187389976516842 refl:, 1.2074072828613354, 0.632258064516129, 49.76423655717566, 0.13846172295217293 refl:, 1.2074072828613354, 0.6222222222222222, 48.70099913863975, 0.14446173011718028 refl:, 1.2074072828613354, 0.6124999999999999, 47.69198665804526, 0.14995965765204525 refl:, 1.2074072828613354, 0.6030769230769231, 46.73231696062462, 0.1549825384470223 refl:, 1.2074072828613354, 0.593939393939394, 45.817766249324485, 0.15958402411708034 refl:, 1.2074072828613354, 0.5850746268656717, 44.94465059915575, 0.16380217089082738 refl:, 1.2074072828613354, 0.5764705882352942, 44.109733785683, 0.1676654538139623 refl:, 1.2074072828613354, 0.5681159420289855, 43.31015457776125, 0.17121662693717593 refl:, 1.2074072828613354, 0.56, 42.543368664412704, 0.1744741445938872 refl:, 1.2074072828613354, 0.552112676056338, 41.80710174663229, 0.17746630440999525 refl:, 1.2074072828613354, 0.5444444444444444, 41.099311260706344, 0.1802095438402015 refl:, 1.2074072828613354, 0.536986301369863, 40.418154855029364, 0.18272173169827843 refl:, 1.2074072828613354, 0.5297297297297298, 39.76196420912173, 0.18502119081455226 refl:, 1.2074072828613354, 0.5226666666666667, 39.12922312098109, 0.1871223337651491 refl:, 1.2074072828613354, 0.5157894736842105, 38.51854903626071, 0.18904214560856375 refl:, 1.2074072828613354, 0.509090909090909, 37.928677376424204, 0.190789263829168 refl:, 1.2074072828613354, 0.5025641025641026, 37.35844816099257, 0.192373003931677 refl:, 1.2074072828613354, 0.4962025316455696, 36.806794523773725, 0.19380521484150348 refl:, 1.2074072828613354, 0.49, 36.272732803336254, 0.19509368153538312 refl:, 1.2074072828613354, 0.4839506172839506, 35.755353950220794, 0.19624859910356954 refl:, 1.2074072828613354, 0.47804878048780486, 35.253816041990945, 0.19727952544337504 refl:, 1.2074072828613354, 0.47228915662650606, 34.76733773550362, 0.19818701041137332 refl:, 1.2074072828613354, 0.4666666666666667, 34.295192516152845, 0.19897646054575882 refl:, 1.2074072828613354, 0.4611764705882353, 33.83670362811667, 0.19965889435875267 refl:, 1.2074072828613354, 0.4558139534883721, 33.391239589169516, 0.20023796589376308 refl:, 1.2074072828613354, 0.4505747126436782, 32.95821020943854, 0.20072106469864034 refl:, 1.2074072828613354, 0.4454545454545454, 32.537063046367535, 0.2011073384968122 refl:, 1.2074072828613354, 0.44044943820224725, 32.12728023870626, 0.2013932866552523 refl:, 1.2074072828613354, 0.43555555555555553, 31.728375671037586, 0.20159478768644698 refl:, 1.2074072828613354, 0.4307692307692308, 31.33989242755132, 0.20171671281423958 refl:, 1.2074072828613354, 0.4260869565217391, 30.96140049976025, 0.20176796051792995 refl:, 1.2074072828613354, 0.421505376344086, 30.592494717856905, 0.20173546160628877 refl:, 1.2074072828613354, 0.41702127659574467, 30.23279287960783, 0.2015919333831556 refl:, 1.2074072828613354, 0.4126315789473684, 29.88193405422147, 0.20140043913698885 refl:, 1.2074072828613354, 0.4083333333333333, 29.539577041619506, 0.20116952947891878 refl:, 1.2074072828613354, 0.4041237113402062, 29.2053989700849, 0.2008832766544487 refl:, 1.2074072828613354, 0.4, 28.879094017427605, 0.200502480384686 ----------- Initializing structure... time for choose_chunkdivision = 0.000110865 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 time for set_epsilon = 0.00754905 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.1392134275610875 / 0.1392134275610875 = 1.0 field decay(t = 100.01): 1.3144200424788948e-11 / 0.1392134275610875 = 9.441761944278838e-11 run 0 finished at t = 100.01 (10001 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 0.000108957 s Working in 3D dimensions. Computational cell is 0.02 x 0.02 x 12 with resolution 50 block, center = (0,0,3) size (1e+20,1e+20,6) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (12.25,12.25,12.25) time for set_epsilon = 0.019413 s ----------- Meep: using complex fields. field decay(t = 50.01): 0.13921342793500532 / 0.13921342793500532 = 1.0 field decay(t = 100.01): 9.871965144118156e-11 / 0.13921342793500532 = 7.091244925544889e-10 run 0 finished at t = 100.01 (10001 timesteps) refl:, 1.23100969126526, 0.8, 79.99999999999994, 0.028667884842985224 refl:, 1.23100969126526, 0.784, 74.82079670263325, 0.0021347453734021486 refl:, 1.23100969126526, 0.7686274509803922, 71.11813585396203, 0.0041164349187047915 refl:, 1.23100969126526, 0.7538461538461539, 68.12392493650243, 0.017504855468406976 refl:, 1.23100969126526, 0.739622641509434, 65.57213417097068, 0.03282494395652445 refl:, 1.23100969126526, 0.7259259259259259, 63.331955818023665, 0.04792423816150501 refl:, 1.23100969126526, 0.7127272727272727, 61.32721657609733, 0.06193140028783135 refl:, 1.23100969126526, 0.7, 59.508764013939604, 0.07458578833476998 refl:, 1.23100969126526, 0.6877192982456141, 57.84259890532987, 0.08617897772571746 refl:, 1.23100969126526, 0.6758620689655173, 56.30398980710504, 0.0966374263474978 refl:, 1.23100969126526, 0.664406779661017, 54.87424701948212, 0.10612404031633497 refl:, 1.23100969126526, 0.6533333333333333, 53.53881754577338, 0.1147736073308378 refl:, 1.23100969126526, 0.6426229508196721, 52.2860932123638, 0.12261878935484406 refl:, 1.23100969126526, 0.632258064516129, 51.10662968356778, 0.1297955968814688 refl:, 1.23100969126526, 0.6222222222222222, 49.992614952243095, 0.1363423963503902 refl:, 1.23100969126526, 0.6124999999999999, 48.937495982303304, 0.14232561295216176 refl:, 1.23100969126526, 0.6030769230769231, 47.93570930843333, 0.1477973393651342 refl:, 1.23100969126526, 0.593939393939394, 46.98248211140486, 0.1527979150330299 refl:, 1.23100969126526, 0.5850746268656717, 46.07368235996428, 0.15737616810594365 refl:, 1.23100969126526, 0.5764705882352942, 45.205703916101704, 0.16157175886307418 refl:, 1.23100969126526, 0.5681159420289855, 44.37537706820699, 0.16542702710683707 refl:, 1.23100969126526, 0.56, 43.57989789517631, 0.16896483740231436 refl:, 1.23100969126526, 0.552112676056338, 42.81677180338089, 0.17222042834178233 refl:, 1.23100969126526, 0.5444444444444444, 42.083767886770914, 0.17520079945615433 refl:, 1.23100969126526, 0.536986301369863, 41.378881661291835, 0.1779352569426369 refl:, 1.23100969126526, 0.5297297297297298, 40.70030435652892, 0.18044129732794095 refl:, 1.23100969126526, 0.5226666666666667, 40.046397397863686, 0.182734107792854 refl:, 1.23100969126526, 0.5157894736842105, 39.41567103836258, 0.18483861166276766 refl:, 1.23100969126526, 0.509090909090909, 38.806766338771496, 0.18675601688212784 refl:, 1.23100969126526, 0.5025641025641026, 38.218439871701385, 0.1884996414001091 refl:, 1.23100969126526, 0.4962025316455696, 37.64955065969264, 0.19008080645459632 refl:, 1.23100969126526, 0.49, 37.099048958377544, 0.1915079168724047 refl:, 1.23100969126526, 0.4839506172839506, 36.56596657389412, 0.19279885897563126 refl:, 1.23100969126526, 0.47804878048780486, 36.049408464086156, 0.19395530230339666 refl:, 1.23100969126526, 0.47228915662650606, 35.54854542021205, 0.19498186511670068 refl:, 1.23100969126526, 0.4666666666666667, 35.062607663064846, 0.19588322966927102 refl:, 1.23100969126526, 0.4611764705882353, 34.59087921692149, 0.19666472127374038 refl:, 1.23100969126526, 0.4558139534883721, 34.13269294834013, 0.19734835806450404 refl:, 1.23100969126526, 0.4505747126436782, 33.6874261758209, 0.197930496667696 refl:, 1.23100969126526, 0.4454545454545454, 33.25449677173539, 0.19840503576945626 refl:, 1.23100969126526, 0.44044943820224725, 32.83335969046958, 0.19877805012361263 refl:, 1.23100969126526, 0.43555555555555553, 32.42350386700291, 0.19904737688450852 refl:, 1.23100969126526, 0.4307692307692308, 32.0244494386133, 0.19925404758581897 refl:, 1.23100969126526, 0.4260869565217391, 31.63574524940934, 0.1993954429097289 refl:, 1.23100969126526, 0.421505376344086, 31.2569666032255, 0.19942422718572894 refl:, 1.23100969126526, 0.41702127659574467, 30.887713235292672, 0.1993498431053417 refl:, 1.23100969126526, 0.4126315789473684, 30.527607477190394, 0.19917542826890433 refl:, 1.23100969126526, 0.4083333333333333, 30.176292593038497, 0.19903411145233238 refl:, 1.23100969126526, 0.4041237113402062, 29.83343126780691, 0.19890224385206828 refl:, 1.23100969126526, 0.4, 29.498704231103652, 0.1984945592268268
Two-dimensional plots of the angular reflectance spectrum based on the simulated data and the analytic Fresnel equations are generated using the script below. The plots are shown in the accompanying figure with four insets. The top left inset shows the simulated and analytic reflectance spectra at a wavelength of 0.6 μm. The top right inset shows the simulated reflectance spectrum as a function of the source wavelength λ and Bloch-periodic wavevector $k_x$: $R(\lambda, k_x)$. The lower left inset is a transformation of $R(\lambda, k_x)$ into $R(\lambda, \theta)$. Note how the range of angles depends on the wavelength. For a particular angle, the reflectance is a constant for all wavelengths due to the dispersionless dielectric. The lower right inset is the analytic reflectance spectrum computed using the Fresnel equations. There is agreement between the simulated and analytic results. The Brewster's angle, where the transmittance is 1 and the reflectance is 0, is tan-1(3.5/1)=74.1°. This is also verified by the simulated results.
In order to generate results for the missing portion of the reflectance spectrum (i.e., the white region), we will need to rerun the simulations for different wavelength spectra.
plt.figure(dpi=200)
plt.pcolormesh(
kxs, wvls, Rmeep, cmap="inferno", shading="gouraud", vmin=0, vmax=Rmeep.max()
)
plt.axis([kxs[0, 0], kxs[0, -1], wvl_min, wvl_max])
plt.yticks([t for t in np.linspace(0.4, 0.8, 5)])
plt.xlabel(r"Bloch-periodic wavevector (k$_x$/2π)")
plt.ylabel("wavelength (μm)")
plt.title("reflectance (meep)")
cbar = plt.colorbar()
cbar.set_ticks([t for t in np.linspace(0, 0.4, 5)])
cbar.set_ticklabels(["{:.1f}".format(t) for t in np.linspace(0, 0.4, 5)])
/home/gcrouch/anaconda3/envs/mp/lib/python3.11/site-packages/numpy/core/getlimits.py:500: UserWarning: The value of the smallest subnormal for <class 'numpy.float64'> type is zero. setattr(self, word, getattr(machar, word).flat[0]) /home/gcrouch/anaconda3/envs/mp/lib/python3.11/site-packages/numpy/core/getlimits.py:89: UserWarning: The value of the smallest subnormal for <class 'numpy.float64'> type is zero. return self._float_to_str(self.smallest_subnormal)
plt.figure(dpi=200)
plt.pcolormesh(
thetas, wvls, Rmeep, cmap="inferno", shading="gouraud", vmin=0, vmax=Rmeep.max()
)
plt.axis([thetas.min(), thetas.max(), wvl_min, wvl_max])
plt.xticks([t for t in range(0, 100, 20)])
plt.yticks([t for t in np.linspace(0.4, 0.8, 5)])
plt.xlabel("angle of incident planewave (degrees)")
plt.ylabel("wavelength (μm)")
plt.title("reflectance (meep)")
cbar = plt.colorbar()
cbar.set_ticks([t for t in np.linspace(0, 0.4, 5)])
cbar.set_ticklabels(["{:.1f}".format(t) for t in np.linspace(0, 0.4, 5)])
n1 = 1
n2 = 3.5
# compute angle of refracted planewave in medium n2
# for incident planewave in medium n1 at angle theta_in
theta_out = lambda theta_in: math.asin(n1 * math.sin(theta_in) / n2)
# compute Fresnel reflectance for P-polarization in medium n2
# for incident planewave in medium n1 at angle theta_in
Rfresnel = (
lambda theta_in: math.fabs(
(n1 * math.cos(theta_out(theta_in)) - n2 * math.cos(theta_in))
/ (n1 * math.cos(theta_out(theta_in)) + n2 * math.cos(theta_in))
)
** 2
)
Ranalytic = np.empty((nfreq, theta_in.size))
for m in range(wvl.size):
for n in range(theta_in.size):
Ranalytic[m, n] = Rfresnel(math.radians(thetas[m, n]))
plt.figure(dpi=200)
plt.pcolormesh(
thetas,
wvls,
Ranalytic,
cmap="inferno",
shading="gouraud",
vmin=0,
vmax=Ranalytic.max(),
)
plt.axis([thetas.min(), thetas.max(), wvl_min, wvl_max])
plt.xticks([t for t in range(0, 100, 20)])
plt.yticks([t for t in np.linspace(0.4, 0.8, 5)])
plt.xlabel("angle of incident planewave (degrees)")
plt.ylabel("wavelength (μm)")
plt.title("reflectance (analytic)")
cbar = plt.colorbar()
cbar.set_ticks([t for t in np.linspace(0, 0.4, 5)])
cbar.set_ticklabels(["{:.1f}".format(t) for t in np.linspace(0, 0.4, 5)])
# note: this plots at 0.8 μm, not 0.6 as in the top left image above, in order to correctly show the Brewster angle
wvl_idx = 0
plt.figure(dpi=200)
plt.plot(theta_in,Rmeep[wvl_idx,:],'-s',theta_in,Ranalytic[wvl_idx,:],'-s')
plt.xlabel("angle of incident planewave (degrees)")
plt.ylabel(f"reflectance at $\lambda$ = {wvl[wvl_idx]:.3f} μm")
plt.legend(["meep","analytic"])
plt.show()