# ifdef ANSI_HEADERS
#   include <cmath>
#   include <cstdlib>
#   include <cstring>
#   include <ctime>
# else
#   include <math.h>
#   include <stdlib.h>
#   include <string.h>
#   include <time.h>
# endif

# include <iostream>
# include <iomanip>
# include <fstream>

using namespace std;

# include "fsu.hpp"

int main ( );
void quality_test01 ( int dim_num, int n, double z[], int ns, int seed_init );
void quality_test02 ( int dim_num, int n, double z[], int ns, int seed_init );
void quality_test03 ( int dim_num, int n, double z[], int ns, int seed_init );
void quality_test04 ( int dim_num, int n, double z[], int ns, int seed_init );

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

int main ( )

//****************************************************************************80
//
//  Purpose:
//
//    FSU_QUALITY_TEST calls the FSU_QUALITY routines.
//
//  Modified:
//
//    10 November 2006
//
//  Author:
//
//    Max Gunzburger
//    John Burkardt
//
//  License:
//
//    Copyright (C) 2004  John Burkardt and Max Gunzburger
//
//    This library is free software; you can redistribute it and/or
//    modify it under the terms of the GNU Lesser General Public
//    License as published by the Free Software Foundation; either
//    version 2.1 of the License, or (at your option) any later version.
//
//    This library 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
//    Lesser General Public License for more details.
//
//    You should have received a copy of the GNU Lesser General Public
//    License along with this library; if not, write to the Free Software
//    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
{
  char input_filename[80];
  int n;
  int dim_num;
  int ns;
  int seed_init;
  double *z;

  timestamp ( );

  cout << "\n";
  cout << "FSU_QUALITY_TEST\n";
  cout << "  Test the C++ FSU_QUALITY library.\n";

  ns = 100000;
  seed_init = 123456789;
  strcpy ( input_filename, "halton_02_00100.txt" );

  dtable_header_read ( input_filename, &dim_num, &n );

  cout << "\n";
  cout << "FSU_QUALITY_PROBLEMS:\n";
  cout << "  Measures of uniform point dispersion.\n";
  cout << "\n";
  cout << "  The pointset was read from \"" << input_filename << "\"\n";
  cout << "\n";
  cout << "  The spatial dimension DIM_NUM =     " << dim_num      << "\n";
  cout << "  The number of points N =         " << n         << "\n";
  cout << "  The number of sample points NS = " << ns        << "\n";
  cout << "  The random number SEED_INIT =    " << seed_init << "\n";
  cout << "\n";

  z = dtable_data_read ( input_filename, dim_num, n );

  r8mat_transpose_print_some ( dim_num, n, z, 1, 1, 5, 5, 
    "  5x5 portion of data read from file:" );

  quality_test01 ( dim_num, n, z, ns, seed_init );
  quality_test02 ( dim_num, n, z, ns, seed_init );
  quality_test03 ( dim_num, n, z, ns, seed_init );
  quality_test04 ( dim_num, n, z, ns, seed_init );

  delete [] z;
//
//  Terminate.
//
  cout << "\n";
  cout << "FSU_QUALITY_TEST\n";
  cout << "  Normal end of execution.\n";
  cout << "\n";
  timestamp ( );

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

void quality_test01 ( int dim_num, int n, double z[], int ns, int seed_init )

//****************************************************************************80
//
//  Purpose:
//
//    QUALITY_TEST01 tests CHI_MEASURE.
//
//  Modified:
//
//    10 November 2006
//
//  Author:
//
//    Max Gunzburger
//    John Burkardt
//
//  License:
//
//    Copyright (C) 2004  John Burkardt and Max Gunzburger
//
//    This library is free software; you can redistribute it and/or
//    modify it under the terms of the GNU Lesser General Public
//    License as published by the Free Software Foundation; either
//    version 2.1 of the License, or (at your option) any later version.
//
//    This library 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
//    Lesser General Public License for more details.
//
//    You should have received a copy of the GNU Lesser General Public
//    License along with this library; if not, write to the Free Software
//    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
{
  cout << "\n";
  cout << "QUALITY_TEST01\n";
  cout << "  CHI_MEASURE computes the CHI measure of quality.\n";
  cout << "  The regularity measure         Chi = "
       << chi_measure ( dim_num, n, z, ns, seed_init ) << "\n";

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

void quality_test02 ( int dim_num, int n, double z[], int ns, int seed_init )

//****************************************************************************80
//
//  Purpose:
//
//    QUALITY_TEST02 tests D_MEASURE.
//
//  Modified:
//
//    10 November 2006
//
//  Author:
//
//    Max Gunzburger
//    John Burkardt
//
//  License:
//
//    Copyright (C) 2004  John Burkardt and Max Gunzburger
//
//    This library is free software; you can redistribute it and/or
//    modify it under the terms of the GNU Lesser General Public
//    License as published by the Free Software Foundation; either
//    version 2.1 of the License, or (at your option) any later version.
//
//    This library 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
//    Lesser General Public License for more details.
//
//    You should have received a copy of the GNU Lesser General Public
//    License along with this library; if not, write to the Free Software
//    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
{
  cout << "\n";
  cout << "QUALITY_TEST02\n";
  cout << "  D_MEASURE computes the D measure of quality.\n";
  cout << "  2nd moment determinant measure   D = "
       << d_measure ( dim_num, n, z, ns, seed_init ) << "\n";

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

void quality_test03 ( int dim_num, int n, double z[], int ns, int seed_init )

//****************************************************************************80
//
//  Purpose:
//
//    QUALITY_TEST03 tests H_MEASURE.
//
//  Modified:
//
//    10 November 2006
//
//  Author:
//
//    Max Gunzburger
//    John Burkardt
//
//  License:
//
//    Copyright (C) 2004  John Burkardt and Max Gunzburger
//
//    This library is free software; you can redistribute it and/or
//    modify it under the terms of the GNU Lesser General Public
//    License as published by the Free Software Foundation; either
//    version 2.1 of the License, or (at your option) any later version.
//
//    This library 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
//    Lesser General Public License for more details.
//
//    You should have received a copy of the GNU Lesser General Public
//    License along with this library; if not, write to the Free Software
//    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
{
  cout << "\n";
  cout << "QUALITY_TEST03\n";
  cout << "  H_MEASURE computes the H measure of quality.\n";
  cout << "  The point distribution norm      H = "
       << h_measure ( dim_num, n, z, ns, seed_init ) << "\n";

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

void quality_test04 ( int dim_num, int n, double z[], int ns, int seed_init )

//****************************************************************************80
//
//  Purpose:
//
//    QUALITY_TEST04 tests TAU_MEASURE.
//
//  Modified:
//
//    10 November 2006
//
//  Author:
//
//    Max Gunzburger
//    John Burkardt
//
//  License:
//
//    Copyright (C) 2004  John Burkardt and Max Gunzburger
//
//    This library is free software; you can redistribute it and/or
//    modify it under the terms of the GNU Lesser General Public
//    License as published by the Free Software Foundation; either
//    version 2.1 of the License, or (at your option) any later version.
//
//    This library 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
//    Lesser General Public License for more details.
//
//    You should have received a copy of the GNU Lesser General Public
//    License along with this library; if not, write to the Free Software
//    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
{
  cout << "\n";
  cout << "QUALITY_TEST04\n";
  cout << "  TAU_MEASURE computes the Tau measure of quality.\n";
  cout << "  2nd moment trace measure       Tau = "
       << tau_measure ( dim_num, n, z, ns, seed_init ) << "\n";

  return;
}