1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88
|
<simulation xmds-version="2">
<testing>
<command_line> mpirun -n 2 ./mpi_forward_plan_bug </command_line>
<xsil_file name="mpi_forward_plan_bug.xsil" expected="mpi_forward_plan_bug_expected.xsil" absolute_tolerance="1e-7" relative_tolerance="1e-5" />
</testing>
<name> mpi_forward_plan_bug </name>
<features>
<auto_vectorise />
<globals>
<![CDATA[
real N = 10.0; // number of atoms
real g = 1.0; // nonlinear coupling
]]>
</globals>
</features>
<geometry>
<propagation_dimension> t </propagation_dimension>
<transverse_dimensions>
<dimension name="x" lattice="256" domain="(-7, 7)" />
<dimension name="y" lattice="256" domain="(-7, 7)" />
</transverse_dimensions>
</geometry>
<driver name="distributed-mpi" />
<vector name="potential" type="real">
<components> V </components>
<initialisation>
<![CDATA[
V = 0.5*(x*x + y*y);
]]>
</initialisation>
</vector>
<vector name="wavefunction" type="complex">
<components> psi </components>
<initialisation>
<![CDATA[
psi = sqrt(N) * pow(M_PI, -0.25) * exp(-(x*x + y*y)/2);
]]>
</initialisation>
</vector>
<sequence>
<integrate algorithm="ARK45" interval="1" tolerance="1e-5">
<samples> 2 </samples>
<operators>
<operator kind="ip" dimensions="x">
<operator_names> Tx </operator_names>
<![CDATA[
Tx = -i * 0.5 * (kx * kx);
]]>
</operator>
<operator kind="ip" dimensions="y">
<operator_names> Ty </operator_names>
<![CDATA[
Ty = -i * 0.5 * (ky * ky);
]]>
</operator>
<integration_vectors>
wavefunction
</integration_vectors>
<dependencies> potential </dependencies>
<![CDATA[
dpsi_dt = Tx[psi] + Ty[psi] - i * (V + g * mod2(psi)) * psi;
]]>
</operators>
</integrate>
</sequence>
<output format="hdf5">
<sampling_group initial_sample="yes" basis="kx(8) ky(8)">
<dependencies> wavefunction </dependencies>
<moments> psireal psiimag </moments>
<![CDATA[
psireal = Re(psi);
psiimag = Im(psi);
]]>
</sampling_group>
</output>
</simulation>
|