funit 0.9.0

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,47 @@
+real :: leftState, rightState, interfaceFlux
+
+beginsetup
+  leftState  = 0
+  rightState = 1
+endsetup
+
+beginTest FluxZero
+  real :: state
+  state = 0
+  IsEqualWithin( 0, Flux(state), 0.00001 )
+endTest
+ 
+beginTest FluxOne
+  real :: state = 1
+  IsEqualWithin( 0.5, Flux(state), 0.00001 )
+endTest
+
+beginTest RoeAvgZero
+  IsRealEqual( 0, RoeAvg(0.0,0.0) )
+  IsFalse( RoeAvg(0.0,0.0)==1 )
+endTest
+
+beginTest RoeAvgKnown
+  IsRealEqual( 0.5, RoeAvg(leftState,rightState) )
+  IsTrue( RoeAvg(leftState,rightState) > 0 )
+endTest
+
+beginTest CentralFluxKnown
+  call CentralFlux( leftState, rightState, interfaceFlux )
+  IsEqualWithin( 0.25, interfaceFlux, 0.001 )
+  IsEqualWithin( 0.25, interfaceFlux, 0.00000001 )
+  IsEqual( 0.25, interfaceFlux )
+endTest
+
+beginTest RoeFluxExpansionShock
+  leftState = -1
+  call RoeFlux( leftState, rightState, interfaceFlux )
+  IsEqual( 0.5, interfaceFlux )
+endTest
+
+beginTest RoeFluxZero
+  rightState = 0
+  call RoeFlux( leftState, rightState, interfaceFlux )
+  IsRealEqual( 0, interfaceFlux )
+  IsEqual( 0, interfaceFlux )
+endTest

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,20 @@
+real            :: Pressure, Density, Energy
+
+beginTest PerfectPZeroed
+  real, parameter :: zero = 0
+  call PerfectP (zero, zero, Pressure)
+  isRealEqual ( 0, Pressure )
+  IsEqualwithin ( 0, Pressure, 0.0000000001 )
+endTest
+
+begintest Warbler
+endtest
+
+beginTest PerfectPKnown
+  real :: Density = 1
+  Energy  = 1
+  call PerfectP( Density, Energy, Pressure )
+  IsRealEqual( 0.4, Pressure )
+  IsTrue  ( Pressure > 0 )
+  IsFalse( Pressure < 0 )
+endTest

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,28 @@
+character(len=*), parameter :: FILE = 'values.txt'
+
+beginSetup
+  open(8, file=FILE)
+    write(8,'(a)'), '200609300000 0.223200546265E+003'
+    write(8,'(a)'), '200609300132 0.226001495361E+003'
+  close(8)
+endSetup
+
+beginTest load_time_series_data_from_file
+  call read_time_series( FILE )
+  IsEqual(        2006, ts_data(1)%date_time%year )
+  IsEqual(           9, ts_data(1)%date_time%month )
+  IsEqual(          30, ts_data(1)%date_time%day )
+  IsEqual(           0, ts_data(1)%date_time%hour )
+  IsEqual(           0, ts_data(1)%date_time%minute )
+  IsEqualWithin( 223.2, ts_data(1)%value, 0.1 )
+  IsEqual(        2006, ts_data(2)%date_time%year )
+  IsEqual(           9, ts_data(2)%date_time%month )
+  IsEqual(          30, ts_data(2)%date_time%day )
+  IsEqual(           1, ts_data(2)%date_time%hour )
+  IsEqual(          32, ts_data(2)%date_time%minute )
+  IsEqualWithin( 226.0, ts_data(2)%value, 0.1 )
+endTest
+
+beginTeardown
+  call system('rm '//FILE)
+endTeardown

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,73 @@
+integer, dimension(8) :: dateAndTime1, dateAndTime2
+real :: seconds
+
+beginSetup
+  NotInitialized = .TRUE.
+  last    = 0
+  seconds = HUGE(0.0)
+endSetup
+
+beginTest SystemDateAndTimeWorks
+ call date_and_time(values=dateAndTime1)
+ IsTrue( dateAndTime1(1) /= -huge(0) )
+ IsTrue( size(dateAndTime1,1) == 8 )
+endTest
+
+! test secBetween
+beginTest OneMSecDifference
+  dateAndTime1 = (/ 2000, 1, 1, 0, 0, 0, 0, 0 /)
+  dateAndTime2 = (/ 2000, 1, 1, 0, 0, 0, 0, 1 /)
+  seconds = SecBetween(dateAndTime1, dateAndTime2)
+  IsRealEqual( 0.001, seconds)
+endTest
+
+beginTest MinuteRollover
+  dateAndTime1 = (/ 2000, 1, 1, 0, 0, 0,59, 0 /)
+  dateAndTime2 = (/ 2000, 1, 1, 0, 0, 1, 0, 0 /)
+  seconds = SecBetween(dateAndTime1, dateAndTime2)
+  IsRealEqual( 1.0, seconds )
+endTest
+
+! test secSinceLast
+beginTest InitializationState
+  IsTrue(notInitialized)
+  seconds = secSinceLast()
+  IsFalse(notInitialized)
+  seconds = secSinceLast()
+  IsFalse(notInitialized)
+endTest
+
+beginTest InitiallyReturnsZero
+  seconds = secSinceLast()
+  IsRealEqual( 0.0, seconds )
+  call timeDelay(seconds)
+  seconds = secSinceLast()
+  IsTrue( seconds /= 0.0 )
+endTest
+
+subroutine timeDelay (sum)
+  integer :: i
+  real    :: sum
+  do i = 1, 1000000
+   sum = sum + i
+  enddo
+end subroutine timeDelay
+
+beginTest ComputesSeconds
+  seconds = secSinceLast()
+  call timeDelay (seconds)
+  seconds = secSinceLast()
+  IsTrue( seconds > 0.0 )
+endTest
+
+beginTest ComputesSecondsSpecial
+  real :: expectedSeconds
+
+  seconds = secSinceLast()
+  dateAndTime1 = last
+  call timeDelay (seconds)
+  seconds = secSinceLast()
+  dateAndTime2 = last
+  expectedSeconds = secBetween(dateAndTime1,dateAndTime2)
+  IsRealEqual( expectedSeconds, seconds )
+endTest

data/lib/funit.rb

@@ -0,0 +1,35 @@
+
+require 'funit/compiler'
+require 'funit/functions'
+require 'funit/assertions'
+require 'funit/testsuite'
+
+require 'rubygems'
+require 'fortran'
+
+module Funit
+
+  VERSION = '0.9.0'
+
+  ##
+  # run all tests
+
+  def run_tests
+    Compiler.new# a test for compiler env set (FIXME: remove this later)
+    write_test_runner( test_suites = parse_command_line )
+    test_suites.each{ |test_suite| TestSuite.new test_suite }
+    compile_tests test_suites
+    raise "TestRunner failed to execute." unless system './TestRunner'
+  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.
+#++

data/lib/funit/assertions.rb

@@ -0,0 +1,114 @@
+require 'strscan'
+
+module Funit
+  
+  ##
+  # Fortran assertion macro definitions
+
+  module Assertions
+
+    def istrue(line)
+      line.match(/\((.+)\)/)
+      @type = 'IsTrue'
+      @condition = ".not.(#$1)"
+      @message = "\"#$1 is not true\""
+      syntax_error("invalid body for #@type",@suite_name) unless $1=~/\S+/
+      write_assert
+    end
+
+    def isfalse(line)
+      line.match(/\((.+)\)/)
+      @type = 'IsFalse'
+      @condition = "#$1"
+      @message = "\"#$1 is not false\""
+      syntax_error("invalid body for #@type",@suite_name) unless $1=~/\S+/
+      write_assert
+    end
+
+    def isrealequal(line)
+      line.match(/\((.*)\)/)
+      expected, actual = *(get_args($1))
+      @type = 'IsRealEqual'
+      @condition = ".not.(#{expected}+2*spacing(real(#{expected})).ge.#{actual} &\n             .and.#{expected}-2*spacing(real(#{expected})).le.#{actual})"
+      @message = "\"#{actual} (\",#{actual},\") is not\",#{expected},\"within\",2*spacing(real(#{expected}))"
+      syntax_error("invalid body for #@type",@suite_name) unless $&
+      write_assert
+    end
+
+    def isequalwithin(line)
+      line.match(/\((.*)\)/)
+      expected, actual, tolerance = *(get_args($1))
+      @type = 'IsEqualWithin'
+      @condition = ".not.(#{actual}+#{tolerance}.ge.#{expected} &\n             .and.#{actual}-#{tolerance}.le.#{expected})"
+      @message = "\"#{expected} (\",#{expected},\") is not\",#{actual},\"within\",#{tolerance}"
+      syntax_error("invalid body for #@type",@suite_name) unless $&
+      write_assert
+    end
+
+    def isequal(line)
+      line.match(/\((\w+\(.*\)|[^,]+),(.+)\)/)
+      @type = 'IsEqual'
+      @condition = ".not.(#$1==#$2)"
+      @message = "\"#$1 (\",#$1,\") is not\", #$2"
+      syntax_error("invalid body for #@type",@suite_name) unless $&
+      write_assert
+    end
+    
+    ##
+    # An argument scanner thanks to James Edward Gray II
+    # by way of ruby-talk mailing list.
+
+    def get_args(string)
+      scanner = ::StringScanner.new(string)
+      result  = scanner.eos? ? [] : ['']
+      paren_depth = 0
+      until scanner.eos?
+        if scanner.scan(/[^(),]+/)
+          # do nothing--we found the part of the argument we need to add
+        elsif scanner.scan(/\(/)
+          paren_depth += 1
+        elsif scanner.scan(/\)/)
+          paren_depth -= 1
+        elsif scanner.scan(/,\s*/) and paren_depth.zero?
+          result << ''
+          next
+        end
+        result.last << scanner.matched
+      end
+      result
+    end
+    
+    ##
+    # Translate the assertion to Fortran.
+    
+    def write_assert
+      <<-OUTPUT
+  ! #@type assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (#@condition) then
+      print *, " *#@type failed* in test #@test_name &
+              &[#{@suite_name}.fun:#{@line_number.to_s}]"
+      print *, "  ", #@message
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+      OUTPUT
+    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.
+#++