diff options
Diffstat (limited to 'useful_functions.m')
| -rw-r--r-- | useful_functions.m | 13 |
1 files changed, 6 insertions, 7 deletions
diff --git a/useful_functions.m b/useful_functions.m index 2ae4ced..1765dd9 100644 --- a/useful_functions.m +++ b/useful_functions.m @@ -59,9 +59,8 @@ function L=Liouville_operator_matrix( modulation_freq, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c ) %------------------------- - useful_constants; L=zeros(N); % NxN matrix - + useful_constants; Nfreq=length(modulation_freq); for p=1:N @@ -70,16 +69,15 @@ function L=Liouville_operator_matrix( k=rhoLiouville_c(p); m=rhoLiouville_r(s); n=rhoLiouville_c(s); - % we garanted to know frequency of final and initial rhoLiouville + % we guaranteed to know frequency of final and initial rhoLiouville w1i=rhoLiouville_w(p); w2i=rhoLiouville_w(s); w_jk=modulation_freq(w1i); w_mn=modulation_freq(w2i); % thus we know L matrix element frequency which we need to match w_l=w_jk-w_mn; - % lets search this wrequency in the list of available frequencyes - % but since we not garanteed to find it lets assign temporary 0 to Liouville matrix element - L(p,s)=0; + % lets search this frequency in the list of available frequencies + % but since we not guaranteed to find it lets assign temporary 0 to Liouville matrix element decay_part=0; Lt=0; for w3i=1:Nfreq @@ -108,7 +106,8 @@ function L=Liouville_operator_matrix( if ((p == s)) Lt+=-im_one*w_jk; endif - L(p,s)=Lt; + Lps=Lt; + L(p,s)=Lps; endfor endfor endfunction |