more descriptive name for the simple MOR directory

1 files changed, 0 insertions, 148 deletions

diff --git a/simple_thin_layer_model/delta_2photon_noise.xmds b/simple_thin_layer_model/delta_2photon_noise.xmds
deleted file mode 100644
index 2cb3a52..0000000
--- a/simple_thin_layer_model/delta_2photon_noise.xmds
+++ /dev/null
@@ -1,148 +0,0 @@
-<?xml version="1.0"?>
-<simulation xmds-version="2">
-
-	<name>delta_2photon_noise</name>
-
-	<author>E. Mikhailov, M. Guidry</author>
-	<description>
-		Solving 3 level atom in lambda field configuration,
-		no field propagation included. Noise on delta_2photon.
-	</description>
-
-	<features>
-		<benchmark />
-		<bing />
-		<fftw plan="patient" />
-		<openmp />
-		<auto_vectorise />
-		<globals>
-			<![CDATA[
-			const complex g=10;
-			const complex gbc=.001;
-			const complex Split = .1;
-			const double tset=5;
-			const double alpha = 0.001;
-
-			complex Gab, Gca, Gcb, GCB;
-			complex Ep,Ed,Epc,Edc;
-			complex noise;
-			double delta_2photon, delta_d;
-			complex rba,rac,rbc;
-			complex rBA,rAC,rBC;
-			complex drive_field0=1, probe_field0=1;
-			]]>
-		</globals>
-	</features>
-
-	<geometry>
-		<propagation_dimension> t </propagation_dimension>
-		<transverse_dimensions>
-			<dimension name="delta_2photon" lattice="256"  domain="(-2, 2)" />
-		</transverse_dimensions>
-	</geometry>
-
-	<driver name="multi-path" paths="16" />
-
-	<noise_vector name="drivingNoise" kind="wiener" type="real" method="dsfmt" >
-    		<components>Eta</components>
-  	</noise_vector>
-
-	<vector name="density_matrix" type="complex" dimensions="delta_2photon">
-		<components>raa rbb rcc rab rca rcb rAB rCA rCB</components>
-		<initialisation>
-			<![CDATA[
-			raa = 0;
-			rbb = 1;
-			rcc = 0;
-			rab = 0;
-			rca = 0;
-			rcb= 0;
-			rAB = 0;
-			rCA = 0;
-			rCB= 0;
-			]]>
-		</initialisation>
-	</vector>
-
-	<sequence>
-		<integrate algorithm="ARK45" interval="13.0" tolerance="1e-7">
-			<samples>100</samples>
-			<operators>
-				<integration_vectors>density_matrix</integration_vectors>
-				<dependencies>drivingNoise</dependencies>
-				<![CDATA[
-				delta_d=0;
-
-				noise = Eta*alpha;
-				Ed=drive_field0;
-				Ep=probe_field0;
-
-				Epc = conj(Ep);
-				Edc = conj(Ed);
-
-				Gab=g+i*(delta_d+delta_2photon+1*noise);
-				Gca=g-i*(delta_d);
-				Gcb=gbc+i*(Split + delta_2photon+0*noise);
-				GCB=gbc+i*(-Split + delta_2photon+0*noise);
-
-
-				rba=conj(rab);
-				rac=conj(rca);
-				rbc=conj(rcb);
-				rBA=conj(rAB);
-				rAC=conj(rCA);
-				rBC=conj(rCB);
-
-				draa_dt = -i*Epc*rab+i*Ep*rba-i*Edc*rac+i*Ed*rca-2*g*raa
-				          -i*Edc*rAB+i*Ed*rBA-i*Epc*rAC+i*Ep*rCA;
-				drbb_dt =  i*Epc*rab-i*Ep*rba+g*raa-gbc*rbb+gbc*rcc
-				          +i*Edc*rAB-i*Ed*rBA;
-				drcc_dt =  i*Edc*rac-i*Ed*rca+g*raa-gbc*rcc+gbc*rbb
-				          +i*Epc*rAC-i*Ep*rCA;
-
-				drab_dt = -Gab*rab+i*Ep*(rbb-raa)+i*Ed*rcb;
-				drca_dt = -Gca*rca+i*Edc*(raa-rcc)-i*Epc*rcb;
-				drcb_dt = -Gcb*rcb-i*Ep*rca+i*Edc*rab;
-
-				drAB_dt = -Gab*rAB+i*Ed*(rbb-raa)+i*Ep*rCB;
-				drCA_dt = -Gca*rCA+i*Epc*(raa-rcc)-i*Edc*rCB;
-				drCB_dt = -GCB*rCB-i*Ed*rCA+i*Epc*rAB;
-					]]>
-				</operators>
-			</integrate>
-		</sequence>
-
-
-
-
-		<!-- The output to generate -->
-		<!--<output format="ascii" filename="delta_2photon_noise_alpha_1.xsil">-->
-		<output>
-			<group>
-				<sampling basis="delta_2photon" initial_sample="yes">
-					<dependencies>density_matrix </dependencies>
-					<moments>eit_signal mor_signal rab_rlOut rab_imOut rca_rlOut rca_imOut raaOut rbbOut rccOut noiseOut PDl PDr</moments>
-					<![CDATA[
-					rab_rlOut = rab.Re();
-					rab_imOut = rab.Im();
-					rca_rlOut = rca.Re();
-					rca_imOut = rca.Im();
-					raaOut = raa.Re();
-					rbbOut = rbb.Re();
-					rccOut = rcc.Re();
-					noiseOut = noise.Re();
-					PDl = norm( 1/sqrt(2)*(rab - rAB) + i/sqrt(2)*(rca + rCA) );
-					PDr = norm( 1/sqrt(2)*(rab - rAB) - i/sqrt(2)*(rca + rCA) );
-					eit_signal = PDl + PDr;
-					mor_signal = PDl - PDr;
-					]]>
-				</sampling>
-			</group>
-		</output>
-
-	</simulation>
-	
-<!-- 
-vim: ts=2 sw=2 foldmethod=indent: 
--->
-