# include <cstdlib>
# include <iostream>
# include <cmath>
# include <ctime>

using namespace std;

# include "bisection_integer.hpp"

int main ( );
void test01 ( );
int f01 ( int x );
void timestamp ( );

//****************************************************************************80

int main ( )

//****************************************************************************80
//
//  Purpose:
//
//    bisection_integer_test() tests bisection_integer().
//
//  Licensing:
//
//    This code is distributed under the MIT license.
//
//  Modified:
//
//    23 August 2012
//
//  Author:
//
//    John Burkardt
//
{
  timestamp ( );
  cout << "\n";
  cout << "BISECTION_INTEGER_TEST():\n";
  cout << "  C++ version\n";
  cout << "  Test BISECTION_INTEGER().\n";

  test01 ( );
//
//  Terminate.
//
  cout << "\n";
  cout << "BISECTION_INTEGER_TEST():\n";
  cout << "  Normal end of execution.\n";
  cout << "\n";
  timestamp ( );

  return 0;
}
//****************************************************************************80

void test01 ( )

//****************************************************************************80
//
//  Purpose:
//
//    TEST01 tests BISECTION_INTEGER;
//
//  Licensing:
//
//    This code is distributed under the MIT license.
//
//  Modified:
//
//    23 August 2012
//
//  Author:
//
//    John Burkardt
//
{
  int a;
  int b;
  int c;
  int fc;

  cout << "\n";
  cout << "TEST01\n";
  cout << "  BISECTION_INTEGER attempts to locate an integer root C\n";
  cout << "  of an equation F(C) = 0.\n";
  cout << "  The user supplies a change of sign interval [A,B].\n";
  cout << "  The function considered here has two real roots\n";
  cout << "  as well as an integer root, so the algorithm can\n";
  cout << "  fail depending on how the change of sign interval is chosen.\n";

  a = 4;
  b = 100;

  cout << "\n";
  cout << "  The initial change of sign interval is:\n";
  cout << "  F(" << a << ") = " << f01 ( a ) << "\n";
  cout << "  F(" << b << ") = " << f01 ( b ) << "\n";

  bisection_integer ( f01, a, b, c, fc );

  if ( fc == 0 )
  {
    cout << "\n";
    cout << "  An exact root was found at C = " << c << "\n";
  }
  else
  {
    cout << "\n";
    cout << "  An exact root was NOT found.\n";
    cout << "  The change of sign interval is now:\n";
    cout << "  F(" << a << ") = " << f01 ( a ) << "\n";
    cout << "  F(" << b << ") = " << f01 ( b ) << "\n";
  }

  a = -10;
  b = 15;

  cout << "\n";
  cout << "  The initial change of sign interval is:\n";
  cout << "  F(" << a << ") = " << f01 ( a ) << "\n";
  cout << "  F(" << b << ") = " << f01 ( b ) << "\n";

  bisection_integer ( f01, a, b, c, fc );

  if ( fc == 0 )
  {
    cout << "\n";
    cout << "  An exact root was found at C = " << c << "\n";
  }
  else
  {
    cout << "\n";
    cout << "  An exact root was NOT found.\n";
    cout << "  The change of sign interval is now:\n";
    cout << "  F(" << a << ") = " << f01 ( a ) << "\n";
    cout << "  F(" << b << ") = " << f01 ( b ) << "\n";
  }

  return;
}
//****************************************************************************80

int f01 ( int n )

//****************************************************************************80
//
//  Purpose:
//
//    F01 is a test function.
//
//  Discussion:
//
//    The polynomial has roots 1/2, 7/2, and 10.
//
//  Licensing:
//
//    This code is distributed under the MIT license.
//
//  Modified:
//
//    23 August 2012
//
//  Author:
//
//    John Burkardt
//
//  Parameters:
//
//    Input, int N, the argument.
//
//    Output, int F01, the function value.
//
{
  int value;

  value = ( 2 * n - 7 ) * ( 2 * n - 1 ) * ( n - 10 );

  return value;
}
//****************************************************************************80

void timestamp ( )

//****************************************************************************80
//
//  Purpose:
//
//    TIMESTAMP prints the current YMDHMS date as a time stamp.
//
//  Example:
//
//    31 May 2001 09:45:54 AM
//
//  Licensing:
//
//    This code is distributed under the MIT license.
//
//  Modified:
//
//    08 July 2009
//
//  Author:
//
//    John Burkardt
//
{
# define TIME_SIZE 40

  static char time_buffer[TIME_SIZE];
  const struct std::tm *tm_ptr;
  std::time_t now;

  now = std::time ( NULL );
  tm_ptr = std::localtime ( &now );

  std::strftime ( time_buffer, TIME_SIZE, "%d %B %Y %I:%M:%S %p", tm_ptr );

  std::cout << time_buffer << "\n";

  return;
# undef TIME_SIZE
}