diff options
Diffstat (limited to 'test.m')
| -rw-r--r-- | test.m | 13 |
1 files changed, 13 insertions, 0 deletions
@@ -136,10 +136,16 @@ E_field.left=Efld; E_field.linear=Efld; % calculate E_field independent properties of athe atom % to be used as sub matrix templates for Liouville operator matrix +t1 = clock (); % we will use this latter to calculate elapsed time [L0m, polarizability_m]=L0_and_polarization_submatrices( ... Nlevels, ... H0, g_decay, g_dephasing, dipole_elements ... ); +elapsed_time = etime (clock (), t1); +fprintf (stderr, "elapsed time for polarazability creation is %.3f sec\n",elapsed_time); +fflush (stderr); + +t1 = clock (); % we will use this latter to calculate elapsed time % Liouville operator matrix construction L=Liouville_operator_matrix( N, @@ -153,10 +159,17 @@ L=Liouville_operator_matrix( [rhoLiouville_dot, L]=constrain_rho_and_match_L( N, L, modulation_freq, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c); +elapsed_time = etime (clock (), t1); +fprintf (stderr, "elapsed time for full L creation is %.3f sec\n",elapsed_time); +fflush (stderr); +t1 = clock (); % we will use this latter to calculate elapsed time %solving for density matrix vector rhoLiouville=L\rhoLiouville_dot; + elapsed_time = etime (clock (), t1); +fprintf (stderr, "elapsed time for rhoL solving is %.3f sec\n",elapsed_time); +fflush (stderr); rho_2=rhoOfFreq(rhoLiouville, 2, Nlevels); rho_1=rhoOfFreq(rhoLiouville, 1, Nlevels); |