diff options
author | Eugeniy Mikhailov <evgmik@gmail.com> | 2013-01-16 14:27:10 -0500 |
---|---|---|
committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2013-01-16 14:27:10 -0500 |
commit | 37d09478a8fca4e1229dcb1246c9f23f1f6118df (patch) | |
tree | 79e2d5f9179987b46b5b743c552c29f94ffb999e /example_beam_trace_diamond.m | |
parent | 28ef703f718ce7f540864ca8f9a301f8d7d3f3e0 (diff) | |
download | wgmr-37d09478a8fca4e1229dcb1246c9f23f1f6118df.tar.gz wgmr-37d09478a8fca4e1229dcb1246c9f23f1f6118df.zip |
move code to directory with speaking name
Diffstat (limited to 'example_beam_trace_diamond.m')
-rw-r--r-- | example_beam_trace_diamond.m | 76 |
1 files changed, 0 insertions, 76 deletions
diff --git a/example_beam_trace_diamond.m b/example_beam_trace_diamond.m deleted file mode 100644 index e3cc30d..0000000 --- a/example_beam_trace_diamond.m +++ /dev/null @@ -1,76 +0,0 @@ -%beam1.k=[-1,0]; -%beam1.origin=[2,0.5]; -beam_initial_travel_angle=(70+30-180)/180*pi; -beam1.k=[cos(beam_initial_travel_angle), sin(beam_initial_travel_angle)]; -beam1.origin=[0.5,2]; -%beam1.k=[-1,-1]; -%beam1.origin=[2,2]; -beam1.face=NA; -beam1.intensity = 1; -beam1.polarization = 1; % 1 for s and 2 for p -beam1.status='incoming'; % could be reflected, refracted, incoming - -beam2=beam1; -beam1.polarization = 2; % 1 for s and 2 for p - -beams={beam1,beam2}; - -n_rutile=[2.7532, 2.4885]; % n_s and n_p -n_diamond=[2.4 , 2.4 ]; % n_s and n_p -n_disk= [1.3751, 1.3868]; % MgF2 book value -n_disk= [1.465 , 1.465 ]; % MgF2 measured -n_disk= [1.39 , 1.39 ]; % CaF2 measured for p-polarization -n_air= [1., 1.]; - -face1.vertex1=[-1,0]; -face1.vertex2=[1,0]; -%face1.n_right = n_air; -face1.n_right = n_disk; -face1.n_left= n_diamond; - -face2.vertex1=[1,0]; -face2.vertex2=[0,sqrt(3)]; -face2.n_right = [1, 1]; -face2.n_left= n_diamond; - -face3.vertex1=[0,sqrt(3)]; -face3.vertex2=[-1,0]; -face3.n_right = [1, 1]; -face3.n_left= n_diamond; - -faces={face1,face2,face3}; - - -border1.vertex1=[-2,-2]; -border1.vertex2=[2,-2]; -border1.n_right = [1, 1]; -border1.n_left= [1, 1]; - -border2.vertex1=[2,-2]; -border2.vertex2=[2,2]; -border2.n_right = [1, 1]; -border2.n_left= [1, 1]; - -border3.vertex1=[2,2]; -border3.vertex2=[-2,2]; -border3.n_right = [1, 1]; -border3.n_left= [1, 1]; - -border4.vertex1=[-2,2]; -border4.vertex2=[-2,-2]; -border4.n_right = [1, 1]; -border4.n_left= [1, 1]; - -borders={border1,border2, border3, border4}; - - -border_limits=[-2,-2, 2,2]; -figure(1); -clf; -[processed_beams] = beam_trace(beams, faces, borders, border_limits ); -text(-1,-1, 'red: s-polarization; blue: p-polarization'); - - -print('beams.eps','-depsc2'); -print('beams.png'); - |