From d0071faa593907567ffafc0fedf63f209ab1293f Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Wed, 9 Dec 2009 19:47:57 +0000
Subject: moved some useful functions to separate file

---
 useful_functions.m | 30 ++++++++++++++++++++++++++++++
 1 file changed, 30 insertions(+)
 create mode 100644 useful_functions.m

(limited to 'useful_functions.m')

diff --git a/useful_functions.m b/useful_functions.m
new file mode 100644
index 0000000..748a677
--- /dev/null
+++ b/useful_functions.m
@@ -0,0 +1,30 @@
+1;
+
+% calculate total decay for particular level taking in account all branches
+function ret=decay_total(g_decay,i)
+	ret=0;
+	for k=1:size(g_decay)(1)
+		ret=ret+g_decay(i,k);
+	endfor
+endfunction
+
+% kroneker delta symbol
+function ret=kron_delta(i,j)
+	if ((i==j))
+		ret=1;
+	else
+		ret=0;
+	endif
+endfunction
+
+% this function create from Liouville density vector 
+% the density matrix with given modulation frequency
+function rho=rhoOfFreq(rhoLiouville, freqIndex, Nlevels, Nfreq)
+	rho=zeros(Nlevels);
+	for r=1:Nlevels
+		for c=1:Nlevels
+			rho(r,c)=rhoLiouville((freqIndex-1)*Nlevels^2+(r-1)*Nlevels+c);
+		endfor
+	endfor
+endfunction
+
-- 
cgit v1.2.3