int64x64.cc

Go to the documentation of this file.

/* -*- Mode:C++; c-file-style:"gnu"; indent-tabs-mode:nil; -*- */
/*
 * Copyright (c) 2010 INRIA
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License version 2 as
 * published by the Free Software Foundation;
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */

#include "int64x64.h"
#include <stdint.h>
#include <iostream>
#include <iomanip>  // showpos
#include <sstream>
#include "assert.h"
#include "log.h"


namespace ns3 {

// Note:  Logging in this file is largely avoided due to the
// number of calls that are made to these functions and the possibility
// of causing recursions leading to stack overflow
NS_LOG_COMPONENT_DEFINE ("int64x64");

std::ostream &operator << (std::ostream &os, const int64x64_t &value)
{
  const bool negative = (value < 0);
  const int64x64_t absVal = (negative ? -value : value);
  
  int64_t hi = absVal.GetHigh ();

  // Save stream format flags
  const std::streamsize precision = os.precision ();
  std::ios_base::fmtflags ff = os.flags ();
  const bool floatfield = os.flags () & std::ios_base::floatfield;
  os << std::setw (1) << std::noshowpos;
  
  os << std::right << (negative ? "-" : "+");

  std::ostringstream oss;
  oss << hi << ".";  // collect the digits here so we can round properly


  int64x64_t low(0, absVal.GetLow ());
  int places = 0;    // Number of decimal places printed so far
  bool more = true;  // Should we print more digits?

#define HEXHILOW(hi, lo) \
  std::hex << std::setfill ('0') << std::right << " (0x"        \
       << std::setw (16) << hi << " "               \
       << std::setw (16) << lo                  \
       << std::dec << std::setfill (' ') << std::left << ")"

  
  NS_LOG_LOGIC (std::endl
        << (floatfield ? " f" : "  ")
        << "[" << precision << "] " << hi << ". "
        << HEXHILOW (hi, low.GetLow ())
        );

  int64_t digit;
  do 
    {
      low *= 10;
      digit = low.GetHigh ();
      NS_ASSERT_MSG ( (0 <= digit) && (digit <= 9),
              "digit " << digit << " out of range [0,9] "
              << " streaming out "
              << HEXHILOW (value.GetHigh (), value.GetLow ()) );
      low -= digit;

      oss << std::setw (1) << digit;

      ++places;
      if (floatfield)
    {
      more = places < precision;
    }
      else  // default
    {
      // Full resolution is 20 decimal digits
      more = low.GetLow () && (places < 20);
    }

      NS_LOG_LOGIC ((more ? "+" : " ")
            << (floatfield ? "f" : " ")
            << "[" << places << "] " << digit
            << HEXHILOW (low.GetHigh (), low.GetLow ())
            << std::dec << std::setfill (' ' ) << std::left);

    } while (more);

  // Check if we need to round the last printed digit,
  // based on the first unprinted digit
  std::string digits = oss.str ();
  low *= 10;
  int64_t nextDigit = low.GetHigh ();
  if ( (nextDigit > 5) || ((nextDigit == 5) && (digit % 2 == 1)) )
    {
      // Walk backwards with the carry
      bool carry = true;
      for (std::string::reverse_iterator rit = digits.rbegin ();
       rit != digits.rend ();
       ++rit)
    {
      if (*rit == '.')  // Skip over the decimal point
        {
          continue ;
        }
      
      ++(*rit);         // Add the carry
      if (*rit <= '9')  // Relies on character order...
        {
          carry = false;
          break ;       // Carry complete
        }
      else
        {
          *rit = '0';     // Continue carry to next higher digit
        }
    }
      if (carry)            // If we still have a carry...
    {
      digits.insert (digits.begin (), '1');
    }
    }
  os << digits;
  
  os.flags (ff);  // Restore stream flags
  return os;
}

static uint64_t ReadHiDigits (std::string str)
{
  const char *buf = str.c_str ();
  uint64_t retval = 0;
  while (*buf != 0)
    {
      retval *= 10;
      retval += *buf - '0';
      buf++;
    }
  return retval;
}

static uint64_t ReadLoDigits (std::string str)
{
  int64x64_t low;
  const int64x64_t round (0, 5);  // Round last place in division

  for (std::string::const_reverse_iterator rit = str.rbegin ();
       rit != str.rend ();
       ++rit)
    {
      int digit = *rit - '0';
      NS_ASSERT_MSG ( (0 <= digit) && (digit <= 9),
              "digit " << digit << " out of range [0,9]"
              << " streaming in low digits \"" << str << "\"");
      low = (low + digit + round) / 10; 
    }
  
  return low.GetLow ();
}
    
std::istream &operator >> (std::istream &is, int64x64_t &value)
{
  std::string str;

  is >> str;
  bool negative;
  // skip heading spaces
  std::string::size_type cur;
  cur = str.find_first_not_of (" ");
  std::string::size_type next;
  // first, remove the sign.
  next = str.find ("-", cur);
  if (next != std::string::npos)
    {
      negative = true;
      next++;
    }
  else
    {
      next = str.find ("+", cur);
      if (next != std::string::npos)
        {
          next++;
        }
      else
        {
          next = cur;
        }
      negative = false;
    }
  cur = next;
  int64_t hi;
  uint64_t lo;
  next = str.find (".", cur);
  if (next != std::string::npos)
    {
      hi = ReadHiDigits (str.substr (cur, next-cur));
      lo = ReadLoDigits (str.substr (next+1, str.size ()-(next+1)));
    }
  else if (cur != std::string::npos)
    {
      hi = ReadHiDigits (str.substr (cur, str.size ()-cur));
      lo = 0;
    }
  else
    {
      hi = 0;
      lo = 0;
    }
  
  value = int64x64_t (hi, lo);
  value = negative ? -value : value;

  return is;
}

} // namespace ns3