TimeDomainCovariance.h

/***************************************************************************
 *  *
 *  *   Copyright (C) 2011 by Stefan Oslowski
 *  *   Licensed under the Academic Free License version 2.1
 *  *
 ****************************************************************************/
 
#ifndef __Pulsar__TimeDomainCovariance_h
#define __Pulsar__TimeDomainCovariance_h
 
#include "Pulsar/ProfileCovariance.h"
 
namespace Pulsar
{
  class TimeDomainCovariance : public ProfileCovariance
  {
  public:
    
    TimeDomainCovariance ();
 
    void choose_bins ( unsigned, unsigned );
 
    std::string get_name () const;
 
    void add_Profile (const Profile* );
    void add_Profile (const Profile*, float wt );
 
    void set_covariance_matrix ( const double * );
 
    void get_covariance_matrix_copy ( double* dest );
      
    double get_covariance_matrix_value ( unsigned row, unsigned col );
 
    void finalize ();
 
    void eigen ();
 
    void set_eigenvectors ( const std::vector<double>& );
    void set_eigenvectors ( const double* );
 
    void get_eigenvectors_copy ( std::vector<double>& );
    void get_eigenvectors_copy ( double* );
    
    double get_eigenvectors_value ( unsigned row, unsigned col );
    const double* get_eigenvectors_pointer ();
 
    void set_eigenvalues ( const std::vector<double>& );
    void set_eigenvalues ( const double* );
 
    void get_eigenvalues_copy ( std::vector<double>& );
    void get_eigenvalues_copy ( double* );
    
    double get_eigenvalue ( unsigned row );
    const double* get_eigenvalues_pointer ();
    
    unsigned get_count ();
 
    void set_count ( unsigned );
 
    unsigned get_rank ();
 
    void set_rank ( unsigned );
 
    void reset ();
 
  protected:
    
    unsigned first_bin;
 
    unsigned last_bin;
 
    bool subtract_mean;
 
    std::vector<double> covariance_matrix;
    
    std::vector<double> mean;
 
    std::vector<double> eigenvectors;
    
    std::vector<double> eigenvalues;
 
    bool eigen_decomposed;
 
  };
}
 
#endif