E.2.9 Statistics Class Unit Test Program

This lightly commented program performs a unit test on the Statistics Class.

program Unit_Test

  use Caesar_Data_Structures_Module
  use Caesar_Statistics_Class
  use Caesar_Numbers_Module, only: one
  implicit none

  type(Statistics_type) :: Variable_Stats
  type(Communication_type) :: Comm

  ! Initialize communications.

  call Initialize (Comm)
  call Output (Comm)
  if (this_is_IO_PE) write (6,*)

  ! Testing statements.

  ! Odd Numbers.

  call Initialize (Variable_Stats, 'Odd Numbers')

  call Add_Value (Variable_Stats, 1.0d0)
  call Add_Value (Variable_Stats, 3.0d0)
  call Add_Value (Variable_Stats, 5.0d0)
  call Add_Value (Variable_Stats, 7.0d0)
  call Add_Value (Variable_Stats, 9.0d0)

  call Output (Variable_Stats, Global=.false., Verbose=.true.)
  if (this_is_IO_PE) write (6,*)
  call Output (Variable_Stats, Global=.true.,  Verbose=.true.)
  if (this_is_IO_PE) write (6,*)

  VERIFY(Valid_State(Variable_Stats),0)

  call Finalize (Variable_Stats)

  ! Powers of 2.

  call Initialize (Variable_Stats, 'Powers of 2')

  call Add_Value (Variable_Stats, 1.0d0)
  call Add_Value (Variable_Stats, 2.0d0)
  call Add_Value (Variable_Stats, 4.0d0)
  call Add_Value (Variable_Stats, 8.0d0)
  call Add_Value (Variable_Stats, 16.0d0)
  call Add_Value (Variable_Stats, 32.0d0)
  call Add_Value (Variable_Stats, 64.0d0)

  call Output (Variable_Stats, Global=.false., Verbose=.true.)
  if (this_is_IO_PE) write (6,*)
  call Output (Variable_Stats, Global=.true.,  Verbose=.true.)
  if (this_is_IO_PE) write (6,*)

  VERIFY(Valid_State(Variable_Stats),0)

  call Finalize (Variable_Stats)

  ! Reciprocals.

  call Initialize (Variable_Stats, 'Reciprocals')

  call Add_Value (Variable_Stats, one/1.0d0)
  call Add_Value (Variable_Stats, one/2.0d0)
  call Add_Value (Variable_Stats, one/3.0d0)
  call Add_Value (Variable_Stats, one/4.0d0)
  call Add_Value (Variable_Stats, one/5.0d0)
  call Add_Value (Variable_Stats, one/6.0d0)
  call Add_Value (Variable_Stats, one/7.0d0)
  call Add_Value (Variable_Stats, one/8.0d0)
  call Add_Value (Variable_Stats, one/9.0d0)
  call Add_Value (Variable_Stats, one/55.0d0)

  call Output (Variable_Stats, Global=.false., Verbose=.true.)
  if (this_is_IO_PE) write (6,*)
  call Output (Variable_Stats, Global=.true.,  Verbose=.true.)
  if (this_is_IO_PE) write (6,*)

  VERIFY(Valid_State(Variable_Stats),0)

  call Finalize (Variable_Stats)

  ! Negatives.

  call Initialize (Variable_Stats, 'Negatives')

  call Add_Value (Variable_Stats, 1.0d0)
  call Add_Value (Variable_Stats, -1.0d0)
  call Add_Value (Variable_Stats, 2.0d0)
  call Add_Value (Variable_Stats, -2.0d0)
  call Add_Value (Variable_Stats, 3.0d0)
  call Add_Value (Variable_Stats, -3.0d0)

  call Output (Variable_Stats, Global=.false., Verbose=.true.)
  if (this_is_IO_PE) write (6,*)
  call Output (Variable_Stats, Global=.true.,  Verbose=.true.)
  if (this_is_IO_PE) write (6,*)

  VERIFY(Valid_State(Variable_Stats),0)

  call Finalize (Variable_Stats)

  ! Parallel test.

  call Initialize (Variable_Stats, 'Parallel Test')

  call Add_Value (Variable_Stats, changetype(real, this_PE))

  call Output (Variable_Stats, Global=.false., Verbose=.true.)
  if (this_is_IO_PE) write (6,*)
  call Output (Variable_Stats, Global=.true.,  Verbose=.true.)

  VERIFY(Valid_State(Variable_Stats),0)

  call Finalize (Variable_Stats)

  ! Finalize communications.

  call Finalize (Comm)

end