# include <cstdlib>
# include <cmath>
# include <ctime>
# include <iostream>
# include <iomanip>
# include <fstream>
# include <string>

using namespace std;

namespace triangle_quality {

  int main ( int argc, char *argv[] );
  void alpha_measure ( int n, double z[], int triangle_order, int triangle_num,
                       int triangle_node[], double *alpha_min, double *alpha_ave, 
                       double *alpha_area );
  double arc_cosine ( double c );
  void area_measure ( int n, double z[], int triangle_order, int triangle_num,
                      int triangle_node[], double &area_min, double &area_max, double &area_ratio,
                      double &area_ave, double &area_std, int &area_negative, int &area_zero );
  void bandwidth_mesh ( int element_order, int element_num, int element_node[],
                        int *ml, int *mu, int *m );
  char ch_cap ( char c );
  bool ch_eqi ( char c1, char c2 );
  int ch_to_digit ( char c );
  int file_column_count ( string filename );
  int file_row_count ( string filename );
  int i4_max ( int i1, int i2 );
  int i4_min ( int i1, int i2 );
  int *i4mat_data_read ( string input_filename, int m, int n );
  void i4mat_header_read ( string input_filename, int *m, int *n );
  void i4mat_transpose_print ( int m, int n, int a[], string title );
  void i4mat_transpose_print_some ( int m, int n, int a[], int ilo, int jlo, 
                                    int ihi, int jhi, string title );
  void mesh_base_zero ( int node_num, int element_order, 
                        int element_num, int element_node[] );
  void q_measure ( int n, const double z[], int triangle_order, int triangle_num,
                   const int triangle_node[], double *q_min, double *q_max, double *q_ave, 
                   double *q_area, double * q_vec = NULL, int offset=1 );
  double r8_abs ( double x );
  double r8_huge ( );
  double r8_max ( double x, double y );
  double r8_min ( double x, double y );
  double *r8mat_data_read ( string input_filename, int m, int n );
  void r8mat_header_read ( string input_filename, int *m, int *n );
  void r8mat_transpose_print ( int m, int n, double a[], string title );
  void r8mat_transpose_print_some ( int m, int n, double a[], int ilo, int jlo, 
                                    int ihi, int jhi, string title );
  int s_len_trim ( string s );
  int s_to_i4 ( string s, int *last, bool *error );
  bool s_to_i4vec ( string s, int n, int ivec[] );
  double s_to_r8 ( string s, int *lchar, bool *error );
  bool s_to_r8vec ( string s, int n, double rvec[] );
  int s_word_count ( string s );
  void timestamp ( );
}