soliton.xmds

<?xml version="1.0"?>
<!--Example Exactly Solvable Coupled NLSE-->

<simulation>

  <name>soliton</name>
  <prop_dim>z</prop_dim>
  <error_check>yes</error_check>
  <stochastic>no</stochastic>

  <globals>
    <![CDATA[
    const double alpha = 1;
    const double beta = 1;
    ]]>
  </globals>

  <field>
    <name> main </name>
    <dimensions> t    </dimensions>
    <lattice>    50   </lattice>
    <domains>  (-5,5) </domains>
    <samples> 1 </samples>
    <vector>
      <name> main </name>
      <type> complex </type>
      <components>u v</components>
      <fourier_space>no</fourier_space>
      <![CDATA[
      u = complex(exp(-t*t/alpha/alpha/4),0);
      v = complex(exp(-t*t/beta/beta/4),0);
      ]]>
    </vector>
  </field>

  <sequence>
    <integrate>
      <algorithm>RK4IP</algorithm>
      <interval>10</interval>
      <lattice>2400</lattice>
      <samples>30</samples>

      <k_operators>
	<constant>yes</constant>
	<operator_names>L</operator_names>
	<![CDATA[
	L = i*(-kt*kt/2);
	]]>
      </k_operators>

      <iterations>3</iterations>
      <![CDATA[
      const double density = u.re*u.re+u.im*u.im+v.re*v.re+v.im*v.im;

      du_dz = L[u] + i*u*density;
      dv_dz = L[v] + i*v*density;
      ]]>
    </integrate>
  </sequence>

  <output>
    <group>
      <sampling>
	<fourier_space> no </fourier_space>
	<lattice>       25  </lattice>
	<moments>ure uim vre vim</moments>
	<![CDATA[
	ure = u;
	uim = -i*u;
	vre = v;
	vim = -i*v;
	]]>
      </sampling>
    </group>
  </output>
</simulation>