From e7f9feea33c3a4f10a4581c64f0077f7588cf76f Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Tue, 15 Dec 2009 06:53:47 +0000 Subject: small modification and removing redundant initialization --- liouville.m | 17 +++++++++++------ 1 file changed, 11 insertions(+), 6 deletions(-) (limited to 'liouville.m') diff --git a/liouville.m b/liouville.m index c6603b4..1225304 100644 --- a/liouville.m +++ b/liouville.m @@ -10,7 +10,9 @@ useful_functions; useful_constants; % load atom energy levels and decay description +%rb87_D1_line; four_levels_with_polarization; +%four_levels_with_polarization; %four_levels; %three_levels; %two_levels; @@ -18,7 +20,7 @@ four_levels_with_polarization; % load EM field description field_description; -Nfreq=length(modulation_freq); +%Nfreq=length(modulation_freq); @@ -32,10 +34,8 @@ kappa_p =zeros(1,N_detun_steps+1); kappa_m =zeros(1,N_detun_steps+1); detun_step=(detuning_p_max-detuning_p_min)/N_detun_steps; -% now we create Liouville indexes list -[N, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c]=unfold_density_matrix(Nlevels,Nfreq); -rhoLiouville=zeros(N,1); - +fprintf (stderr, "calculating atom properties\n"); +fflush (stderr); % calculate E_field independent properties of the atom % to be used as sub matrix templates for Liouville operator matrix [L0m, polarizability_m]=L0_and_polarization_submatrices( ... @@ -48,8 +48,11 @@ atom_properties.L0m=L0m; atom_properties.polarizability_m=polarizability_m; atom_properties.dipole_elements=dipole_elements; +fprintf (stderr, "tuning laser in forloop to set conditions vs detuning\n"); +fflush (stderr); for detuning_p_cntr=1:N_detun_steps+1; - wp0=w12; + wp0=w_pf1; + wd=w_pf1-w_hpf_ground; detuning_p=detuning_p_min+detun_step*(detuning_p_cntr-1); wp=wp0+detuning_p; wm=wd-(wp-wd); @@ -72,6 +75,8 @@ for detuning_p_cntr=1:N_detun_steps+1; detuning_freq(detuning_p_cntr)=detuning_p; endfor +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); -- cgit v1.2.3