εzz

Author: Jon Schlipf

src/Common/Common.jl

@@ -8,7 +8,7 @@ include("models.jl")
 include("grids.jl")
 export Eigenmodes,Halfspace
 export eigenmodes,halfspace
-export a2e2d,a2p,e2p,slicehalf,getfields
+export a2e2d,a2p,e2p,slicehalf,getfields,get_εzz
 """
     Eigenmodes(V,W,X,q)
 
@@ -175,6 +175,36 @@ function eigenmodes(g::RCWAGrid,λ,l::Array{Layer,1})
     end
     return rt
 end
+function get_εzz(dnx,dny,λ,l::PatternedLayer)
+    #get the base permittivity
+    εzz=get_permittivity(l.materials[1],λ,1)*I
+    #add the permittivity for all inclusions
+    if minimum([typeof(m)<:Isotropic for m in l.materials])
+        #all isotropic
+        εzz=get_permittivity(l.materials[1],λ)*I
+        for ct=1:length(l.geometries)
+            rec=reciprocal(l.geometries[ct],dnx,dny)
+            εzz+=rec*(get_permittivity(l.materials[ct+1],λ)-get_permittivity(l.materials[ct],λ))
+        end
+    else
+        #anisotropic
+        εzz=get_permittivity(l.materials[1],λ,5)*I
+        for ct=1:length(l.geometries)
+            rec=reciprocal(l.geometries[ct],dnx,dny)
+            εzz+=rec*(get_permittivity(l.materials[ct+1],λ,5)-get_permittivity(l.materials[ct],λ,5))
+        end
+    end	 	
+    return εzz
+end
+
+function get_εzz(dnx,dny,λ,l:SimpleLayer)
+    return get_permittivity(l.material,λ)*I
+end
+
+function get_εzz(dnx,dny,λ,l:AnisotropicLayer)
+    return get_permittivity(l.material,λ,5)*I
+end
+
 """
     halfspace(Kx,Ky,material,λ)
 
@@ -305,7 +335,7 @@ computes the electric and magnetic fields within a layer
 * `efield` : 4D tensor for the electric field (dimensions are x, y, z, and the component (E_x or E_y or E_z)
 * `hfield` : 4D tensor for the magnetic field (dimensions are x, y, z, and the component (E_x or E_y or E_z)
 """
-function getfields(ain,bout,em::Eigenmodes,grd::RCWAGrid,xypoints,zpoints,λ,window="Hann",padding=[0,0])
+function getfields(ain,bout,em::Eigenmodes,grd::RCWAGrid,xypoints,zpoints,λ,εzz,window="Hann",padding=[0,0])
     Nx=Int(floor(xypoints[1]/2))
     Ny=Int(floor(xypoints[2]/2))
     nx=[r for r in -Nx:Nx-1, c in -Ny:Ny-1]
@@ -327,7 +357,7 @@ function getfields(ain,bout,em::Eigenmodes,grd::RCWAGrid,xypoints,zpoints,λ,win
         #convert amplitude vectors to electric fields
         ex,ey=a2e2d(a+b,em.W)
         hx,hy=a2e2d(a-b,em.V)
-        ez=-1im*(grd.Kx*hy-grd.Ky*hx)
+        ez=-1im*εzz\(grd.Kx*hy-grd.Ky*hx)
         hz=-1im*(grd.Kx*ey-grd.Ky*ex)
         #convert from reciprocal lattice vectors to real space distribution
         efield[:,:,zind,1]=recipvec2real(Array(grd.nx),Array(grd.ny),Array(ex),nx,ny,windowfunction)

src/ETM/ETM.jl

@@ -188,7 +188,7 @@ function etm_getfields_stack(mdl::RCWAModel,grd::RCWAGrid,xypoints,zpoints,λ,a,
     H=zeros(xypoints[1],xypoints[2],length(zpoints),3)*1im
     for i in eachindex(mdl.layers)
         zpoints_separated=zpoints[zstack.<=zpoints.<zstack+mdl.layers[i].thickness]
-        Ei,Hi=getfields(a[i],b[i],em[i],grd,xypoints,zpoints_separated.-zstack,λ,window,padding)
+        Ei,Hi=getfields(a[i],b[i],em[i],grd,xypoints,zpoints_separated.-zstack,λ,get_εzz(grd.dnx,grd.dny,λ,mdl.layers[i]),window,padding)
         E[:,:,zind:zind+length(zpoints_separated)-1,:]=Ei
         H[:,:,zind:zind+length(zpoints_separated)-1,:]=Hi
         zstack=zstack+mdl.layers[i].thickness