diff options
Diffstat (limited to 'psr/psr.m')
| -rw-r--r-- | psr/psr.m | 9 |
1 files changed, 3 insertions, 6 deletions
@@ -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; |