r8psp_tolsum.f

C  utility routines for (r8)pspltsub
 
      subroutine r8psp_tolsum(refval,tolval,sum)
 
C  for R4 sum.eq.refval expected, for R8 difference will be detectable
C  refval =~ 1.0; tolval =~1.0e-6
 
      IMPLICIT NONE
      INTEGER, PARAMETER :: R8=selected_real_kind(12,100)
      real*8 refval, tolval  ! input
      real*8 sum             ! output
 
      sum = refval + tolval*tolval
 
      return
      end