program main

!*****************************************************************************80
!
!! MAIN is the main program for LATIN_COVER_TEST.
!
!  Discussion:
!
!    LATIN_COVER_TEST tests the LATIN_COVER library.
!
!  Licensing:
!
!    This code is distributed under the GNU LGPL license.
!
!  Modified:
!
!    24 June 2012
!
!  Author:
!
!    John Burkardt
!
  implicit none

  call timestamp ( )
  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'LATIN_COVER_TEST:'
  write ( *, '(a)' ) '  FORTRAN90 version'
  write ( *, '(a)' ) '  Test the LATIN_COVER library.'

  call test01 ( )
  call test02 ( )
  call test03 ( )
!
!  Terminate.
!
  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'LATIN_COVER_TEST:'
  write ( *, '(a)' ) '  Normal end of execution.'
  write ( *, '(a)' ) ' '
  call timestamp ( )

  stop 0
end
subroutine test01 ( )

!*****************************************************************************80
!
!! TEST01 tests LATIN_COVER.
!
!  Licensing:
!
!    This code is distributed under the GNU LGPL license.
!
!  Modified:
!
!    24 June 2012
!
!  Author:
!
!    John Burkardt
!
  implicit none

  integer ( kind = 4 ), allocatable :: a(:,:)
  integer ( kind = 4 ), parameter :: base = 1
  integer ( kind = 4 ) n
  integer ( kind = 4 ), allocatable :: p(:)
  integer ( kind = 4 ) seed
  integer ( kind = 4 ) test

  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'TEST01'
  write ( *, '(a)' ) '  LATIN_COVER:'

  do n = 3, 9, 2

    allocate ( p(1:n) )
    allocate ( a(1:n,1:n) )

    seed = 123456789

    do test = 1, 3

      call perm_uniform ( n, base, seed, p )
 
      call perm_print ( n, p, '  Permutation' )

      call latin_cover ( n, p, a )

      call i4mat_print ( n, n, a, '  Latin cover' )

    end do

    deallocate ( a )
    deallocate ( p )

  end do

  return
end
subroutine test02 ( )

!*****************************************************************************80
!
!! TEST02 tests LATIN_COVER_2D.
!
!  Licensing:
!
!    This code is distributed under the GNU LGPL license.
!
!  Modified:
!
!    24 June 2012
!
!  Author:
!
!    John Burkardt
!
  implicit none

  integer ( kind = 4 ), allocatable :: a(:,:)
  integer ( kind = 4 ), parameter :: base = 1
  integer ( kind = 4 ) n
  integer ( kind = 4 ), allocatable :: p(:,:)
  integer ( kind = 4 ) seed
  integer ( kind = 4 ) test

  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'TEST02'
  write ( *, '(a)' ) '  LATIN_COVER_2D:'

  do n = 3, 9, 2

    allocate ( p(2,1:n) )
    allocate ( a(1:n,1:n) )

    seed = 123456789

    do test = 1, 3

      call perm_uniform ( n, base, seed, p(1,1:n) )
 
      call perm_print ( n, p(1,1:n), '  Permutation 1' )

      call perm_uniform ( n, base, seed, p(2,1:n) ) 
 
      call perm_print ( n, p(2,1:n), '  Permutation 2' )
      call latin_cover_2d ( n, p, a )

      call i4mat_print ( n, n, a, '  Latin cover' )

    end do

    deallocate ( a )
    deallocate ( p )

  end do

  return
end
subroutine test03 ( )

!*****************************************************************************80
!
!! TEST03 tests LATIN_COVER_3D.
!
!  Licensing:
!
!    This code is distributed under the GNU LGPL license.
!
!  Modified:
!
!    24 June 2012
!
!  Author:
!
!    John Burkardt
!
  implicit none

  integer ( kind = 4 ), allocatable :: a(:,:,:)
  integer ( kind = 4 ), parameter :: base = 1
  integer ( kind = 4 ) n
  integer ( kind = 4 ), allocatable :: p(:,:)
  integer ( kind = 4 ) seed
  integer ( kind = 4 ) test

  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'TEST03'
  write ( *, '(a)' ) '  LATIN_COVER_3D:'

  do n = 3, 9, 2

    allocate ( p(3,1:n) )
    allocate ( a(1:n,1:n,1:n) )

    seed = 123456789

    do test = 1, 3

      call perm_uniform ( n, base, seed, p(1,1:n) )
 
      call perm_print ( n, p(1,1:n), '  Permutation 1' )

      call perm_uniform ( n, base, seed, p(2,1:n) ) 
 
      call perm_print ( n, p(2,1:n), '  Permutation 2' )

      call perm_uniform ( n, base, seed, p(3,1:n) ) 
 
      call perm_print ( n, p(3,1:n), '  Permutation 3' )

      call latin_cover_3d ( n, p, a )

      call i4block_print ( n, n, n, a, '  Latin cover' )

    end do

    deallocate ( a )
    deallocate ( p )

  end do

  return
end