diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-12 22:46:42 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-12 22:46:42 +0000 |
| commit | b9a5a0b8349bdebd5dd1fec3460d002d964d5100 (patch) | |
| tree | e818b80b454ec23c5bdbc7583032acc8d0be1d2f /useful_functions.m | |
| parent | c46ae2c7bbe38393813a9cca11872f2550898f8a (diff) | |
| download | multi_mode_eit-b9a5a0b8349bdebd5dd1fec3460d002d964d5100.tar.gz multi_mode_eit-b9a5a0b8349bdebd5dd1fec3460d002d964d5100.zip | |
matlab compatible long splitted strings
Diffstat (limited to 'useful_functions.m')
| -rw-r--r-- | useful_functions.m | 22 |
1 files changed, 13 insertions, 9 deletions
diff --git a/useful_functions.m b/useful_functions.m index a1135f7..606fe49 100644 --- a/useful_functions.m +++ b/useful_functions.m @@ -82,8 +82,12 @@ function [L0m, polarizability_m]=L0_and_polarization_submatrices( ... % calculate unperturbed part (Hamiltonian without EM field) L0m(p,s)=H0(j,m)*kron_delta(k,n)-H0(n,k)*kron_delta(j,m); - decay_part_m(p,s)=\ - ( decay_total(g_decay,k)/2 + decay_total(g_decay,j)/2 + g_dephasing(j,k) )* kron_delta(j,m)*kron_delta(k,n) \ + decay_part_m(p,s)= ... + ( ... + decay_total(g_decay,k)/2 ... + + decay_total(g_decay,j)/2 ... + + g_dephasing(j,k) ... + )* kron_delta(j,m)*kron_delta(k,n) ... - kron_delta(m,n)*kron_delta(j,k)*g_decay(m,j) ; polarizability_m(p,s)= ( dipole_elements(j,m)*kron_delta(k,n)-dipole_elements(n,k)*kron_delta(j,m) ); endfor @@ -93,11 +97,11 @@ endfunction % Liouville operator matrix construction % based on recipe from Eugeniy Mikhailov thesis -function L=Liouville_operator_matrix( - N, - L0m, polarizability_m, - E_field, - modulation_freq, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c +function L=Liouville_operator_matrix( ... + N, ... + L0m, polarizability_m, ... + E_field, ... + modulation_freq, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c ... ) %------------------------- useful_constants; @@ -163,8 +167,8 @@ endfunction % sum(rho_ii)=1 (sum of all populations (with zero modulation frequency) scales to 1 % we will replace first row of Liouville operator with this condition % thus rhoLiouville_dot(1)=1 -function [rhoLiouville_dot, L]=constrain_rho_and_match_L( - N, L, +function [rhoLiouville_dot, L]=constrain_rho_and_match_L( ... + N, L, ... modulation_freq, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c) for i=1:N w2i=rhoLiouville_w(i); |