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
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
import math
import unittest
import numpy as np
import meep as mp
class TestGetPoint(unittest.TestCase):
def test_get_point(self):
sxy = 6 # cell size
dpml = 1 # thickness of PML
def sinusoid(p):
r = (p.x**2 + p.y**2) ** 0.5
return mp.Medium(index=1.0 + math.sin(2 * math.pi * r) ** 2)
geometry = [
mp.Block(center=mp.Vector3(), size=mp.Vector3(sxy, sxy), material=sinusoid)
]
src = [
mp.Source(
mp.GaussianSource(1.0, fwidth=0.1), component=mp.Ez, center=mp.Vector3()
)
]
sim = mp.Simulation(
cell_size=mp.Vector3(sxy, sxy),
geometry=geometry,
sources=src,
k_point=mp.Vector3(),
resolution=20,
symmetries=[mp.Mirror(mp.X), mp.Mirror(mp.Y)],
boundary_layers=[mp.PML(dpml)],
)
sim.run(until_after_sources=100)
## reference values for Ez and epsilon from serial run
ez_ref = [
-0.0002065983,
-0.0001954795,
-0.0000453570,
0.0000311267,
-0.0000121473,
-0.0000410032,
-0.0000341301,
-0.0000275021,
-0.0000397990,
-0.0000351730,
0.0000079602,
0.0000227437,
-0.0001092821,
-0.0002202751,
-0.0001408186,
0.0006325076,
0.0024890489,
0.0027476069,
0.0014815873,
0.0004714913,
-0.0004332029,
-0.0007101315,
-0.0003818581,
-0.0000748507,
0.0001408819,
0.0001119776,
0.0000395008,
0.0000078844,
-0.0000010431,
]
eps_ref = [
1.6458346134,
1.2752837068,
1.0974010956,
1.0398089537,
1.0465784716,
1.0779924737,
1.1059439286,
1.1135579291,
1.0971979186,
1.0653178566,
1.0391657283,
1.0513779677,
1.1466009312,
1.3882154483,
1.8496939317,
2.5617731415,
3.3788212533,
3.9019494270,
3.6743431894,
2.7285622651,
1.6635165033,
1.0891237010,
1.1485969863,
1.9498398061,
3.3100416367,
3.9038800599,
2.8471862395,
1.4742605488,
1.0370162714,
]
x = np.linspace(-0.865692, 2.692867, 29)
places = 5 if mp.is_single_precision() else 10
for j in range(x.size):
self.assertAlmostEqual(
np.real(sim.get_field_point(mp.Ez, mp.Vector3(x[j], -0.394862))),
ez_ref[j],
places=places,
)
self.assertAlmostEqual(
sim.get_epsilon_point(mp.Vector3(x[j], 2.967158)),
eps_ref[j],
places=places,
)
if __name__ == "__main__":
unittest.main()
|
