cgaugenoisethetaphi.xmds

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

<simulation>

  <name>cgaugenoisethetaphi</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 kappa = 1;
    const double sqrtkappa = sqrt(kappa);
    const double impsi = -1.6;
    const double expimpsi = exp(impsi);
    const double expmimpsi = exp(-impsi);
  ]]>
  </globals>

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

  <sequence>
    <integrate>
      <algorithm>SIEX</algorithm>
      <interval>1.5</interval>
      <lattice>15000</lattice>
      <samples>1500</samples>
      <iterations>3</iterations>
      <![CDATA[
        dtheta_dt = 0.5*sqrtkappa*expimpsi*(n_1-i*n_2);
        dphi_dt = -kappa*exp(2*real(theta))*cos(2*imag(theta))+0.5*kappa-0.5*sqrtkappa*expmimpsi*(n_1+i*n_2);
        domega_dt = omega*(0.5*kappa*expimpsi*expimpsi*c_exp(2*conj(theta))+sqrtkappa*expimpsi*exp(2*real(theta))*sin(2*imag(theta))*(i*n_1+n_2));
      ]]>
    </integrate>
  </sequence>

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