timing output is added

1 files changed, 13 insertions, 0 deletions

diff --git a/test.m b/test.m
index 17a1622..4df7c26 100644
--- a/test.m
+++ b/test.m
@@ -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);