RCTGLPBC used to simulate infinite slab
  infinite slab , m=(1.33,0.01), thickness = 0.32um, wave=0.50um
  using 16x1x1 TUC (M,N)=(0,0)

a = 0.32um = slab thickness
b = X/16 = 0.02um
c = X/16 = 0.02um
V = 0.32*0.02*0.02 um^3
aeff = (3*V/4*pi)^{1/3} = 0.031264 um


s1 = 16	= X/d
s2 = 1	= Y/d
s3 = 1  = Z/d
s4 = 1  = P_y/d
s5 = 1  = P_z/d

wave = 0.5um
allowed (M,N)
M_max = int(L_y/wave)=int(0.02/0.5)=0

m=1.33+0.01i
results with gamma=1.e-2
              1x1      2x2        4x4       8x8      16x16     32x32     40x40
S11(0,0) = 8.7691e-1 8.7683e-1 8.7683e-1 8.7682e-1 8.7682e-1 8.7684e-1 8.7681e-1
	 = 4.8522e-2 4.8522e-1 4.8518e-2 4.8519e-2 4.8519e-2 4.8519e-2 4.8519e-2
S11(1,0) =                                                    < 1e-11   < 1e-11
                                                              < 1e-11   < 1e-11
S11(0,1)						      < 1e-11   < 1e-11
                                                              < 1e-11   < 1e-11
S11(1,1)                                                                < 1e-11
                                                                        < 1e-11
R+T =    = 0.92543   0.92535   0.92535
1-R+T	 = 0.07457   0.07465   0.07465
calc A   = 0.07465   0.07465   0.07465   0.074654  0.074654  0.074654  0.74655

=========================================================================
2x2
a = 0.32um
b = X/8 = 0.04
c = X/8 = 0.04
V = 0.32*0.04*0.04 um^3
aeff = (3*V/4*pi)^{1/3} = 0.049628

M_max = int(L_y/wave)=int(0.04/0.5)=0

=========================================================================
4x4
a = 0.32um
b = X/4 = 0.08
c = X/4 = 0.08
V = 0.32*0.08*0.08 um^3
aeff = (3*V/4*pi)^{1/3} = 0.078780

M_max = int(L_y/wave)=int(0.08/0.5)=0
=========================================================================
8x8
a = 0.32um
b = X/2 = 0.16
c = X/2 = 0.16
V = 0.32*0.16*0.16 um^3
aeff = (3*V/4*pi)^{1/3} = 0.125055

M_max = int(L_y/wave)=int(0.16/0.5)=0
=========================================================================
16x16  (4096)
a = 0.32um
b = X = 0.32
c = X = 0.32
V = 0.32*0.32*0.32 um^3
aeff = (3*V/4*pi)^{1/3} = 0.198512

M_max = int(L_y/wave)=int(0.32/0.5)=0
=========================================================================
32x32  (16384)
a = 0.32um
b = 2*X = 0.64
c = 2*X = 0.64
V = 0.32*0.64*0.64 um^3
aeff = (3*V/4*pi)^{1/3} = 0.315118

M_max = int(L_y/wave)=int(0.64/0.5)=1
=========================================================================
40x40
a = 0.32um
b = 0.80
c = 0.80
V = 0.32*0.80*0.80 um^3
aeff = (3*V/4*pi)^{1/3} = 0.365663

M_max = int(L_y/wave)=int(.80/0.5)=1
N_max(M=1)=int(sqrt(1.6**2-1)=int(1.249)
thus have (0,0),(1,0),(1,1) allowed
=========================================================================
