unit_test_r2 Subroutine

private subroutine unit_test_r2(this, res, expected, tol, msg)

Type Bound

unit_test

Arguments

Type IntentOptional Attributes Name
class(unit_test), intent(inout) :: this
real(kind=rk), intent(in), dimension(:,:) :: res
real(kind=rk), intent(in), dimension(:,:) :: expected
real(kind=rk), intent(in), optional :: tol
character(len=*), intent(in), optional :: msg

Calls

proc~~unit_test_r2~~CallsGraph proc~unit_test_r2 unit_test%unit_test_r2 proc~print_msg unit_test%print_msg proc~unit_test_r2->proc~print_msg colorize colorize proc~print_msg->colorize

Called by

proc~~unit_test_r2~~CalledByGraph proc~unit_test_r2 unit_test%unit_test_r2 none~check unit_test%check none~check->proc~unit_test_r2 program~demo demo program~demo->none~check program~test test program~test->none~check

Source Code

   subroutine unit_test_r2(this, res, expected, tol, msg)
      class(unit_test), intent(inout) :: this
      real(rk), intent(in), dimension(:,:) :: res, expected
      real(rk), intent(in), optional :: tol
      character(*), intent(in), optional :: msg
      logical :: condition
      real(rk) :: rel_err

      if (present(msg)) then
         this%msg = msg
      else
         this%msg = 'forunittest'
      end if

      if (norm2(expected)<tiny(0.0_rk)) then
         rel_err = norm2(res-expected)
      else
         rel_err = norm2(res-expected)/norm2(expected)
      end if

      if (present(tol)) then
         condition = rel_err < tol
      else
         ! condition = rel_err == 0.0_rk
         condition = abs(rel_err) < 2.0_rk*epsilon(0.0_rk)
      end if

      call this%print_msg(condition)
   end subroutine unit_test_r2