From c3154f49eba151466cbf933ae65a37eeaead7c16 Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Wed, 11 Jan 2012 23:11:48 -0500
Subject: 4wm xmds perturbative: length of the cell can be set

use --Lcel  at runtime
---
 .../Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds          | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds
index 4fc9086..43ba51a 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds
@@ -83,6 +83,8 @@
       <argument name="delta1"  type="real" default_value="0.0" />
       <argument name="delta2"  type="real" default_value="0.0" />
       <argument name="delta3"  type="real" default_value="0.0" />
+			<!--Cell length-->
+      <argument name="Lcell"  type="real" default_value="1.5e-2" />
     </arguments>
 		<bing />
 		<fftw plan="patient" />
@@ -202,7 +204,11 @@
 
 	<sequence>
 		<!--For this set of conditions ARK45 is faster than ARK89-->
-		<integrate algorithm="ARK45" tolerance="1e-5" interval="1.5e-2">
+		<!--ARK45 is good for small detuning when all frequency like term are close to zero-->
+		<integrate algorithm="ARK45" tolerance="1e-5" interval="Lcell"> 
+		<!--RK4 is good for large detunings when frequency like term are big, it does not try to be too smart about adaptive step which ARK seems to make too small-->
+		<!--When ARK45 works it about 3 times faster then RK4 with 1000 steps-->
+		<!--<integrate algorithm="RK4" steps="100" interval="1.5e-2">-->
 		<!--SIC algorithm seems to be much slower and needs fine 'z'  step tuning and much finer time grid-->
 		<!--For example I had to quadruple the time grid from 1000 to 4000 when increased z distance from 0.02 to 0.04-->
 
-- 
cgit v1.2.3