indii/fmri/neural/MultiSeedNeuralModel.hpp

#ifndef INDII_ML_ODE_MULTISEEDNEURALMODEL_HPP
#define INDII_ML_ODE_MULTISEEDNEURALMODEL_HPP

#include "indii/ml/aux/vector.hpp"
#include "indii/ml/aux/matrix.hpp"
#include "indii/ml/ode/DifferentialModel.hpp"

namespace indii {
  namespace fmri {
    namespace neural {

/**
 * Multi-seed differential %model.
 *
 * @author Lawrence Murray <lawrence@indii.org>
 * @version $Rev: 285 $
 * @date $Date: 2007-07-20 17:25:40 +0100 (Fri, 20 Jul 2007) $
 *
 * Neural activities of \f$N\f$ seeds \f$\mathbf{x}\f$ are related
 * according to the differential system:
 *
 * \f[
 *   \frac{d\mathbf{x}}{dt} = A\mathbf{x}
 * \f]
 *
 * where \f$A\f$ is an \f$N \times N\f$ matrix of efficacies between
 * seeds.
 */
class MultiSeedNeuralModel : public indii::ml::ode::DifferentialModel {
public:
  /**
   * Construct new %model with given number of seeds. The efficacy
   * matrix is uninitialised.
   *
   * @param N \f$N\f$; the number of seeds.
   */
  MultiSeedNeuralModel(const unsigned int N);

  /**
   * Construct new %model with given efficacy matrix.
   *
   * @param A \f$A\f$; efficacies between seeds, assumed to be \f$N
   * \times N\f$ in size.
   */
  MultiSeedNeuralModel(const indii::ml::aux::matrix& A);

  /**
   * Destructor.
   */
  virtual ~MultiSeedNeuralModel();

  /**
   * Set all efficacies.
   *
   * @param A \f$A\f$; efficacies between seeds, must be \f$N \times
   * N\f$ in size.
   */
  void setEfficacies(const indii::ml::aux::matrix& A);

  /**
   * Set the efficacy between two seeds.
   *
   * @param i \f$i\f$; index of the affected seed.
   * @param j \f$j\f$; index of the affecting seed.
   * @param A_ij \f$A_{ij}\f$; efficacy of seed \f$j\f$ on seed
   * \f$i\f$.
   */
  void setEfficacy(const unsigned int i, const unsigned int j,
      const double A_ij);

  /**
   * @see indii::ml::ode::DifferentialModel
   */
  virtual void calculateDerivatives(double t, const double y[], double dydt[]);

private:
  /**
   * \f$N\f$n; the number of seeds in the model.
   */
  const unsigned int N;

  /**
   * \f$A\f$; the efficacy matrix.
   */
  indii::ml::aux::matrix A;

};

    }
  }
}

#endif

---