From 45f970c010cd2981f86eb4b7ce3cd25f9947e96f Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Wed, 16 Jun 2010 19:35:56 +0000 Subject: small change in initial conditions --- psr/psr.m | 9 +++------ 1 file changed, 3 insertions(+), 6 deletions(-) (limited to 'psr/psr.m') diff --git a/psr/psr.m b/psr/psr.m index 7c8f604..6e3f0c9 100644 --- a/psr/psr.m +++ b/psr/psr.m @@ -29,7 +29,7 @@ field_description; detuning_p=0; N_detun_steps=100; %detuning_p_min=-B_field*gmg*4; % span +/-4 Zeeman splitting -detuning_p_min=-10040.0; +detuning_p_min=-40.0; detuning_p_max=-detuning_p_min; detuning_freq=zeros(1,N_detun_steps+1); kappa_p =zeros(1,N_detun_steps+1); @@ -91,7 +91,7 @@ phi=pi*2/8; % theta is angle between lab z axis (light propagation direction) and magnetic field axis (z') theta=0; % psi_el is the ellipticity parameter (phase difference between left and right polarization) -psi_el=-3/180*pi; +psi_el=-2/180*pi; % we define light as linearly polarized % where phi is angle between light polarization and axis x @@ -133,11 +133,8 @@ save '/tmp/problem_definition.mat' problems_cell_array atom_properties detuning fprintf (stderr, "now really hard calculations begin\n"); fflush (stderr); % once we define all problems the main job is done here -%kappa_p=cellfun( @susceptibility_steady_state_at_freq, problems_cell_array); -%kappa_p=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array); [xi_linear, xi_left, xi_right]=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array); -%relative_transmission_vs_detuning=parcellfun(2, @total_relative_transmission, problems_cell_array); -%relative_transmission_vs_detuning=cellfun(@total_relative_transmission, problems_cell_array); +%[xi_linear, xi_left, xi_right]=cellfun( @susceptibility_steady_state_at_freq, problems_cell_array); save '/tmp/xi_vs_detuning.mat' detuning_freq xi_linear xi_left xi_right ; %save '/tmp/relative_transmission_vs_detuning.mat' detuning_freq relative_transmission_vs_detuning; -- cgit v1.2.3