Diffraction Spectrum of Liquid-Crystal Polarization Gratings

As a final demonstration of mode decomposition, we compute the diffraction spectrum of a liquid-crystal polarization grating. These types of beam splitters use birefringence to produce diffraction orders which are circularly polarized. We will investigate two kinds of polarization gratings: (1) a homogeneous uniaxial grating (commonly known as a circular-polarization grating), and (2) a twisted-nematic bilayer grating as described in Optics Letters, Vol. 33, No. 20, pp. 2287-9, 2008 (pdf). The homogeneous uniaxial grating is just a special case of the twisted-nematic grating with a nematic director rotation angle of φ=0°.

A schematic of the grating geometry is shown below. The grating is a 2d slab in the xy-plane with two parameters: birefringence (Δn) and thickness (d). The twisted-nematic grating consists of two layers of thickness d each with equal and opposite rotation angles of φ=70° for the nematic director. Both gratings contain only three diffraction orders: m=0, ±1. The m=0 order is linearly polarized and the m=±1 orders are circularly polarized with opposite chirality. For the uniaxial grating, the diffraction efficiencies for a mode with wavelength λ can be computed analytically: η0=cos2(πΔnd/λ), η±1=0.5sin2(πΔnd/λ). The derivation of these formulas is presented in Optics Letters, Vol. 24, No. 9, pp. 584-6, 1999. We will verify these analytic results and also demonstrate that the twisted-nematic grating produces a broader bandwidth response for the ±1 orders than the homogeneous uniaxial grating. An important property of these polarization gratings for e.g. display applications is that for a circular-polarized input planewave and phase delay (Δnd/λ) of nearly 0.5, there is only a single diffraction order (+1 or -1) with opposite chiraity to that of the input. This is also demonstrated below.

In this example, the input is a linear-polarized planewave pulse at normal incidence with center wavelength of λ=0.54 μm. The linear polarization is in the yz-plane with a rotation angle of 45° counter clockwise around the x axis. Two sets of mode coefficients are computed in the air region adjacent to the grating for each orthogonal polarization: ODD_Z+EVEN_Y and EVEN_Z+ODD_Y, which correspond to +ky + -ky (cosine) and +ky - -ky (sine) modes. From these coefficients for linear-polarized modes, the power in the circular-polarized modes can be computed: |ODD_Z+EVEN_Y|2+|EVEN_Z+ODD_Y|2. The power is identical for the two circular-polarized modes with opposite chiralities since the input is linearly polarized and at normal incidence. The transmittance for the diffraction orders are computed from the mode coefficients. As usual, this requires a separate normalization run to compute the power of the input planewave.

The main part of the script is the function pol_grating which computes the mode coefficients for a grating with thickness d, twisted-nematic rotation angle ph, and periodicity gp. The anisotropic permittivity of the grating is specified using the material function lc_mat which involves a position-dependent rotation of the diagonal ε tensor about the x-axis. For φ=0°, the nematic director is oriented along the z-axis: Ez has a larger permittivity than Ey where the birefringence (Δn) is 0.159. The grating has a periodicity of Λ=6.5 μm in the y direction.

In [1]:
import meep as mp
import math
import numpy as np
import matplotlib.pyplot as plt

resolution = 30        # pixels/μm

dpml = 1.0             # PML thickness
dsub = 1.0             # substrate thickness
dpad = 1.0             # padding thickness

k_point = mp.Vector3(0,0,0)

pml_layers = [mp.PML(thickness=dpml,direction=mp.X)]

n_0 = 1.55
delta_n = 0.159
epsilon_diag = mp.Matrix(mp.Vector3(n_0**2,0,0),mp.Vector3(0,n_0**2,0),mp.Vector3(0,0,(n_0+delta_n)**2))

wvl = 0.54             # center wavelength
fcen = 1/wvl           # center frequency

def pol_grating(d,ph,gp,nmode):
    sx = dpml+dsub+d+d+dpad+dpml
    sy = gp

    cell_size = mp.Vector3(sx,sy,0)

    # twist angle of nematic director; from equation 1b
    def phi(p):
        xx  = p.x-(-0.5*sx+dpml+dsub)
        if (xx >= 0) and (xx <= d):
            return math.pi*p.y/gp + ph*xx/d
        else:
            return math.pi*p.y/gp - ph*xx/d + 2*ph

    # return the anisotropic permittivity tensor for a uniaxial, twisted nematic liquid crystal
    def lc_mat(p):
        # rotation matrix for rotation around x axis
        Rx = mp.Matrix(mp.Vector3(1,0,0),mp.Vector3(0,math.cos(phi(p)),math.sin(phi(p))),mp.Vector3(0,-math.sin(phi(p)),math.cos(phi(p))))
        lc_epsilon = Rx * epsilon_diag * Rx.transpose()
        lc_epsilon_diag = mp.Vector3(lc_epsilon[0].x,lc_epsilon[1].y,lc_epsilon[2].z)
        lc_epsilon_offdiag = mp.Vector3(lc_epsilon[1].x,lc_epsilon[2].x,lc_epsilon[2].y)
        return mp.Medium(epsilon_diag=lc_epsilon_diag,epsilon_offdiag=lc_epsilon_offdiag)

    geometry = [mp.Block(center=mp.Vector3(-0.5*sx+0.5*(dpml+dsub)),size=mp.Vector3(dpml+dsub,mp.inf,mp.inf),material=mp.Medium(index=n_0)),
                mp.Block(center=mp.Vector3(-0.5*sx+dpml+dsub+d),size=mp.Vector3(2*d,mp.inf,mp.inf),material=lc_mat)]

    # linear-polarized planewave pulse source
    src_pt = mp.Vector3(-0.5*sx+dpml+0.3*dsub,0,0)
    sources = [mp.Source(mp.GaussianSource(fcen,fwidth=0.05*fcen), component=mp.Ez, center=src_pt, size=mp.Vector3(0,sy,0)),
               mp.Source(mp.GaussianSource(fcen,fwidth=0.05*fcen), component=mp.Ey, center=src_pt, size=mp.Vector3(0,sy,0))]

    sim = mp.Simulation(resolution=resolution,
                        cell_size=cell_size,
                        boundary_layers=pml_layers,
                        k_point=k_point,
                        sources=sources,
                        default_material=mp.Medium(index=n_0))

    tran_pt = mp.Vector3(0.5*sx-dpml-0.5*dpad,0,0)
    tran_flux = sim.add_flux(fcen, 0, 1, mp.FluxRegion(center=tran_pt, size=mp.Vector3(0,sy,0)))

    sim.run(until_after_sources=100)

    input_flux = mp.get_fluxes(tran_flux)
    input_flux_data = sim.get_flux_data(tran_flux)

    sim.reset_meep()

    sim = mp.Simulation(resolution=resolution,
                        cell_size=cell_size,
                        boundary_layers=pml_layers,
                        k_point=k_point,
                        sources=sources,
                        geometry=geometry)

    tran_flux = sim.add_flux(fcen, 0, 1, mp.FluxRegion(center=tran_pt, size=mp.Vector3(0,sy,0)))

    sim.run(until_after_sources=300)

    res1 = sim.get_eigenmode_coefficients(tran_flux, range(1,nmode+1), eig_parity=mp.ODD_Z+mp.EVEN_Y)
    res2 = sim.get_eigenmode_coefficients(tran_flux, range(1,nmode+1), eig_parity=mp.EVEN_Z+mp.ODD_Y)
    angles = [math.degrees(math.acos(kdom.x/fcen)) for kdom in res1.kdom]

    return input_flux[0], angles, res1.alpha[:,0,0], res2.alpha[:,0,0];

The properties of the two gratings are computed over the thickness range of 0.1 to 3.4 μm corresponding to phase delays (Δnd/λ) of approximately 0 to 1.

In [2]:
ph_uniaxial = 0               # chiral layer twist angle for uniaxial grating
ph_twisted = 70               # chiral layer twist angle for bilayer grating
gp = 6.5                      # grating period
nmode = 5                     # number of mode coefficients to compute
dd = np.arange(0.2,3.5,0.2)   # chiral layer thickness

m0_uniaxial = np.zeros(dd.size)
m1_uniaxial = np.zeros(dd.size)
ang_uniaxial = np.zeros(dd.size)

m0_twisted = np.zeros(dd.size)
m1_twisted = np.zeros(dd.size)
ang_twisted = np.zeros(dd.size)

for k in range(len(dd)):
    input_flux, angles, coeffs1, coeffs2 = pol_grating(0.5*dd[k],math.radians(ph_uniaxial),gp,nmode)
    tran = (abs(coeffs1)**2+abs(coeffs2)**2)/input_flux
    for m in range(nmode):
        print("tran (uniaxial):, {}, {:.2f}, {:.5f}".format(m,angles[m],tran[m]))
    m0_uniaxial[k] = tran[0]
    m1_uniaxial[k] = tran[1]
    ang_uniaxial[k] = angles[1]

    input_flux, angles, coeffs1, coeffs2 = pol_grating(dd[k],math.radians(ph_twisted),gp,nmode)
    tran = (abs(coeffs1)**2+abs(coeffs2)**2)/input_flux
    for m in range(nmode):
        print("tran (twisted):, {}, {:.2f}, {:.5f}".format(m,angles[m],tran[m]))
    m0_twisted[k] = tran[0]
    m1_twisted[k] = tran[1]
    ang_twisted[k] = angles[1]
-----------
Initializing structure...
time for choose_chunkdivision = 5.79357e-05 s
Working in 2D dimensions.
Computational cell is 4.2 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.0339479 s
-----------
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 2.71797e-05 s
Working in 2D dimensions.
Computational cell is 4.2 x 6.5 x 0 with resolution 30
     block, center = (-1.1,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (-8.32667e-17,0,0)
          size (0.2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 0.758735 s
-----------
Meep progress: 224.88333333333333/408.0 = 55.1% done in 4.0s, 3.3s to go
on time step 13493 (time=224.883), 0.000296471 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 65 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 37 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 47 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 42 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.89247
tran (uniaxial):, 1, 4.77, 0.01875
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 4.41074e-05 s
Working in 2D dimensions.
Computational cell is 4.4 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.0652409 s
-----------
Meep progress: 193.48333333333332/208.0 = 93.0% done in 4.0s, 0.3s to go
on time step 11609 (time=193.483), 0.000344592 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.00815e-05 s
Working in 2D dimensions.
Computational cell is 4.4 x 6.5 x 0 with resolution 30
     block, center = (-1.2,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (-1.66533e-16,0,0)
          size (0.4,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 3.35973 s
-----------
Meep progress: 100.36666666666666/408.0 = 24.6% done in 4.0s, 12.3s to go
on time step 6022 (time=100.367), 0.000664325 s/step
Meep progress: 200.35/408.0 = 49.1% done in 8.0s, 8.3s to go
on time step 12022 (time=200.367), 0.000666721 s/step
Meep progress: 301.4166666666667/408.0 = 73.9% done in 12.0s, 4.2s to go
on time step 18088 (time=301.467), 0.000659487 s/step
Meep progress: 400.9/408.0 = 98.3% done in 16.0s, 0.3s to go
on time step 24058 (time=400.967), 0.000670094 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 29 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 55 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.84732
tran (twisted):, 1, 4.77, 0.03707
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 8.70228e-05 s
Working in 2D dimensions.
Computational cell is 4.4 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.0973799 s
-----------
Meep progress: 104.83333333333333/208.0 = 50.4% done in 4.0s, 3.9s to go
on time step 6290 (time=104.833), 0.00063597 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 5.91278e-05 s
Working in 2D dimensions.
Computational cell is 4.4 x 6.5 x 0 with resolution 30
     block, center = (-1.2,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (-1.66533e-16,0,0)
          size (0.4,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 3.38153 s
-----------
Meep progress: 100.73333333333333/408.0 = 24.7% done in 4.0s, 12.2s to go
on time step 6044 (time=100.733), 0.00066186 s/step
Meep progress: 201.35/408.0 = 49.4% done in 8.0s, 8.2s to go
on time step 12083 (time=201.383), 0.000662443 s/step
Meep progress: 301.6333333333333/408.0 = 73.9% done in 12.0s, 4.2s to go
on time step 18102 (time=301.7), 0.000664572 s/step
Meep progress: 403.3333333333333/408.0 = 98.9% done in 16.0s, 0.2s to go
on time step 24205 (time=403.417), 0.000655435 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 50 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 37 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 66 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.79050
tran (uniaxial):, 1, 4.77, 0.06529
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 6.79493e-05 s
Working in 2D dimensions.
Computational cell is 4.8 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.0874732 s
-----------
Meep progress: 102.21666666666667/208.0 = 49.1% done in 4.0s, 4.1s to go
on time step 6133 (time=102.217), 0.000652292 s/step
Meep progress: 205.86666666666667/208.0 = 99.0% done in 8.0s, 0.1s to go
on time step 12353 (time=205.883), 0.000643173 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.41346e-05 s
Working in 2D dimensions.
Computational cell is 4.8 x 6.5 x 0 with resolution 30
     block, center = (-1.4,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (1.11022e-16,0,0)
          size (0.8,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 6.61127 s
-----------
Meep progress: 97.56666666666666/408.0 = 23.9% done in 4.0s, 12.7s to go
on time step 5854 (time=97.5667), 0.000683333 s/step
Meep progress: 195.35/408.0 = 47.9% done in 8.0s, 8.7s to go
on time step 11722 (time=195.367), 0.000681668 s/step
Meep progress: 293.75/408.0 = 72.0% done in 12.0s, 4.7s to go
on time step 17627 (time=293.783), 0.000677395 s/step
Meep progress: 391.4/408.0 = 95.9% done in 16.0s, 0.7s to go
on time step 23487 (time=391.45), 0.000682659 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 38 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 36 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 34 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.64693
tran (twisted):, 1, 4.77, 0.14670
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 7.70092e-05 s
Working in 2D dimensions.
Computational cell is 4.6 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.0854559 s
-----------
Meep progress: 102.86666666666666/208.0 = 49.5% done in 4.0s, 4.1s to go
on time step 6172 (time=102.867), 0.000648176 s/step
Meep progress: 207.0/208.0 = 99.5% done in 8.0s, 0.0s to go
on time step 12421 (time=207.017), 0.000640153 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.69956e-05 s
Working in 2D dimensions.
Computational cell is 4.6 x 6.5 x 0 with resolution 30
     block, center = (-1.3,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (2.22045e-16,0,0)
          size (0.6,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 4.98919 s
-----------
Meep progress: 98.01666666666667/408.0 = 24.0% done in 4.0s, 12.7s to go
on time step 5881 (time=98.0167), 0.000680194 s/step
Meep progress: 197.98333333333332/408.0 = 48.5% done in 8.0s, 8.5s to go
on time step 11880 (time=198), 0.00066678 s/step
Meep progress: 297.0833333333333/408.0 = 72.8% done in 12.0s, 4.5s to go
on time step 17827 (time=297.117), 0.000672654 s/step
Meep progress: 395.51666666666665/408.0 = 96.9% done in 16.0s, 0.5s to go
on time step 23735 (time=395.583), 0.000677137 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 70 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 31 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 70 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 50 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 44 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 33 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.66119
tran (uniaxial):, 1, 4.77, 0.14065
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 2.69413e-05 s
Working in 2D dimensions.
Computational cell is 5.2 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.090482 s
-----------
Meep progress: 97.23333333333333/208.0 = 46.7% done in 4.0s, 4.6s to go
on time step 5834 (time=97.2333), 0.000685701 s/step
Meep progress: 196.66666666666666/208.0 = 94.6% done in 8.0s, 0.5s to go
on time step 11801 (time=196.683), 0.000670365 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.22272e-05 s
Working in 2D dimensions.
Computational cell is 5.2 x 6.5 x 0 with resolution 30
     block, center = (-1.6,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (1.2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 10.1716 s
-----------
Meep progress: 89.91666666666667/408.0 = 22.0% done in 4.0s, 14.2s to go
on time step 5395 (time=89.9167), 0.000741585 s/step
Meep progress: 183.15/408.0 = 44.9% done in 8.0s, 9.8s to go
on time step 10990 (time=183.167), 0.000715001 s/step
Meep progress: 274.8333333333333/408.0 = 67.4% done in 12.0s, 5.8s to go
on time step 16493 (time=274.883), 0.000726979 s/step
Meep progress: 367.5333333333333/408.0 = 90.1% done in 16.0s, 1.8s to go
on time step 22056 (time=367.6), 0.000719077 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 36 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 27 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 50 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 31 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.39057
tran (twisted):, 1, 4.77, 0.27182
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.000102043 s
Working in 2D dimensions.
Computational cell is 4.8 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.0766759 s
-----------
Meep progress: 102.18333333333334/208.0 = 49.1% done in 4.0s, 4.1s to go
on time step 6131 (time=102.183), 0.000652499 s/step
Meep progress: 205.46666666666667/208.0 = 98.8% done in 8.0s, 0.1s to go
on time step 12330 (time=205.5), 0.000645285 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.38962e-05 s
Working in 2D dimensions.
Computational cell is 4.8 x 6.5 x 0 with resolution 30
     block, center = (-1.4,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (1.11022e-16,0,0)
          size (0.8,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 6.60763 s
-----------
Meep progress: 93.03333333333333/408.0 = 22.8% done in 4.0s, 13.5s to go
on time step 5582 (time=93.0333), 0.00071674 s/step
Meep progress: 188.96666666666667/408.0 = 46.3% done in 8.0s, 9.3s to go
on time step 11340 (time=189), 0.000694738 s/step
Meep progress: 285.1666666666667/408.0 = 69.9% done in 12.0s, 5.2s to go
on time step 17113 (time=285.217), 0.0006929 s/step
Meep progress: 381.0/408.0 = 93.4% done in 16.0s, 1.1s to go
on time step 22865 (time=381.083), 0.000695525 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 58 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 54 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.47140
tran (uniaxial):, 1, 4.77, 0.23506
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 8.98838e-05 s
Working in 2D dimensions.
Computational cell is 5.6 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.100205 s
-----------
Meep progress: 92.38333333333333/208.0 = 44.4% done in 4.0s, 5.0s to go
on time step 5543 (time=92.3833), 0.000721698 s/step
Meep progress: 186.2/208.0 = 89.5% done in 8.0s, 0.9s to go
on time step 11174 (time=186.233), 0.000710356 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.50883e-05 s
Working in 2D dimensions.
Computational cell is 5.6 x 6.5 x 0 with resolution 30
     block, center = (-1.8,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (2.22045e-16,0,0)
          size (1.6,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 71.5192% done, 1.70329 s remaining
subpixel-averaging is 71.5192% done, 1.68457 s remaining
subpixel-averaging is 71.5192% done, 1.68515 s remaining
time for set_epsilon = 13.1647 s
-----------
Meep progress: 87.73333333333333/408.0 = 21.5% done in 4.0s, 14.6s to go
on time step 5264 (time=87.7333), 0.00076 s/step
Meep progress: 176.35/408.0 = 43.2% done in 8.0s, 10.5s to go
on time step 10583 (time=176.383), 0.000752144 s/step
Meep progress: 265.5833333333333/408.0 = 65.1% done in 12.0s, 6.4s to go
on time step 15938 (time=265.633), 0.000747004 s/step
Meep progress: 353.1166666666667/408.0 = 86.5% done in 16.0s, 2.5s to go
on time step 21191 (time=353.183), 0.000761606 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 51 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 34 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 64 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 43 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 31 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 33 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.17500
tran (twisted):, 1, 4.77, 0.38032
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 7.10487e-05 s
Working in 2D dimensions.
Computational cell is 5 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.0600371 s
-----------
Meep progress: 100.0/208.0 = 48.1% done in 4.0s, 4.3s to go
on time step 6000 (time=100), 0.000666774 s/step
Meep progress: 201.83333333333334/208.0 = 97.0% done in 8.0s, 0.2s to go
on time step 12113 (time=201.883), 0.000654451 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.10352e-05 s
Working in 2D dimensions.
Computational cell is 5 x 6.5 x 0 with resolution 30
     block, center = (-1.5,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (1,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 8.54403 s
-----------
Meep progress: 93.58333333333333/408.0 = 22.9% done in 4.0s, 13.4s to go
on time step 5615 (time=93.5833), 0.000712491 s/step
Meep progress: 189.16666666666666/408.0 = 46.4% done in 8.0s, 9.3s to go
on time step 11352 (time=189.2), 0.000697323 s/step
Meep progress: 283.1666666666667/408.0 = 69.4% done in 12.0s, 5.3s to go
on time step 16993 (time=283.217), 0.000709178 s/step
Meep progress: 378.1166666666667/408.0 = 92.7% done in 16.0s, 1.3s to go
on time step 22691 (time=378.183), 0.000702061 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 30 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 35 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 61 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 44 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.27041
tran (uniaxial):, 1, 4.77, 0.32724
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 3.60012e-05 s
Working in 2D dimensions.
Computational cell is 6 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.131748 s
-----------
Meep progress: 87.66666666666667/208.0 = 42.1% done in 4.0s, 5.5s to go
on time step 5260 (time=87.6667), 0.000760564 s/step
Meep progress: 179.5/208.0 = 86.3% done in 8.0s, 1.3s to go
on time step 10772 (time=179.533), 0.000725698 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.29425e-05 s
Working in 2D dimensions.
Computational cell is 6 x 6.5 x 0 with resolution 30
     block, center = (-2,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 68.5946% done, 1.84034 s remaining
subpixel-averaging is 64.3072% done, 2.3664 s remaining
subpixel-averaging is 64.3072% done, 2.43912 s remaining
time for set_epsilon = 15.5353 s
-----------
Meep progress: 83.8/408.0 = 20.5% done in 4.0s, 15.5s to go
on time step 5028 (time=83.8), 0.000795612 s/step
Meep progress: 168.46666666666667/408.0 = 41.3% done in 8.0s, 11.4s to go
on time step 10110 (time=168.5), 0.000787236 s/step
Meep progress: 253.21666666666667/408.0 = 62.1% done in 12.0s, 7.3s to go
on time step 15197 (time=253.283), 0.000786469 s/step
Meep progress: 337.46666666666664/408.0 = 82.7% done in 16.0s, 3.3s to go
on time step 20253 (time=337.55), 0.00079122 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 38 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 46 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 28 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.04972
tran (twisted):, 1, 4.77, 0.44671
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.000128031 s
Working in 2D dimensions.
Computational cell is 5.2 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.082103 s
-----------
Meep progress: 98.0/208.0 = 47.1% done in 4.0s, 4.5s to go
on time step 5880 (time=98), 0.000680371 s/step
Meep progress: 197.83333333333334/208.0 = 95.1% done in 8.0s, 0.4s to go
on time step 11871 (time=197.85), 0.000667708 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.10352e-05 s
Working in 2D dimensions.
Computational cell is 5.2 x 6.5 x 0 with resolution 30
     block, center = (-1.6,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (-1.11022e-16,0,0)
          size (1.2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 9.882 s
-----------
Meep progress: 92.15/408.0 = 22.6% done in 4.0s, 13.7s to go
on time step 5529 (time=92.15), 0.000723491 s/step
Meep progress: 185.61666666666667/408.0 = 45.5% done in 8.0s, 9.6s to go
on time step 11139 (time=185.65), 0.00071314 s/step
Meep progress: 278.5833333333333/408.0 = 68.3% done in 12.0s, 5.6s to go
on time step 16719 (time=278.65), 0.000716963 s/step
Meep progress: 371.01666666666665/408.0 = 90.9% done in 16.0s, 1.6s to go
on time step 22267 (time=371.117), 0.000721017 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 42 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 44 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 33 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 62 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 37 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.14816
tran (uniaxial):, 1, 4.77, 0.39120
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.000138044 s
Working in 2D dimensions.
Computational cell is 6.4 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.109352 s
-----------
Meep progress: 85.88333333333333/208.0 = 41.3% done in 4.0s, 5.7s to go
on time step 5153 (time=85.8833), 0.00077639 s/step
Meep progress: 172.66666666666666/208.0 = 83.0% done in 8.0s, 1.6s to go
on time step 10361 (time=172.683), 0.000768107 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.69956e-05 s
Working in 2D dimensions.
Computational cell is 6.4 x 6.5 x 0 with resolution 30
     block, center = (-2.2,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (-2.22045e-16,0,0)
          size (2.4,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 58.4164% done, 3.10351 s remaining
subpixel-averaging is 58.4164% done, 3.12547 s remaining
subpixel-averaging is 58.4164% done, 3.0991 s remaining
time for set_epsilon = 19.9794 s
-----------
Meep progress: 80.01666666666667/408.0 = 19.6% done in 4.0s, 16.4s to go
on time step 4801 (time=80.0167), 0.000833327 s/step
Meep progress: 160.8/408.0 = 39.4% done in 8.0s, 12.3s to go
on time step 9649 (time=160.817), 0.000825088 s/step
Meep progress: 241.7/408.0 = 59.2% done in 12.0s, 8.3s to go
on time step 14504 (time=241.733), 0.000823907 s/step
Meep progress: 320.8333333333333/408.0 = 78.6% done in 16.0s, 4.3s to go
on time step 19253 (time=320.883), 0.000842394 s/step
Meep progress: 401.3/408.0 = 98.4% done in 20.0s, 0.3s to go
on time step 24084 (time=401.4), 0.000828145 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 65 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 44 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 29 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 42 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 27 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.00536
tran (twisted):, 1, 4.77, 0.46488
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 6.10352e-05 s
Working in 2D dimensions.
Computational cell is 5.4 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.0821059 s
-----------
Meep progress: 93.78333333333333/208.0 = 45.1% done in 4.0s, 4.9s to go
on time step 5627 (time=93.7833), 0.000710922 s/step
Meep progress: 189.56666666666666/208.0 = 91.1% done in 8.0s, 0.8s to go
on time step 11376 (time=189.6), 0.000695889 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.29425e-05 s
Working in 2D dimensions.
Computational cell is 5.4 x 6.5 x 0 with resolution 30
     block, center = (-1.7,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (-1.11022e-16,0,0)
          size (1.4,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
time for set_epsilon = 11.8445 s
-----------
Meep progress: 89.43333333333334/408.0 = 21.9% done in 4.0s, 14.2s to go
on time step 5366 (time=89.4333), 0.000745554 s/step
Meep progress: 178.88333333333333/408.0 = 43.8% done in 8.0s, 10.2s to go
on time step 10735 (time=178.917), 0.000745129 s/step
Meep progress: 269.46666666666664/408.0 = 66.0% done in 12.0s, 6.2s to go
on time step 16172 (time=269.533), 0.000735702 s/step
Meep progress: 359.71666666666664/408.0 = 88.2% done in 16.0s, 2.1s to go
on time step 21588 (time=359.8), 0.000738555 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 27 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 63 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 42 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.02627
tran (uniaxial):, 1, 4.77, 0.45435
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 7.79629e-05 s
Working in 2D dimensions.
Computational cell is 6.8 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.105287 s
-----------
Meep progress: 82.33333333333333/208.0 = 39.6% done in 4.0s, 6.1s to go
on time step 4940 (time=82.3333), 0.000809818 s/step
Meep progress: 165.56666666666666/208.0 = 79.6% done in 8.0s, 2.1s to go
on time step 9936 (time=165.6), 0.000800751 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 7.00951e-05 s
Working in 2D dimensions.
Computational cell is 6.8 x 6.5 x 0 with resolution 30
     block, center = (-2.4,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (-2.22045e-16,0,0)
          size (2.8,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 53.5143% done, 3.8263 s remaining
subpixel-averaging is 53.5143% done, 3.81257 s remaining
subpixel-averaging is 53.5143% done, 3.81274 s remaining
time for set_epsilon = 23.647 s
-----------
Meep progress: 74.46666666666667/408.0 = 18.3% done in 4.0s, 17.9s to go
on time step 4468 (time=74.4667), 0.000895362 s/step
Meep progress: 150.38333333333333/408.0 = 36.9% done in 8.0s, 13.7s to go
on time step 9024 (time=150.4), 0.000878078 s/step
Meep progress: 226.41666666666666/408.0 = 55.5% done in 12.0s, 9.6s to go
on time step 13587 (time=226.45), 0.000876645 s/step
Meep progress: 302.23333333333335/408.0 = 74.1% done in 16.0s, 5.6s to go
on time step 18137 (time=302.283), 0.000879289 s/step
Meep progress: 377.98333333333335/408.0 = 92.6% done in 20.0s, 1.6s to go
on time step 22683 (time=378.05), 0.000879927 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 65 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 40 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 37 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.00085
tran (twisted):, 1, 4.77, 0.46955
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.000133991 s
Working in 2D dimensions.
Computational cell is 5.6 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.096647 s
-----------
Meep progress: 93.56666666666666/208.0 = 45.0% done in 4.0s, 4.9s to go
on time step 5614 (time=93.5667), 0.000712629 s/step
Meep progress: 188.45/208.0 = 90.6% done in 8.0s, 0.8s to go
on time step 11308 (time=188.467), 0.000702542 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.29425e-05 s
Working in 2D dimensions.
Computational cell is 5.6 x 6.5 x 0 with resolution 30
     block, center = (-1.8,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (2.22045e-16,0,0)
          size (1.6,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 71.5192% done, 1.72158 s remaining
subpixel-averaging is 71.5192% done, 1.70206 s remaining
subpixel-averaging is 71.5192% done, 1.70269 s remaining
time for set_epsilon = 13.3156 s
-----------
Meep progress: 87.53333333333333/408.0 = 21.5% done in 4.0s, 14.6s to go
on time step 5252 (time=87.5333), 0.00076165 s/step
Meep progress: 176.3/408.0 = 43.2% done in 8.0s, 10.5s to go
on time step 10580 (time=176.333), 0.000750865 s/step
Meep progress: 263.75/408.0 = 64.6% done in 12.0s, 6.6s to go
on time step 15829 (time=263.817), 0.000762183 s/step
Meep progress: 352.4/408.0 = 86.4% done in 16.0s, 2.5s to go
on time step 21149 (time=352.483), 0.000751941 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 40 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 42 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 37 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 46 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.00005
tran (uniaxial):, 1, 4.77, 0.46429
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 6.69956e-05 s
Working in 2D dimensions.
Computational cell is 7.2 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.0993609 s
-----------
Meep progress: 78.7/208.0 = 37.8% done in 4.0s, 6.6s to go
on time step 4722 (time=78.7), 0.000847187 s/step
Meep progress: 160.36666666666667/208.0 = 77.1% done in 8.0s, 2.4s to go
on time step 9624 (time=160.4), 0.000816089 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.81877e-05 s
Working in 2D dimensions.
Computational cell is 7.2 x 6.5 x 0 with resolution 30
     block, center = (-2.6,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (3.2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 49.3713% done, 4.3752 s remaining
subpixel-averaging is 78.9961% done, 1.1067 s remaining
subpixel-averaging is 49.3713% done, 4.42276 s remaining
subpixel-averaging is 78.9961% done, 1.10731 s remaining
subpixel-averaging is 49.3713% done, 4.43155 s remaining
subpixel-averaging is 78.9961% done, 1.11122 s remaining
time for set_epsilon = 26.3941 s
-----------
Meep progress: 72.93333333333334/408.0 = 17.9% done in 4.0s, 18.4s to go
on time step 4376 (time=72.9333), 0.000914306 s/step
Meep progress: 148.9/408.0 = 36.5% done in 8.0s, 13.9s to go
on time step 8935 (time=148.917), 0.000877482 s/step
Meep progress: 224.29999999999998/408.0 = 55.0% done in 12.0s, 9.8s to go
on time step 13460 (time=224.333), 0.000884124 s/step
Meep progress: 297.3333333333333/408.0 = 72.9% done in 16.0s, 6.0s to go
on time step 17843 (time=297.383), 0.00091272 s/step
Meep progress: 372.43333333333334/408.0 = 91.3% done in 20.0s, 1.9s to go
on time step 22351 (time=372.517), 0.000887407 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 58 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 52 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 52 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 42 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 27 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.00061
tran (twisted):, 1, 4.77, 0.46431
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.000139952 s
Working in 2D dimensions.
Computational cell is 5.8 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.103642 s
-----------
Meep progress: 92.38333333333333/208.0 = 44.4% done in 4.0s, 5.0s to go
on time step 5543 (time=92.3833), 0.000721725 s/step
Meep progress: 186.01666666666665/208.0 = 89.4% done in 8.0s, 0.9s to go
on time step 11163 (time=186.05), 0.000711761 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.69956e-05 s
Working in 2D dimensions.
Computational cell is 5.8 x 6.5 x 0 with resolution 30
     block, center = (-1.9,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (1.11022e-16,0,0)
          size (1.8,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 67.7217% done, 2.0545 s remaining
subpixel-averaging is 67.7217% done, 2.04229 s remaining
subpixel-averaging is 67.7217% done, 2.04423 s remaining
time for set_epsilon = 14.95 s
-----------
Meep progress: 86.31666666666666/408.0 = 21.2% done in 4.0s, 14.9s to go
on time step 5179 (time=86.3167), 0.00077244 s/step
Meep progress: 173.38333333333333/408.0 = 42.5% done in 8.0s, 10.8s to go
on time step 10404 (time=173.4), 0.000765589 s/step
Meep progress: 261.1/408.0 = 64.0% done in 12.0s, 6.8s to go
on time step 15668 (time=261.133), 0.000759921 s/step
Meep progress: 347.5333333333333/408.0 = 85.2% done in 16.0s, 2.8s to go
on time step 20857 (time=347.617), 0.000770998 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 52 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 40 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 62 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.02864
tran (uniaxial):, 1, 4.77, 0.44626
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 5.00679e-05 s
Working in 2D dimensions.
Computational cell is 7.6 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.137346 s
-----------
Meep progress: 78.08333333333333/208.0 = 37.5% done in 4.0s, 6.7s to go
on time step 4685 (time=78.0833), 0.000853827 s/step
Meep progress: 156.88333333333333/208.0 = 75.4% done in 8.0s, 2.6s to go
on time step 9414 (time=156.9), 0.000845896 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.69956e-05 s
Working in 2D dimensions.
Computational cell is 7.6 x 6.5 x 0 with resolution 30
     block, center = (-2.8,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (2.22045e-16,0,0)
          size (3.6,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 45.8237% done, 5.1262 s remaining
subpixel-averaging is 73.3197% done, 1.54776 s remaining
subpixel-averaging is 45.8237% done, 5.0838 s remaining
subpixel-averaging is 73.3197% done, 1.54977 s remaining
subpixel-averaging is 45.8237% done, 5.10381 s remaining
subpixel-averaging is 73.3197% done, 1.55772 s remaining
time for set_epsilon = 30.0956 s
-----------
Meep progress: 70.91666666666667/408.0 = 17.4% done in 4.0s, 19.0s to go
on time step 4255 (time=70.9167), 0.000940257 s/step
Meep progress: 142.71666666666667/408.0 = 35.0% done in 8.0s, 14.9s to go
on time step 8564 (time=142.733), 0.000928472 s/step
Meep progress: 213.79999999999998/408.0 = 52.4% done in 12.0s, 10.9s to go
on time step 12830 (time=213.833), 0.000937849 s/step
Meep progress: 284.56666666666666/408.0 = 69.7% done in 16.0s, 6.9s to go
on time step 17078 (time=284.633), 0.000941825 s/step
Meep progress: 355.2833333333333/408.0 = 87.1% done in 20.0s, 3.0s to go
on time step 21322 (time=355.367), 0.000942543 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 31 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 28 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.00230
tran (twisted):, 1, 4.77, 0.45960
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 5.38826e-05 s
Working in 2D dimensions.
Computational cell is 6 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.0905509 s
-----------
Meep progress: 90.25/208.0 = 43.4% done in 4.0s, 5.2s to go
on time step 5415 (time=90.25), 0.000738743 s/step
Meep progress: 181.95/208.0 = 87.5% done in 8.0s, 1.1s to go
on time step 10918 (time=181.967), 0.000726975 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.41346e-05 s
Working in 2D dimensions.
Computational cell is 6 x 6.5 x 0 with resolution 30
     block, center = (-2,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 64.3072% done, 2.38699 s remaining
subpixel-averaging is 64.3072% done, 2.40457 s remaining
subpixel-averaging is 64.3072% done, 2.40796 s remaining
time for set_epsilon = 16.6695 s
-----------
Meep progress: 82.75/408.0 = 20.3% done in 4.0s, 15.7s to go
on time step 4965 (time=82.75), 0.000805785 s/step
Meep progress: 175.51666666666665/408.0 = 43.0% done in 8.0s, 10.6s to go
on time step 10533 (time=175.55), 0.000718519 s/step
Meep progress: 259.15/408.0 = 63.5% done in 12.0s, 6.9s to go
on time step 15553 (time=259.217), 0.000796936 s/step
Meep progress: 342.56666666666666/408.0 = 84.0% done in 16.0s, 3.1s to go
on time step 20559 (time=342.65), 0.00079916 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 68 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 30 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 34 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 55 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 42 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 34 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.15199
tran (uniaxial):, 1, 4.77, 0.39654
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 8.79765e-05 s
Working in 2D dimensions.
Computational cell is 8 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.125251 s
-----------
Meep progress: 75.01666666666667/208.0 = 36.1% done in 4.0s, 7.1s to go
on time step 4501 (time=75.0167), 0.000888818 s/step
Meep progress: 150.11666666666667/208.0 = 72.2% done in 8.0s, 3.1s to go
on time step 9009 (time=150.15), 0.000887502 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 7.39098e-05 s
Working in 2D dimensions.
Computational cell is 8 x 6.5 x 0 with resolution 30
     block, center = (-3,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (4,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 42.7517% done, 5.69995 s remaining
subpixel-averaging is 68.4044% done, 1.92679 s remaining
subpixel-averaging is 42.7517% done, 5.7258 s remaining
subpixel-averaging is 68.4044% done, 1.91778 s remaining
subpixel-averaging is 42.7517% done, 5.72558 s remaining
subpixel-averaging is 68.4044% done, 1.9235 s remaining
time for set_epsilon = 32.8454 s
-----------
Meep progress: 67.6/408.0 = 16.6% done in 4.0s, 20.1s to go
on time step 4056 (time=67.6), 0.000986321 s/step
Meep progress: 136.23333333333332/408.0 = 33.4% done in 8.0s, 16.0s to go
on time step 8175 (time=136.25), 0.00097132 s/step
Meep progress: 204.93333333333334/408.0 = 50.2% done in 12.0s, 11.9s to go
on time step 12298 (time=204.967), 0.000970339 s/step
Meep progress: 273.5833333333333/408.0 = 67.1% done in 16.0s, 7.9s to go
on time step 16418 (time=273.633), 0.000971036 s/step
Meep progress: 341.8666666666667/408.0 = 83.8% done in 20.0s, 3.9s to go
on time step 20515 (time=341.917), 0.000976364 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 42 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 37 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 30 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.05159
tran (twisted):, 1, 4.77, 0.43938
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 5.29289e-05 s
Working in 2D dimensions.
Computational cell is 6.2 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.104505 s
-----------
Meep progress: 86.15/208.0 = 41.4% done in 4.0s, 5.7s to go
on time step 5169 (time=86.15), 0.000773898 s/step
Meep progress: 174.43333333333334/208.0 = 83.9% done in 8.0s, 1.5s to go
on time step 10468 (time=174.467), 0.000754955 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.50883e-05 s
Working in 2D dimensions.
Computational cell is 6.2 x 6.5 x 0 with resolution 30
     block, center = (-2.1,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (2.2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 61.2204% done, 2.79341 s remaining
subpixel-averaging is 61.2204% done, 2.64891 s remaining
subpixel-averaging is 61.2204% done, 2.75618 s remaining
time for set_epsilon = 18.3206 s
-----------
Meep progress: 81.01666666666667/408.0 = 19.9% done in 4.0s, 16.1s to go
on time step 4861 (time=81.0167), 0.000823039 s/step
Meep progress: 163.05/408.0 = 40.0% done in 8.0s, 12.0s to go
on time step 9784 (time=163.067), 0.00081255 s/step
Meep progress: 246.23333333333332/408.0 = 60.4% done in 12.0s, 7.9s to go
on time step 14776 (time=246.267), 0.000801294 s/step
Meep progress: 327.5333333333333/408.0 = 80.3% done in 16.0s, 3.9s to go
on time step 19655 (time=327.583), 0.000819844 s/step
Meep progress: 407.8666666666667/408.0 = 100.0% done in 20.0s, 0.0s to go
on time step 24477 (time=407.95), 0.000829674 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 51 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 49 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 43 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 36 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 46 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.30148
tran (uniaxial):, 1, 4.77, 0.31081
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 9.58443e-05 s
Working in 2D dimensions.
Computational cell is 8.4 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.128375 s
-----------
Meep progress: 72.65/208.0 = 34.9% done in 4.0s, 7.5s to go
on time step 4359 (time=72.65), 0.000917816 s/step
Meep progress: 144.25/208.0 = 69.4% done in 8.0s, 3.5s to go
on time step 8656 (time=144.267), 0.000930968 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.69956e-05 s
Working in 2D dimensions.
Computational cell is 8.4 x 6.5 x 0 with resolution 30
     block, center = (-3.2,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (4.4,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 40.0657% done, 6.53195 s remaining
subpixel-averaging is 64.1067% done, 2.36242 s remaining
subpixel-averaging is 40.0657% done, 6.48991 s remaining
subpixel-averaging is 64.1067% done, 2.36408 s remaining
subpixel-averaging is 40.0657% done, 6.52898 s remaining
subpixel-averaging is 64.1067% done, 2.35834 s remaining
time for set_epsilon = 36.7265 s
-----------
Meep progress: 65.65/408.0 = 16.1% done in 4.0s, 20.9s to go
on time step 3939 (time=65.65), 0.00101563 s/step
Meep progress: 132.1/408.0 = 32.4% done in 8.0s, 16.7s to go
on time step 7927 (time=132.117), 0.00100305 s/step
Meep progress: 198.91666666666666/408.0 = 48.8% done in 12.0s, 12.6s to go
on time step 11937 (time=198.95), 0.000997598 s/step
Meep progress: 265.51666666666665/408.0 = 65.1% done in 16.0s, 8.6s to go
on time step 15934 (time=265.567), 0.00100085 s/step
Meep progress: 331.93333333333334/408.0 = 81.4% done in 20.0s, 4.6s to go
on time step 19920 (time=332), 0.00100371 s/step
Meep progress: 397.23333333333335/408.0 = 97.4% done in 24.0s, 0.7s to go
on time step 23839 (time=397.317), 0.00102074 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 47 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 43 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 50 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.17455
tran (twisted):, 1, 4.77, 0.38170
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 4.50611e-05 s
Working in 2D dimensions.
Computational cell is 6.4 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.11933 s
-----------
Meep progress: 86.23333333333333/208.0 = 41.5% done in 4.0s, 5.6s to go
on time step 5174 (time=86.2333), 0.000773265 s/step
Meep progress: 173.26666666666665/208.0 = 83.3% done in 8.0s, 1.6s to go
on time step 10398 (time=173.3), 0.000765817 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.50883e-05 s
Working in 2D dimensions.
Computational cell is 6.4 x 6.5 x 0 with resolution 30
     block, center = (-2.2,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (2.4,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 58.4164% done, 3.09306 s remaining
subpixel-averaging is 58.4164% done, 2.89371 s remaining
subpixel-averaging is 58.4164% done, 3.09871 s remaining
time for set_epsilon = 19.6879 s
-----------
Meep progress: 80.91666666666667/408.0 = 19.8% done in 4.0s, 16.2s to go
on time step 4855 (time=80.9167), 0.000824082 s/step
Meep progress: 162.7/408.0 = 39.9% done in 8.0s, 12.1s to go
on time step 9763 (time=162.717), 0.000815005 s/step
Meep progress: 244.65/408.0 = 60.0% done in 12.0s, 8.0s to go
on time step 14681 (time=244.683), 0.000813452 s/step
Meep progress: 326.3833333333333/408.0 = 80.0% done in 16.0s, 4.0s to go
on time step 19586 (time=326.433), 0.00081554 s/step
Meep progress: 407.9166666666667/408.0 = 100.0% done in 20.0s, 0.0s to go
on time step 24479 (time=407.983), 0.000817639 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 69 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 49 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 37 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 51 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 36 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 28 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.46820
tran (uniaxial):, 1, 4.77, 0.22878
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.000147104 s
Working in 2D dimensions.
Computational cell is 8.8 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.139945 s
-----------
Meep progress: 69.91666666666667/208.0 = 33.6% done in 4.0s, 7.9s to go
on time step 4195 (time=69.9167), 0.000953672 s/step
Meep progress: 140.83333333333334/208.0 = 67.7% done in 8.0s, 3.8s to go
on time step 8450 (time=140.833), 0.000940074 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 7.10487e-05 s
Working in 2D dimensions.
Computational cell is 8.8 x 6.5 x 0 with resolution 30
     block, center = (-3.4,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (4.8,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 37.6973% done, 7.13631 s remaining
subpixel-averaging is 60.3172% done, 2.75956 s remaining
subpixel-averaging is 82.9371% done, 0.865693 s remaining
subpixel-averaging is 37.6973% done, 7.11049 s remaining
subpixel-averaging is 60.3172% done, 2.75036 s remaining
subpixel-averaging is 82.9371% done, 0.858843 s remaining
subpixel-averaging is 37.6973% done, 7.12251 s remaining
subpixel-averaging is 60.3172% done, 2.75345 s remaining
subpixel-averaging is 82.9371% done, 0.866829 s remaining
time for set_epsilon = 39.6605 s
-----------
Meep progress: 61.016666666666666/408.0 = 15.0% done in 4.0s, 22.7s to go
on time step 3661 (time=61.0167), 0.00109264 s/step
Meep progress: 122.33333333333333/408.0 = 30.0% done in 8.0s, 18.7s to go
on time step 7342 (time=122.367), 0.00108692 s/step
Meep progress: 185.16666666666666/408.0 = 45.4% done in 12.0s, 14.4s to go
on time step 11113 (time=185.217), 0.00106086 s/step
Meep progress: 248.03333333333333/408.0 = 60.8% done in 16.0s, 10.3s to go
on time step 14886 (time=248.1), 0.00106035 s/step
Meep progress: 310.8/408.0 = 76.2% done in 20.0s, 6.3s to go
on time step 18653 (time=310.883), 0.00106216 s/step
Meep progress: 373.48333333333335/408.0 = 91.5% done in 24.0s, 2.2s to go
on time step 22416 (time=373.6), 0.0010633 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 37 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 27 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 51 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.40540
tran (twisted):, 1, 4.77, 0.26104
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 7.20024e-05 s
Working in 2D dimensions.
Computational cell is 6.6 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.098484 s
-----------
Meep progress: 84.7/208.0 = 40.7% done in 4.0s, 5.8s to go
on time step 5082 (time=84.7), 0.000787265 s/step
Meep progress: 168.76666666666665/208.0 = 81.1% done in 8.0s, 1.9s to go
on time step 10127 (time=168.783), 0.000792943 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.69956e-05 s
Working in 2D dimensions.
Computational cell is 6.6 x 6.5 x 0 with resolution 30
     block, center = (-2.3,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (2.6,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 55.858% done, 3.41444 s remaining
subpixel-averaging is 55.858% done, 3.40405 s remaining
subpixel-averaging is 55.858% done, 3.39686 s remaining
time for set_epsilon = 21.5757 s
-----------
Meep progress: 78.68333333333334/408.0 = 19.3% done in 4.0s, 16.7s to go
on time step 4721 (time=78.6833), 0.000847361 s/step
Meep progress: 157.51666666666665/408.0 = 38.6% done in 8.0s, 12.7s to go
on time step 9452 (time=157.533), 0.000845538 s/step
Meep progress: 237.29999999999998/408.0 = 58.2% done in 12.0s, 8.6s to go
on time step 14240 (time=237.333), 0.000835482 s/step
Meep progress: 315.93333333333334/408.0 = 77.4% done in 16.0s, 4.7s to go
on time step 18959 (time=315.983), 0.000847783 s/step
Meep progress: 394.68333333333334/408.0 = 96.7% done in 20.0s, 0.7s to go
on time step 23686 (time=394.767), 0.000846355 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 53 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 38 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 70 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 47 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.66179
tran (uniaxial):, 1, 4.77, 0.14253
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.000156164 s
Working in 2D dimensions.
Computational cell is 9.2 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.172037 s
-----------
Meep progress: 68.06666666666666/208.0 = 32.7% done in 4.0s, 8.2s to go
on time step 4084 (time=68.0667), 0.000979654 s/step
Meep progress: 137.1/208.0 = 65.9% done in 8.0s, 4.1s to go
on time step 8226 (time=137.1), 0.000965725 s/step
Meep progress: 205.7/208.0 = 98.9% done in 12.0s, 0.1s to go
on time step 12343 (time=205.717), 0.000971618 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.98566e-05 s
Working in 2D dimensions.
Computational cell is 9.2 x 6.5 x 0 with resolution 30
     block, center = (-3.6,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (5.2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 35.5933% done, 7.99426 s remaining
subpixel-averaging is 56.9506% done, 3.19121 s remaining
subpixel-averaging is 78.308% done, 1.1792 s remaining
subpixel-averaging is 35.5933% done, 7.82132 s remaining
subpixel-averaging is 56.9506% done, 3.19462 s remaining
subpixel-averaging is 78.308% done, 1.17005 s remaining
subpixel-averaging is 35.5933% done, 7.79931 s remaining
subpixel-averaging is 56.9506% done, 3.1716 s remaining
subpixel-averaging is 78.308% done, 1.18297 s remaining
time for set_epsilon = 43.3607 s
-----------
Meep progress: 61.21666666666667/408.0 = 15.0% done in 4.0s, 22.7s to go
on time step 3673 (time=61.2167), 0.00108931 s/step
Meep progress: 122.93333333333334/408.0 = 30.1% done in 8.0s, 18.6s to go
on time step 7377 (time=122.95), 0.00108006 s/step
Meep progress: 183.58333333333334/408.0 = 45.0% done in 12.0s, 14.7s to go
on time step 11017 (time=183.617), 0.00109908 s/step
Meep progress: 244.83333333333334/408.0 = 60.0% done in 16.0s, 10.7s to go
on time step 14693 (time=244.883), 0.00108818 s/step
Meep progress: 306.65/408.0 = 75.2% done in 20.0s, 6.6s to go
on time step 18403 (time=306.717), 0.00107845 s/step
Meep progress: 367.8/408.0 = 90.1% done in 24.0s, 2.6s to go
on time step 22073 (time=367.883), 0.00109002 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 47 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.67195
tran (twisted):, 1, 4.77, 0.13686
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 9.08375e-05 s
Working in 2D dimensions.
Computational cell is 6.8 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.099133 s
-----------
Meep progress: 83.68333333333334/208.0 = 40.2% done in 4.0s, 5.9s to go
on time step 5021 (time=83.6833), 0.000796734 s/step
Meep progress: 168.46666666666667/208.0 = 81.0% done in 8.0s, 1.9s to go
on time step 10110 (time=168.5), 0.000786127 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.69956e-05 s
Working in 2D dimensions.
Computational cell is 6.8 x 6.5 x 0 with resolution 30
     block, center = (-2.4,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (-2.22045e-16,0,0)
          size (2.8,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 53.5143% done, 3.79773 s remaining
subpixel-averaging is 53.5143% done, 3.81162 s remaining
subpixel-averaging is 53.5143% done, 3.8089 s remaining
time for set_epsilon = 23.541 s
-----------
Meep progress: 74.93333333333334/408.0 = 18.4% done in 4.0s, 17.8s to go
on time step 4496 (time=74.9333), 0.000889767 s/step
Meep progress: 151.51666666666665/408.0 = 37.1% done in 8.0s, 13.5s to go
on time step 9092 (time=151.533), 0.000870396 s/step
Meep progress: 229.31666666666666/408.0 = 56.2% done in 12.0s, 9.4s to go
on time step 13761 (time=229.35), 0.000856824 s/step
Meep progress: 307.1666666666667/408.0 = 75.3% done in 16.0s, 5.3s to go
on time step 18434 (time=307.233), 0.00085611 s/step
Meep progress: 383.71666666666664/408.0 = 94.0% done in 20.0s, 1.3s to go
on time step 23028 (time=383.8), 0.000870738 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 69 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 34 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.80825
tran (uniaxial):, 1, 4.77, 0.06307
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.000152111 s
Working in 2D dimensions.
Computational cell is 9.6 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.157447 s
-----------
Meep progress: 64.9/208.0 = 31.2% done in 4.0s, 8.8s to go
on time step 3894 (time=64.9), 0.00102741 s/step
Meep progress: 130.88333333333333/208.0 = 62.9% done in 8.0s, 4.7s to go
on time step 7854 (time=130.9), 0.00101027 s/step
Meep progress: 196.55/208.0 = 94.5% done in 12.0s, 0.7s to go
on time step 11795 (time=196.583), 0.00101506 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 7.20024e-05 s
Working in 2D dimensions.
Computational cell is 9.6 x 6.5 x 0 with resolution 30
     block, center = (-3.8,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (-4.44089e-16,0,0)
          size (5.6,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 33.7117% done, 8.53515 s remaining
subpixel-averaging is 53.94% done, 3.58281 s remaining
subpixel-averaging is 74.1684% done, 1.47227 s remaining
subpixel-averaging is 33.7117% done, 8.64079 s remaining
subpixel-averaging is 53.94% done, 3.58974 s remaining
subpixel-averaging is 74.1684% done, 1.46986 s remaining
subpixel-averaging is 33.7117% done, 8.56235 s remaining
subpixel-averaging is 53.94% done, 3.61465 s remaining
subpixel-averaging is 74.1684% done, 1.47245 s remaining
time for set_epsilon = 46.6436 s
-----------
Meep progress: 58.0/408.0 = 14.2% done in 4.0s, 24.1s to go
on time step 3480 (time=58), 0.0011497 s/step
Meep progress: 116.86666666666666/408.0 = 28.6% done in 8.0s, 19.9s to go
on time step 7013 (time=116.883), 0.00113224 s/step
Meep progress: 175.05/408.0 = 42.9% done in 12.0s, 16.0s to go
on time step 10505 (time=175.083), 0.00114558 s/step
Meep progress: 233.93333333333334/408.0 = 57.3% done in 16.0s, 11.9s to go
on time step 14039 (time=233.983), 0.00113194 s/step
Meep progress: 293.5/408.0 = 71.9% done in 20.0s, 7.8s to go
on time step 17613 (time=293.55), 0.00111922 s/step
Meep progress: 352.01666666666665/408.0 = 86.3% done in 24.0s, 3.8s to go
on time step 21125 (time=352.083), 0.00113908 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 47 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 51 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 31 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.86580
tran (twisted):, 1, 4.77, 0.02804
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.000117064 s
Working in 2D dimensions.
Computational cell is 7 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.103678 s
-----------
Meep progress: 80.11666666666666/208.0 = 38.5% done in 4.0s, 6.4s to go
on time step 4807 (time=80.1167), 0.000832301 s/step
Meep progress: 161.7/208.0 = 77.7% done in 8.0s, 2.3s to go
on time step 9704 (time=161.733), 0.000816961 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 7.39098e-05 s
Working in 2D dimensions.
Computational cell is 7 x 6.5 x 0 with resolution 30
     block, center = (-2.5,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (2.22045e-16,0,0)
          size (3,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 51.3594% done, 4.1118 s remaining
subpixel-averaging is 51.3594% done, 4.09325 s remaining
subpixel-averaging is 51.3594% done, 4.09493 s remaining
time for set_epsilon = 24.965 s
-----------
Meep progress: 75.21666666666667/408.0 = 18.4% done in 4.0s, 17.7s to go
on time step 4513 (time=75.2167), 0.000886419 s/step
Meep progress: 150.08333333333334/408.0 = 36.8% done in 8.0s, 13.7s to go
on time step 9006 (time=150.1), 0.00089035 s/step
Meep progress: 222.83333333333334/408.0 = 54.6% done in 12.0s, 10.0s to go
on time step 13373 (time=222.883), 0.000916091 s/step
Meep progress: 297.3666666666667/408.0 = 72.9% done in 16.0s, 6.0s to go
on time step 17846 (time=297.433), 0.000894295 s/step
Meep progress: 371.5/408.0 = 91.1% done in 20.0s, 2.0s to go
on time step 22295 (time=371.583), 0.000899168 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 64 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 45 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.88675
tran (uniaxial):, 1, 4.77, 0.01655
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 7.20024e-05 s
Working in 2D dimensions.
Computational cell is 10 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.155606 s
-----------
Meep progress: 63.06666666666666/208.0 = 30.3% done in 4.0s, 9.2s to go
on time step 3784 (time=63.0667), 0.0010573 s/step
Meep progress: 127.85/208.0 = 61.5% done in 8.0s, 5.0s to go
on time step 7672 (time=127.867), 0.001029 s/step
Meep progress: 191.65/208.0 = 92.1% done in 12.0s, 1.0s to go
on time step 11501 (time=191.683), 0.00104491 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 7.00951e-05 s
Working in 2D dimensions.
Computational cell is 10 x 6.5 x 0 with resolution 30
     block, center = (-4,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (4.44089e-16,0,0)
          size (6,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 32.019% done, 9.22186 s remaining
subpixel-averaging is 51.2317% done, 3.98466 s remaining
subpixel-averaging is 70.4445% done, 1.77788 s remaining
subpixel-averaging is 32.019% done, 9.30279 s remaining
subpixel-averaging is 51.2317% done, 4.03405 s remaining
subpixel-averaging is 70.4445% done, 1.79632 s remaining
subpixel-averaging is 32.019% done, 9.17258 s remaining
subpixel-averaging is 51.2317% done, 4.02283 s remaining
subpixel-averaging is 70.4445% done, 1.78612 s remaining
time for set_epsilon = 50.0774 s
-----------
Meep progress: 55.21666666666667/408.0 = 13.5% done in 4.0s, 25.6s to go
on time step 3313 (time=55.2167), 0.00120767 s/step
Meep progress: 112.33333333333333/408.0 = 27.5% done in 8.0s, 21.1s to go
on time step 6741 (time=112.35), 0.00116717 s/step
Meep progress: 168.81666666666666/408.0 = 41.4% done in 12.0s, 17.0s to go
on time step 10131 (time=168.85), 0.00118012 s/step
Meep progress: 225.75/408.0 = 55.3% done in 16.0s, 12.9s to go
on time step 13549 (time=225.817), 0.00117055 s/step
Meep progress: 281.48333333333335/408.0 = 69.0% done in 20.0s, 9.0s to go
on time step 16894 (time=281.567), 0.00119605 s/step
Meep progress: 338.5333333333333/408.0 = 83.0% done in 24.0s, 4.9s to go
on time step 20317 (time=338.617), 0.00116858 s/step
Meep progress: 395.55/408.0 = 96.9% done in 28.0s, 0.9s to go
on time step 23739 (time=395.65), 0.00116892 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 65 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 45 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.94211
tran (twisted):, 1, 4.77, 0.00255
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.00011301 s
Working in 2D dimensions.
Computational cell is 7.2 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.111866 s
-----------
Meep progress: 80.61666666666666/208.0 = 38.8% done in 4.0s, 6.3s to go
on time step 4837 (time=80.6167), 0.000827066 s/step
Meep progress: 162.21666666666667/208.0 = 78.0% done in 8.0s, 2.3s to go
on time step 9734 (time=162.233), 0.000816884 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.91414e-05 s
Working in 2D dimensions.
Computational cell is 7.2 x 6.5 x 0 with resolution 30
     block, center = (-2.6,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (3.2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 49.3713% done, 4.4399 s remaining
subpixel-averaging is 78.9961% done, 1.12404 s remaining
subpixel-averaging is 49.3713% done, 4.44133 s remaining
subpixel-averaging is 78.9961% done, 1.12126 s remaining
subpixel-averaging is 49.3713% done, 4.42219 s remaining
subpixel-averaging is 78.9961% done, 1.12413 s remaining
time for set_epsilon = 26.6699 s
-----------
Meep progress: 73.56666666666666/408.0 = 18.0% done in 4.0s, 18.2s to go
on time step 4414 (time=73.5667), 0.00090629 s/step
Meep progress: 147.21666666666667/408.0 = 36.1% done in 8.0s, 14.2s to go
on time step 8834 (time=147.233), 0.000905054 s/step
Meep progress: 221.16666666666666/408.0 = 54.2% done in 12.0s, 10.1s to go
on time step 13272 (time=221.2), 0.000901458 s/step
Meep progress: 294.68333333333334/408.0 = 72.2% done in 16.0s, 6.2s to go
on time step 17684 (time=294.733), 0.000906649 s/step
Meep progress: 368.93333333333334/408.0 = 90.4% done in 20.0s, 2.1s to go
on time step 22140 (time=369), 0.000897794 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 45 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 28 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 65 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 37 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.93860
tran (uniaxial):, 1, 4.77, 0.00001
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 0.000154018 s
Working in 2D dimensions.
Computational cell is 10.4 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.161312 s
-----------
Meep progress: 60.9/208.0 = 29.3% done in 4.0s, 9.7s to go
on time step 3654 (time=60.9), 0.00109495 s/step
Meep progress: 122.28333333333333/208.0 = 58.8% done in 8.0s, 5.6s to go
on time step 7338 (time=122.3), 0.00108598 s/step
Meep progress: 183.23333333333332/208.0 = 88.1% done in 12.0s, 1.6s to go
on time step 10996 (time=183.267), 0.00109353 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 7.10487e-05 s
Working in 2D dimensions.
Computational cell is 10.4 x 6.5 x 0 with resolution 30
     block, center = (-4.2,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (6.4,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 30.4883% done, 9.98227 s remaining
subpixel-averaging is 48.7824% done, 4.42182 s remaining
subpixel-averaging is 67.0766% done, 2.07037 s remaining
subpixel-averaging is 85.3708% done, 0.733933 s remaining
subpixel-averaging is 30.4883% done, 9.99655 s remaining
subpixel-averaging is 48.7824% done, 4.44853 s remaining
subpixel-averaging is 67.0766% done, 2.11095 s remaining
subpixel-averaging is 85.3708% done, 0.733249 s remaining
subpixel-averaging is 30.4883% done, 9.9553 s remaining
subpixel-averaging is 48.7824% done, 4.42917 s remaining
subpixel-averaging is 67.0766% done, 2.08722 s remaining
subpixel-averaging is 85.3708% done, 0.729907 s remaining
time for set_epsilon = 53.5962 s
-----------
Meep progress: 54.233333333333334/408.0 = 13.3% done in 4.0s, 26.1s to go
on time step 3254 (time=54.2333), 0.00122939 s/step
Meep progress: 108.7/408.0 = 26.6% done in 8.0s, 22.0s to go
on time step 6523 (time=108.717), 0.0012238 s/step
Meep progress: 163.1/408.0 = 40.0% done in 12.0s, 18.0s to go
on time step 9787 (time=163.117), 0.00122558 s/step
Meep progress: 217.76666666666665/408.0 = 53.4% done in 16.0s, 14.0s to go
on time step 13068 (time=217.8), 0.00121939 s/step
Meep progress: 272.51666666666665/408.0 = 66.8% done in 20.0s, 9.9s to go
on time step 16354 (time=272.567), 0.0012174 s/step
Meep progress: 327.21666666666664/408.0 = 80.2% done in 24.0s, 5.9s to go
on time step 19637 (time=327.283), 0.00121851 s/step
Meep progress: 380.9166666666667/408.0 = 93.4% done in 28.0s, 2.0s to go
on time step 22860 (time=381), 0.00124152 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 62 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 46 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 62 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 42 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.77362
tran (twisted):, 1, 4.77, 0.07215
tran (twisted):, 2, 9.56, 0.00000
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 3.60012e-05 s
Working in 2D dimensions.
Computational cell is 7.4 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.107628 s
-----------
Meep progress: 78.06666666666666/208.0 = 37.5% done in 4.0s, 6.7s to go
on time step 4684 (time=78.0667), 0.000854101 s/step
Meep progress: 157.66666666666666/208.0 = 75.8% done in 8.0s, 2.6s to go
on time step 9461 (time=157.683), 0.000837429 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 6.69956e-05 s
Working in 2D dimensions.
Computational cell is 7.4 x 6.5 x 0 with resolution 30
     block, center = (-2.7,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (3.4,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 47.5314% done, 4.80929 s remaining
subpixel-averaging is 76.0521% done, 1.33452 s remaining
subpixel-averaging is 47.5314% done, 4.83031 s remaining
subpixel-averaging is 76.0521% done, 1.34121 s remaining
subpixel-averaging is 47.5314% done, 4.85467 s remaining
subpixel-averaging is 76.0521% done, 1.34843 s remaining
time for set_epsilon = 28.6162 s
-----------
Meep progress: 69.88333333333333/408.0 = 17.1% done in 4.0s, 19.4s to go
on time step 4193 (time=69.8833), 0.000954222 s/step
Meep progress: 141.2/408.0 = 34.6% done in 8.0s, 15.1s to go
on time step 8473 (time=141.217), 0.000934654 s/step
Meep progress: 212.45/408.0 = 52.1% done in 12.0s, 11.0s to go
on time step 12749 (time=212.483), 0.000935672 s/step
Meep progress: 283.9166666666667/408.0 = 69.6% done in 16.0s, 7.0s to go
on time step 17038 (time=283.967), 0.000932807 s/step
Meep progress: 355.21666666666664/408.0 = 87.1% done in 20.0s, 3.0s to go
on time step 21317 (time=355.283), 0.000934988 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 50 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 33 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (uniaxial):, 0, 0.00, 0.90581
tran (uniaxial):, 1, 4.77, 0.01950
tran (uniaxial):, 2, 9.56, 0.00000
tran (uniaxial):, 3, 14.43, 0.00000
tran (uniaxial):, 4, 19.41, 0.00000
-----------
Initializing structure...
time for choose_chunkdivision = 6.91414e-05 s
Working in 2D dimensions.
Computational cell is 10.8 x 6.5 x 0 with resolution 30
time for set_epsilon = 0.160437 s
-----------
Meep progress: 59.15/208.0 = 28.4% done in 4.0s, 10.1s to go
on time step 3549 (time=59.15), 0.00112713 s/step
Meep progress: 119.0/208.0 = 57.2% done in 8.0s, 6.0s to go
on time step 7141 (time=119.017), 0.00111385 s/step
Meep progress: 179.25/208.0 = 86.2% done in 12.0s, 1.9s to go
on time step 10757 (time=179.283), 0.00110643 s/step
run 0 finished at t = 208.0 (12480 timesteps)
-----------
Initializing structure...
time for choose_chunkdivision = 7.39098e-05 s
Working in 2D dimensions.
Computational cell is 10.8 x 6.5 x 0 with resolution 30
     block, center = (-4.4,0,0)
          size (2,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025)
     block, center = (0,0,0)
          size (6.8,1e+20,1e+20)
          axes (1,0,0), (0,1,0), (0,0,1)
subpixel-averaging is 29.0972% done, 10.5851 s remaining
subpixel-averaging is 46.5566% done, 4.83 s remaining
subpixel-averaging is 64.0161% done, 2.40549 s remaining
subpixel-averaging is 81.4755% done, 0.964889 s remaining
subpixel-averaging is 29.0972% done, 10.5591 s remaining
subpixel-averaging is 46.5566% done, 4.83987 s remaining
subpixel-averaging is 64.0161% done, 2.39451 s remaining
subpixel-averaging is 81.4755% done, 0.972159 s remaining
subpixel-averaging is 29.0972% done, 10.6613 s remaining
subpixel-averaging is 46.5566% done, 4.82302 s remaining
subpixel-averaging is 64.0161% done, 2.39088 s remaining
subpixel-averaging is 81.4755% done, 0.967715 s remaining
time for set_epsilon = 56.7605 s
-----------
Meep progress: 51.55/408.0 = 12.6% done in 4.0s, 27.7s to go
on time step 3093 (time=51.55), 0.00129338 s/step
Meep progress: 103.8/408.0 = 25.4% done in 8.0s, 23.4s to go
on time step 6229 (time=103.817), 0.00127578 s/step
Meep progress: 157.21666666666667/408.0 = 38.5% done in 12.0s, 19.1s to go
on time step 9435 (time=157.25), 0.0012479 s/step
Meep progress: 210.46666666666667/408.0 = 51.6% done in 16.0s, 15.0s to go
on time step 12631 (time=210.517), 0.00125168 s/step
Meep progress: 263.68333333333334/408.0 = 64.6% done in 20.0s, 10.9s to go
on time step 15825 (time=263.75), 0.0012525 s/step
Meep progress: 316.7833333333333/408.0 = 77.6% done in 24.0s, 6.9s to go
on time step 19011 (time=316.85), 0.0012555 s/step
Meep progress: 370.1166666666667/408.0 = 90.7% done in 28.0s, 2.9s to go
on time step 22212 (time=370.2), 0.00124999 s/step
run 0 finished at t = 408.0 (24480 timesteps)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 43 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
MPB solved for omega_1(1.85185,0,0) = 1.85185 after 56 iters
Dominant planewave for band 1: (1.851852,-0.000000,0.000000)
MPB solved for omega_2(1.85185,0,0) = 1.85823 after 51 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters
Dominant planewave for band 2: (1.845450,-0.153846,0.000000)
MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters
MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters
MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters
Dominant planewave for band 3: (1.826111,-0.307692,0.000000)
MPB solved for omega_4(1.85185,0,0) = 1.9085 after 30 iters
MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters
Dominant planewave for band 4: (1.793415,-0.461538,0.000000)
MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters
MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters
Dominant planewave for band 5: (1.746613,-0.615385,0.000000)
tran (twisted):, 0, 0.00, 0.49695
tran (twisted):, 1, 4.77, 0.22184
tran (twisted):, 2, 9.56, 0.00001
tran (twisted):, 3, 14.43, 0.00000
tran (twisted):, 4, 19.41, 0.00000

The diffraction spectra is plotted using the script below and shown in the accompanying figure.

In [3]:
cos_angles = [math.cos(math.radians(t)) for t in ang_uniaxial]
tran = m0_uniaxial+2*m1_uniaxial
eff_m0 = m0_uniaxial/tran
eff_m1 = (2*m1_uniaxial/tran)/cos_angles

phase = delta_n*dd/wvl
eff_m0_analytic = [math.cos(math.pi*p)**2 for p in phase]
eff_m1_analytic = [math.sin(math.pi*p)**2 for p in phase]
In [4]:
plt.figure(dpi=150)
plt.subplot(1,2,1)
plt.plot(phase,eff_m0,'bo-',clip_on=False,label='0th order (meep)')
plt.plot(phase,eff_m0_analytic,'b--',clip_on=False,label='0th order (analytic)')
plt.plot(phase,eff_m1,'ro-',clip_on=False,label='±1 orders (meep)')
plt.plot(phase,eff_m1_analytic,'r--',clip_on=False,label='±1 orders (analytic)')
plt.axis([0, 1.0, 0, 1])
plt.xticks([t for t in np.arange(0,1.2,0.2)])
plt.xlabel("phase delay Δnd/λ")
plt.ylabel("diffraction efficiency @ λ = 0.54 μm")
plt.legend(loc='center')
plt.title("homogeneous uniaxial grating")

cos_angles = [math.cos(math.radians(t)) for t in ang_twisted]
tran = m0_twisted+2*m1_twisted
eff_m0 = m0_twisted/tran
eff_m1 = (2*m1_twisted/tran)/cos_angles

plt.subplot(1,2,2)
plt.plot(phase,eff_m0,'bo-',clip_on=False,label='0th order (meep)')
plt.plot(phase,eff_m1,'ro-',clip_on=False,label='±1 orders (meep)')
plt.axis([0, 1.0, 0, 1])
plt.xticks([t for t in np.arange(0,1.2,0.2)])
plt.xlabel("phase delay Δnd/λ")
plt.ylabel("diffraction efficiency @ λ = 0.54 μm")
plt.legend(loc='center')
plt.title("bilayer twisted-nematic grating")

plt.tight_layout()
plt.show()

The left figure shows good agreement between the simulation results and analytic theory for the homogeneous uniaxial grating. Approximately 6% of the power in the input planewave is lost due to reflection from the grating. This value is an average over all phase delays. The total transmittance is therefore around 94%. The twisted-nematic grating, with results shown in the right figure, produces ±1 diffraction orders with nearly-constant peak transmittance over a broader bandwidth around Δnd/λ=0.5 than the homogeneous uniaxial polarization grating. This is consistent with results from the reference. The average reflectance and transmittance for the twisted-nematic grating are similar to those for the homogeneous uniaxial grating.

Finally, we demonstrate that when Δnd/λ=0.5 a circular-polarized planewave input produces just a single ±1 diffraction order. To specify a Ez+iEy circular-polarized planewave requires setting the amplitude of the Ey source to an imaginary number (from its default of 1):

sources = [mp.Source(mp.GaussianSource(fcen,fwidth=0.05*fcen), component=mp.Ez, center=src_pt, size=mp.Vector3(0,sy,0)),
           mp.Source(mp.GaussianSource(fcen,fwidth=0.05*fcen), component=mp.Ey, center=src_pt, size=mp.Vector3(0,sy,0), amplitude=1j)]

Note that even though the Jy current amplitude is complex in this example, only its real part is used and the resulting fields are therefore still real (the default).

The figure below shows a snapshot of Ez within the cell for four different cases: phase delays (Δnd/λ) of 0.5 and 1.0, and planewave circular polarization of Ez+iEy and Ez-iEy. The empty regions on the cell sides are PMLs. The thin solid black line denotes the boundary between the grating (on the left) and air. As expected, for Δnd/λ=0.5 there is just a single ±1 diffraction order which depends on the chirality of the input planewave (this is not the case for a linear-polarized planewave). The angle of this diffracted order (±4.8°) agrees with the analytic result. Snapshots of Ey are similar.