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

using namespace std;

# include "gmsh_io.hpp"

int main ( );
void test01 ( );
void test02 ( );
void timestamp ( );

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

int main ( )

//****************************************************************************80
//
//  Purpose:
//
//    gmsh_io_test() tests gmsh_io().
//
//  Licensing:
//
//    This code is distributed under the MIT license. 
//
//  Modified:
//
//    21 October 2014
//
//  Author:
//
//    John Burkardt
//
{
  timestamp ( );
  cout << "\n";
  cout << "gmsh_io_test():\n";
  cout << "  C++ version\n";
  cout << "  Test GMSH_IO().\n";

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

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

void test01 ( )

//****************************************************************************80
//
//  Purpose:
//
//    TEST01 gets the example 2D data and writes it to a file.
//
//  Licensing:
//
//    This code is distributed under the MIT license. 
//
//  Modified:
//
//    20 October 2014
//
//  Author:
//
//    John Burkardt
//
{
  int *element_node;
  int element_num;
  int element_order;
  string gmsh_filename = "example_2d.msh";
  int m;
  int node_num;
  double *node_x;

  cout << "\n";
  cout << "TEST01:\n";
  cout << "  Get example 2D data, write to a file.\n";
//
//  Get sizes.
//
  gmsh_mesh2d_node_size_example ( node_num, m );

  gmsh_mesh2d_element_size_example ( element_num, element_order );
//
//  Print the sizes.
//
  cout << "\n";
  cout << "  Number of nodes = " << node_num << "\n";
  cout << "  Spatial dimension = " << m << "\n";
  cout << "  Number of elements = " << element_num << "\n";
  cout << "  Order of elements = " << element_order << "\n";
//
//  Get the data.
//
  node_x = gmsh_mesh2d_node_data_example ( node_num, m );

  element_node = gmsh_mesh2d_element_data_example ( element_num, element_order );
//
//  Print some of the data.
//
  r8mat_transpose_print_some ( m, node_num, node_x, 
    1, 1, m, 10, "  Coordinates for first 10 nodes:" );

  i4mat_transpose_print_some ( element_order, element_num, element_node, 
    1, 1, element_order, 10, "  Node connectivity of first 10 elements:" );
//
//  Write the GMSH file.
//
  gmsh_mesh2d_write ( gmsh_filename, m, node_num, node_x, 
    element_order, element_num, element_node );

  cout << "\n";
  cout << "  Wrote example data to file \"" << gmsh_filename << "\"\n";
//
//  Clean up.
//
  delete [] element_node;
  delete [] node_x;

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

void test02 ( )

//****************************************************************************80
//
//  Purpose:
//
//    TEST02 reads the example data from a file.
//
//  Licensing:
//
//    This code is distributed under the MIT license. 
//
//  Modified:
//
//    21 October 2014
//
//  Author:
//
//   John Burkardt
//
{
  int *element_node;
  int element_num;
  int element_order;
  string gmsh_filename = "example_2d.msh";
  int m;
  int node_num;
  double *node_x;

  cout << "\n";
  cout << "TEST02:\n";
  cout << "  Read data from a file.\n";
//
//  Get the data size.
//
  gmsh_size_read ( gmsh_filename, node_num, m, element_num, 
    element_order );
//
//  Print the sizes.
//
  cout << "\n";
  cout << "  Node data read from file \"" << gmsh_filename << "\"\n";
  cout << "\n";
  cout << "  Number of nodes = " << node_num << "\n";
  cout << "  Spatial dimension = " << m << "\n";
  cout << "  Number of elements = " << element_num << "\n";
  cout << "  Element order = " << element_order << "\n";
//
//  Allocate memory.
//
  node_x = new double[m * node_num];
  element_node = new int[ element_order * element_num];
//
//  Get the data.
//
  gmsh_data_read ( gmsh_filename, m, node_num, node_x, 
    element_order, element_num, element_node );
//
//  Print some of the data.
//
  r8mat_transpose_print_some ( m, node_num, node_x, 
    1, 1, m, 10, "  Coordinates for first 10 nodes:" );

  i4mat_transpose_print_some ( element_order, element_num, element_node, 
    1, 1, element_order, 10, "  Connectivity for first 10 elements:" );
//
//  Clean up.
//
  delete [] element_node;
  delete [] node_x;

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

void timestamp ( )

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

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

  now = time ( NULL );
  tm = localtime ( &now );

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

  cout << time_buffer << "\n";

  return;
# undef TIME_SIZE
}