uvfits.cc File Reference

#include "purify/config.h"
#include "purify/types.h"
#include "catch2/catch_all.hpp"
#include "purify/logging.h"
#include <iostream>
#include "purify/directories.h"
#include "purify/uvfits.h"

Include dependency graph for uvfits.cc:

Go to the source code of this file.

Functions
	TEST_CASE ("readfile")

Function Documentation

TEST_CASE()

TEST_CASE

(

"readfile"

)

Definition at line 12 of file uvfits.cc.

                      {
  purify::logging::set_level("debug");
  const std::string filename = atca_filename("0332-391");
  SECTION("filter") {
    const auto uvfits = pfitsio::read_uvfits(filename + ".uvfits", true, stokes::I);
    const auto vis = utilities::read_visibility(filename + ".vis", true);
    REQUIRE(299792458. == constant::c);
    REQUIRE(1483999958.029309988 == Approx(uvfits.frequencies(0)).epsilon(1e-12));
    REQUIRE(13 == uvfits.frequencies.size());
    CHECK(uvfits.size() / 22675 == 10);
    CHECK(vis.size() / 22675 == 10);
    CHECK(uvfits.size() == vis.size());
    for (int i = 0; i < uvfits.size(); i++) {
      CAPTURE(i);
      CAPTURE(vis.u(i));
      CAPTURE(uvfits.u(i));
      CAPTURE(vis.v(i));
      CAPTURE(uvfits.v(i));
      CAPTURE(vis.w(i));
      CAPTURE(uvfits.w(i));
      CAPTURE(vis.vis(i));
      CAPTURE(uvfits.vis(i));
      CAPTURE(vis.weights(i));
      CAPTURE(uvfits.weights(i));
      REQUIRE(std::abs((uvfits.u(i) - vis.u(i)) / vis.u(i)) < 1e-4);
      REQUIRE(std::abs((uvfits.v(i) - vis.v(i)) / vis.v(i)) < 1e-4);
      REQUIRE(std::abs((uvfits.w(i) - vis.w(i)) / vis.w(i)) < 1e-4);
      if (std::abs(vis.weights(i)) > 0)
        REQUIRE(std::abs((uvfits.weights(i) - vis.weights(i)) / vis.weights(i)) < 1e-4);
      if (std::abs(uvfits.weights(i)) > 0)
        REQUIRE(std::abs((uvfits.weights(i) - vis.weights(i)) / uvfits.weights(i)) < 1e-4);
      if (std::abs(vis.vis(i)) > 0)
        REQUIRE(std::abs((uvfits.vis(i) - vis.vis(i)) / vis.vis(i)) < 1e-4);
      if (std::abs(uvfits.vis(i)) > 0)
        REQUIRE(std::abs((uvfits.vis(i) - vis.vis(i)) / uvfits.vis(i)) < 1e-4);
    }
  }
  SECTION("no_filter") {
    const auto uvfits = pfitsio::read_uvfits(filename + ".uvfits", false, stokes::I);
    const auto vis = utilities::read_visibility(filename + "_nofilter.vis", true);
    REQUIRE(299792458. == constant::c);
    REQUIRE(1483999958.029309988 == Approx(uvfits.frequencies(0)).epsilon(1e-12));
    REQUIRE(uvfits.size() == vis.size());
    for (int i = 0; i < uvfits.size(); i++) {
      CAPTURE(i);
      CAPTURE(vis.u(i));
      CAPTURE(uvfits.u(i));
      CAPTURE(vis.v(i));
      CAPTURE(uvfits.v(i));
      CAPTURE(vis.w(i));
      CAPTURE(uvfits.w(i));
      CAPTURE(vis.vis(i));
      CAPTURE(uvfits.vis(i));
      CAPTURE(vis.weights(i));
      CAPTURE(uvfits.weights(i));
      REQUIRE(uvfits.u(i) == Approx(vis.u(i)).epsilon(1e-4));
      REQUIRE(uvfits.v(i) == Approx(vis.v(i)).epsilon(1e-4));
      REQUIRE(uvfits.w(i) == Approx(vis.w(i)).epsilon(1e-4));
    }
  }
}

References purify::atca_filename(), purify::constant::c, CHECK, purify::I, purify::pfitsio::read_uvfits(), purify::utilities::read_visibility(), and purify::logging::set_level().