From 94eb84b73e3a2e90000da5fff2ae83bcc09411c4 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 10:12:23 -0500 Subject: Added TODO for xmds --- xmds2/TODO | 8 ++++++++ 1 file changed, 8 insertions(+) create mode 100644 xmds2/TODO diff --git a/xmds2/TODO b/xmds2/TODO new file mode 100644 index 0000000..daee897 --- /dev/null +++ b/xmds2/TODO @@ -0,0 +1,8 @@ +Fix Ndens argument treatment: Right now code is oblivius to --Ndens +parameter, since eta ~1/Ndens sits in global parameter section and already +defined. So it is not recalculated after the parameter assignment. + +Ugly hack would be to make variable eta_times_Ndens instead of eta, and +recalculate eta=eta_times_Ndens/Ndens in the ODE part around where it +needed. It will have performance penalty but probably not too large. + -- cgit v1.2.3 From bc0bf295b05c838a6b016d506c9bce279265e504 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 10:28:55 -0500 Subject: plot and postprocessing files are symlinked --- .../plot_fields_propagation_I2.gp | 16 +--- .../plot_fields_propagation_I4.gp | 1 + .../fast_to_slow_switch/pp_I2.m | 92 +--------------------- .../plot_fields_propagation_I2.gp | 15 ++++ .../plot_fields_propagation_I4.gp | 15 ++++ .../pp_I2.m | 63 ++++++++++----- 6 files changed, 75 insertions(+), 127 deletions(-) mode change 100644 => 120000 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp create mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I4.gp mode change 100644 => 120000 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I2.gp create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I4.gp diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp deleted file mode 100644 index 77b4870..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp +++ /dev/null @@ -1,15 +0,0 @@ -set terminal postscript portrait enhanced color solid size 5,3.5 -set output 'fields_propagation_I2.eps' -set dgrid3d 100,100 qnorm 4 -set pm3d map -#set contour -set hidden3d -set palette rgb 10,13,31 negative - - -set xlabel "z (cm)" -set ylabel "t ({/Symbol m}S)" -set zlabel "I_2 (1/S)" -set nokey -#set view map -splot [0:][-0.2:0.2] 'I2.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp new file mode 120000 index 0000000..e11c569 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp @@ -0,0 +1 @@ +../plot_fields_propagation_I2.gp \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I4.gp new file mode 120000 index 0000000..134a380 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I4.gp @@ -0,0 +1 @@ +../plot_fields_propagation_I4.gp \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m deleted file mode 100644 index c4b6333..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m +++ /dev/null @@ -1,91 +0,0 @@ -Nlevels_no_dopler_with_z_4wm; - -%% field I2 propagation -z_1=z_1*100; % z in cm -t_1=t_1*1e6; % time now measured in uS -figure(1) -%set(gca,'fontsize',20); -imagesc(z_1, t_1, I2_out_1); colorbar -tmin=-0.4; -tmax= 0.4; -ylim([tmin,tmax],'manual'); -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_2') -title('I_2') -print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps') - -desired_x_size=200; -desired_y_size=200; -xskip=ceil(length(z_1)/desired_x_size); -yskip=ceil(length(t_1)/desired_y_size); -map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); - - -%% field I4 propagation -figure(2) - -imagesc(z_1, t_1, I4_out_1); colorbar -tmin=-0.4; -tmax= 0.4; -ylim([tmin,tmax],'manual'); -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_2') -title('I_2') -print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I4.eps') - - - - - -%% fields before and after the cell -figure(3) -%set(gca,'fontsize',30); -plot( ... - t_1,I2_out_1(:,1),'.-;I_2 before;', "linewidth", 4 ... - ,t_1,I2_out_1(:,end), '-;I_2 after;', "linewidth", 4 ... - ,t_1,I4_out_1(:,end), '-;I_4 after;', "linewidth", 4 ... - ) -xlabel('t (uS)') -ylabel('I (1/s)^2') -title('Fields before and after cell') -legend('location', 'northeast'); - -[b, a]=butter(3, 0.05); -I2_out_after=I2_out_1(:,end); -I2_out_after_filtered=filtfilt(b,a,I2_out_after); -settling_time=0.01; %uS -t_good_indx=t_1> min(t_1 + settling_time); -[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); -delay_time=t_1(max_pos_after)-t_1(max_pos_before); -printf('Second field delay time = %f uS\n',delay_time); - -%set(gca,'fontsize',40); -%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell.eps') - -figure(4) -I2_max_in=max(I2_out_1(t_good_indx,1)); -I2_max_out=max(I2_out_1(t_good_indx,end)); -I4_max_out=max(I4_out_1(t_good_indx,end)); -I2_in_norm=(I2_out_1(:,1))/I2_max_in; -I2_out_norm=(I2_out_1(:,end))/I2_max_out; -I4_out_norm=(I4_out_1(:,end))/I4_max_out; -tmin=-.05; -tmax=.05; -indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region -plot( ... - t_1(indx),I2_in_norm(indx),'.-;I_2 before;', "linewidth", 4, ... - t_1(indx),I2_out_norm(indx), '-;I_2 after;', "linewidth", 4 ... - ,t_1(indx),I4_out_norm(indx), '-;I_4 after;', "linewidth", 4 ... - ) -legend('location', 'southeast'); -xlim([tmin,tmax],'manual'); -ylim([0.5,1],'manual'); -xlabel('t (uS)') -ylabel('I normalized') -title('Fields before and after cell normalized') -%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_normalized.eps') - diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m new file mode 120000 index 0000000..4a1ed9b --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m @@ -0,0 +1 @@ +../pp_I2.m \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I2.gp new file mode 100644 index 0000000..77b4870 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I2.gp @@ -0,0 +1,15 @@ +set terminal postscript portrait enhanced color solid size 5,3.5 +set output 'fields_propagation_I2.eps' +set dgrid3d 100,100 qnorm 4 +set pm3d map +#set contour +set hidden3d +set palette rgb 10,13,31 negative + + +set xlabel "z (cm)" +set ylabel "t ({/Symbol m}S)" +set zlabel "I_2 (1/S)" +set nokey +#set view map +splot [0:][-0.2:0.2] 'I2.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I4.gp new file mode 100644 index 0000000..3849e2f --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I4.gp @@ -0,0 +1,15 @@ +set terminal postscript portrait enhanced color solid size 5,3.5 +set output 'fields_propagation_I4.eps' +set dgrid3d 100,100 qnorm 4 +set pm3d map +#set contour +set hidden3d +set palette rgb 10,13,31 negative + + +set xlabel "z (cm)" +set ylabel "t ({/Symbol m}S)" +set zlabel "I_2 (1/S)" +set nokey +#set view map +splot [0:][-0.2:0.2] 'I4.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m index 45d4913..f632d0a 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m @@ -1,6 +1,6 @@ -Nlevels_no_dopler_with_z_4wm +Nlevels_with_doppler_with_z_4wm -%% field propagation +%% field I2 propagation z_1=z_1*100; % z in cm t_1=t_1*1e6; % time now measured in uS figure(1) @@ -13,33 +13,50 @@ xlabel('z (cm)') ylabel('t (uS)') zlabel('I_2') title('I_2') +print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps') -xskip=1; -yskip=10; -%map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); +desired_x_size=200; +desired_y_size=200; +xskip=ceil(length(z_1)/desired_x_size); +yskip=ceil(length(t_1)/desired_y_size); +map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); +%% field I4 propagation +figure(2) + +imagesc(z_1, t_1, I4_out_1); colorbar +tmin=-0.4; +tmax= 0.4; +ylim([tmin,tmax],'manual'); +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_2') +title('I_2') +print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I4.eps') +map2dat('I4.dat',z_1,t_1, I4_out_1, xskip, yskip); + -print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps') %% fields before and after the cell -figure(2) +figure(3) %set(gca,'fontsize',30); plot( ... - t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4, ... - t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ... + t_1,I2_out_1(:,1),'.-;I_2 before;', "linewidth", 4 ... + ,t_1,I2_out_1(:,end), '-;I_2 after;', "linewidth", 4 ... + ,t_1,I4_out_1(:,end), '-;I_4 after;', "linewidth", 4 ... ) xlabel('t (uS)') -ylabel('I_2 (1/s)^2') -title('I_2 before and after cell') -legend('location', 'southwest'); +ylabel('I (1/s)^2') +title('Fields before and after cell') +legend('location', 'northeast'); [b, a]=butter(3, 0.05); I2_out_after=I2_out_1(:,end); I2_out_after_filtered=filtfilt(b,a,I2_out_after); -settling_time=0.8; %uS +settling_time=0.01; %uS t_good_indx=t_1> min(t_1 + settling_time); [m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); delay_time=t_1(max_pos_after)-t_1(max_pos_before); @@ -47,25 +64,29 @@ printf('Second field delay time = %f uS\n',delay_time); %set(gca,'fontsize',40); %set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_I2.eps') +print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell.eps') figure(4) I2_max_in=max(I2_out_1(t_good_indx,1)); I2_max_out=max(I2_out_1(t_good_indx,end)); +I4_max_out=max(I4_out_1(t_good_indx,end)); I2_in_norm=(I2_out_1(:,1))/I2_max_in; I2_out_norm=(I2_out_1(:,end))/I2_max_out; -tmin=-0.05; -tmax=0.05; +I4_out_norm=(I4_out_1(:,end))/I4_max_out; +tmin=-.05; +tmax=.05; indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region plot( ... - t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4, ... - t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ... + t_1(indx),I2_in_norm(indx),'.-;I_2 before;', "linewidth", 4, ... + t_1(indx),I2_out_norm(indx), '-;I_2 after;', "linewidth", 4 ... + ,t_1(indx),I4_out_norm(indx), '-;I_4 after;', "linewidth", 4 ... ) legend('location', 'southeast'); xlim([tmin,tmax],'manual'); +ylim([0.5,1],'manual'); xlabel('t (uS)') -ylabel('I_2') -title('I_2 before and after cell normalized') +ylabel('I normalized') +title('Fields before and after cell normalized') %set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','-depsc2', '-tight','-S200,120', 'probe_before_after_cell_I2_normalized.eps') +print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_normalized.eps') -- cgit v1.2.3 From 5173b4afe998c3584e81ffd3f9bde24189fd8518 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 10:31:41 -0500 Subject: Added new parameters to makefile I.e. Temperature, cell length and density --- .../fast_to_slow_switch/Makefile | 9 ++++++++- .../Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/Makefile | 9 ++++++++- 2 files changed, 16 insertions(+), 2 deletions(-) diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile index 7581c04..d67ae2c 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile @@ -9,7 +9,14 @@ GNUPLOT_FILES = $(wildcard *.gp) XSIL2GRAPHICS = xsil2graphics # fast light -PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 +PARAMS = \ + --Ndens=1e15 \ + --Lcell=1.5e-2 \ + --Temperature=1e-9 \ + --delta1=0 --delta2=0 --delta3=0 \ + --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 + + # slow light EIT #PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0 diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/Makefile b/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/Makefile index 11c0320..12efac9 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/Makefile +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/Makefile @@ -10,7 +10,14 @@ XSIL2GRAPHICS = xsil2graphics # fast light #PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=1e1 -PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=0.1e7 --E2o=1e4 --E3o=.3e7 --E4o=0 +#PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=0.1e7 --E2o=1e4 --E3o=.3e7 --E4o=0 --Lcell=1.5e-2 --Temperature=.0001 +PARAMS = \ + --Ndens=1e15 \ + --Lcell=1.5e-2 \ + --Temperature=1e-9 \ + --delta1=0 --delta2=0 --delta3=0 \ + --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 + # slow light EIT #PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0 -- cgit v1.2.3 From ef0e5222bc255dcc0a4dceeaac8b526311759d09 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 10:33:47 -0500 Subject: TODO: more items --- xmds2/TODO | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/xmds2/TODO b/xmds2/TODO index daee897..c152093 100644 --- a/xmds2/TODO +++ b/xmds2/TODO @@ -1,4 +1,6 @@ -Fix Ndens argument treatment: Right now code is oblivius to --Ndens +- TODO Fix Ndens argument treatment: + +Right now code is oblivius to --Ndens parameter, since eta ~1/Ndens sits in global parameter section and already defined. So it is not recalculated after the parameter assignment. @@ -6,3 +8,7 @@ Ugly hack would be to make variable eta_times_Ndens instead of eta, and recalculate eta=eta_times_Ndens/Ndens in the ODE part around where it needed. It will have performance penalty but probably not too large. + +- TODO Fix Temperature argument treatment: + +Same problem as above -- cgit v1.2.3 From 77bd2b73b7fd57bf79507c9a7535b28a8b920d2f Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 10:44:18 -0500 Subject: bug fix: wrong data loading script --- xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m index f632d0a..ea75b38 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m @@ -1,4 +1,4 @@ -Nlevels_with_doppler_with_z_4wm +Nlevels_no_dopler_with_z_4wm; %% field I2 propagation z_1=z_1*100; % z in cm -- cgit v1.2.3 From 9a48d51435fccbe66579b4d6e465205cb7ef8e61 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 10:44:58 -0500 Subject: Readme and Makefile use new run time parameters --- .../fast_to_slow_switch/Makefile | 2 +- .../fast_to_slow_switch/README | 7 ++++++- 2 files changed, 7 insertions(+), 2 deletions(-) diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile index d67ae2c..9c65fd6 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile @@ -11,7 +11,7 @@ XSIL2GRAPHICS = xsil2graphics # fast light PARAMS = \ --Ndens=1e15 \ - --Lcell=1.5e-2 \ + --Lcell=10.0e-2 \ --Temperature=1e-9 \ --delta1=0 --delta2=0 --delta3=0 \ --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README index 70f7262..ec7d200 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README @@ -1,5 +1,10 @@ fast to slow transition -PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4 +PARAMS = \ + --Ndens=1e15 \ + --Lcell=10.0e-2 \ + --Temperature=1e-9 \ + --delta1=0 --delta2=0 --delta3=0 \ + --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 1.5 cm seems to give largest advance 2.5 cm gives zero delay -- cgit v1.2.3 From 4318279d6329bee3db6484e1c8771b87c7d8ab63 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 13:10:58 -0500 Subject: moved experimental files Nlevels_no_dopler_with_z_4wm_for_irina_pqe -> Nlevels_no_dopler_with_z_4wm --- .../Nlevels_no_dopler_with_z_4wm.xmds | 4 +- .../fast_to_slow_switch/Makefile | 11 +- .../fast_to_slow_switch/README | 18 ++ .../plot_fields_propagation_I4.gp | 1 + .../fast_to_slow_switch/pp_I2.m | 68 +--- .../plot_fields_propagation_I2.gp | 2 +- .../plot_fields_propagation_I4.gp | 15 + xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m | 63 ++-- .../pulse_split/Makefile | 56 ++++ .../pulse_split/map2dat.m | 1 + .../pulse_split/plot_fields_propagation_I2.gp | 15 + .../pulse_split/pp_I2.m | 68 ++++ .../Makefile | 54 ---- .../Nlevels_no_dopler_with_z_4wm.xmds | 349 --------------------- .../fast_to_slow_switch/Makefile | 61 ---- .../fast_to_slow_switch/README | 18 -- .../fast_to_slow_switch/map2dat.m | 1 - .../plot_fields_propagation_I2.gp | 1 - .../plot_fields_propagation_I4.gp | 1 - .../fast_to_slow_switch/pp_I2.m | 1 - .../map2dat.m | 1 - .../plot_fields_propagation_I2.gp | 15 - .../plot_fields_propagation_I4.gp | 15 - .../pp.m | 190 ----------- .../pp_I2.m | 92 ------ .../pulse_split/Makefile | 56 ---- .../pulse_split/map2dat.m | 1 - .../pulse_split/plot_fields_propagation_I2.gp | 15 - .../pulse_split/pp_I2.m | 68 ---- 29 files changed, 228 insertions(+), 1033 deletions(-) create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/README create mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/plot_fields_propagation_I4.gp mode change 100644 => 120000 xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I4.gp create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/Makefile create mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/map2dat.m create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp create mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/Makefile delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/Nlevels_no_dopler_with_z_4wm.xmds delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README delete mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m delete mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp delete mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I4.gp delete mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m delete mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/map2dat.m delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I2.gp delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I4.gp delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp.m delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/Makefile delete mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/plot_fields_propagation_I2.gp delete mode 100644 xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/pp_I2.m diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds index 44b13f7..5102c46 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds @@ -112,7 +112,7 @@ z - + @@ -310,7 +310,7 @@ - + density_matrix r11_out r22_out r33_out r44_out diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile index 697fd27..9c65fd6 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile @@ -9,11 +9,16 @@ GNUPLOT_FILES = $(wildcard *.gp) XSIL2GRAPHICS = xsil2graphics # fast light -# PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4 +PARAMS = \ + --Ndens=1e15 \ + --Lcell=10.0e-2 \ + --Temperature=1e-9 \ + --delta1=0 --delta2=0 --delta3=0 \ + --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 + + # slow light EIT #PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0 -#Fast light to Slow light switch -PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=2e7 --E2o=3e3 --E3o=6e6 --E4o=3e3 all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/README b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/README new file mode 100644 index 0000000..ec7d200 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/README @@ -0,0 +1,18 @@ +fast to slow transition +PARAMS = \ + --Ndens=1e15 \ + --Lcell=10.0e-2 \ + --Temperature=1e-9 \ + --delta1=0 --delta2=0 --delta3=0 \ + --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 + +1.5 cm seems to give largest advance +2.5 cm gives zero delay +10cm slow light again + +fast light with decay in 1.5 cm cell +PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 + + +PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=3.0e7 --E2o=3.1e5 --E3o=3.0e7 --E4o=6.3e3 +10 cm first fast light then around 8 cm transition to slow light diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/plot_fields_propagation_I4.gp new file mode 120000 index 0000000..134a380 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/plot_fields_propagation_I4.gp @@ -0,0 +1 @@ +../plot_fields_propagation_I4.gp \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m deleted file mode 100644 index 73487e0..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m +++ /dev/null @@ -1,67 +0,0 @@ -Nlevels_no_dopler_with_z_4wm - -%% field propagation -z_1=z_1*100; % z in cm -t_1=t_1*1e6; % time now measured in uS -figure(1) -set(gca,'fontsize',20); -imagesc(z_1, t_1, I2_out_1); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_2') -title('I_2') - -xskip=1; -yskip=10; -map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); - - - -print('-color','fields_propagation_I2.eps') - - - -%% fields before and after the cell -figure(2) -%set(gca,'fontsize',30); -plot( ... - t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4, ... - t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ... - ) -xlabel('t (uS)') -ylabel('I_2 (1/s)^2') -title('I_2 before and after cell') -legend('location', 'southwest'); - -[b, a]=butter(3, 0.05); -I2_out_after=I2_out_1(:,end); -I2_out_after_filtered=filtfilt(b,a,I2_out_after); -settling_time=0.8; %uS -t_good_indx=t_1> min(t_1 + settling_time); -[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); -delay_time=t_1(max_pos_after)-t_1(max_pos_before); -printf('Second field delay time = %f uS\n',delay_time); - -%set(gca,'fontsize',40); -set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','fields_before_after_cell_I2.eps') - -figure(4) -I2_max_in=max(I2_out_1(t_good_indx,1)); -I2_max_out=max(I2_out_1(t_good_indx,end)); -I2_in_norm=(I2_out_1(:,1))/I2_max_in; -I2_out_norm=(I2_out_1(:,end))/I2_max_out; -tmin=-0.05; -tmax=0.05; -indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region -plot( ... - t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4, ... - t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ... - ) -xlim([tmin,tmax],'manual'); -xlabel('t (uS)') -ylabel('I_2') -title('I_2 before and after cell normalized') -set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','probe_before_after_cell_I2_normalized.eps') - diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m new file mode 120000 index 0000000..4a1ed9b --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m @@ -0,0 +1 @@ +../pp_I2.m \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I2.gp index b6721d4..77b4870 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I2.gp +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I2.gp @@ -12,4 +12,4 @@ set ylabel "t ({/Symbol m}S)" set zlabel "I_2 (1/S)" set nokey #set view map -splot [0:1.5][-0.4:0.4] 'I2.dat' +splot [0:][-0.2:0.2] 'I2.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I4.gp new file mode 100644 index 0000000..3849e2f --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I4.gp @@ -0,0 +1,15 @@ +set terminal postscript portrait enhanced color solid size 5,3.5 +set output 'fields_propagation_I4.eps' +set dgrid3d 100,100 qnorm 4 +set pm3d map +#set contour +set hidden3d +set palette rgb 10,13,31 negative + + +set xlabel "z (cm)" +set ylabel "t ({/Symbol m}S)" +set zlabel "I_2 (1/S)" +set nokey +#set view map +splot [0:][-0.2:0.2] 'I4.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m index 45d4913..ea75b38 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m @@ -1,6 +1,6 @@ -Nlevels_no_dopler_with_z_4wm +Nlevels_no_dopler_with_z_4wm; -%% field propagation +%% field I2 propagation z_1=z_1*100; % z in cm t_1=t_1*1e6; % time now measured in uS figure(1) @@ -13,33 +13,50 @@ xlabel('z (cm)') ylabel('t (uS)') zlabel('I_2') title('I_2') +print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps') -xskip=1; -yskip=10; -%map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); +desired_x_size=200; +desired_y_size=200; +xskip=ceil(length(z_1)/desired_x_size); +yskip=ceil(length(t_1)/desired_y_size); +map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); +%% field I4 propagation +figure(2) + +imagesc(z_1, t_1, I4_out_1); colorbar +tmin=-0.4; +tmax= 0.4; +ylim([tmin,tmax],'manual'); +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_2') +title('I_2') +print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I4.eps') +map2dat('I4.dat',z_1,t_1, I4_out_1, xskip, yskip); + -print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps') %% fields before and after the cell -figure(2) +figure(3) %set(gca,'fontsize',30); plot( ... - t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4, ... - t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ... + t_1,I2_out_1(:,1),'.-;I_2 before;', "linewidth", 4 ... + ,t_1,I2_out_1(:,end), '-;I_2 after;', "linewidth", 4 ... + ,t_1,I4_out_1(:,end), '-;I_4 after;', "linewidth", 4 ... ) xlabel('t (uS)') -ylabel('I_2 (1/s)^2') -title('I_2 before and after cell') -legend('location', 'southwest'); +ylabel('I (1/s)^2') +title('Fields before and after cell') +legend('location', 'northeast'); [b, a]=butter(3, 0.05); I2_out_after=I2_out_1(:,end); I2_out_after_filtered=filtfilt(b,a,I2_out_after); -settling_time=0.8; %uS +settling_time=0.01; %uS t_good_indx=t_1> min(t_1 + settling_time); [m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); delay_time=t_1(max_pos_after)-t_1(max_pos_before); @@ -47,25 +64,29 @@ printf('Second field delay time = %f uS\n',delay_time); %set(gca,'fontsize',40); %set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_I2.eps') +print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell.eps') figure(4) I2_max_in=max(I2_out_1(t_good_indx,1)); I2_max_out=max(I2_out_1(t_good_indx,end)); +I4_max_out=max(I4_out_1(t_good_indx,end)); I2_in_norm=(I2_out_1(:,1))/I2_max_in; I2_out_norm=(I2_out_1(:,end))/I2_max_out; -tmin=-0.05; -tmax=0.05; +I4_out_norm=(I4_out_1(:,end))/I4_max_out; +tmin=-.05; +tmax=.05; indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region plot( ... - t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4, ... - t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ... + t_1(indx),I2_in_norm(indx),'.-;I_2 before;', "linewidth", 4, ... + t_1(indx),I2_out_norm(indx), '-;I_2 after;', "linewidth", 4 ... + ,t_1(indx),I4_out_norm(indx), '-;I_4 after;', "linewidth", 4 ... ) legend('location', 'southeast'); xlim([tmin,tmax],'manual'); +ylim([0.5,1],'manual'); xlabel('t (uS)') -ylabel('I_2') -title('I_2 before and after cell normalized') +ylabel('I normalized') +title('Fields before and after cell normalized') %set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','-depsc2', '-tight','-S200,120', 'probe_before_after_cell_I2_normalized.eps') +print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_normalized.eps') diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/Makefile new file mode 100644 index 0000000..95ebb0a --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/Makefile @@ -0,0 +1,56 @@ +### -*- make -*- +### This makefile can be used to build and run the XMDS examples + + +XSIL_FILES = Nlevels_no_dopler_with_z_4wm.xsil +M_FILES = $(patsubst %.xsil,%.m,$(XSIL_FILES)) +GNUPLOT_FILES = $(wildcard *.gp) + +XSIL2GRAPHICS = xsil2graphics + +# fast light +# PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4 +# slow light EIT +#PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0 +#Fast light to Slow light switch +PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=2e7 --E2o=3e3 --E3o=6e6 --E4o=3e2 + +all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png + +Nlevels_no_dopler_with_z_4wm.xsil: ../Nlevels_no_dopler_with_z_4wm.run + $< $(PARAMS) | grep "Time elapsed for simulation is:" > exact_analysis_execution_time.txt + +%.m: %.xsil + $(XSIL2GRAPHICS) $< + +plot: $(M_FILES) $(GNUPLOT_FILES) + octave pp_I2.m + gnuplot plot_fields_propagation_I2.gp + +clean: + rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf + rm -f $(pdf_targets) + rm -f $(eps_targets) + +real_clean: clean + rm -f $(png_targets) + +eps_targets = $(wildcard *.eps) +pdf_targets = $(eps_targets:%.eps=%.pdf) +png_targets = $(pdf_targets:%.pdf=%.png) + +pdf: $(pdf_targets) + +$(pdf_targets): %.pdf : %.eps + cat $< | ps2eps -B > __tt.eps + epspdf __tt.eps $@ + rm -f __tt.eps + #ps2eps -B $< | epspdf $< $@ + +png: pdf $(png_targets) + +$(png_targets): %.png : %.pdf + convert -density 300 $< $@ + +.PRECIOUS: %.run %.xsil %.m +.PHONY: all clean diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/map2dat.m new file mode 120000 index 0000000..14fae30 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/map2dat.m @@ -0,0 +1 @@ +../map2dat.m \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp new file mode 100644 index 0000000..a0cd31f --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp @@ -0,0 +1,15 @@ +set terminal postscript portrait enhanced color solid size 5,3.5 +set output 'fields_propagation_I2.eps' +set dgrid3d 100,100 qnorm 4 +set pm3d map +#set contour +set hidden3d +set palette rgb 10,13,31 negative + + +set xlabel "z (cm)" +set ylabel "t ({/Symbol m}S)" +set zlabel "I_2 (1/S)" +set nokey +#set view map +splot [0:][-0.2:0.4] 'I2.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m new file mode 100644 index 0000000..ecfdbbb --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m @@ -0,0 +1,68 @@ +Nlevels_no_dopler_with_z_4wm + +%% field propagation +z_1=z_1*100; % z in cm +t_1=t_1*1e6; % time now measured in uS +figure(1) +set(gca,'fontsize',20); +imagesc(z_1, t_1, I2_out_1); colorbar +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_2') +title('I_2') + +xskip=1; +yskip=10; +map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); + + + +print('-color','fields_propagation_I2.eps') + + + +%% fields before and after the cell +figure(2) +%set(gca,'fontsize',30); +plot( ... + t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4, ... + t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ... + ) +xlabel('t (uS)') +ylabel('I_2 (1/s)^2') +title('I_2 before and after cell') +legend('location', 'northeast'); + +[b, a]=butter(3, 0.05); +I2_out_after=I2_out_1(:,end); +I2_out_after_filtered=filtfilt(b,a,I2_out_after); +settling_time=0.01; %uS +t_good_indx=t_1> min(t_1 + settling_time); +[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); +delay_time=t_1(max_pos_after)-t_1(max_pos_before); +printf('Second field delay time = %f uS\n',delay_time); + +%set(gca,'fontsize',40); +set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts +print('-color','fields_before_after_cell_I2.eps') + +figure(4) +I2_max_in=max(I2_out_1(t_good_indx,1)); +I2_max_out=max(I2_out_1(t_good_indx,end)); +I2_in_norm=(I2_out_1(:,1))/I2_max_in; +I2_out_norm=(I2_out_1(:,end))/I2_max_out; +tmin=-0.05; +tmax=0.05; +indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region +plot( ... + t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4, ... + t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ... + ) +xlim([tmin,tmax],'manual'); +xlabel('t (uS)') +ylabel('I_2') +title('I_2 before and after cell normalized') +legend('location', 'northeast'); +set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts +print('-color','probe_before_after_cell_I2_normalized.eps') + diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/Makefile deleted file mode 100644 index 578f2b0..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/Makefile +++ /dev/null @@ -1,54 +0,0 @@ -### -*- make -*- -### This file is part of the Debian xmds package -### Copyright (C) 2006 Rafael Laboissiere -### This file is relased under the GNU General Public License -### NO WARRANTIES! - -### This makefile can be used to build and run the XMDS examples - -XMDS_FILES = $(shell ls *.xmds) -RUN_FILES = $(patsubst %.xmds,%.run,$(XMDS_FILES)) -CC_FILES = $(patsubst %.xmds,%.cc,$(XMDS_FILES)) -XSIL_FILES = $(patsubst %.xmds,%.xsil,$(XMDS_FILES)) -M_FILES = $(patsubst %.xmds,%.m,$(XMDS_FILES)) - -XMDS = xmds2 -XSIL2GRAPHICS = xsil2graphics - -all: $(M_FILES) - -%.run: %.xmds - $(XMDS) $< - mv $(patsubst %.xmds,%,$<) $@ - -%.xsil: %.run - ./$< - -%.m: %.xsil - $(XSIL2GRAPHICS) $< - -plot: $(M_FILES) - octave pp.m - -clean: - rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf - rm -f $(png_targets) - -eps_targets = $(wildcard *.eps) -pdf_targets = $(eps_targets:%.eps=%.pdf) -png_targets = $(pdf_targets:%.pdf=%.png) - -png: pdf $(png_targets) - -$(png_targets): %.png : %.pdf - convert -density 300 $< $@ - -pdf: $(pdf_targets) - -$(pdf_targets): %.pdf : %.eps - cat $< | ps2eps -B > __tt.eps - epspdf __tt.eps $@ - rm -f __tt.eps - #ps2eps -B $< | epspdf $< $@ -.PRECIOUS: %.run %.xsil %.m -.PHONY: all clean diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/Nlevels_no_dopler_with_z_4wm.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/Nlevels_no_dopler_with_z_4wm.xmds deleted file mode 100644 index 5102c46..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/Nlevels_no_dopler_with_z_4wm.xmds +++ /dev/null @@ -1,349 +0,0 @@ - - - - Nlevels_no_dopler_with_z_4wm - - Eugeniy Mikhailov - - License GPL. - - Solving 4 level atom in N-field configuration, - with field propagation along spatial axis Z - no Doppler broadening - - We assume four-wave mixing condition when w3-w4=w2-w1 i.e. fields E3 and E4 drive the same - resonance as fields E2 and E1. - - - * --------------- |4> - * \ \ - * \ E3 \ -------- |3> - * \ E4 \ / \ - * \ \ / E2 \ - * \ / \ E1 - * |2> -------------- \ - * \ \ - * \ \ - * ------------- |1> - * - - - We are solving - dE/dz+(1/c)*dE/dt=i*eta*rho_ij, where j level is higher then i. - Note that E is actually a Rabi frequency of electromagnetic field not the EM field - in xmds terms it looks like - dE_dz = i*eta*rhoij - 1/c*L[E], here we moved t dependence to Fourier space - - VERY IMPORTANT: all Rabi frequency should be given in [1/s], if you want to - normalize it to something else look drho/dt equation. - No need to renormalizes eta as long as its express through i - the upper level decay rate in the same units as Rabi frequency. - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - z - - - - - - - - E1 E2 E3 E4 - - - - - - - - - r11 r22 r33 r12 r13 r14 r23 r24 r34 r44 - - - - - - - E_field - - - - - - - - - - - - - - - - 200 200 - - - density_matrix - E_field - - - - - - - - Lt - - - E_field - density_matrix - - - - - - - - - - - - - E_field - I1_out I2_out I3_out I4_out - - - - - - - density_matrix - - r11_out r22_out r33_out r44_out - r12_re_out r12_im_out r13_re_out r13_im_out r14_re_out r14_im_out - r23_re_out r23_im_out r24_re_out r24_im_out - r34_re_out r34_im_out - - - - - - - - - diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile deleted file mode 100644 index 9c65fd6..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile +++ /dev/null @@ -1,61 +0,0 @@ -### -*- make -*- -### This makefile can be used to build and run the XMDS examples - - -XSIL_FILES = Nlevels_no_dopler_with_z_4wm.xsil -M_FILES = $(patsubst %.xsil,%.m,$(XSIL_FILES)) -GNUPLOT_FILES = $(wildcard *.gp) - -XSIL2GRAPHICS = xsil2graphics - -# fast light -PARAMS = \ - --Ndens=1e15 \ - --Lcell=10.0e-2 \ - --Temperature=1e-9 \ - --delta1=0 --delta2=0 --delta3=0 \ - --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 - - -# slow light EIT -#PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0 - -all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png - -Nlevels_no_dopler_with_z_4wm.xsil: ../Nlevels_no_dopler_with_z_4wm.run - $< $(PARAMS) | grep "Time elapsed for simulation is:" > exact_analysis_execution_time.txt - -%.m: %.xsil - $(XSIL2GRAPHICS) $< - -plot: $(M_FILES) $(GNUPLOT_FILES) - octave pp_I2.m - gnuplot plot_fields_propagation_I2.gp - -clean: - rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf - rm -f $(pdf_targets) - rm -f $(eps_targets) - -real_clean: clean - rm -f $(png_targets) - -eps_targets = $(wildcard *.eps) -pdf_targets = $(eps_targets:%.eps=%.pdf) -png_targets = $(pdf_targets:%.pdf=%.png) - -pdf: $(pdf_targets) - -$(pdf_targets): %.pdf : %.eps - cat $< | ps2eps -B > __tt.eps - epspdf __tt.eps $@ - rm -f __tt.eps - #ps2eps -B $< | epspdf $< $@ - -png: pdf $(png_targets) - -$(png_targets): %.png : %.pdf - convert -density 300 $< $@ - -.PRECIOUS: %.run %.xsil %.m -.PHONY: all clean diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README deleted file mode 100644 index ec7d200..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README +++ /dev/null @@ -1,18 +0,0 @@ -fast to slow transition -PARAMS = \ - --Ndens=1e15 \ - --Lcell=10.0e-2 \ - --Temperature=1e-9 \ - --delta1=0 --delta2=0 --delta3=0 \ - --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 - -1.5 cm seems to give largest advance -2.5 cm gives zero delay -10cm slow light again - -fast light with decay in 1.5 cm cell -PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 - - -PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=3.0e7 --E2o=3.1e5 --E3o=3.0e7 --E4o=6.3e3 -10 cm first fast light then around 8 cm transition to slow light diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m deleted file mode 120000 index 14fae30..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m +++ /dev/null @@ -1 +0,0 @@ -../map2dat.m \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp deleted file mode 120000 index e11c569..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp +++ /dev/null @@ -1 +0,0 @@ -../plot_fields_propagation_I2.gp \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I4.gp deleted file mode 120000 index 134a380..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I4.gp +++ /dev/null @@ -1 +0,0 @@ -../plot_fields_propagation_I4.gp \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m deleted file mode 120000 index 4a1ed9b..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m +++ /dev/null @@ -1 +0,0 @@ -../pp_I2.m \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/map2dat.m deleted file mode 120000 index 90fc33d..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/map2dat.m +++ /dev/null @@ -1 +0,0 @@ -../../matlab_helper_files/map2dat.m \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I2.gp deleted file mode 100644 index 77b4870..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I2.gp +++ /dev/null @@ -1,15 +0,0 @@ -set terminal postscript portrait enhanced color solid size 5,3.5 -set output 'fields_propagation_I2.eps' -set dgrid3d 100,100 qnorm 4 -set pm3d map -#set contour -set hidden3d -set palette rgb 10,13,31 negative - - -set xlabel "z (cm)" -set ylabel "t ({/Symbol m}S)" -set zlabel "I_2 (1/S)" -set nokey -#set view map -splot [0:][-0.2:0.2] 'I2.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I4.gp deleted file mode 100644 index 3849e2f..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/plot_fields_propagation_I4.gp +++ /dev/null @@ -1,15 +0,0 @@ -set terminal postscript portrait enhanced color solid size 5,3.5 -set output 'fields_propagation_I4.eps' -set dgrid3d 100,100 qnorm 4 -set pm3d map -#set contour -set hidden3d -set palette rgb 10,13,31 negative - - -set xlabel "z (cm)" -set ylabel "t ({/Symbol m}S)" -set zlabel "I_2 (1/S)" -set nokey -#set view map -splot [0:][-0.2:0.2] 'I4.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp.m deleted file mode 100644 index b88fa6e..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp.m +++ /dev/null @@ -1,190 +0,0 @@ -Nlevels_no_dopler_with_z_4wm - -%% field propagation -z_1=z_1*100; % z in cm -t_1=t_1*1e6; % time now measured in uS -figure(1) -subplot(2,2,1); imagesc(z_1, t_1, I1_out_1); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_1') -title('I_1') -subplot(2,2,2); imagesc(z_1, t_1, I2_out_1); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_2') -title('I_2') -subplot(2,2,3); imagesc(z_1, t_1, I3_out_1); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_3') -title('I_3') -subplot(2,2,4); imagesc(z_1, t_1, I4_out_1); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_4') -title('I_4') - - -print('-color','fields_propagation.eps') - - - -%% fields before and after the cell -figure(2) -subplot(2,2,1); -plot( ... - t_1,I1_out_1(:,1),'-;before;', "linewidth", 4, ... - t_1,I1_out_1(:,end), '-;after;', "linewidth", 4 ... - ) -xlabel('t (uS)') -ylabel('I_1 (1/s)^2') -title('I_1 before and after cell') -subplot(2,2,2); -plot( ... - t_1,I2_out_1(:,1),'-;before;', "linewidth", 4, ... - t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ... - ) -xlabel('t (uS)') -ylabel('I_2 (1/s)^2') -title('I_2 before and after cell') -subplot(2,2,3); -plot( ... - t_1,I3_out_1(:,1),'-;before;', "linewidth", 4, ... - t_1,I3_out_1(:,end), '-;after;', "linewidth", 4 ... - ) -xlabel('t (uS)') -ylabel('I_3 (1/s)^2') -title('I_3 before and after cell') - -[b, a]=butter(3, 0.05); -I2_out_after=I2_out_1(:,end); -I2_out_after_filtered=filtfilt(b,a,I2_out_after); -settling_time=0.8; %uS -t_good_indx=t_1> min(t_1 + settling_time); -[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); -delay_time=t_1(max_pos_after)-t_1(max_pos_before); -printf('Second field delay time = %f uS\n',delay_time); - -print('-color','fields_before_after_cell.eps') - -subplot(2,2,4); -plot( ... - t_1,I4_out_1(:,1),'-;before;', "linewidth", 4, ... - t_1,I4_out_1(:,end), '-;after;', "linewidth", 4 ... - ) -xlabel('t (uS)') -ylabel('I_3 (1/s)^2') -title('I_3 before and after cell') - -figure(4) -I2_max_in=max(I2_out_1(t_good_indx,1)); -I2_max_out=max(I2_out_1(t_good_indx,end)); -I2_in_norm=(I2_out_1(:,1))/I2_max_in; -I2_out_norm=(I2_out_1(:,end))/I2_max_out; -tmin=-0.05; -tmax=0.05; -indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region -plot( ... - t_1(indx),I2_in_norm(indx),'-;before;', "linewidth", 4, ... - t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ... - ) -xlim([tmin,tmax],'manual'); -xlabel('t (uS)') -ylabel('I_2') -title('I_2 before and after cell normalized') -print('-color','probe_before_after_cell.eps') - -return; - -%% all density matrix elements in one plot -% diagonal populations, -% upper triangle real part of coherences, -% lower diagonal imaginary part of coherences -z_2=z_2*100; % z in cm -t_2=t_2*1e6; % time now measured in uS -figure(3) -subplot(4,4,1); imagesc (z_2, t_2, r11_out_2); caxis([0,1]); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('rho_{11}') -title('rho_{11}') -subplot(4,4,6); imagesc (z_2, t_2, r22_out_2); caxis([0,1]); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('rho_{22}') -title('rho_{22}') -subplot(4,4,11); imagesc (z_2, t_2, r33_out_2); caxis([0,1]); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('rho_{33}') -title('rho_{33}') -subplot(4,4,16); imagesc (z_2, t_2, r44_out_2); caxis([0,1]); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('rho_{44}') -title('rho_{44}') -% real parts of coherences -subplot(4,4,2); imagesc(z_2, t_2, r12_re_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Real(rho_{12})') -title('Real(rho_{12})') -subplot(4,4,3); imagesc(z_2, t_2, r13_re_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Real(rho_{13})') -title('Real(rho_{13})') -subplot(4,4,4); imagesc(z_2, t_2, r14_re_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Real(rho_{14})') -title('Real(rho_{14})') -subplot(4,4,7); imagesc(z_2, t_2, r23_re_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Real(rho_{23})') -title('Real(rho_{23})') -subplot(4,4,8); imagesc(z_2, t_2, r24_re_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Real(rho_{24})') -title('Real(rho_{24})') -subplot(4,4,12); imagesc(z_2, t_2, r34_re_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Real(rho_{34})') -title('Real(rho_{34})') -% imaginary parts of coherences -subplot(4,4,5); imagesc(z_2, t_2, r12_im_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Imag(rho_{12})') -title('Imag(rho_{12})') -subplot(4,4,9); imagesc(z_2, t_2, r13_im_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Imag(rho_{13})') -title('Imag(rho_{13})') -subplot(4,4,10); imagesc(z_2, t_2, r23_im_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Imag(rho_{23})') -title('Imag(rho_{23})') -subplot(4,4,13); imagesc(z_2, t_2, r14_im_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Imag(rho_{14})') -title('Imag(rho_{14})') -subplot(4,4,14); imagesc(z_2, t_2, r24_im_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Imag(rho_{24})') -title('Imag(rho_{24})') -subplot(4,4,15); imagesc(z_2, t_2, r34_im_out_2); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('Imag(rho_{34})') -title('Imag(rho_{34})') - - diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m deleted file mode 100644 index ea75b38..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pp_I2.m +++ /dev/null @@ -1,92 +0,0 @@ -Nlevels_no_dopler_with_z_4wm; - -%% field I2 propagation -z_1=z_1*100; % z in cm -t_1=t_1*1e6; % time now measured in uS -figure(1) -%set(gca,'fontsize',20); -imagesc(z_1, t_1, I2_out_1); colorbar -tmin=-0.4; -tmax= 0.4; -ylim([tmin,tmax],'manual'); -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_2') -title('I_2') -print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps') - -desired_x_size=200; -desired_y_size=200; -xskip=ceil(length(z_1)/desired_x_size); -yskip=ceil(length(t_1)/desired_y_size); -map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); - - -%% field I4 propagation -figure(2) - -imagesc(z_1, t_1, I4_out_1); colorbar -tmin=-0.4; -tmax= 0.4; -ylim([tmin,tmax],'manual'); -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_2') -title('I_2') -print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I4.eps') -map2dat('I4.dat',z_1,t_1, I4_out_1, xskip, yskip); - - - - - -%% fields before and after the cell -figure(3) -%set(gca,'fontsize',30); -plot( ... - t_1,I2_out_1(:,1),'.-;I_2 before;', "linewidth", 4 ... - ,t_1,I2_out_1(:,end), '-;I_2 after;', "linewidth", 4 ... - ,t_1,I4_out_1(:,end), '-;I_4 after;', "linewidth", 4 ... - ) -xlabel('t (uS)') -ylabel('I (1/s)^2') -title('Fields before and after cell') -legend('location', 'northeast'); - -[b, a]=butter(3, 0.05); -I2_out_after=I2_out_1(:,end); -I2_out_after_filtered=filtfilt(b,a,I2_out_after); -settling_time=0.01; %uS -t_good_indx=t_1> min(t_1 + settling_time); -[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); -delay_time=t_1(max_pos_after)-t_1(max_pos_before); -printf('Second field delay time = %f uS\n',delay_time); - -%set(gca,'fontsize',40); -%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell.eps') - -figure(4) -I2_max_in=max(I2_out_1(t_good_indx,1)); -I2_max_out=max(I2_out_1(t_good_indx,end)); -I4_max_out=max(I4_out_1(t_good_indx,end)); -I2_in_norm=(I2_out_1(:,1))/I2_max_in; -I2_out_norm=(I2_out_1(:,end))/I2_max_out; -I4_out_norm=(I4_out_1(:,end))/I4_max_out; -tmin=-.05; -tmax=.05; -indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region -plot( ... - t_1(indx),I2_in_norm(indx),'.-;I_2 before;', "linewidth", 4, ... - t_1(indx),I2_out_norm(indx), '-;I_2 after;', "linewidth", 4 ... - ,t_1(indx),I4_out_norm(indx), '-;I_4 after;', "linewidth", 4 ... - ) -legend('location', 'southeast'); -xlim([tmin,tmax],'manual'); -ylim([0.5,1],'manual'); -xlabel('t (uS)') -ylabel('I normalized') -title('Fields before and after cell normalized') -%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_normalized.eps') - diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/Makefile deleted file mode 100644 index 95ebb0a..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/Makefile +++ /dev/null @@ -1,56 +0,0 @@ -### -*- make -*- -### This makefile can be used to build and run the XMDS examples - - -XSIL_FILES = Nlevels_no_dopler_with_z_4wm.xsil -M_FILES = $(patsubst %.xsil,%.m,$(XSIL_FILES)) -GNUPLOT_FILES = $(wildcard *.gp) - -XSIL2GRAPHICS = xsil2graphics - -# fast light -# PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4 -# slow light EIT -#PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0 -#Fast light to Slow light switch -PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=2e7 --E2o=3e3 --E3o=6e6 --E4o=3e2 - -all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png - -Nlevels_no_dopler_with_z_4wm.xsil: ../Nlevels_no_dopler_with_z_4wm.run - $< $(PARAMS) | grep "Time elapsed for simulation is:" > exact_analysis_execution_time.txt - -%.m: %.xsil - $(XSIL2GRAPHICS) $< - -plot: $(M_FILES) $(GNUPLOT_FILES) - octave pp_I2.m - gnuplot plot_fields_propagation_I2.gp - -clean: - rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf - rm -f $(pdf_targets) - rm -f $(eps_targets) - -real_clean: clean - rm -f $(png_targets) - -eps_targets = $(wildcard *.eps) -pdf_targets = $(eps_targets:%.eps=%.pdf) -png_targets = $(pdf_targets:%.pdf=%.png) - -pdf: $(pdf_targets) - -$(pdf_targets): %.pdf : %.eps - cat $< | ps2eps -B > __tt.eps - epspdf __tt.eps $@ - rm -f __tt.eps - #ps2eps -B $< | epspdf $< $@ - -png: pdf $(png_targets) - -$(png_targets): %.png : %.pdf - convert -density 300 $< $@ - -.PRECIOUS: %.run %.xsil %.m -.PHONY: all clean diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m deleted file mode 120000 index 14fae30..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m +++ /dev/null @@ -1 +0,0 @@ -../map2dat.m \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/plot_fields_propagation_I2.gp deleted file mode 100644 index a0cd31f..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/plot_fields_propagation_I2.gp +++ /dev/null @@ -1,15 +0,0 @@ -set terminal postscript portrait enhanced color solid size 5,3.5 -set output 'fields_propagation_I2.eps' -set dgrid3d 100,100 qnorm 4 -set pm3d map -#set contour -set hidden3d -set palette rgb 10,13,31 negative - - -set xlabel "z (cm)" -set ylabel "t ({/Symbol m}S)" -set zlabel "I_2 (1/S)" -set nokey -#set view map -splot [0:][-0.2:0.4] 'I2.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/pp_I2.m deleted file mode 100644 index ecfdbbb..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/pp_I2.m +++ /dev/null @@ -1,68 +0,0 @@ -Nlevels_no_dopler_with_z_4wm - -%% field propagation -z_1=z_1*100; % z in cm -t_1=t_1*1e6; % time now measured in uS -figure(1) -set(gca,'fontsize',20); -imagesc(z_1, t_1, I2_out_1); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_2') -title('I_2') - -xskip=1; -yskip=10; -map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); - - - -print('-color','fields_propagation_I2.eps') - - - -%% fields before and after the cell -figure(2) -%set(gca,'fontsize',30); -plot( ... - t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4, ... - t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ... - ) -xlabel('t (uS)') -ylabel('I_2 (1/s)^2') -title('I_2 before and after cell') -legend('location', 'northeast'); - -[b, a]=butter(3, 0.05); -I2_out_after=I2_out_1(:,end); -I2_out_after_filtered=filtfilt(b,a,I2_out_after); -settling_time=0.01; %uS -t_good_indx=t_1> min(t_1 + settling_time); -[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); -delay_time=t_1(max_pos_after)-t_1(max_pos_before); -printf('Second field delay time = %f uS\n',delay_time); - -%set(gca,'fontsize',40); -set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','fields_before_after_cell_I2.eps') - -figure(4) -I2_max_in=max(I2_out_1(t_good_indx,1)); -I2_max_out=max(I2_out_1(t_good_indx,end)); -I2_in_norm=(I2_out_1(:,1))/I2_max_in; -I2_out_norm=(I2_out_1(:,end))/I2_max_out; -tmin=-0.05; -tmax=0.05; -indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region -plot( ... - t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4, ... - t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ... - ) -xlim([tmin,tmax],'manual'); -xlabel('t (uS)') -ylabel('I_2') -title('I_2 before and after cell normalized') -legend('location', 'northeast'); -set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','probe_before_after_cell_I2_normalized.eps') - -- cgit v1.2.3 From 3186c94147bc3d1d4cb7c26fd842520edc7d4aa7 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 13:14:09 -0500 Subject: Nlevels_with_doppler_with_z_4wm: similar scripts symlinked --- .../fast_light/plot_fields_propagation_I4.gp | 1 + .../pulse_split/plot_fields_propagation_I2.gp | 16 +---- .../pulse_split/plot_fields_propagation_I4.gp | 1 + .../pulse_split/pp_I2.m | 69 +--------------------- .../slow_light/plot_fields_propagation_I4.gp | 1 + 5 files changed, 5 insertions(+), 83 deletions(-) create mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm/fast_light/plot_fields_propagation_I4.gp mode change 100644 => 120000 xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp create mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I4.gp mode change 100644 => 120000 xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m create mode 120000 xmds2/Nlevels_no_dopler_with_z_4wm/slow_light/plot_fields_propagation_I4.gp diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_light/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_light/plot_fields_propagation_I4.gp new file mode 120000 index 0000000..134a380 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_light/plot_fields_propagation_I4.gp @@ -0,0 +1 @@ +../plot_fields_propagation_I4.gp \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp deleted file mode 100644 index a0cd31f..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp +++ /dev/null @@ -1,15 +0,0 @@ -set terminal postscript portrait enhanced color solid size 5,3.5 -set output 'fields_propagation_I2.eps' -set dgrid3d 100,100 qnorm 4 -set pm3d map -#set contour -set hidden3d -set palette rgb 10,13,31 negative - - -set xlabel "z (cm)" -set ylabel "t ({/Symbol m}S)" -set zlabel "I_2 (1/S)" -set nokey -#set view map -splot [0:][-0.2:0.4] 'I2.dat' diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp new file mode 120000 index 0000000..e11c569 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp @@ -0,0 +1 @@ +../plot_fields_propagation_I2.gp \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I4.gp new file mode 120000 index 0000000..134a380 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I4.gp @@ -0,0 +1 @@ +../plot_fields_propagation_I4.gp \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m deleted file mode 100644 index ecfdbbb..0000000 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m +++ /dev/null @@ -1,68 +0,0 @@ -Nlevels_no_dopler_with_z_4wm - -%% field propagation -z_1=z_1*100; % z in cm -t_1=t_1*1e6; % time now measured in uS -figure(1) -set(gca,'fontsize',20); -imagesc(z_1, t_1, I2_out_1); colorbar -xlabel('z (cm)') -ylabel('t (uS)') -zlabel('I_2') -title('I_2') - -xskip=1; -yskip=10; -map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip); - - - -print('-color','fields_propagation_I2.eps') - - - -%% fields before and after the cell -figure(2) -%set(gca,'fontsize',30); -plot( ... - t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4, ... - t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ... - ) -xlabel('t (uS)') -ylabel('I_2 (1/s)^2') -title('I_2 before and after cell') -legend('location', 'northeast'); - -[b, a]=butter(3, 0.05); -I2_out_after=I2_out_1(:,end); -I2_out_after_filtered=filtfilt(b,a,I2_out_after); -settling_time=0.01; %uS -t_good_indx=t_1> min(t_1 + settling_time); -[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx)); -delay_time=t_1(max_pos_after)-t_1(max_pos_before); -printf('Second field delay time = %f uS\n',delay_time); - -%set(gca,'fontsize',40); -set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','fields_before_after_cell_I2.eps') - -figure(4) -I2_max_in=max(I2_out_1(t_good_indx,1)); -I2_max_out=max(I2_out_1(t_good_indx,end)); -I2_in_norm=(I2_out_1(:,1))/I2_max_in; -I2_out_norm=(I2_out_1(:,end))/I2_max_out; -tmin=-0.05; -tmax=0.05; -indx=(t_1>=tmin & t_1<=tmax); % soom in in time to this region -plot( ... - t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4, ... - t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ... - ) -xlim([tmin,tmax],'manual'); -xlabel('t (uS)') -ylabel('I_2') -title('I_2 before and after cell normalized') -legend('location', 'northeast'); -set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts -print('-color','probe_before_after_cell_I2_normalized.eps') - diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m new file mode 120000 index 0000000..4a1ed9b --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m @@ -0,0 +1 @@ +../pp_I2.m \ No newline at end of file diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/slow_light/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/slow_light/plot_fields_propagation_I4.gp new file mode 120000 index 0000000..134a380 --- /dev/null +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/slow_light/plot_fields_propagation_I4.gp @@ -0,0 +1 @@ +../plot_fields_propagation_I4.gp \ No newline at end of file -- cgit v1.2.3 From da6d030b4364e9016899cf86bbbb7e19af4d3d18 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 14:54:19 -0500 Subject: plot title name fix --- xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m | 4 ++-- xmds2/Nlevels_with_doppler_with_z_4wm/pp_I2.m | 4 ++-- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m index ea75b38..3e5c174 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m @@ -31,8 +31,8 @@ tmax= 0.4; ylim([tmin,tmax],'manual'); xlabel('z (cm)') ylabel('t (uS)') -zlabel('I_2') -title('I_2') +zlabel('I_4') +title('I_4') print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I4.eps') map2dat('I4.dat',z_1,t_1, I4_out_1, xskip, yskip); diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/pp_I2.m b/xmds2/Nlevels_with_doppler_with_z_4wm/pp_I2.m index f632d0a..6b90568 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/pp_I2.m +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/pp_I2.m @@ -31,8 +31,8 @@ tmax= 0.4; ylim([tmin,tmax],'manual'); xlabel('z (cm)') ylabel('t (uS)') -zlabel('I_2') -title('I_2') +zlabel('I_4') +title('I_4') print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I4.eps') map2dat('I4.dat',z_1,t_1, I4_out_1, xskip, yskip); -- cgit v1.2.3 From 4c2c6743f88505c0a1f9cf22e2d59ff1ce69c95c Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 14:56:03 -0500 Subject: Makefile pulse width parameter addded --- .../Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile | 10 +++++++--- 1 file changed, 7 insertions(+), 3 deletions(-) diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile index 9c65fd6..fcc817c 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile @@ -13,8 +13,9 @@ PARAMS = \ --Ndens=1e15 \ --Lcell=10.0e-2 \ --Temperature=1e-9 \ + --Pwidth=0.1e-6 \ --delta1=0 --delta2=0 --delta3=0 \ - --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 + --E1o=1.9e7 --E2o=1e2 --E3o=3.8e7 --E4o=0 # slow light EIT @@ -28,9 +29,12 @@ Nlevels_no_dopler_with_z_4wm.xsil: ../Nlevels_no_dopler_with_z_4wm.run %.m: %.xsil $(XSIL2GRAPHICS) $< -plot: $(M_FILES) $(GNUPLOT_FILES) - octave pp_I2.m +pretty_plots: plot $(GNUPLOT_FILES) gnuplot plot_fields_propagation_I2.gp + gnuplot plot_fields_propagation_I4.gp + +plot: $(M_FILES) + octave pp_I2.m clean: rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf -- cgit v1.2.3 From 644c58d69b8be68cec5fceb97e8890b90b082fc2 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 14:57:15 -0500 Subject: Compile only for default target --- xmds2/Nlevels_no_dopler_with_z_4wm/Makefile | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm/Makefile index 578f2b0..f8b3559 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/Makefile +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/Makefile @@ -15,7 +15,7 @@ M_FILES = $(patsubst %.xmds,%.m,$(XMDS_FILES)) XMDS = xmds2 XSIL2GRAPHICS = xsil2graphics -all: $(M_FILES) +all: $(RUN_FILES) %.run: %.xmds $(XMDS) $< -- cgit v1.2.3 From 416702a44bdc76cf49e9cd3b980fdf73f50d8caf Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 15:29:44 -0500 Subject: output is more finite --- .../Nlevels_with_doppler_with_z_4wm.xmds | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds index d63704f..c63870e 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds @@ -117,7 +117,7 @@ z - + @@ -368,7 +368,7 @@ - + E_field_avgd I1_out I2_out I3_out I4_out Date: Thu, 12 Jan 2012 15:30:56 -0500 Subject: Doppler feature added: --Ndens now is taken in account --- .../Nlevels_with_doppler_with_z_4wm.xmds | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds index c63870e..d8fd6a4 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds @@ -50,7 +50,8 @@ const double lambda=794.7e-9; //wavelength in m const double Kvec = 2*M_PI/lambda; // k-vector const double Gamma_super=6*(2*M_PI*1e6); // characteristic decay rate of upper level used for eta calculations expressed in [1/s] - const double eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] + // eta will be calculated in the sequence -> filter element to use --Ndens parameter + double eta = 0; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] // --------- Atom and cell properties ------------------------- // range of Maxwell distribution atomic velocities @@ -287,6 +288,11 @@ + + + -- cgit v1.2.3 From 76b91d4697a673ced8a474982e5f8fca79a75c77 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 15:56:12 -0500 Subject: no Doppler feature added: --Ndens now taken in account --- .../Nlevels_no_dopler_with_z_4wm.xmds | 10 ++++++++-- .../Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds | 10 ++++++++-- 2 files changed, 16 insertions(+), 4 deletions(-) diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds index 5102c46..045adac 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds @@ -36,7 +36,7 @@ VERY IMPORTANT: all Rabi frequency should be given in [1/s], if you want to normalize it to something else look drho/dt equation. - No need to renormalizes eta as long as its express through i + No need to renormalizes eta as long as its express through the upper level decay rate in the same units as Rabi frequency. @@ -49,7 +49,8 @@ const double lambda=794.7e-9; //wavelength in m const double Kvec = 2*M_PI/lambda; // k-vector const double Gamma_super=6*(2*M_PI*1e6); // characteristic decay rate of upper level used for eta calculations expressed in [1/s] - const double eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] + // eta will be calculated in the sequence -> filter element to use --Ndens parameter + double eta = 0; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] // --------- Atom and cell properties ------------------------- // range of Maxwell distribution atomic velocities @@ -218,6 +219,11 @@ + + + diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds index e0440ef..d662779 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds @@ -41,7 +41,7 @@ VERY IMPORTANT: all Rabi frequency should be given in [1/s], if you want to normalize it to something else look drho/dt equation. - No need to renormalizes eta as long as its express through i + No need to renormalizes eta as long as its express through the upper level decay rate in the same units as Rabi frequency. @@ -54,7 +54,8 @@ const double lambda=794.7e-9; //wavelength in m const double Kvec = 2*M_PI/lambda; // k-vector const double Gamma_super=6*(2*M_PI*1e6); // characteristic decay rate of upper level used for eta calculations expressed in [1/s] - const double eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] + // eta will be calculated in the sequence -> filter element to use --Ndens parameter + double eta = 0; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] // --------- Atom and cell properties ------------------------- // range of Maxwell distribution atomic velocities @@ -225,6 +226,11 @@ + + + -- cgit v1.2.3 From b11e6e19fcf6f4c02a5efb82831bdecd39c00d84 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 16:05:04 -0500 Subject: mat2dat: timing output commented out --- matlab_helper_files/map2dat.m | 21 ++++++++++----------- 1 file changed, 10 insertions(+), 11 deletions(-) diff --git a/matlab_helper_files/map2dat.m b/matlab_helper_files/map2dat.m index 4f4ca52..853ddd8 100644 --- a/matlab_helper_files/map2dat.m +++ b/matlab_helper_files/map2dat.m @@ -12,8 +12,8 @@ Ny=length(y); Nxs=ceil(Nx/xskip); Nys=ceil(Ny/yskip); points=zeros(1,3*Nxs*Nys); -%points=[]; -tic; + +%tic; for i=1:Nxs for k=1:Nys %points=[points x(i*xskip) y(k*yskip) z(k*yskip,i*xskip)]; @@ -22,12 +22,11 @@ for i=1:Nxs points((i-1)*(Nys)*3+3*(k-1)+3) = z(1+(k-1)*yskip,1+(i-1)*xskip); end end -disp('=== points formation complete ==='); -toc -tic; -%points -fd = fopen(outfile, "wt"); -fprintf (fd, "%g %g %g\n", points); -fclose(fd); -disp('=== points saving complete ==='); -toc; +%disp('=== points formation complete ==='); +%toc +%tic; +%fd = fopen(outfile, "wt"); +%fprintf (fd, "%g %g %g\n", points); +%fclose(fd); +%disp('=== points saving complete ==='); +%toc; -- cgit v1.2.3 From f1d0fb6eafd1261b7141e1357bc663a0de41478c Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 16:05:47 -0500 Subject: white space fix --- .../Nlevels_with_doppler_with_z_4wm.xmds | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds index d8fd6a4..7f13652 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds @@ -37,7 +37,7 @@ VERY IMPORTANT: all Rabi frequency should be given in [1/s], if you want to normalize it to something else look drho/dt equation. - No need to renormalizes eta as long as its express through + No need to renormalizes eta as long as its express through the upper level decay rate in the same units as Rabi frequency. -- cgit v1.2.3 From 49fdac80f40267b5decf1c82b6c1913b3975116a Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 16:15:29 -0500 Subject: some updates in attempt to find fast light with Doppler --- .../fast_to_slow_switch/Makefile | 13 ++++++++----- 1 file changed, 8 insertions(+), 5 deletions(-) diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/Makefile b/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/Makefile index 12efac9..7dbc3f4 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/Makefile +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/Makefile @@ -13,10 +13,11 @@ XSIL2GRAPHICS = xsil2graphics #PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=0.1e7 --E2o=1e4 --E3o=.3e7 --E4o=0 --Lcell=1.5e-2 --Temperature=.0001 PARAMS = \ --Ndens=1e15 \ - --Lcell=1.5e-2 \ - --Temperature=1e-9 \ + --Lcell=10.0e-2 \ + --Temperature=5 \ + --Pwidth=0.4e-6 \ --delta1=0 --delta2=0 --delta3=0 \ - --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 + --E1o=2e7 --E2o=1e2 --E3o=4e7 --E4o=1e0 # slow light EIT @@ -30,11 +31,13 @@ Nlevels_with_doppler_with_z_4wm.xsil: ../Nlevels_with_doppler_with_z_4wm.run %.m: %.xsil $(XSIL2GRAPHICS) $< -plot: $(M_FILES) $(GNUPLOT_FILES) - octave pp_I2.m +pretty_plots: plot $(GNUPLOT_FILES) gnuplot plot_fields_propagation_I2.gp gnuplot plot_fields_propagation_I4.gp +plot: $(M_FILES) + octave pp_I2.m + clean: rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf rm -f $(pdf_targets) -- cgit v1.2.3 From dfa371fd839e92434503628457da13df3e6ba0b4 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 12 Jan 2012 16:18:46 -0500 Subject: Doppler: documented fast light --- .../fast_to_slow_switch/README | 11 ++++++++++- 1 file changed, 10 insertions(+), 1 deletion(-) diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/README b/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/README index 035178c..34ec5be 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/README +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/fast_to_slow_switch/README @@ -1,4 +1,13 @@ -Temperature = 1e-2K +Fast light with advancement 0.006000 uS it is barely visible but there +PARAMS = \ + --Ndens=1e15 \ + --Lcell=10.0e-2 \ + --Temperature=5 \ + --Pwidth=0.4e-6 \ + --delta1=0 --delta2=0 --delta3=0 \ + --E1o=2e7 --E2o=1e2 --E3o=4e7 --E4o=1e0 + + Fast light with decay in 1.5 cm cell PARAMS = --delta1=0 --delta2=0 --delta3=0 --E1o=0.1e7 --E2o=1e4 --E3o=.3e7 --E4o=0 -- cgit v1.2.3 From 8a56ad5f67b534aca62b6e6c561aaa9fac0893f6 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Fri, 13 Jan 2012 00:15:25 -0500 Subject: Doppler: Temperature and Ndens argument properly treated --- .../Nlevels_with_doppler_with_z_4wm.xmds | 26 +++++++++++++++------- xmds2/TODO | 15 ++----------- 2 files changed, 20 insertions(+), 21 deletions(-) diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds index 7f13652..bcb6b39 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds @@ -50,7 +50,7 @@ const double lambda=794.7e-9; //wavelength in m const double Kvec = 2*M_PI/lambda; // k-vector const double Gamma_super=6*(2*M_PI*1e6); // characteristic decay rate of upper level used for eta calculations expressed in [1/s] - // eta will be calculated in the sequence -> filter element to use --Ndens parameter + // eta will be calculated in the section double eta = 0; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] // --------- Atom and cell properties ------------------------- @@ -59,9 +59,10 @@ // above mass expression is written as (expression is isotopic_mass * atomic_mass_unit) // Average sqrt(v^2) in Maxwell distribution for one dimension - const double v_thermal_averaged=sqrt(k_boltzmann*Temperature/mass); + // Maxwell related parameters will be calculated in section + double v_thermal_averaged=0; // Maxwell distribution velocities range to take in account in [m/s] - const double V_maxwell_min = -4*v_thermal_averaged, V_maxwell_max = -V_maxwell_min; // there is almost zero probability for higher velocity p(4*v_av) = 3.3e-04 * p(0) + double V_maxwell_min = 0, V_maxwell_max = 0; // repopulation rate (atoms flying in/out the laser beam) in [1/s] const double gt=0.01 *(2*M_PI*1e6); @@ -107,6 +108,20 @@ I am guessing detunings are too large and thus it became a stiff equation--> + + + 0 to provide range for Maxwell velocity distribution\n"); + v_thermal_averaged=sqrt(k_boltzmann*Temperature/mass); + // Maxwell distribution velocities range to take in account in [m/s] + // there is almost zero probability for higher velocity p(4*v_av) = 3.3e-04 * p(0) + V_maxwell_min = -4*v_thermal_averaged; V_maxwell_max = -V_maxwell_min; + + // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] + eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI; + ]]> @@ -288,11 +303,6 @@ - - - diff --git a/xmds2/TODO b/xmds2/TODO index c152093..c77e332 100644 --- a/xmds2/TODO +++ b/xmds2/TODO @@ -1,14 +1,3 @@ -- TODO Fix Ndens argument treatment: +- TODO Fix Temperature > 10 misbehavior + Suspected stiff equation at this condition -Right now code is oblivius to --Ndens -parameter, since eta ~1/Ndens sits in global parameter section and already -defined. So it is not recalculated after the parameter assignment. - -Ugly hack would be to make variable eta_times_Ndens instead of eta, and -recalculate eta=eta_times_Ndens/Ndens in the ODE part around where it -needed. It will have performance penalty but probably not too large. - - -- TODO Fix Temperature argument treatment: - -Same problem as above -- cgit v1.2.3 From 3c5b44a1d51edfa24dfa22190265f13285c8d421 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Fri, 13 Jan 2012 00:25:59 -0500 Subject: no Doppler: Temperature and Ndens argument properly treated --- .../Nlevels_no_dopler_with_z_4wm.xmds | 24 ++++++++++++++-------- ...ls_no_dopler_with_z_4wm_with_perturbations.xmds | 24 ++++++++++++++-------- 2 files changed, 32 insertions(+), 16 deletions(-) diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds index 045adac..bb0feb2 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds +++ b/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds @@ -49,7 +49,7 @@ const double lambda=794.7e-9; //wavelength in m const double Kvec = 2*M_PI/lambda; // k-vector const double Gamma_super=6*(2*M_PI*1e6); // characteristic decay rate of upper level used for eta calculations expressed in [1/s] - // eta will be calculated in the sequence -> filter element to use --Ndens parameter + // eta will be calculated in the section double eta = 0; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] // --------- Atom and cell properties ------------------------- @@ -58,9 +58,10 @@ // above mass expression is written as (expression is isotopic_mass * atomic_mass_unit) // Average sqrt(v^2) in Maxwell distribution for one dimension - const double v_thermal_averaged=sqrt(k_boltzmann*Temperature/mass); + // Maxwell related parameters will be calculated in section + double v_thermal_averaged=0; // Maxwell distribution velocities range to take in account in [m/s] - const double V_maxwell_min = -4*v_thermal_averaged, V_maxwell_max = -V_maxwell_min; // there is almost zero probability for higher velocity p(4*v_av) = 3.3e-04 * p(0) + double V_maxwell_min = 0, V_maxwell_max = 0; // repopulation rate (atoms flying in/out the laser beam) in [1/s] const double gt=0.01 *(2*M_PI*1e6); @@ -102,6 +103,18 @@ I am guessing detunings are too large and thus it became a stiff equation--> + + + @@ -219,11 +232,6 @@ - - - diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds index d662779..e336db7 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds @@ -54,7 +54,7 @@ const double lambda=794.7e-9; //wavelength in m const double Kvec = 2*M_PI/lambda; // k-vector const double Gamma_super=6*(2*M_PI*1e6); // characteristic decay rate of upper level used for eta calculations expressed in [1/s] - // eta will be calculated in the sequence -> filter element to use --Ndens parameter + // eta will be calculated in the section double eta = 0; // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)] // --------- Atom and cell properties ------------------------- @@ -63,9 +63,10 @@ // above mass expression is written as (expression is isotopic_mass * atomic_mass_unit) // Average sqrt(v^2) in Maxwell distribution for one dimension - const double v_thermal_averaged=sqrt(k_boltzmann*Temperature/mass); + // Maxwell related parameters will be calculated in section + double v_thermal_averaged=0; // Maxwell distribution velocities range to take in account in [m/s] - const double V_maxwell_min = -4*v_thermal_averaged, V_maxwell_max = -V_maxwell_min; // there is almost zero probability for higher velocity p(4*v_av) = 3.3e-04 * p(0) + double V_maxwell_min = 0, V_maxwell_max = 0; // repopulation rate (atoms flying in/out the laser beam) in [1/s] const double gt=0.01 *(2*M_PI*1e6); @@ -108,6 +109,18 @@ I am guessing detunings are too large and thus it became a stiff equation--> + + + @@ -226,11 +239,6 @@ - - - -- cgit v1.2.3