# include <math.h>
# include <stdio.h>
# include <stdlib.h>
# include <string.h>

# include "log_normal.h"

int main ( );

void log_normal_cdf_test ( );
void log_normal_sample_test ( );

/******************************************************************************/

int main ( )

/******************************************************************************/
/*
  Purpose:

    log_normal_test() tests log_normal().

  Licensing:

    This code is distributed under the MIT license.

  Modified:

    22 March 2016

  Author:

    John Burkardt
*/
{
  timestamp ( );
  printf ( "\n" );
  printf ( "LOG_NORMAL_TEST\n" );
  printf ( "  C version\n" );
  printf ( "  Test the LOG_NORMAL library.\n" );

  log_normal_cdf_test ( );
  log_normal_sample_test ( );
/*
  Terminate.
*/
  printf ( "\n" );
  printf ( "LOG_NORMAL_TEST\n" );
  printf ( "  Normal end of execution.\n" );
  printf ( "\n" );
  timestamp ( );

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

void log_normal_cdf_test ( )

/******************************************************************************/
/*
  Purpose:

    LOG_NORMAL_CDF_TEST tests LOG_NORMAL_CDF.

  Licensing:

    This code is distributed under the MIT license.

  Modified:

    27 March 2016

  Author:

    John Burkardt
*/
{
  double cdf;
  int i;
  double mu;
  double pdf;
  int seed = 123456789;
  double sigma;
  double x;
  double x2;

  printf ( "\n" );
  printf ( "LOG_NORMAL_CDF_TEST\n" );
  printf ( "  LOG_NORMAL_CDF evaluates the Log Normal CDF;\n" );
  printf ( "  LOG_NORMAL_CDF_INV inverts the Log Normal CDF.\n" );
  printf ( "  LOG_NORMAL_PDF evaluates the Log Normal PDF;\n" );

  mu = 10.0;
  sigma = 2.25;

  printf ( "\n" );
  printf ( "  PDF parameter MU =     %g\n", mu );
  printf ( "  PDF parameter SIGMA =  %g\n", sigma );

  if ( ! log_normal_check ( mu, sigma ) )
  {
    printf ( "\n" );
    printf ( "LOG_NORMAL_CDF_TEST - Fatal error!\n" );
    printf ( "  The parameters are not legal.\n" );
    return;
  }

  printf ( "\n" );
  printf ( "       X            PDF           CDF            CDF_INV\n" );
  printf ( "\n" );

  for ( i = 1; i <= 10; i++ )
  {
    x = log_normal_sample ( mu, sigma, &seed );
    pdf = log_normal_pdf ( x, mu, sigma );
    cdf = log_normal_cdf ( x, mu, sigma );
    x2 = log_normal_cdf_inv ( cdf, mu, sigma );

    printf ( "  %12g  %12g  %12g  %12g\n", x, pdf, cdf, x2 );
  }

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

void log_normal_sample_test ( )

/******************************************************************************/
/*
  Purpose:

    LOG_NORMAL_SAMPLE_TEST tests LOG_NORMAL_SAMPLE.

  Licensing:

    This code is distributed under the MIT license.

  Modified:

    15 September 2013

  Author:

    John Burkardt
*/
{
# define SAMPLE_NUM 1000

  int i;
  double mean;
  double mu;
  int seed = 123456789;
  double sigma;
  double variance;
  double x[SAMPLE_NUM];
  double xmax;
  double xmin;

  printf ( "\n" );
  printf ( "LOG_NORMAL_SAMPLE_TEST\n" );
  printf ( "  LOG_NORMAL_MEAN computes the Log Normal mean;\n" );
  printf ( "  LOG_NORMAL_SAMPLE samples the Log Normal distribution;\n" );
  printf ( "  LOG_NORMAL_VARIANCE computes the Log Normal variance;\n" );

  mu = 1.0;
  sigma = 2.0;

  printf ( "\n" );
  printf ( "  PDF parameter MU =      %g\n", mu );
  printf ( "  PDF parameter SIGMA =   %g\n", sigma );

  if ( ! log_normal_check ( mu, sigma ) )
  {
    printf ( "\n" );
    printf ( "LOG_NORMAL_SAMPLE_TEST - Fatal error!\n" );
    printf ( "  The parameters are not legal.\n" );
    return;
  }

  mean = log_normal_mean ( mu, sigma);
  variance = log_normal_variance ( mu, sigma );

  printf ( "\n" );
  printf ( "  PDF mean =     %g\n", mean     );
  printf ( "  PDF variance = %g\n", variance );

  for ( i = 0; i < SAMPLE_NUM; i++ )
  {
    x[i] = log_normal_sample ( mu, sigma, &seed );
  }

  mean = r8vec_mean ( SAMPLE_NUM, x );
  variance = r8vec_variance ( SAMPLE_NUM, x );
  xmax = r8vec_max ( SAMPLE_NUM, x );
  xmin = r8vec_min ( SAMPLE_NUM, x );

  printf ( "\n" );
  printf ( "  Sample size =     %d\n", SAMPLE_NUM  );
  printf ( "  Sample mean =     %g\n", mean     );
  printf ( "  Sample variance = %g\n", variance );
  printf ( "  Sample maximum =  %g\n", xmax     );
  printf ( "  Sample minimum =  %g\n", xmin     );

  return;
# undef SAMPLE_NUM
}