summaryrefslogtreecommitdiff
path: root/example_beam_trace_diamond.m
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2013-01-16 14:27:10 -0500
committerEugeniy Mikhailov <evgmik@gmail.com>2013-01-16 14:27:10 -0500
commit37d09478a8fca4e1229dcb1246c9f23f1f6118df (patch)
tree79e2d5f9179987b46b5b743c552c29f94ffb999e /example_beam_trace_diamond.m
parent28ef703f718ce7f540864ca8f9a301f8d7d3f3e0 (diff)
downloadwgmr-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.m76
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');
-