melica.h

/*  MELODIC - Multivariate exploratory linear optimized decomposition into 
              independent components
    
    melica.h - ICA estimation

    Christian F. Beckmann, FMRIB Image Analysis Group
    
    Copyright (C) 1999-2007 University of Oxford */

/*  CCOPYRIGHT */

#ifndef __MELODICICA_h
#define __MELODICICA_h
#include "newimage/newimageall.h"
#include "utils/log.h"
#include "melpca.h"
#include "meloptions.h"
#include "meldata.h"
#include "melodic.h"
#include "newmatap.h"
#include "newmatio.h"
#include "melreport.h"
#include "ggmix.h"

using namespace Utilities;
using namespace NEWIMAGE;

namespace Melodic{
  
  class MelodicICA{
    public:
      MelodicICA(MelodicData &pmelodat, MelodicOptions &popts, Log &plogger, MelodicReport &preport):  
				melodat(pmelodat),
				opts(popts),
				logger(plogger),
				report(preport){} 
      
      void perf_ica(const Matrix &Data);
      
      bool no_convergence;

    private:
      MelodicData &melodat;
      MelodicOptions &opts;
      Log &logger;
      MelodicReport &report;

      int dim;
      int samples; 
      Matrix redUMM;

      void ica_fastica_symm(const Matrix &Data);
      void ica_fastica_defl(const Matrix &Data);
      void ica_maxent(const Matrix &Data);
      void ica_jade(const Matrix &Data);
      //void tica();
      //void tica_concat();
      //void parafac();
      Matrix randm(const int dim1, const int dim2);
      void sort();
      Matrix sign(const Matrix &Inp);

  };   
}

#endif