atfa_pol.hpp

Go to the documentation of this file.

#ifndef EXPLICIT_MAPMAKING_ATFA_POL_HPP
#define EXPLICIT_MAPMAKING_ATFA_POL_HPP

#include "atfa.hpp"

namespace mapmaking {
class AtFAPol : public AtFA {
 public:
  AtFAPol(const Matrix& matrix) : AtFA(matrix) {}
  AtFAPol(Grid* grid, DataManager* data_manager) : 
    AtFA(grid, data_manager, data_manager->n_pixel()*2) {}

 private:
  inline int NStokes() { return 2; }
  inline 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) {
    //QQ
    afa_part[index] += contrib * ces.cos(t1) * ces.cos(t2);
    //UQ
    index++;
    afa_part[index] += contrib * ces.sin(t1) * ces.cos(t2);
    //UU
    index += n_row_this_block;
    afa_part[index] += contrib * ces.sin(t1) * ces.sin(t2);
    //QU
    index--;
    afa_part[index] += contrib * ces.cos(t1) * ces.sin(t2);
  }
};
}

#endif  // EXPLICIT_MAPMAKING_ATFA_POL_HPP