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