diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-15 21:25:54 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-15 21:25:54 +0000 |
| commit | e5b9a01f36acaa186b462feae2cff4da62d5491a (patch) | |
| tree | 16519209e6523f63cab0eba6d7429ddd1f5bf470 /useful_functions.m | |
| parent | e25827cbfc12e4064c6ad5596369fd7e37e1389a (diff) | |
| download | multi_mode_eit-e5b9a01f36acaa186b462feae2cff4da62d5491a.tar.gz multi_mode_eit-e5b9a01f36acaa186b462feae2cff4da62d5491a.zip | |
call to kroneker delta function replaced to array element of eye matrix
Diffstat (limited to 'useful_functions.m')
| -rw-r--r-- | useful_functions.m | 13 |
1 files changed, 7 insertions, 6 deletions
diff --git a/useful_functions.m b/useful_functions.m index 4b239a5..913b505 100644 --- a/useful_functions.m +++ b/useful_functions.m @@ -66,6 +66,7 @@ function [L0m, polarizability_m]=L0_and_polarization_submatrices( ... % now we create Liouville indexes list [Ndummy, rhoLiouville_w_notused, rhoLiouville_r, rhoLiouville_c]=unfold_density_matrix(Nlevels,1); + kron_delta_m=eye(Nlevels); % note that L0 and decay parts depend only on combination of indexes % jk,mn but repeats itself for every frequency L0m=zeros(rho_size); % (Nlevels^2)x(Nlevels^2) matrix @@ -86,17 +87,17 @@ function [L0m, polarizability_m]=L0_and_polarization_submatrices( ... n=rhoLiouville_c(s); % calculate unperturbed part (Hamiltonian without EM field) - L0m(p,s)=H0(j,m)*kron_delta(k,n)-H0(n,k)*kron_delta(j,m); + L0m(p,s)=H0(j,m)*kron_delta_m(k,n)-H0(n,k)*kron_delta_m(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) ... - - kron_delta(m,n)*kron_delta(j,k)*g_decay(m,j) ; - polarizability_m.linear(p,s)= ( dipole_elements.linear(j,m)*kron_delta(k,n)-dipole_elements.linear(n,k)*kron_delta(j,m) ); - polarizability_m.left(p,s)= ( dipole_elements.left(j,m)*kron_delta(k,n)-dipole_elements.left(n,k)*kron_delta(j,m) ); - polarizability_m.right(p,s)= ( dipole_elements.right(j,m)*kron_delta(k,n)-dipole_elements.right(n,k)*kron_delta(j,m) ); + )* kron_delta_m(j,m)*kron_delta_m(k,n) ... + - kron_delta_m(m,n)*kron_delta_m(j,k)*g_decay(m,j) ; + polarizability_m.linear(p,s)= ( dipole_elements.linear(j,m)*kron_delta_m(k,n)-dipole_elements.linear(n,k)*kron_delta_m(j,m) ); + polarizability_m.left(p,s)= ( dipole_elements.left(j,m)*kron_delta_m(k,n)-dipole_elements.left(n,k)*kron_delta_m(j,m) ); + polarizability_m.right(p,s)= ( dipole_elements.right(j,m)*kron_delta_m(k,n)-dipole_elements.right(n,k)*kron_delta_m(j,m) ); endfor endfor L0m=-im_one/hbar*L0m - decay_part_m; |