purify_fitsio.cc File Reference

#include "catch2/catch_all.hpp"
#include "purify/config.h"
#include "purify/types.h"
#include "purify/logging.h"
#include "purify/directories.h"
#include "purify/pfitsio.h"

Include dependency graph for purify_fitsio.cc:

Go to the source code of this file.

Functions
	TEST_CASE ("header")

Function Documentation

TEST_CASE()

TEST_CASE

(

"header"

)

Definition at line 90 of file purify_fitsio.cc.

                    {
  const std::string fits_name = "test_image.fits";
  const t_real mean_frequency = 123;  // in MHz
  const t_real cell_x = 5;            // in arcseconds
  const t_real cell_y = 2;            // in arcseconds
  const t_real ra = 12;               // in radians, converted to decimal degrees before write
  const t_real dec = 98;              // in radians, converted to decimal degrees before write
  const std::string pix_units = "Jy/PIXEL";
  const t_real channels_total = 2;
  const t_real channel_width = 11;  // in MHz
  const t_real polarisation = stokes_int.at(stokes::I);
  const int niters = 10;           // number of iterations
  const bool hasconverged = true;  // stating if model has converged
  const t_real relative_variation = 1e-3;
  const t_real residual_convergence = 1e-4;
  const t_real epsilon = 10;
  pfitsio::header_params header;
  header.fits_name = fits_name;
  header.mean_frequency = mean_frequency;
  header.cell_x = cell_x;
  header.cell_y = cell_y;
  header.ra = ra;
  header.dec = dec;
  header.pix_units = pix_units;
  header.channels_total = channels_total;
  header.channel_width = channel_width;
  header.polarisation = polarisation;
  header.niters = niters;
  header.hasconverged = hasconverged;
  header.relative_variation = relative_variation;
  header.residual_convergence = residual_convergence;
  header.epsilon = epsilon;
  const auto header_test = pfitsio::header_params(
      fits_name, pix_units, channels_total, ra, dec, stokes::I, cell_x, cell_y, mean_frequency,
      channel_width, niters, hasconverged, relative_variation, residual_convergence, epsilon);
  CHECK(header_test == header);
  CHECK(header_test == pfitsio::header_params(fits_name, pix_units, channels_total, ra, dec,
                                              stokes_string.at("p"), cell_x, cell_y, mean_frequency,
                                              channel_width, niters, hasconverged,
                                              relative_variation, residual_convergence, epsilon));
  CHECK_THROWS(header_test == pfitsio::header_params(fits_name, pix_units, channels_total, ra, dec,
                                                     stokes_string.at("z"), cell_x, cell_y,
                                                     mean_frequency, channel_width, niters,
                                                     hasconverged, relative_variation,
                                                     residual_convergence, epsilon));
  CHECK(header_test != pfitsio::header_params(fits_name, pix_units, channels_total, ra, dec,
                                              stokes::Q, cell_x, cell_y, mean_frequency,
                                              channel_width, niters, hasconverged,
                                              relative_variation, residual_convergence, epsilon));
}

References purify::pfitsio::header_params::cell_x, purify::pfitsio::header_params::cell_y, purify::pfitsio::header_params::channel_width, purify::pfitsio::header_params::channels_total, CHECK, CHECK_THROWS, purify::pfitsio::header_params::dec, purify::pfitsio::header_params::epsilon, purify::pfitsio::header_params::fits_name, purify::pfitsio::header_params::hasconverged, purify::I, purify::pfitsio::header_params::mean_frequency, purify::pfitsio::header_params::niters, purify::pfitsio::header_params::pix_units, purify::pfitsio::header_params::polarisation, purify::Q, purify::pfitsio::header_params::ra, purify::pfitsio::header_params::relative_variation, purify::pfitsio::header_params::residual_convergence, purify::stokes_int, and purify::stokes_string.