BlockDesignHarness.cpp

Go to the documentation of this file.

#include "indii/fmri/experimental/BlockDesignModel.hpp"

#include <iostream>
#include <fstream>

using namespace std;
using namespace indii::fmri::experimental;

/**
 * @file BlockDesignHarness.cpp
 *
 * Basic test of BlockDesignModel.
 *
 * This test sets up a BlockDesignModel and simulates it through a set
 * period. Results are output to the files
 * <tt>results/BlockDesignHarness_activityfirst.out</tt> and
 * <tt>results/BlockDesignHarness_restfirst.out</tt>, tab delimited,
 * with columns representing:
 *
 * \li \f$t\f$; time
 * \li \f$f(t)\f$; the function value
 *
 * Results are as follows:
 *
 * \image html BlockDesignHarness_activityfirst.png "Results, activity first"
 * \image latex BlockDesignHarness_activityfirst.eps "Results, activity first"
 *
 * \image html BlockDesignHarness_restfirst.png "Results, rest first"
 * \image latex BlockDesignHarness_restfirst.eps "Results, rest first"
 */

const unsigned int BLOCK_LENGTH = 10;
const double REST_LEVEL = 0.0;
const double ACTIVITY_LEVEL = 1.0;

const double LENGTH = 60.0;
const double STEP = 2.0;

/**
 * Run tests.
 */
int main(int argc, const char* argv[]) {
  BlockDesignModel restFirst(BLOCK_LENGTH, true, REST_LEVEL, ACTIVITY_LEVEL);
  BlockDesignModel activityFirst(BLOCK_LENGTH, false, REST_LEVEL,
      ACTIVITY_LEVEL);
  double t;

  ofstream fActivityFirst("results/BlockDesignHarness_activityfirst.out");
  ofstream fRestFirst("results/BlockDesignHarness_restfirst.out");

  for (t = 0.0; t <= LENGTH; t += STEP) {
    fActivityFirst << t << '\t' << activityFirst.evaluate(t, NULL) << endl;
    fRestFirst << t << '\t' << restFirst.evaluate(t, NULL) << endl;
  }

  fActivityFirst.close();
  fRestFirst.close();

  return 0;
}

---