update sources to v12.2.4

[ceph.git] / ceph / src / boost / libs / sort / example / double.cpp

// spreadsort double sorting example.\r
//\r
//  Copyright Steven Ross 2009-2014.\r
//\r
// Distributed under the Boost Software License, Version 1.0.\r
//    (See accompanying file LICENSE_1_0.txt or copy at\r
//          http://www.boost.org/LICENSE_1_0.txt)\r
\r
//  See http://www.boost.org/libs/sort for library home page.\r
\r
#include <boost/sort/spreadsort/spreadsort.hpp>\r
#include <time.h>\r
#include <stdio.h>\r
#include <stdlib.h>\r
#include <algorithm>\r
#include <vector>\r
#include <string>\r
#include <fstream>\r
#include <iostream>\r
using namespace boost::sort::spreadsort;\r
\r
#define DATA_TYPE double\r
#define CAST_TYPE boost::int64_t\r
\r
//Pass in an argument to test std::sort\r
//Note that this converts NaNs and -0.0 to 0.0, so that sorting results are\r
//identical every time\r
int main(int argc, const char ** argv) {\r
  size_t uCount,uSize=sizeof(DATA_TYPE);\r
  bool stdSort = false;\r
  unsigned loopCount = 1;\r
  for (int u = 1; u < argc; ++u) {\r
    if (std::string(argv[u]) == "-std")\r
      stdSort = true;\r
    else\r
      loopCount = atoi(argv[u]);\r
  }\r
  std::ifstream input("input.txt", std::ios_base::in | std::ios_base::binary);\r
  if (input.fail()) {\r
    printf("input.txt could not be opened\n");\r
    return 1;\r
  }\r
  double total = 0.0;\r
  std::vector<DATA_TYPE> array;\r
  input.seekg (0, std::ios_base::end);\r
  size_t length = input.tellg();\r
  uCount = length/uSize;\r
  //Using this to support compilers that don't support 64-bit constants\r
  CAST_TYPE exponent_mask = 0x7ff00000;\r
  exponent_mask <<= 32;\r
  CAST_TYPE top_exponent_bit = 0x40000000;\r
  top_exponent_bit <<= 32;\r
  //Run multiple loops, if requested\r
  for (unsigned u = 0; u < loopCount; ++u) {\r
    input.seekg (0, std::ios_base::beg);\r
    //Conversion to a vector\r
    array.resize(uCount);\r
    unsigned v = 0;\r
    while (input.good() && v < uCount) {\r
      input.read(reinterpret_cast<char *>(&(array[v])), uSize );\r
     //Checking for denormalized numbers\r
     if (!(float_mem_cast<DATA_TYPE, CAST_TYPE>(array[v]) & exponent_mask)) {\r
       //Make the top exponent bit high\r
       CAST_TYPE temp = top_exponent_bit |\r
         float_mem_cast<DATA_TYPE, CAST_TYPE>(array[v]);\r
       memcpy(&(array[v]), &temp, sizeof(DATA_TYPE));\r
     }\r
     //Testcase doesn't sort NaNs; they just cause confusion\r
     if (!(array[v] < 0.0) && !(0.0 < array[v]))\r
      array[v] = 0.0;\r
     ++v;\r
    }\r
    clock_t start, end;\r
    double elapsed;\r
    start = clock();\r
    if (stdSort)\r
      //std::sort(&(array[0]), &(array[0]) + uCount);\r
      std::sort(array.begin(), array.end());\r
    else\r
      //float_sort(&(array[0]), &(array[0]) + uCount);\r
      float_sort(array.begin(), array.end());\r
    end = clock();\r
    elapsed = static_cast<double>(end - start) ;\r
    std::ofstream ofile;\r
    if (stdSort)\r
      ofile.open("standard_sort_out.txt", std::ios_base::out |\r
                 std::ios_base::binary | std::ios_base::trunc);\r
    else\r
      ofile.open("boost_sort_out.txt", std::ios_base::out |\r
                 std::ios_base::binary | std::ios_base::trunc);\r
    if (ofile.good()) {\r
      for (unsigned v = 0; v < array.size(); ++v) {\r
        ofile.write(reinterpret_cast<char *>(&(array[v])), sizeof(array[v]) );\r
      }\r
      ofile.close();\r
    }\r
    total += elapsed;\r
    array.clear();\r
  }\r
  if (stdSort)\r
    printf("std::sort elapsed time %f\n", total / CLOCKS_PER_SEC);\r
  else\r
    printf("spreadsort elapsed time %f\n", total / CLOCKS_PER_SEC);\r
  return 0;\r
}\r