Simulate.h

Go to the documentation of this file.

// Simulate.h: Simulate visibilities for a patch of sources.
//
// Copyright (C) 2023 ASTRON (Netherlands Institute for Radio Astronomy)
// SPDX-License-Identifier: GPL-3.0-or-later
 
 
#ifndef DP3_BASE_SIMULATE_H_
#define DP3_BASE_SIMULATE_H_
 
#include <memory>
 
#include <xtensor/containers/xtensor.hpp>
 
#include <aocommon/xt/span.h>
 
#include "Baseline.h"
#include "Cursor.h"
 
#include "model/Patch.h"
 
#include "Direction.h"
#include "DPBuffer.h"
 
namespace dp3 {
namespace base {
 
std::vector<int> SetupUvwSplitting(unsigned int nant,
                                   const std::vector<int>& ant1,
                                   const std::vector<int>& ant2);
 
std::vector<size_t> GetBaselinesSortedByLength(
    const std::vector<int>& antennas1, const std::vector<int>& antennas2,
    const std::vector<std::array<double, 3>>& antenna_positions);
 
std::vector<size_t> GetBaselineSelection(
    size_t n_antennas, const std::vector<size_t>& bl_idx_sorted,
    const std::vector<int>& antennas1, const std::vector<int>& antennas2);
 
std::vector<int> SetupUvwSplitting(
    unsigned int nant, const std::vector<int>& ant1,
    const std::vector<int>& ant2,
    const std::vector<std::array<double, 3>>& antenna_positions);
 
void SplitUvw(const std::vector<int>& baseline_indices,
              const std::vector<Baseline>& baselines,
              const DPBuffer::UvwType& uvw_bl, xt::xtensor<double, 2>& uvw_ant);
 
void rotateUVW(const Direction& from, const Direction& to, size_t nUVW,
               double* uvw);
 
void simulate(const Direction& reference,
              const std::shared_ptr<const model::Patch>& patch, size_t nStation,
              size_t nBaseline, size_t nChannel,
              const_cursor<Baseline> baselines, const_cursor<double> freq,
              const_cursor<double> uvw, cursor<std::complex<double>> buffer);
 
}  // namespace base
}  // namespace dp3
 
#endif