-
Notifications
You must be signed in to change notification settings - Fork 0
/
verification.log
88 lines (83 loc) · 5 KB
/
verification.log
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
./verification.out ipw.txt multilayer.txt
-- incident plane wave --
wavelength in vacuum : 0.6328
refractive index : 1.515
p-polarization amplitude : 1.000000e-03 +0.000000e+00 I
s-polarization amplitude : 0.000000e+00 +0.000000e+00 I
incident angle (theta) [rad] : 0.75049157835756
parameter for incident direction (phi) [rad] : 0
-- multilayer --
number of layers : 3
layer 0, refractive index n_00 = 1.515000 +0.000000 I
z_01 ------------------------------------------------ : 0
layer 1, refractive index n_01 = 0.066600 +4.045000 I
z_02 ------------------------------------------------ : 0.054
layer 2, refractive index n_02 = 1.000000 +0.000000 I
-- verification using boundary conditions --
boundary 1, z=0.00000000000000
electromagnetic field at r=(7.8, 5.4, 0), layer id=0
E_x=-1.22710656443026e-04 -7.91224193484585e-04 I
E_y= 0.00000000000000e+00 -0.00000000000000e+00 I
E_z= 7.64096515291166e-06 +6.20694909070772e-04 I
H_x= 0.00000000000000e+00 -0.00000000000000e+00 I
H_y=-1.69737390652968e-05 -1.37881971909147e-03 I
H_z= 0.00000000000000e+00 -0.00000000000000e+00 I
electromagnetic field at r=(7.8, 5.4, 0), layer id=1
E_x=-1.22710656443026e-04 -7.91224193484585e-04 I
E_y= 0.00000000000000e+00 -0.00000000000000e+00 I
E_z= 1.79462255986847e-06 -8.70340673377630e-05 I
H_x= 0.00000000000000e+00 -0.00000000000000e+00 I
H_y=-1.69737390652971e-05 -1.37881971909148e-03 I
H_z= 0.00000000000000e+00 -0.00000000000000e+00 I
boundary condition : n dot D_0 = n dot D_1, n dot B_0 = n dot B_1, n=(0,0,1)
n dot D_0 = 0.00001753773424 +0.00142463447267 I
n dot D_1 = 0.00001753773424 +0.00142463447267 I
n dot B_0 = 0.00000000000000 -0.00000000000000 I
n dot B_1 = 0.00000000000000 -0.00000000000000 I
boundary condition : n times E_0 = n times E_1, n times H_0 = n times H_1
n times E_0 = (-0.00000000000000 +0.00000000000000 I, -0.00012271065644 -0.00079122419348 I, 0)
n times E_1 = (-0.00000000000000 +0.00000000000000 I, -0.00012271065644 -0.00079122419348 I, 0)
n times H_0 = ( 0.00001697373907 +0.00137881971909 I, 0.00000000000000 -0.00000000000000 I, 0)
n times H_1 = ( 0.00001697373907 +0.00137881971909 I, 0.00000000000000 -0.00000000000000 I, 0)
boundary 2, z=0.05400000000000
electromagnetic field at r=(7.8, 5.4, 0.054), layer id=1
E_x= 9.21683161083001e-04 -3.77991423084813e-03 I
E_y= 0.00000000000000e+00 -0.00000000000000e+00 I
E_z=-9.10214145972415e-04 -2.53964749868545e-04 I
H_x= 0.00000000000000e+00 -0.00000000000000e+00 I
H_y=-1.45425317855029e-02 -3.54600814929202e-03 I
H_z= 0.00000000000000e+00 -0.00000000000000e+00 I
electromagnetic field at r=(7.8, 5.4, 0.054), layer id=2
E_x= 9.21683161083000e-04 -3.77991423084813e-03 I
E_y= 0.00000000000000e+00 -0.00000000000000e+00 I
E_z= 1.50257439857377e-02 +3.66383319001689e-03 I
H_x= 0.00000000000000e+00 -0.00000000000000e+00 I
H_y=-1.45425317855029e-02 -3.54600814929202e-03 I
H_z= 0.00000000000000e+00 -0.00000000000000e+00 I
boundary condition : n dot D_0 = n dot D_1, n dot B_0 = n dot B_1, n=(0,0,1)
n dot D_0 = 0.01502574398574 +0.00366383319002 I
n dot D_1 = 0.01502574398574 +0.00366383319002 I
n dot B_0 = 0.00000000000000 -0.00000000000000 I
n dot B_1 = 0.00000000000000 -0.00000000000000 I
boundary condition : n times E_0 = n times E_1, n times H_0 = n times H_1
n times E_0 = (-0.00000000000000 +0.00000000000000 I, 0.00092168316108 -0.00377991423085 I, 0)
n times E_1 = (-0.00000000000000 +0.00000000000000 I, 0.00092168316108 -0.00377991423085 I, 0)
n times H_0 = ( 0.01454253178550 +0.00354600814929 I, 0.00000000000000 -0.00000000000000 I, 0)
n times H_1 = ( 0.01454253178550 +0.00354600814929 I, 0.00000000000000 -0.00000000000000 I, 0)
-- verification using Maxwell's equations --
electromagnetic field at r=(5.2, -9.6, 3.6), layer id=2
E_x=-3.97786629661388e-07 -1.09680281017993e-07 I
E_y=-0.00000000000000e+00 +0.00000000000000e+00 I
E_z= 4.35996089384902e-07 -1.58126340790022e-06 I
H_x=-0.00000000000000e+00 +0.00000000000000e+00 I
H_y=-4.21974911475483e-07 +1.53041163169484e-06 I
H_z=-0.00000000000000e+00 +0.00000000000000e+00 I
Faraday equation : rot E = -dB/dt
rot E = ( 0.00000000000000 +0.00000000000000 I, -0.00001519573306 -0.00000418986499 I, -0.00000000000000 +0.00000000000000 I)
-dB/dt = (-0.00000000000000 +0.00000000000000 I, -0.00001519573306 -0.00000418986499 I, -0.00000000000000 +0.00000000000000 I)
Ampere equation : rot H = dD/dt
rot H = (-0.00000108903529 +0.00000394969518 I, 0.00000000000000 +0.00000000000000 I, -0.00001570064951 -0.00000432908379 I)
dD/dt = (-0.00000108903529 +0.00000394969518 I, 0.00000000000000 +0.00000000000000 I, -0.00001570064951 -0.00000432908379 I)
Gauss's law in this system : div D = 0, div B = 0
div D = -0.00000000000000 +0.00000000000000 I
div B = -0.00000000000000 +0.00000000000000 I