esputr_test.f90

Go to the documentation of this file.

!    Copyright (c) 2016  Forschungszentrum Juelich GmbH
!                        Markus Brenneis, Vladislav Kotov 
!
!    This file is part of ESPUTR.
!
!    ESPUTR is free software: you can redistribute it and/or modify
!    it under the terms of the GNU General Public License as published by
!    the Free Software Foundation, either version 3 of the License, or
!    (at your option) any later version.
!
!    ESPUTR is distributed in the hope that it will be useful,
!    but WITHOUT ANY WARRANTY; without even the implied warranty of
!    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!    GNU General Public License for more details.
!
!    You should have received a copy of the GNU General Public License
!    along with ESPUTR.  If not, see <http://www.gnu.org/licenses/>.
!

module esputr_test
    use esputr
    implicit none
    public ispassed, printfailed
    
    real(kind=ESPUTR_DP), parameter, public :: esputr_tol = 2e-7
    
    contains
    
    logical elemental function ispassed(expected, calculated)
        real(kind=ESPUTR_DP), intent(in) :: expected, calculated
        intrinsic tiny
        real(kind=ESPUTR_DP) :: eps
        eps=10.*tiny(expected)
        ispassed = abs(calculated - expected) < esputr_tol*expected + eps
    end function ispassed
    
    subroutine printfailed(expected, calculated, passed)
        real(kind=ESPUTR_DP), dimension(:), intent(in) :: expected, calculated
        logical, dimension(size(expected)), intent(in) :: passed
        integer :: i
        do i = 1, size(expected)
            if(.not. passed(i)) then
                print *, "  Expected", expected(i), ", got", calculated(i)
                print *, "  Relative Error", abs(calculated(i) - expected(i))*100/expected(i), "%"
            end if
        end do
    end subroutine printfailed
    
end module esputr_test