-
Notifications
You must be signed in to change notification settings - Fork 587
/
mpb_sq_rods.py
49 lines (36 loc) · 1.13 KB
/
mpb_sq_rods.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
from __future__ import division
import time
import meep as mp
from meep import mpb
# Compute band structure for a square lattice of dielectric rods
# in air.
# Define various parameters with define_param so that they are
# settable from the command_line (with mpb <param>=<value>):
r = 0.2 # radius of the rods
eps = 11.56 # dielectric constant
k_interp = 4 # number of k points to interpolate
GaAs = mp.Medium(epsilon=eps)
geometry_lattice = mp.Lattice(size=mp.Vector3(1, 1)) # 2d cell
geometry = [mp.Cylinder(r, material=GaAs)]
Gamma = mp.Vector3()
X = mp.Vector3(0.5, 0)
M = mp.Vector3(0.5, 0.5)
k_points = mp.interpolate(k_interp, [Gamma, X, M, Gamma])
resolution = 32
num_bands = 8
ms = mpb.ModeSolver(
geometry_lattice=geometry_lattice,
geometry=geometry,
k_points=k_points,
resolution=resolution,
num_bands=num_bands
)
def main():
# Compute the TE and TM bands and report the total elapsed time:
t0 = time.time()
ms.run_te()
ms.run_tm()
print("total time for both TE and TM bands: {:.2f} seconds".format(time.time() - t0))
ms.display_eigensolver_stats()
if __name__ == '__main__':
main()