full3d.xmds

<?xml version="1.0"?>
<!-- Simulates evolution of coupled atomic-molecular BECs on full 3D lattice -->

<simulation>

  <name>full3d</name>

  <author>Tim Vaughan</author>
  <description>
    Simulates evolution of coupled atomic-molecular BECs on full 3D
    lattice
  </description>

  <prop_dim>t</prop_dim>
  <error_check>no</error_check>
  <stochastic>no</stochastic>

  <globals>
  <![CDATA[
    const double alpha_11 = 0.1;
    const double gam = 1;

    const double A = 3.7919;
    const double B = 4.4037;
    const double a = 0.2756;
    const double b = 0.2901;

    const double R = 10;
    const double AG = 1;
    const double PG = log(0.0363206335777221)/(log(sin((1-5/R)*M_PI/2)));
  ]]>
  </globals>

  <field>
    <name>main</name>
    <dimensions>x y z</dimensions>
    <lattice>100 100 100</lattice>
    <domains>(-10,10) (-10,10) (-10,10)</domains>
    <samples>1 1</samples>
    <vector>
      <name>main</name>
      <type>complex</type>
      <components>u v</components>
      <fourier_space>no no no</fourier_space>
      <![CDATA[
        double r2 = x*x + y*y + z*z;
        u = complex(A*exp(-a*r2),0);
        v = complex(B*exp(-b*r2),0);
    ]]>
    </vector>
    <vector>
      <name>vc1</name>
      <type>complex</type>
      <components>damp</components>
      <fourier_space>no no no</fourier_space>
      <![CDATA[
        double r = sqrt(x*x + y*y + z*z);
        damp = complex(AG*pow(sin(r/R*0.5*M_PI),PG),0);
      ]]>
    </vector>
  </field>
 
 <sequence>
    <integrate>
      <algorithm>SIIP</algorithm>
      <interval>10</interval>
      <lattice>100</lattice>
      <samples>100 100</samples>
      <k_operators>
        <constant>yes</constant>
        <operator_names>L1 L2</operator_names>
        <![CDATA[
          L1 = -i*(kx*kx + ky*ky + kz*kz);
          L2 = -i*0.5*(kx*kx + ky*ky + kz*kz);
        ]]>
      </k_operators>
      <iterations>3</iterations>
      <vectors>main vc1</vectors>
      <![CDATA[
        du_dt = L1[u] - i*(u - conj(u)*v + alpha_11*conj(u)*u*u) - u*damp;
        dv_dt = L2[v] - 0.5*i*(gam*v - u*u) - 0.5*v*damp;
      ]]>
    </integrate>
  </sequence>

  <output>
    <group>
      <sampling>
        <fourier_space>no no no</fourier_space>
        <lattice>100 1 1</lattice>
        <moments>Iu Iv</moments>
        <![CDATA[
          Iu = mod2(u);
          Iv = mod2(v);
        ]]>
      </sampling>
    </group>
    <group>
      <sampling>
        <fourier_space>no no no</fourier_space>
        <lattice>0 0 0</lattice>
        <moments>Na Nb N</moments>
        <![CDATA[
          Na = mod2(u);
          Nb = mod2(v);
          N = Na + 2*Nb;
        ]]>
      </sampling>
    </group>
  </output>
</simulation>