funit-12 0.12.1

data/bin/funit-12

@@ -0,0 +1,72 @@
+#! /usr/bin/env ruby
+
+# Main fUnit routine
+
+begin require 'rubygems'; rescue LoadError; end
+
+$:.push "#{File.dirname $0}/../lib"
+require 'funit'
+require 'optparse'
+require 'ostruct'
+
+include Funit
+
+funit_data = OpenStruct.new
+funit_data.prog_source_dirs = ['.']
+funit_data.c_code = Funit::C_compile.new
+  
+OptionParser.new do |opts|
+  opts.banner = <<-END_OF_HELP
+        To use fUnit, type:
+          funit [-options] [test_file_name(s)]
+        The argument(s) is optional. If no argument is given, then all the .fun files inside the working directory will be used.
+      END_OF_HELP
+      
+  opts.separator ""
+  opts.separator "Specific options:"
+  
+  opts.on("-s","--source <dir>",String,"Specify a directory for the non-test source") do |dir|
+    funit_data.prog_source_dirs << dir.chomp("/")
+  end
+  
+  opts.on("-l","--link_c_code <dir>",String,"Specify a directory of C/C++ code that must be compiled/linked") do |dir|
+    funit_data.c_code.compile_library(dir)
+  end
+  
+  opts.on_tail("--clean","To remove the files generated by funit") do 
+    funit_data.c_code.clean_c_code
+    Funit::clean_genFiles
+    exit
+  end
+  
+  opts.separator ""
+  opts.separator "Common options:"
+  
+  opts.on_tail("-v","--version","Print version") do 
+    puts "funit version #{Funit::VERSION}"
+    exit
+  end
+  
+ opts.on_tail("-h", "--help", "Detailed overview of funit options") do
+   puts opts
+   exit
+ end
+
+end.parse!
+
+funit_data.prog_source_dirs.flatten!
+funit_data.prog_source_dirs.uniq!
+
+Funit::check_for_FSFLAG
+
+Funit::run_tests(funit_data)
+
+#--
+# Copyright 2006-2007 United States Government as represented by
+# NASA Langley Research Center. No copyright is claimed in
+# the United States under Title 17, U.S. Code. All Other Rights
+# Reserved.
+#
+# This file is governed by the NASA Open Source Agreement.
+# See License.txt for details.
+#++

data/examples/CFD/FluxFunctions.f90

@@ -0,0 +1,34 @@
+module FluxFunctions
+
+  implicit none
+
+contains
+
+  subroutine CentralFlux( LeftState, RightState, InterfaceFlux )
+    real, intent(in)  :: leftState
+    real, intent(in)  :: rightState
+    real, intent(out) :: interfaceFlux
+    interfaceFlux = 0.5*(Flux(leftState)+Flux(rightState))
+  end subroutine CentralFlux
+
+  subroutine RoeFlux( LeftState, RightState, InterfaceFlux )
+    real, intent(in)  :: leftState
+    real, intent(in)  :: rightState
+    real, intent(out) :: interfaceFlux
+    interfaceFlux = 0.5*(Flux(leftState)+Flux(rightState)) &
+                  - 0.5*RoeAvg(leftState,rightState)*(rightState-leftState)
+  end subroutine RoeFlux
+
+  function Flux( state )
+    real             :: Flux
+    real, intent(in) :: state
+    Flux = 0.5*state**2
+  end function Flux
+
+  function RoeAvg( leftState, rightState )
+    real             :: RoeAvg
+    real, intent(in) :: leftState, rightState
+    RoeAvg = 0.5*(leftState+rightState)
+  end function RoeAvg
+
+end module FluxFunctions

data/examples/CFD/FluxFunctions.fun

@@ -0,0 +1,51 @@
+test_suite FluxFunctions
+
+real :: leftState, rightState, interfaceFlux
+
+setup
+  leftState  = 0
+  rightState = 1
+end setup
+
+test FluxZero
+  real :: state
+  state = 0
+  Assert_Equal_Within( 0, Flux(state), 0.00001 )
+end test
+ 
+test FluxOne
+  real :: state = 1
+  Assert_Equal_Within( 0.5, Flux(state), 0.00001 )
+end test
+
+test RoeAvgZero
+  Assert_Real_Equal( 0, RoeAvg(0.0,0.0) )
+  Assert_False( RoeAvg(0.0,0.0)==1 )
+end test
+
+test RoeAvgKnown
+  Assert_Real_Equal( 0.5, RoeAvg(leftState,rightState) )
+  Assert_True( RoeAvg(leftState,rightState) > 0 )
+end test
+
+test CentralFluxKnown
+  call CentralFlux( leftState, rightState, interfaceFlux )
+  Assert_Equal_Within( 0.25, interfaceFlux, 0.001 )
+  Assert_Equal_Within( 0.25, interfaceFlux, 0.00000001 )
+  Assert_Equal( 0.25, interfaceFlux )
+end test
+
+test RoeFluxExpansionShock
+  leftState = -1
+  call RoeFlux( leftState, rightState, interfaceFlux )
+  Assert_Equal( 0.5, interfaceFlux )
+end test
+
+test RoeFluxZero
+  rightState = 0
+  call RoeFlux( leftState, rightState, interfaceFlux )
+  Assert_Real_Equal( 0, interfaceFlux )
+  Assert_Equal( 0, interfaceFlux )
+end test
+
+end test_suite

data/examples/CFD/Gammas.f90

@@ -0,0 +1,7 @@
+module Gammas
+
+  implicit none
+  
+  real, parameter :: Gamma = 1.4
+
+end module Gammas

data/examples/CFD/GasModel.f90

@@ -0,0 +1,16 @@
+module GasModel
+
+ use Gammas
+
+  implicit none
+
+contains
+
+  subroutine PerfectP (Density, Energy, Pressure)
+    real, intent(in)  :: Density 
+    real, intent(in)  :: Energy
+    real, intent(out) :: Pressure
+    Pressure = Density * Energy * ( Gamma - 1.0 )
+  end subroutine PerfectP
+
+end module GasModel

data/examples/CFD/GasModel.fun

@@ -0,0 +1,24 @@
+test_suite GasModel
+
+  real            :: Pressure, Density, Energy
+
+test PerfectPZeroed
+  real, parameter :: zero = 0
+  call PerfectP (zero, zero, Pressure)
+  Assert_Real_Equal( 0, Pressure )
+  Assert_Equal_within( 0, Pressure, 0.0000000001 )
+end test
+
+test Warbler
+end test
+
+test PerfectPKnown
+  real :: Density = 1
+  Energy  = 1
+  call PerfectP( Density, Energy, Pressure )
+  Assert_Real_Equal( 0.4, Pressure )
+  Assert_True( Pressure > 0 )
+  Assert_False( Pressure < 0 )
+end test
+
+end test_suite

data/examples/ReadData/time_series_data.f90

@@ -0,0 +1,27 @@
+module time_series_data
+  implicit none
+  integer, parameter :: MAX_POINTS = 100000 ! or use allocate
+  type date_time
+    integer :: year, month, day, hour, minute
+  end type date_time
+  type time_series
+    type(date_time) :: date_time
+    real            :: value      = 0.0
+  end type time_series
+  type(time_series), dimension(MAX_POINTS), save :: ts_data
+contains
+  subroutine read_time_series( filename )
+    character(len=*), intent(in) :: filename
+    integer :: i, ios
+    open(10, file=filename, iostat=ios)
+    if (ios/=0) then
+      print *, 'failed to open file: >>', filename, '<<'
+    else
+      do i = 1, MAX_POINTS
+        read( 10, fmt='(i4,i2,i2,i2,i2,e20.12)', end=20 ) ts_data(i)
+      end do
+      print *, 'quit reading data after ', MAX_POINTS, ' points'
+20    continue
+    end if
+  end subroutine read_time_series
+end module time_series_data

data/examples/ReadData/time_series_data.fun

@@ -0,0 +1,33 @@
+test_suite time_series_data
+
+  character(len=10), parameter :: FILE = 'values.txt'
+
+setup
+  open(8, file=FILE)
+    write(8,'(a)'), '200609300000 0.223200546265E+003'
+    write(8,'(a)'), '200609300132 0.226001495361E+003'
+  close(8)
+end setup
+
+test load_time_series_data_from_file
+  call read_time_series( FILE )
+  Assert_Equal(         2006, ts_data(1)%date_time%year )
+  Assert_Equal(            9, ts_data(1)%date_time%month )
+  Assert_Equal(           30, ts_data(1)%date_time%day )
+  Assert_Equal(            0, ts_data(1)%date_time%hour )
+  Assert_Equal(            0, ts_data(1)%date_time%minute )
+  Assert_Equal_Within( 223.2, ts_data(1)%value, 0.1 )
+  Assert_Equal(         2006, ts_data(2)%date_time%year )
+  Assert_Equal(            9, ts_data(2)%date_time%month )
+  Assert_Equal(           30, ts_data(2)%date_time%day )
+  Assert_Equal(            1, ts_data(2)%date_time%hour )
+  Assert_Equal(           32, ts_data(2)%date_time%minute )
+  Assert_Equal_Within( 226.0, ts_data(2)%value, 0.1 )
+end test
+
+teardown
+  call system('/bin/rm '//FILE)
+end teardown
+
+end test_suite
+

data/examples/StopWatch/StopWatch.f90

@@ -0,0 +1,50 @@
+module StopWatch
+
+  implicit none
+
+  public
+
+  logical :: notInitialized = .TRUE.
+
+  integer, dimension(8) :: last
+
+  contains
+
+  real function SecBetween(beginDnT, endDnT)
+    integer, dimension(8), intent(in)  :: beginDnT, endDnT
+    real :: days, hours, minutes, seconds
+    integer, parameter :: yr=1, mo=2, day=3, utc=4, hr=5, mn=6, s=7, ms=8
+    
+    continue
+
+    if ( endDnT(day) == beginDnT(day) ) then
+       days = 0
+    else
+       days = 1 ! note: assuming one day
+    endif
+
+    hours   = endDnT(hr) - beginDnT(hr) + 24*days
+    minutes = endDnT(mn) - beginDnT(mn) + 60*hours
+    seconds = endDnT(s)  - beginDnT(s)  + 60*minutes
+
+    SecBetween = seconds + ( endDnT(ms) - beginDnT(ms) ) / 1000.
+
+  end function SecBetween
+    
+  real function secSinceLast()
+
+    integer, dimension(8) :: now
+
+    if (notInitialized) then
+      notInitialized = .FALSE.
+      secSinceLast = 0.0
+      call date_and_time(values=last)
+    else
+      call date_and_time(values=now)
+      secSinceLast = secBetween(last, now)
+      last = now
+    endif
+
+  end function secSinceLast
+
+end module StopWatch

data/examples/StopWatch/StopWatch.fun

@@ -0,0 +1,77 @@
+test_suite StopWatch
+
+integer, dimension(8) :: dateAndTime1, dateAndTime2
+real :: seconds
+
+setup
+  NotInitialized = .TRUE.
+  last    = 0
+  seconds = HUGE(0.0)
+end setup
+
+test SystemDateAndTimeWorks
+ call date_and_time(values=dateAndTime1)
+ Assert_True( dateAndTime1(1) /= -huge(0) )
+ Assert_True( size(dateAndTime1,1) == 8 )
+end test
+
+! test secBetween
+test OneMSecDifference
+  dateAndTime1 = (/ 2000, 1, 1, 0, 0, 0, 0, 0 /)
+  dateAndTime2 = (/ 2000, 1, 1, 0, 0, 0, 0, 1 /)
+  seconds = SecBetween(dateAndTime1, dateAndTime2)
+  Assert_Real_Equal( 0.001, seconds)
+end test
+
+test MinuteRollover
+  dateAndTime1 = (/ 2000, 1, 1, 0, 0, 0,59, 0 /)
+  dateAndTime2 = (/ 2000, 1, 1, 0, 0, 1, 0, 0 /)
+  seconds = SecBetween(dateAndTime1, dateAndTime2)
+  Assert_Real_Equal( 1.0, seconds )
+end test
+
+! test secSinceLast
+test InitializationState
+  Assert_True(notInitialized)
+  seconds = secSinceLast()
+  Assert_False(notInitialized)
+  seconds = secSinceLast()
+  Assert_False(notInitialized)
+end test
+
+test InitiallyReturnsZero
+  seconds = secSinceLast()
+  Assert_Real_Equal( 0.0, seconds )
+  call timeDelay(seconds)
+  seconds = secSinceLast()
+  Assert_True( seconds /= 0.0 )
+end test
+
+subroutine timeDelay (sum)
+  integer :: i
+  real    :: sum
+  do i = 1, 1000000
+   sum = sum + i
+  enddo
+end subroutine timeDelay
+
+test ComputesSeconds
+  seconds = secSinceLast()
+  call timeDelay (seconds)
+  seconds = secSinceLast()
+  Assert_True( seconds > 0.0 )
+end test
+
+test ComputesSecondsSpecial
+  real :: expectedSeconds
+
+  seconds = secSinceLast()
+  dateAndTime1 = last
+  call timeDelay (seconds)
+  seconds = secSinceLast()
+  dateAndTime2 = last
+  expectedSeconds = secBetween(dateAndTime1,dateAndTime2)
+  Assert_Real_Equal( expectedSeconds, seconds )
+end test
+
+end test_suite

data/lib/funit.rb

@@ -0,0 +1,87 @@
+begin
+  require 'rubygems'
+  require 'fortran_dependencies'
+rescue LoadError
+  STDERR.puts "gem install fortran"
+  exit 1
+end
+
+require 'funit/compiler'
+require 'funit/functions'
+require 'funit/assertions'
+require 'funit/testsuite'
+require 'funit/c_tools'
+require 'fileutils'
+
+module Funit
+
+  VERSION = '0.12.1'
+
+  ##
+  # run all tests
+
+  def run_tests(funit_data)
+    Compiler.new# a test for compiler env set (FIXME: remove this later)
+    write_test_runner( test_files = parse_command_line )
+    test_suites = []
+    test_files.each{ |test_file|
+      tf_content = IO.read(test_file+'.fun')
+      tf_content.scan(/test_suite\s+(\w+)(.*?)end\s+test_suite/m).each{|ts|
+        ts_name = $1
+        ts_content = $2
+        if((!File.exist?(ts_name+"_fun.f90")) || File.mtime(ts_name+"_fun.f90") < File.mtime(test_file+".fun")) then
+          if ( File.read(ts_name+'.f90').match(/\s*module\s+#{ts_name}/i) ) then
+            TestSuite.new(ts_name, ts_content, false)
+          else
+            TestSuite.new(ts_name, ts_content, true)
+          end
+        end
+        test_suites.push(ts_name)
+      }
+    }
+    compile_tests(test_suites,funit_data)
+    exit 1 unless system "env PATH='.' TestRunner"
+  end
+
+  ##
+  # remove files generated by fUnit
+
+  def clean_genFiles
+    module_names = Dir["**/*.fun"].map{|mn| mn.chomp(".fun")}
+    tbCancelled = module_names.map{|mn| mn+"_fun."} + ["TestRunner."]
+    tbCancelled = tbCancelled.map{|tbc| [tbc+"f90",tbc+"o",tbc+"MOD"]}.flatten
+    tbCancelled += Dir["**/TestRunner"]
+    tbCancelled += Dir["**/makeTestRunner"]
+    tbCancelled = (tbCancelled+tbCancelled.map{|tbc| tbc.downcase}).uniq
+    FileUtils.rm_f(tbCancelled)
+  end
+
+  def check_for_FSFLAG
+    if(ENV['FSFLAG'].nil?) then
+      puts <<-EOF
+No environment variable FSFLAG set.
+
+For example for most compilers such as gfortran you will need: -I
+      sh: export FSFLAG=-I
+     csh: setenv FSFLAG -I
+ windows: set FSFLAG=-I
+
+but for some such as Sun's f95 compiler you will need: -M
+      sh: export FSFLAG=-M
+     csh: setenv FSFLAG -M
+ windows: set FSFLAG=-M
+    EOF
+    exit
+    end
+  end
+end
+
+#--
+# Copyright 2006-2007 United States Government as represented by
+# NASA Langley Research Center. No copyright is claimed in
+# the United States under Title 17, U.S. Code. All Other Rights
+# Reserved.
+#
+# This file is governed by the NASA Open Source Agreement.
+# See License.txt for details.
+#++