program main

c*********************************************************************72
c
cc test_interp_nd_test() tests test_interp_nd().
c
c  Licensing:
c
c    This code is distributed under the MIT license.
c
c  Modified:
c
c    30 August 2012
c
c  Author:
c
c    John Burkardt
c
      implicit none

      integer m
      integer n

      call timestamp (  )
      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'test_interp_nd_test():'
      write ( *, '(a)' ) '  FORTRAN77 version'
      write ( *, '(a)' ) '  Test test_interp_nd().'
      write ( *, '(a)' ) '  The R8LIB library is also needed.'

      call test01 ( )

      n = 10
      do m = 2, 4
        call test02 ( m, n )
      end do

      n = 2
      do m = 2, 4
        call test03 ( m, n )
      end do

      n = 10000
      m = 4
      call test04 ( m, n )
c
c  Terminate.
c
      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'test_interp_nd_test():'
      write ( *, '(a)' ) '  Normal end of execution.'
      write ( *, '(a)' ) ' '
      call timestamp ( )

      stop
      end
      subroutine test01 ( )

c*********************************************************************72
c
cc TEST01 prints the title of each test function.
c
c  Licensing:
c
c    This code is distributed under the MIT license.
c
c  Modified:
c
c    28 August 2012
c
c  Author:
c
c    John Burkardt
c
      implicit none

      integer prob
      integer prob_num
      character * ( 80 ) title

      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'TEST01'
      write ( *, '(a)' ) '  P00_TITLE returns the problem title.'

      call p00_prob_num ( prob_num )
      write ( *, '(a)' ) ' '
      write ( *, '(a,i2,a)' ) 
     &  '  There are a total of ', prob_num, ' problems.'
      write ( *, '(a)' ) ' '

      do prob = 1, prob_num
        call p00_title ( prob, title )
        write ( *, '(2x,i2,2x,a)' ) prob, '"' // trim ( title ) // '"'
      end do

      return
      end
      subroutine test02 ( m, n )

c*********************************************************************72
c
cc TEST02 samples each function in M dimensions, at N points.
c
c  Licensing:
c
c    This code is distributed under the MIT license.
c
c  Modified:
c
c    28 August 2012
c
c  Author:
c
c    John Burkardt
c
c  Parameters:
c
c    Input, integer M, the spatial dimension.
c
c    Input, integer N, the number of evaluation points.
c
      implicit none

      integer m
      integer n

      double precision c(m)
      double precision f(n)
      integer j
      integer prob
      integer prob_num
      integer seed
      double precision w(m)
      double precision x(m,n)

      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'TEST02'
      write ( *, '(a)' ) '  P00_F evaluates the function.'
      write ( *, '(a,i4)' ) '  Here, we use spatial dimension M = ', m
      write ( *, '(a,i6)' ) '  The number of points is N = ', n

      seed = 123456789
      call r8mat_uniform_01 ( m, n, seed, x )

      call p00_prob_num ( prob_num )

      do prob = 1, prob_num

        write ( *, '(a)' ) ' '
        write ( *, '(a,i8)' ) '  Problem  ', prob

        call p00_c ( prob, m, seed, c )
        call r8vec_print ( m, c, '  C parameters:' )

        call p00_w ( prob, m, seed, w )
        call r8vec_print ( m, w, '  W parameters:' )

        write ( *, '(a)' ) ' '
        write ( *, '(a)' ) '      F(X)              X(1)      X(2) ...'
        write ( *, '(a)' ) ' '

        call p00_f ( prob, m, c, w, n, x, f )

        do j = 1, n
          write ( *, '(2x,g14.6,2x,10f10.4)' ) f(j), x(1:m,j)
        end do

      end do

      return
      end
      subroutine test03 ( m, n )

c*********************************************************************72
c
cc TEST03 samples each derivative component in M dimensions, at N points.
c
c  Licensing:
c
c    This code is distributed under the MIT license.
c
c  Modified:
c
c    28 August 2012
c
c  Author:
c
c    John Burkardt
c
c  Parameters:
c
c    Input, integer M, the spatial dimension.
c
c    Input, integer N, the number of evaluation points.
c
      implicit none

      integer m
      integer n

      double precision c(m)
      double precision d(m,n)
      double precision f(n)
      integer id
      integer j
      integer prob
      integer prob_num
      integer seed
      double precision w(m)
      double precision x(m,n)

      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'TEST03'
      write ( *, '(a)' ) '  P00_D evaluates derivative components.'
      write ( *, '(a,i4)' ) '  Here, we use spatial dimension M = ', m
      write ( *, '(a,i6)' ) '  The number of points is N = ', n

      seed = 123456789
      call r8mat_uniform_01 ( m, n, seed, x )

      call p00_prob_num ( prob_num )

      do prob = 1, prob_num

        write ( *, '(a)' ) ' '
        write ( *, '(a,i8)' ) '  Problem  ', prob

        call p00_c ( prob, m, seed, c )
        call r8vec_print ( m, c, '  C parameters:' )

        call p00_w ( prob, m, seed, w )
        call r8vec_print ( m, w, '  W parameters:' )

        write ( *, '(a)' ) ' '
        write ( *, '(a)' ) '                         X(1)      X(2) ...'
        write ( *, '(a)' ) '      F(X)            dFdX(1)   dFdX(2) ...'

        call p00_f ( prob, m, c, w, n, x, f )

        do id = 1, m
          call p00_d ( prob, m, id, c, w, n, x, d(id,1:n) )
        end do

        do j = 1, n
          write ( *, '(a)' ) ' '
          write ( *, '(2x,14x,2x,10f10.4)' )       x(1:m,j)
          write ( *, '(2x,g14.6,2x,10f10.4)' ) f(j), d(1:m,j)
        end do

      end do

      return
      end
      subroutine test04 ( m, n )

c*********************************************************************72
c
cc TEST04 estimates integrals in M dimensions, using N points.
c
c  Licensing:
c
c    This code is distributed under the MIT license.
c
c  Modified:
c
c    30 August 2012
c
c  Author:
c
c    John Burkardt
c
c  Parameters:
c
c    Input, integer M, the spatial dimension.
c
c    Input, integer N, the number of evaluation points.
c
      implicit none

      integer m
      integer n

      double precision c(m)
      double precision f(n)
      integer prob
      integer prob_num
      double precision q1
      double precision q2
      double precision r8vec_sum
      integer seed
      double precision w(m)
      double precision x(m,n)

      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'TEST04'
      write ( *, '(a)' ) 
     &  '  P00_Q returns the integral of F over [0,1]^m.'
      write ( *, '(a,i4)' ) '  Here, we use spatial dimension M = ', m
      write ( *, '(a,i6)' ) '  The number of sample points is N = ', n

      seed = 123456789
      call r8mat_uniform_01 ( m, n, seed, x )

      call p00_prob_num ( prob_num )

      do prob = 1, prob_num

        write ( *, '(a)' ) ' '
        write ( *, '(a,i8)' ) '  Problem  ', prob

        call p00_c ( prob, m, seed, c )
        call r8vec_print ( m, c, '  C parameters:' )

        call p00_w ( prob, m, seed, w )
        call r8vec_print ( m, w, '  W parameters:' )

        write ( *, '(a)' ) ' '
        write ( *, '(a)' ) '      Exact Integral     Q'
        write ( *, '(a)' ) ' '

        call p00_q ( prob, m, c, w, q1 )

        call p00_f ( prob, m, c, w, n, x, f )
        q2 = r8vec_sum ( n, f ) / dble ( n )

        write ( *, '(2x,g14.6,2x,g14.6)' ) q1, q2

      end do

      return
      end
      subroutine timestamp ( )

c*********************************************************************72
c
cc timestamp() prints the YMDHMS date as a timestamp.
c
c  Licensing:
c
c    This code is distributed under the MIT license.
c
c  Modified:
c
c    12 June 2014
c
c  Author:
c
c    John Burkardt
c
      implicit none

      character * ( 8 ) ampm
      integer d
      character * ( 8 ) date
      integer h
      integer m
      integer mm
      character * ( 9 ) month(12)
      integer n
      integer s
      character * ( 10 ) time
      integer y

      save month

      data month /
     &  'January  ', 'February ', 'March    ', 'April    ', 
     &  'May      ', 'June     ', 'July     ', 'August   ', 
     &  'September', 'October  ', 'November ', 'December ' /

      call date_and_time ( date, time )

      read ( date, '(i4,i2,i2)' ) y, m, d
      read ( time, '(i2,i2,i2,1x,i3)' ) h, n, s, mm

      if ( h .lt. 12 ) then
        ampm = 'AM'
      else if ( h .eq. 12 ) then
        if ( n .eq. 0 .and. s .eq. 0 ) then
          ampm = 'Noon'
        else
          ampm = 'PM'
        end if
      else
        h = h - 12
        if ( h .lt. 12 ) then
          ampm = 'PM'
        else if ( h .eq. 12 ) then
          if ( n .eq. 0 .and. s .eq. 0 ) then
            ampm = 'Midnight'
          else
            ampm = 'AM'
          end if
        end if
      end if

      write ( *, 
     &  '(i2,1x,a,1x,i4,2x,i2,a1,i2.2,a1,i2.2,a1,i3.3,1x,a)' ) 
     &  d, trim ( month(m) ), y, h, ':', n, ':', s, '.', mm, 
     &  trim ( ampm )

      return
      end