test4.cpp

Go to the documentation of this file.

#include "indii/fmri/hemodynamic/NeuralBalloonModel.hpp"
#include "indii/fmri/hemodynamic/FlowBalloonModel.hpp"
#include "indii/fmri/hemodynamic/BOLDCalculator.hpp"
#include "indii/ml/ode/AdaptiveRungeKutta.hpp"

#include <iostream>
#include <string.h>
#include <math.h>

using namespace std;
using namespace indii::fmri::hemodynamic;
using namespace indii::ml::ode;

/**
 * @file test4.cpp
 *
 * Recurrent test of NeuralBalloonModel.
 *
 * This test uses a default NeuralBalloonModel with a driving neural
 * activity boxcar function.
 *
 * \li \f$t\f$; time
 * \li \f$f_{in}(t)\f$
 * \li \f$f_{out}(t)\f$
 * \li \f$u(t)\f$
 * \li \f$q\f$
 * \li \f$v\f$
 * \li \f$f\f$
 * \li \f$s\f$
 * \li \f$y\f$; BOLD response
 *
 * Results are as follows:
 *
 * \image html test4.png "Results"
 * \image latex test4.eps "Results"
 *
 * Note in particular that the curve for the state variable \f$f\f$
 * should match that for the function \f$f_{in}(t)\f$.
 */


/**
 * Duration of simulation (s).
 */
static const double END = 60.0;

/**
 * \f$u(t)\f$; 1 for 1 s, then 0.
 */
double U(double t, const double y[], void* o) {
  double u;

  int block = (int)floor(t / 10.0);
  if (block >> 1 << 1 == block) {
    // even number
    u = 1.0;
  } else {
    u = 0.0;
  }
  return u;
}

/**
 * Run tests.
 */
int main(int argc, const char* argv[]) {
  NeuralBalloonModel balloonModel;
  BOLDCalculator boldCalculator(&balloonModel);
  AdaptiveRungeKutta ode(&balloonModel, balloonModel.suggestInitialState());

  /* set up balloon model */
  balloonModel.setFunction(NeuralBalloonModel::U, U);

  double t = 0.0;
  while (t < END) {
    cout << t << "\t";
    cout << balloonModel.getFunction(FlowBalloonModel::F_IN) << "\t";
    cout << balloonModel.getFunction(FlowBalloonModel::F_OUT) << "\t";
    cout << balloonModel.getFunction(NeuralBalloonModel::U) << "\t";
    cout << ode.getVariable(FlowBalloonModel::Q) << "\t";
    cout << ode.getVariable(FlowBalloonModel::V) << "\t";
    cout << ode.getVariable(NeuralBalloonModel::F) << "\t";
    cout << ode.getVariable(NeuralBalloonModel::S) << "\t";
    cout << boldCalculator.calculate(ode.getState()) * 100.0 << "\t";
    cout << endl;
    t = ode.step(END);
  }

  return 0;
}

---