diff options
Diffstat (limited to 'xmds2')
| -rw-r--r-- | xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z.xmds | 16 |
1 files changed, 11 insertions, 5 deletions
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z.xmds index a4d2020..2182fe0 100644 --- a/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z.xmds +++ b/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z.xmds @@ -68,7 +68,7 @@ complex E1c, E2c, E3c, E4c; // Complex conjugated Rabi frequencies - complex r11, r13, r22, r24, r33, r44; // pre-calculated density matrix elements + //complex r11, r13, r22, r24, r33, r44; // pre-calculated density matrix elements complex r21, r31, r41, r32, r42, r43; // density matrix elements ]]> </globals> @@ -113,12 +113,17 @@ </initialisation> </vector> - <vector name="density_matrix_zero" type="complex" initial_space="t"> + <vector name="zero_approximation_rho" type="complex" initial_space="t"> <!--Zero order density matrix time and z independent--> <components>r11 r22 r33 r13 r24 r44</components> <initialisation> <dependencies>E_field</dependencies> <![CDATA[ + E1c = conj(E1); + E2c = conj(E2); + E3c = conj(E3); + E4c = conj(E4); + r11 = (gt*(4*mod2(E1) + (G3 + gt)*(G3 + 2*gt))*((G4 + 2*gt)*(4*mod2(E3) + gt*(G4 + gt)) + 4*mod2(E3)*G4*(R41 - R42)))/(2.*(-16*mod2(E1)*mod2(E3)*G3*G4*R32*R41 + (G3 + @@ -202,7 +207,7 @@ <operators> <operator kind="cross_propagation" algorithm="SI" propagation_dimension="t"> <integration_vectors>density_matrix</integration_vectors> - <dependencies>E_field density_matrix_zero</dependencies> + <dependencies>E_field zero_approximation_rho</dependencies> <boundary_condition kind="left"> <!-- <![CDATA[ @@ -241,12 +246,13 @@ ]]> </operator> <integration_vectors>E_field</integration_vectors> - <dependencies>density_matrix</dependencies> + <dependencies>density_matrix zero_approximation_rho</dependencies> <![CDATA[ dE1_dz = i*eta*conj(r13) -Lt[E1] ; dE2_dz = i*eta*conj(r23) -Lt[E2] ; dE3_dz = i*eta*conj(r24) -Lt[E3] ; dE4_dz = i*eta*conj(r14) -Lt[E4] ; + // printf("z= %g, t= %g, r13 = %g\n", z, t, mod2(r13)); //dbg ]]> </operators> </integrate> @@ -272,7 +278,7 @@ <group> <sampling basis="t(100)" initial_sample="yes"> - <dependencies>density_matrix</dependencies> + <dependencies>density_matrix zero_approximation_rho</dependencies> <moments> r11_out r22_out r33_out r44_out r12_re_out r12_im_out r13_re_out r13_im_out r14_re_out r14_im_out |