thetaphinoise.xmds

<?xml version="1.0"?>
<!--diffusion and complex drift gauge single mode bose field simulation-->

<simulation>

  <name>thetaphinoise</name>

  <author>Mark Dowling</author>
  <description>
    diffusion and complex drift gauge single mode bose field simulation
  </description>

  <prop_dim>t</prop_dim>
  <error_check>yes</error_check>
  <stochastic>yes</stochastic>
  <paths>1000</paths>
  <use_mpi>no</use_mpi>
  <seed>1 2</seed>
  <noises>2</noises>
  
  <globals>
  <![CDATA[
    const double pi = 3.141592653589;
    const double kappa = 2*pi;
    const double sqrtkappa = sqrt(kappa);
    const double impsi = -1;
    const double expimpsi = exp(impsi);
    const double expmimpsi = exp(-impsi);
    const double lambda = expimpsi;
  ]]>
  </globals>

  <field>
    <name>main</name>
    <samples>1</samples>
    <vector>
      <name>main</name>
      <type>complex</type>
      <components>phi theta </components>
      <![CDATA[
        theta = 0.5*log(1.0);
        phi = 0;
      ]]>
    </vector>
  </field>

  <sequence>
    <integrate>
      <algorithm>SIEX</algorithm>
      <interval>1</interval>
      <lattice>10000</lattice>
      <samples>1000</samples>
      <iterations>3</iterations>
      <![CDATA[
        dtheta_dt = 0.5*sqrtkappa*expimpsi*(n_1-i*n_2);
        dphi_dt = -kappa*c_exp(2*theta)+0.5*kappa-0.5*sqrtkappa*expmimpsi*(n_1+i*n_2);
      ]]>
    </integrate>
  </sequence>

  <output>
    <group>
      <sampling>
        <moments>no var X Y rea2 ima2 retheta imtheta rephi imphi retheta2 imtheta2 rephi2 imphi2 </moments>
        <![CDATA[
          no = c_exp(2*theta);
          var = c_exp(4*theta);
          X = real(c_exp(theta+i*phi) + conj(c_exp(theta-i*phi)))/2;
          Y = imag(c_exp(theta+i*phi) + conj(c_exp(theta-i*phi)))/2;
          rea2 = real(c_exp(2*(theta+i*phi))+conj(c_exp(2*(theta-i*phi))))/2; 
          ima2 = imag(c_exp(2*(theta+i*phi))+conj(c_exp(2*(theta-i*phi))))/2; 
          retheta = real(theta);
          imtheta = imag(theta);
          rephi = real(phi);
          imphi = imag(phi);
          retheta2 = real(theta)*real(theta);
          imtheta2 = imag(theta)*real(theta);
          rephi2 = real(phi)*real(phi);
          imphi2 = imag(phi)*imag(phi);
        ]]>
      </sampling>
    </group>
  </output>
</simulation>