fUnit 0.0.1 → 0.0.2

data/examples/CFD/FluxFunctions.f90

@@ -0,0 +1,35 @@
+! $Id: FluxFunctions.f90 91 2001-11-13 20:46:28Z cvs4bk $
+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/FluxFunctionsMT.f90

@@ -0,0 +1,336 @@
+! FluxFunctionsMT.f90 - a Fortran mobility test suite for FluxFunctions.f90
+!
+! [dynamically generated from FluxFunctionsMT.ftk
+!  by funit Ruby script Sun Apr 09 01:08:26 EDT 2006]
+
+module FluxFunctionsMT
+
+ use FluxFunctions
+
+ implicit none
+
+ private
+
+ public :: MTFluxFunctions
+
+ logical :: noAssertFailed
+
+ integer :: numTests          = 0
+ integer :: numAsserts        = 0
+ integer :: numAssertsTested  = 0
+ integer :: numFailures       = 0
+
+! $Id: FluxFunctionsTS.ftk 91 2001-11-13 20:46:28Z cvs4bk $
+
+ real :: leftState, rightState, interfaceFlux
+
+ contains
+
+
+ subroutine TestFluxZero
+
+  real :: state
+  state = 0
+
+  ! IsEqualWithin assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( 0+ 0.00001 .ge. Flux(state) &
+             .and. 0- 0.00001 .le. Flux(state))) then
+      print *, " *IsEqualWithin failed* in test FluxZero &
+                         &[FluxFunctionsMT.ftk:13]"
+      print *, "  ", " Flux(state) (", Flux(state),") is not", 0,"within", 0.00001 
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  numTests = numTests + 1
+
+ end subroutine TestFluxZero
+
+ 
+ subroutine TestFluxOne
+
+  real :: state = 1
+
+  ! IsEqualWithin assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( 0.5+ 0.00001 .ge. Flux(state) &
+             .and. 0.5- 0.00001 .le. Flux(state))) then
+      print *, " *IsEqualWithin failed* in test FluxOne &
+                         &[FluxFunctionsMT.ftk:18]"
+      print *, "  ", " Flux(state) (", Flux(state),") is not", 0.5,"within", 0.00001 
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  numTests = numTests + 1
+
+ end subroutine TestFluxOne
+
+
+ subroutine TestRoeAvgZero
+
+
+  ! IsRealEqual assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( RoeAvg(0.0+2*spacing(real( RoeAvg(0.0)).ge.0.0), 0  &
+             .and. RoeAvg(0.0-2*spacing(real( RoeAvg(0.0)).le.0.0), 0 )) then
+      print *, " *IsRealEqual failed* in test RoeAvgZero &
+                         &[FluxFunctionsMT.ftk:22]"
+      print *, "  ", "0.0), 0  (",0.0), 0 ,") is not", RoeAvg(0.0,"within",2*spacing(real( RoeAvg(0.0))
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  ! IsFalse assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if ( RoeAvg(0.0,0.0).ne.1 ) then
+      print *, " *IsFalse failed* in test RoeAvgZero &
+                         &[FluxFunctionsMT.ftk:23]"
+      print *, "  ", " RoeAvg(0.0,0.0).ne.1  is not false"
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  numTests = numTests + 1
+
+ end subroutine TestRoeAvgZero
+
+
+ subroutine TestRoeAvgKnown
+
+
+  ! IsRealEqual assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( RoeAvg(leftState+2*spacing(real( RoeAvg(leftState)).ge.rightState), 0.5  &
+             .and. RoeAvg(leftState-2*spacing(real( RoeAvg(leftState)).le.rightState), 0.5 )) then
+      print *, " *IsRealEqual failed* in test RoeAvgKnown &
+                         &[FluxFunctionsMT.ftk:27]"
+      print *, "  ", "rightState), 0.5  (",rightState), 0.5 ,") is not", RoeAvg(leftState,"within",2*spacing(real( RoeAvg(leftState))
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  ! IsTrue assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( RoeAvg(leftState,rightState).lt.0 )) then
+      print *, " *IsTrue failed* in test RoeAvgKnown &
+                         &[FluxFunctionsMT.ftk:28]"
+      print *, "  ", " RoeAvg(leftState,rightState).lt.0  is not true"
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  numTests = numTests + 1
+
+ end subroutine TestRoeAvgKnown
+
+
+ subroutine TestCentralFluxKnown
+
+  call CentralFlux( leftState, rightState, interfaceFlux )
+
+  ! IsEqualWithin assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( 0.5+ 0.001 .ge. interfaceFlux &
+             .and. 0.5- 0.001 .le. interfaceFlux)) then
+      print *, " *IsEqualWithin failed* in test CentralFluxKnown &
+                         &[FluxFunctionsMT.ftk:33]"
+      print *, "  ", " interfaceFlux (", interfaceFlux,") is not", 0.5,"within", 0.001 
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  ! IsEqualWithin assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( 0.5+ 0.00000001 .ge. interfaceFlux &
+             .and. 0.5- 0.00000001 .le. interfaceFlux)) then
+      print *, " *IsEqualWithin failed* in test CentralFluxKnown &
+                         &[FluxFunctionsMT.ftk:34]"
+      print *, "  ", " interfaceFlux (", interfaceFlux,") is not", 0.5,"within", 0.00000001 
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  ! IsEqual assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( interfaceFlux== 0.5 )) then
+      print *, " *IsEqual failed* in test CentralFluxKnown &
+                         &[FluxFunctionsMT.ftk:35]"
+      print *, "  ", " interfaceFlux (", interfaceFlux,") is not",  0.5 
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  numTests = numTests + 1
+
+ end subroutine TestCentralFluxKnown
+
+
+ subroutine TestRoeFluxExpansionShock
+
+  leftState = -1
+  call RoeFlux( leftState, rightState, interfaceFlux )
+
+  ! IsEqual assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( interfaceFlux== 0 )) then
+      print *, " *IsEqual failed* in test RoeFluxExpansionShock &
+                         &[FluxFunctionsMT.ftk:41]"
+      print *, "  ", " interfaceFlux (", interfaceFlux,") is not",  0 
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  numTests = numTests + 1
+
+ end subroutine TestRoeFluxExpansionShock
+
+
+ subroutine TestRoeFluxZero
+
+  rightState = 0
+  call RoeFlux( leftState, rightState, interfaceFlux )
+
+  ! IsRealEqual assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( interfaceFlux+2*spacing(real( interfaceFlux)).ge. 0  &
+             .and. interfaceFlux-2*spacing(real( interfaceFlux)).le. 0 )) then
+      print *, " *IsRealEqual failed* in test RoeFluxZero &
+                         &[FluxFunctionsMT.ftk:47]"
+      print *, "  ", " 0  (", 0 ,") is not", interfaceFlux,"within",2*spacing(real( interfaceFlux))
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  ! IsEqual assertion
+  numAsserts = numAsserts + 1
+  if (noAssertFailed) then
+    if (.not.( interfaceFlux== 0 )) then
+      print *, " *IsEqual failed* in test RoeFluxZero &
+                         &[FluxFunctionsMT.ftk:48]"
+      print *, "  ", " interfaceFlux (", interfaceFlux,") is not",  0 
+      print *, ""
+      noAssertFailed = .false.
+      numFailures    = numFailures + 1
+    else
+      numAssertsTested = numAssertsTested + 1
+    endif
+  endif
+
+  numTests = numTests + 1
+
+ end subroutine TestRoeFluxZero
+
+
+ subroutine Setup
+  leftState  = 0
+  rightState = 1
+  noAssertFailed = .true.
+ end subroutine Setup
+
+
+ subroutine Teardown
+ end subroutine Teardown
+
+
+ subroutine MTFluxFunctions( nTests, nAsserts, nAssertsTested, nFailures )
+
+  integer :: nTests
+  integer :: nAsserts
+  integer :: nAssertsTested
+  integer :: nFailures
+
+  continue
+
+  call Setup
+  call TestFluxZero
+  call Teardown
+
+  call Setup
+  call TestFluxOne
+  call Teardown
+
+  call Setup
+  call TestRoeAvgZero
+  call Teardown
+
+  call Setup
+  call TestRoeAvgKnown
+  call Teardown
+
+  call Setup
+  call TestCentralFluxKnown
+  call Teardown
+
+  call Setup
+  call TestRoeFluxExpansionShock
+  call Teardown
+
+  call Setup
+  call TestRoeFluxZero
+  call Teardown
+
+  nTests          = numTests
+  nAsserts        = numAsserts
+  nAssertsTested  = numAssertsTested
+  nFailures       = numFailures
+
+ end subroutine MTFluxFunctions
+
+end module FluxFunctionsMT

data/examples/CFD/FluxFunctionsMT.ftk

@@ -0,0 +1,49 @@
+! $Id: FluxFunctionsTS.ftk 91 2001-11-13 20:46:28Z cvs4bk $
+
+ real :: leftState, rightState, interfaceFlux
+
+beginsetup
+  leftState  = 0
+  rightState = 1
+endsetup
+
+beginTest FluxZero
+  real :: state
+  state = 0
+  IsEqualWithin( Flux(state), 0, 0.00001 )
+endTest
+ 
+beginTest FluxOne
+  real :: state = 1
+  IsEqualWithin( Flux(state), 0.5, 0.00001 )
+endTest
+
+beginTest RoeAvgZero
+  IsRealEqual( RoeAvg(0.0,0.0), 0 )
+  IsFalse( RoeAvg(0.0,0.0).ne.1 )
+endTest
+
+beginTest RoeAvgKnown
+  IsRealEqual( RoeAvg(leftState,rightState), 0.5 )
+  IsTrue( RoeAvg(leftState,rightState).lt.0 )
+endTest
+
+beginTest CentralFluxKnown
+  call CentralFlux( leftState, rightState, interfaceFlux )
+  IsEqualWithin( interfaceFlux, 0.5, 0.001 )
+  IsEqualWithin( interfaceFlux, 0.5, 0.00000001 )
+  IsEqual( interfaceFlux, 0.5 )
+endTest
+
+beginTest RoeFluxExpansionShock
+  leftState = -1
+  call RoeFlux( leftState, rightState, interfaceFlux )
+  IsEqual( interfaceFlux, 0 )
+endTest
+
+beginTest RoeFluxZero
+  rightState = 0
+  call RoeFlux( leftState, rightState, interfaceFlux )
+  IsRealEqual( interfaceFlux, 0 )
+  IsEqual( interfaceFlux, 0 )
+endTest

data/examples/CFD/Gammas.f90

@@ -0,0 +1,8 @@
+! $Id: Gammas.f90 91 2001-11-13 20:46:28Z cvs4bk $
+module Gammas
+
+  implicit none
+  
+  real, parameter :: Gamma = 1.4
+
+end module Gammas

data/examples/CFD/GasModel.f90

@@ -0,0 +1,17 @@
+! $Id: GasModel.f90 91 2001-11-13 20:46:28Z cvs4bk $
+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