diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2012-01-11 23:11:48 -0500 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2012-01-11 23:11:48 -0500 |
| commit | c3154f49eba151466cbf933ae65a37eeaead7c16 (patch) | |
| tree | 91f55421f30aed5f95368adf4819786b2ae96f25 | |
| parent | daf766c19865d93d5e74a1e46f99b37a2f10ab86 (diff) | |
| download | Nresonances-c3154f49eba151466cbf933ae65a37eeaead7c16.tar.gz Nresonances-c3154f49eba151466cbf933ae65a37eeaead7c16.zip | |
4wm xmds perturbative: length of the cell can be set
use --Lcel at runtime
| -rw-r--r-- | xmds2/Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.xmds | 8 |
1 files changed, 7 insertions, 1 deletions
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--> |