test5.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 test5.cpp
 *
 * Rest test of NeuralBalloonModel.
 *
 * This test uses a completely vanilla NeuralBalloonModel to test the
 * behaviour of the model with no neural activity.
 *
 * \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 test5.png "Results"
 * \image latex test5.eps "Results"
 *
 * As the system is never perturbed by a change in neural activity all
 * properties should remain constant at their initial values.
 */


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

/**
 * Write out current state as row of results.
 */
void outputState(AdaptiveRungeKutta* ode, BalloonModel* bm,
    BOLDCalculator* bc) {
  cout << ode->getTime() << "\t";
  cout << bm->getFunction(FlowBalloonModel::F_IN) << "\t";
  cout << bm->getFunction(FlowBalloonModel::F_OUT) << "\t";
  cout << bm->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 << bc->calculate(ode->getState()) * 100.0 << "\t";
  cout << endl;
}

/**
 * Run tests.
 */
int main(int argc, const char* argv[]) {
  NeuralBalloonModel bm;
  BOLDCalculator bc(&bm);
  AdaptiveRungeKutta ode(&bm, bm.suggestInitialState());
  double t = 0.0;

  outputState(&ode, &bm, &bc);
  while (t < END) {
    t = ode.step(END);
    outputState(&ode, &bm, &bc);
  }

  return 0;
}

---