text: /*
      Copyright 2013 Christian Henning
      Use, modification and distribution are subject to the Boost Software License,
      Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
      http://www.boost.org/LICENSE_1_0.txt).
  */
  
  /// \file xyz_test.cpp
  /// \brief Unit test for XYZ Colorspace
  /// \author Davide Anastasia <davideanastasia@users.sourceforge.net>
  
  #include <cstdint>
  #include <iostream>
  #include <limits>
  #include <boost/gil.hpp>
  #include <boost/gil/extension/toolbox/color_spaces/xyz.hpp>
  
  #include <boost/test/unit_test.hpp>
  
  using namespace boost;
  using namespace std;
  
  const float SKEW = 0.0001f;
  
  BOOST_AUTO_TEST_SUITE( toolbox_tests )
  
  BOOST_AUTO_TEST_CASE(rgb32f_xyz32f_1)
  {
      gil::xyz32f_pixel_t xyz32f;
      gil::rgb32f_pixel_t p32f(.934351f, 0.785446f, .105858f), p32f_b;
      gil::color_convert(p32f, xyz32f);
      gil::color_convert(xyz32f, p32f_b);
  
      BOOST_TEST_MESSAGE( p32f[0] << " "
                          << p32f[1] << " "
                          << p32f[2] << " -> "
                          << xyz32f[0] << " "
                          << xyz32f[1] << " "
                          << xyz32f[2] << " -> "
                          << p32f_b[0] << " "
                          << p32f_b[1] << " "
                          << p32f_b[2] );
  
      BOOST_CHECK( abs(p32f[0] - p32f_b[0]) < SKEW );
      BOOST_CHECK( abs(p32f[1] - p32f_b[1]) < SKEW );
      BOOST_CHECK( abs(p32f[2] - p32f_b[2]) < SKEW );
      BOOST_CHECK( abs( xyz32f[0] - 0.562669) < SKEW );
      BOOST_CHECK( abs( xyz32f[1] - 0.597462) < SKEW );
      BOOST_CHECK( abs( xyz32f[2] - 0.096050) < SKEW );
  }
  
  BOOST_AUTO_TEST_CASE(rgb32f_xyz32f_2)
  {
      gil::xyz32f_pixel_t xyz32f;
      gil::rgb32f_pixel_t p32f(.694617f, 0.173810f, 0.218710f), p32f_b;
      gil::color_convert(p32f, xyz32f);
      gil::color_convert(xyz32f, p32f_b);
  
      BOOST_TEST_MESSAGE( p32f[0] << " "
                          << p32f[1] << " "
                          << p32f[2] << " -> "
                          << xyz32f[0] << " "
                          << xyz32f[1] << " "
                          << xyz32f[2] << " -> "
                          << p32f_b[0] << " "
                          << p32f_b[1] << " "
                          << p32f_b[2] );
  
      BOOST_CHECK( abs(p32f[0] - p32f_b[0]) < SKEW );
      BOOST_CHECK( abs(p32f[1] - p32f_b[1]) < SKEW );
      BOOST_CHECK( abs(p32f[2] - p32f_b[2]) < SKEW );
      BOOST_CHECK( abs( xyz32f[0] - 0.197823) < SKEW );
      BOOST_CHECK( abs( xyz32f[1] - 0.114731) < SKEW );
      BOOST_CHECK( abs( xyz32f[2] - 0.048848) < SKEW );
  }
  
  BOOST_AUTO_TEST_CASE(xyz32f_rgb32f_1)
  {
      gil::xyz32f_pixel_t xyz32f(.332634f, .436288f, .109853f), xyz32f_b;
      gil::rgb32f_pixel_t p32f;
      gil::color_convert(xyz32f, p32f);
      gil::color_convert(p32f, xyz32f_b);
  
      BOOST_TEST_MESSAGE( xyz32f[0] << " "
                          << xyz32f[1] << " "
                          << xyz32f[2] << " -> "
                          << p32f[0] << " "
                          << p32f[1] << " "
                          << p32f[2] << " -> "
                          << xyz32f_b[0] << " "
                          << xyz32f_b[1] << " "
                          << xyz32f_b[2] );
  
      BOOST_CHECK( abs(xyz32f_b[0] - xyz32f[0]) < SKEW );
      BOOST_CHECK( abs(xyz32f_b[1] - xyz32f[1]) < SKEW );
      BOOST_CHECK( abs(xyz32f_b[2] - xyz32f[2]) < SKEW );
      BOOST_CHECK( abs( p32f[0] - 0.628242) < SKEW );
      BOOST_CHECK( abs( p32f[1] - 0.735771) < SKEW );
      BOOST_CHECK( abs( p32f[2] - 0.236473) < SKEW );
  }
  
  BOOST_AUTO_TEST_CASE(xyz32f_rgb32f_2)
  {
      gil::xyz32f_pixel_t xyz32f(.375155f, .352705f, .260025f), xyz32f_b;
      gil::rgb32f_pixel_t p32f;
      gil::color_convert(xyz32f, p32f);
      gil::color_convert(p32f, xyz32f_b);
  
      BOOST_TEST_MESSAGE( xyz32f[0] << " "
                          << xyz32f[1] << " "
                          << xyz32f[2] << " -> "
                          << p32f[0] << " "
                          << p32f[1] << " "
                          << p32f[2] << " -> "
                          << xyz32f_b[0] << " "
                          << xyz32f_b[1] << " "
                          << xyz32f_b[2] );
  
      BOOST_CHECK( abs(xyz32f_b[0] - xyz32f[0]) < SKEW );
      BOOST_CHECK( abs(xyz32f_b[1] - xyz32f[1]) < SKEW );
      BOOST_CHECK( abs(xyz32f_b[2] - xyz32f[2]) < SKEW );
      BOOST_CHECK( abs( p32f[0] - 0.763580) < SKEW );
      BOOST_CHECK( abs( p32f[1] - 0.591622) < SKEW );
      BOOST_CHECK( abs( p32f[2] - 0.510392) < SKEW );
  }
  
  BOOST_AUTO_TEST_CASE(rgb8u_xyz32f_1)
  {
      gil::xyz32f_pixel_t xyz32f;
      gil::rgb8_pixel_t p8u(177, 44, 56), p8u_b;
      gil::color_convert(p8u, xyz32f);
      gil::color_convert(xyz32f, p8u_b);
  
      BOOST_TEST_MESSAGE( static_cast<int>(p8u[0]) << " "
                     << static_cast<int>(p8u[1]) << " "
                     << static_cast<int>(p8u[2]) << " -> "
                     << xyz32f[0] << " "
                     << xyz32f[1] << " "
                     << xyz32f[2] << " -> "
                     << static_cast<int>(p8u_b[0]) << " "
                     << static_cast<int>(p8u_b[1]) << " "
                     << static_cast<int>(p8u_b[2]) );
  
      BOOST_CHECK(p8u[0] == p8u_b[0]);
      BOOST_CHECK(p8u[1] == p8u_b[1]);
      BOOST_CHECK(p8u[2] == p8u_b[2]);
  }
  
  BOOST_AUTO_TEST_CASE(rgb8u_xyz32f_2)
  {
      gil::xyz32f_pixel_t xyz32f;
      gil::rgb8_pixel_t p8u(72, 90, 165), p8u_b;
      gil::color_convert(p8u, xyz32f);
      gil::color_convert(xyz32f, p8u_b);
  
      BOOST_TEST_MESSAGE(
                        static_cast<int>(p8u[0]) << " "
                     << static_cast<int>(p8u[1]) << " "
                     << static_cast<int>(p8u[2]) << " -> "
                     << xyz32f[0] << " "
                     << xyz32f[1] << " "
                     << xyz32f[2] << " -> "
                     << static_cast<int>(p8u_b[0]) << " "
                     << static_cast<int>(p8u_b[1]) << " "
                     << static_cast<int>(p8u_b[2]) );
  
      BOOST_CHECK(p8u[0] == p8u_b[0]);
      BOOST_CHECK(p8u[1] == p8u_b[1]);
      BOOST_CHECK(p8u[2] == p8u_b[2]);
  }
  
  BOOST_AUTO_TEST_CASE(rgb16u_xyz32f_1)
  {
      gil::xyz32f_pixel_t xyz32f;
      gil::rgb16_pixel_t p16u(12564, 20657, 200), p16u_b;
      gil::color_convert(p16u, xyz32f);
      gil::color_convert(xyz32f, p16u_b);
  
      BOOST_TEST_MESSAGE(
                        static_cast<int>(p16u[0]) << " "
                     << static_cast<int>(p16u[1]) << " "
                     << static_cast<int>(p16u[2]) << " -> "
                     << xyz32f[0] << " "
                     << xyz32f[1] << " "
                     << xyz32f[2] << " -> "
                     << static_cast<int>(p16u_b[0]) << " "
                     << static_cast<int>(p16u_b[1]) << " "
                     << static_cast<int>(p16u_b[2]) );
  
      BOOST_CHECK(p16u[0] == p16u_b[0]);
      BOOST_CHECK(p16u[1] == p16u_b[1]);
      BOOST_CHECK(p16u[2] == p16u_b[2]);
  }
  
  BOOST_AUTO_TEST_SUITE_END()
  