with Standard_Natural_Numbers;           use Standard_Natural_Numbers;
with Standard_Integer_Numbers;           use Standard_Integer_Numbers;
with Standard_Floating_Numbers;          use Standard_Floating_Numbers;
with Standard_Natural_Vectors;

package Bits_of_Doubles is

-- DESCRIPTION :
--   Provides functions to manipulate the bits of the fraction
--   of a double precision floating point number.

  procedure expand_52bits
              ( bits : out Standard_Natural_Vectors.Vector;
                nbr : in integer64 );

  -- DESCRIPTION :
  --   Given in bits space for 52 bits, in the vector bits of range 0..51,
  --   fills the bits with the binary expansin of the number nbr.

  procedure write_52bits
              ( bits : in Standard_Natural_Vectors.Vector );

  -- DESCRIPTION :
  --   Writes 52 bits in groups of 4.

  procedure Fraction_Exponent
              ( x : in double_float;
                f : out integer64; e : out integer32 );

  -- DESCRIPTION :
  --   Given a double float number x, returns the fraction f
  --   as a 64-bit integer from the fraction of x, multiplied by 2^52,
  --   and returns the exponent e as a 32-bit integer.

  procedure write_52bits_expo ( x : in double_float );

  -- DESCRIPTION :
  --   Shows the bits sequence of the fraction of x, in groups of 4,
  --   and the decimal representation of the exponent of x.

  function Bit_Equal ( x,y : double_float ) return boolean;

  -- DESCRIPTION :
  --   Returns true if x and y have the same sign, same exponent,
  --   and the same bits in their fraction.

  procedure write_fraction_bits ( nbr : in double_float );

  -- DESCRIPTION :
  --   Writes the bits of the fraction of nbr in binary,
  --   without the use of an auxiliary vector of natural numbers.

  function value_52bits
             ( bits : Standard_Natural_Vectors.Vector ) return integer64;

  -- DESCRIPTION :
  --   Given 52 bits in the vector bits, or range 0..51,
  --   returns its value as a 64-bit integer.

  function chop_last_bits
             ( nbr : double_float; lastbits : natural32 )
             return double_float;

  -- DESCRIPTION :
  --   Returns the number nbr with the last bits, equal as lastbits,
  --   chopped off, that is: set to zero.

  -- REQUIRED : lastbits < 52.

  procedure chop_last_bits
             ( nbr : in out double_float; lastbits : in natural32;
               headbits : out Standard_Natural_Vectors.Vector;
               tailbits : out Standard_Natural_Vectors.Vector );

  -- DESCRIPTION :
  --   Sets the last bits of the number nbr to zero.

  -- REQUIRED : lastbits < 52.

  -- ON ENTRY :
  --   nbr     a 64-bit double floating-point number;
  --   lastbits is the number of last bits to set to zero.

  -- ON RETURN :
  --   nbr     the number with the last bits removed;
  --   headbits are the leading bits which remain unchanged;
  --   tailbits are the last bits of the original number nbr.

  procedure insert_first_bits
              ( bits : in out Standard_Natural_Vectors.Vector;
                firstbits : in natural32;
                headbits : in Standard_Natural_Vectors.Vector );

  -- DESCRIPTION :
  --   Inserts the first bits, as many as the value of firstbits,
  --   of the headbits, into the bits.

  -- REQUIRED : firstbits < 52.

  -- ON ENTRY :
  --   bits     a sequence of 52 bits, of range 0..51;
  --   firstbits is the number of bits to be inserted;
  --   headbits contains the first bits to be inserted.

  -- ON RETURN :
  --   bits     contains the first bits of headbits,
  --            with the original bits shifted to the end.

  procedure insert_first_bits
              ( nbr : in out double_float;
                firstbits : in natural32;
                headbits : in Standard_Natural_Vectors.Vector );

  -- DESCRIPTION :
  --   Inserts the first bits, as many as the value of firstbits,
  --   of the headbits, into the fraction of the number nbr.

  -- REQUIRED : firstbits < 52.

  -- ON ENTRY :
  --   nbr      a 64-bit float;
  --   firstbits is the number of bits to be inserted;
  --   headbits contains the first bits to be inserted.

  -- ON RETURN :
  --   nbr      its fraction contains the first bits of headbits,
  --            with the original bits shifted to the end.

  function insert_first_bits
             ( nbr : double_float;
               firstbits : natural32;
               headbits : Standard_Natural_Vectors.Vector )
             return double_float;

  -- DESCRIPTION :
  --   Inserts the first bits, as many as the value of firstbits,
  --   of the headbits, into the fraction of the number nbr.
  --   Returns the number with the original bits of the fraction
  --   of nbr shifted to the end and the first bits inserted
  --   from headbits.

  -- REQUIRED : firstbits < 52.

  -- ON ENTRY :
  --   nbr      a 64-bit float;
  --   firstbits is the number of bits to be inserted;
  --   headbits contains the first bits to be inserted.

  procedure Mod_Split ( x : in double_float;
                        xhi,xlo : out double_float;
                        verbose : in boolean := true );

  -- DESCRIPTION :
  --   Splits the fraction of the double x in two equal halves,
  --   using modular arithmetic, without vectors of natural numbers.

  procedure Vec_Split ( x : in double_float;
                        xhi,xlo : out double_float );

  -- DESCRIPTION :
  --   Splits the fraction of the double x in two equal halves,
  --   using vectors of natural numbers.

  procedure Split ( x : in double_float;
                    x0,x1,x2,x3 : out double_float );

  -- DESCRIPTION :
  --   Splits the 52 bits in the fraction of x in four equal parts,
  --   returned in x0, x1, x2, x3, with x0 > x1 > x2 > x3.
  --   On return: Bit_Equal(x,x0+x1+x2+x3) is true
  --   and x - (x0 + x1 + x2 + x3) is exactly zero.

  procedure Split ( x : in double_float;
                    x0,x1,x2,x3,x4,x5,x6,x7 : out double_float );

  -- DESCRIPTION :
  --   Splits the 52 bits in the fraction of x in eight parts,
  --   the first four x0, x1, x2, x3 have each 7 bits of x in their fraction,
  --   the next four x4, x5, x6, x7 have each 8 bits of x in their fraction.

  function Last_Zero_Count ( x : integer64 ) return natural32;

  -- DESCRIPTION :
  --   Returns the number of last zeros in x.

  function Last_Zero_Count ( x : double_float ) return natural32;

  -- DESCRIPTION :
  --   Returns the number of last zeros in the fraction of x.

end Bits_of_Doubles;
