eigenstructure.hpp

Go to the documentation of this file.

#ifndef EXPLICIT_MAPMAKING_EIGENSTRUCTURE_HPP_
#define EXPLICIT_MAPMAKING_EIGENSTRUCTURE_HPP_

#include "grid.hpp"
#include "matrix.hpp"
#include <stdio.h>
#include <boost/filesystem.hpp>

namespace fs = boost::filesystem;

namespace mapmaking {
class Eigenstructure {

 public:
  Eigenstructure(Grid* grid, fs::path folder);
  Eigenstructure(Grid* grid, int n_row, int n_col) : 
      grid_(grid), vectors_(grid, n_row, n_col), values_(grid, n_row, 1) {}
  Matrix ModesAmplitude(Matrix& map);
  void WeightModes(Matrix* map, std::function<double(double mode_amplitude, double eigenvalue)> weighting_function);
  void ThresholdModes(Matrix* map, double relative_threshold);
  void InvertEigenvalues();
  Matrix WhiteNoiseMap(double relative_threshold);
  void Save(fs::path folder, const std::vector<int>& obspix);
  void Save(fs::path folder, fs::path obspix_file);

 protected:
  Matrix vectors_; 
  Matrix values_; 
  Grid* grid_; 
  double max_eigenvalue_ = 0;
  std::vector<double>* vectors() { return vectors_.local_part(); }
  void SetEigenvalues(std::vector<double>& values, int root_proc);
  void SetMaxEigenvalue();

 private:
  friend class AtFA;
};
}

#endif  // EXPLICIT_MAPMAKING_EIGENSTRUCTURE_HPP_