diff options
Diffstat (limited to 'MOR_simple_thin_layer_model')
| -rw-r--r-- | MOR_simple_thin_layer_model/Makefile | 37 | ||||
| -rw-r--r-- | MOR_simple_thin_layer_model/delta_2photon_noise.xmds | 148 | ||||
| -rw-r--r-- | MOR_simple_thin_layer_model/pp.m | 21 |
3 files changed, 206 insertions, 0 deletions
diff --git a/MOR_simple_thin_layer_model/Makefile b/MOR_simple_thin_layer_model/Makefile new file mode 100644 index 0000000..83877ee --- /dev/null +++ b/MOR_simple_thin_layer_model/Makefile @@ -0,0 +1,37 @@ +### -*- make -*- +### This file is part of the Debian xmds package +### Copyright (C) 2006 Rafael Laboissiere +### This file is relased under the GNU General Public License +### NO WARRANTIES! + +### This makefile can be used to build and run the XMDS examples + +XMDS_FILES = $(shell ls *.xmds) +RUN_FILES = $(patsubst %.xmds,%.run,$(XMDS_FILES)) +CC_FILES = $(patsubst %.xmds,%.cc,$(XMDS_FILES)) +XSIL_FILES = delta_2photon_noise.xsil +M_FILES = $(patsubst %.xsil,%.m,$(XSIL_FILES)) + +XMDS = xmds2 +XSIL2GRAPHICS = xsil2graphics2 + +all: $(M_FILES) + +%.run: %.xmds + $(XMDS) $< + mv $(patsubst %.xmds,%,$<) $@ + +$(XSIL_FILES): %.xsil: $(RUN_FILES) + ./$< + +$(M_FILES): %.m : %.xsil + $(XSIL2GRAPHICS) $< + +plot: $(M_FILES) + octave pp.m + +clean: + rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf + +.PRECIOUS: %.run %.xsil %.m +.PHONY: all clean diff --git a/MOR_simple_thin_layer_model/delta_2photon_noise.xmds b/MOR_simple_thin_layer_model/delta_2photon_noise.xmds new file mode 100644 index 0000000..2cb3a52 --- /dev/null +++ b/MOR_simple_thin_layer_model/delta_2photon_noise.xmds @@ -0,0 +1,148 @@ +<?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: +--> + diff --git a/MOR_simple_thin_layer_model/pp.m b/MOR_simple_thin_layer_model/pp.m new file mode 100644 index 0000000..f59f19c --- /dev/null +++ b/MOR_simple_thin_layer_model/pp.m @@ -0,0 +1,21 @@ +delta_2photon_noise + +figure(1) +plot(delta_2photon_1, mean_mor_signal_1(end,:) ) +hold all +plot(delta_2photon_1, mean_eit_signal_1(end,:) ) + +legend( {'mor', 'eit'} ) + +figure(2) +plot(delta_2photon_1, stderr_mor_signal_1(end,:) ) +hold all +plot(delta_2photon_1, stderr_eit_signal_1(end,:) ) + +legend( {'mor', 'eit'} ) + +figure(3) +plot(delta_2photon_1, stderr_eit_signal_1(end,:)./stderr_mor_signal_1(end,:) ) +hold all + +legend( {'noise eit/mor'} ) |