fibre.xmds

<?xml version="1.0"?>
<!--Example fibre noise simulation-->

<simulation>

  <name>fibre</name>
  <prop_dim>t</prop_dim>
  <error_check>yes</error_check>
  <stochastic>yes</stochastic>
  <use_mpi>no</use_mpi>
  <paths>1</paths>
  <seed>1 2</seed>
  <noises>2</noises>

  <globals>
    <![CDATA[
    const double ggamma = 1;
    const double beta = sqrt(2*2*M_PI*ggamma/10);
    ]]>
  </globals>

  <field>
    <name>main</name>
    <dimensions>   x   </dimensions>
    <lattice>     50  </lattice>
    <domains>  (-5,5)  </domains>
    <samples>1</samples>

    <vector>
      <name>main</name>
      <type>complex</type>
      <components>phi</components>
      <fourier_space>no</fourier_space>
      <![CDATA[
      phi=0;
      ]]>
    </vector>
  </field>

  <sequence>
    <integrate>
      <algorithm>SIIP</algorithm>
      <interval>2.5</interval>
      <lattice>5000</lattice>
      <samples>50</samples>
      <k_operators>
	<constant>yes</constant>
	<operator_names>L</operator_names>
	<![CDATA[
	L = i*(-kx*kx);
	]]>
      </k_operators>
      <iterations>3</iterations>
      <![CDATA[
      dphi_dt = L[phi] - ggamma*phi + beta/sqrt(2)*complex(n_1,n_2);
      ]]>
    </integrate>
  </sequence>

  <output>
    <group>
      <sampling>
	<fourier_space> yes </fourier_space>
	<lattice>       50  </lattice>
	<moments>pow_dens</moments>
	<![CDATA[
	pow_dens = conj(phi)*phi;
	]]>
      </sampling>
    </group>
  </output>
</simulation>