DiagonalSolver.h

Go to the documentation of this file.

// Copyright (C) 2021 ASTRON (Netherlands Institute for Radio Astronomy)
// SPDX-License-Identifier: GPL-3.0-or-later
 
#ifndef DDECAL_DIAGONAL_SOLVER_H
#define DDECAL_DIAGONAL_SOLVER_H
 
#include "SolverBase.h"
#include "SolveData.h"
 
namespace dp3 {
namespace ddecal {
 
class DiagonalSolver final : public SolverBase {
 public:
  DiagonalSolver() : SolverBase() {}
 
  SolveResult Solve(const FullSolveData& data,
                    std::vector<std::vector<DComplex>>& solutions,
                    double time) override;
 
  size_t NSolutionPolarizations() const override { return 2; }
 
 private:
  void PerformIteration(size_t ch_block,
                        const FullSolveData::ChannelBlockData& cb_data,
                        std::vector<Matrix>& g_times_cs,
                        std::vector<std::vector<Complex>>& vs,
                        const std::vector<DComplex>& solutions,
                        SolutionTensor& next_solutions);
 
  void InitializeModelMatrix(
      const FullSolveData::ChannelBlockData& channel_block_data,
      std::vector<Matrix>& g_times_cs,
      std::vector<std::vector<Complex>>& vs) const;
};
 
}  // namespace ddecal
}  // namespace dp3
 
#endif  // DDECAL_DIAGONAL_SOLVER_H