atfa.hpp

Go to the documentation of this file.

#ifndef EXPLICIT_MAPMAKING_ATFA_HPP
#define EXPLICIT_MAPMAKING_ATFA_HPP

#include "matrix.hpp"
#include "eigenstructure.hpp"
#include "data_manager.hpp"

//TODO AtFa should inherit from a SquaredMatrix class and Invert should be moved
//therein

namespace mapmaking {
class AtFA : public Matrix {
 public:
  AtFA() {}
  AtFA(const Matrix& matrix) : Matrix(matrix) {}
  AtFA(Grid* grid, DataManager* data_manager, int matrix_side) : data_manager_(data_manager),
    Matrix(grid, matrix_side, matrix_side) {}

  void ComputeAtfA();
  Eigenstructure Invert(double cn_threshold=1e-6);
  //AtFA& operator=(const Matrix& matrix);

 private:
  //TODO getter with NULL check or (better) move invert in a SquaredMatrix class
  DataManager* data_manager_ = NULL;  
  virtual int NStokes() = 0;
  virtual void UpdateAtFAPart(const double& contrib, const int& n_row_this_block,
                              const CES& ces, const int t1, const int t2,
                              int& index, std::vector<double>& afa_part) = 0;
};
}

#endif  // EXPLICIT_MAPMAKING_ATFA_HPP