twosum.f

Go to the documentation of this file.

!                   ******************
                    SUBROUTINE twosum
!                   ******************
!
     &(a,b,x,y)
!
!***********************************************************************
! BIEF   V7P3                                  24/02/2016
!***********************************************************************
!
!brief    SCALAR SUM AND CALCULATE ROUNDING ERROR.
!
!history  R.NHEILI (Univerte de Perpignan, DALI)
!+        24/02/2016
!+        V7P3
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| A,B          |-->| DOUBLE PRECISION SCALAR
!| X            |<--| DOUBLE PRECISION SUM RESULT
!| Y            |<--| DOUBLE PRECISION ROUNDING ERROR
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION, INTENT(IN)  :: A,B
      DOUBLE PRECISION, INTENT(OUT) :: X,Y
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION Z, XOUT
!
!-----------------------------------------------------------------------
!
        xout=a+b
        z=xout-a
        y=(a-(xout-z))+(b-z)
        x = xout
!
!
!-----------------------------------------------------------------------
!
      END