generate_fftw_wisdom.fpp Source File

!> A program that calls init_transforms, to be used for generating fftw wisdom 
!! for a given input file.
!!
!! This is free software released under the MIT licence
!!   Written by: Edmund Highcock (edmundhighcock@users.sourceforge.net)
program generate_fftw_wisdom
  !use unit_tests
  use mp, only: init_mp, finish_mp, proc0, broadcast, nproc
  use file_utils, only: init_file_utils, run_name
  use nonlinear_terms, only: init_nonlinear_terms, finish_nonlinear_terms
  use job_manage, only: time_message
  use benchmarks, only: benchmark_identifier
  use constants, only: run_name_size
  implicit none
  character(run_name_size), target :: cbuff
  logical :: dummy=.false.
  real :: time_taken(2) = 0.0

  ! Set up depenencies
  call init_mp
  if (proc0) call init_file_utils(dummy, name="gs")
  if (proc0) then
     cbuff = trim(run_name)
  end if

  call broadcast (cbuff)
  if (.not. proc0) run_name = cbuff

  call time_message(.false., time_taken, "Init transforms")

  call init_nonlinear_terms

  call time_message(.false., time_taken, "Init transforms")

  if (proc0) then
     write(*, '(" Time taken to calculate fftw plans ",I6," procs: ",F3.1," s")') nproc, time_taken(1)
     write(*,*)
     !    call append_output_file(timing_unit, &
     !      benchmark_identifier())
     !    write(timing_unit, '(I6,"   ",F9.3)') nproc, time_taken(1)
     !    call close_output_file(timing_unit)
  end if

  !call finish_nonlinear_terms


  call finish_mp

end program generate_fftw_wisdom