lissage.f

Go to the documentation of this file.

!                       ******************
                        SUBROUTINE lissage
!                       ******************
!
     &(dim,entree,sortie)
!
!***********************************************************************
! ARTEMIS   V7P0                                   July 2014
!***********************************************************************
!
!brief    SMOOTHES A FUNCTION USING LEAST SQUARE METHOD
!         BASED ON A PARABOLA FIT TO 5 EQUALLY SPACED POINTS.
!
!history  S.BOURBAN (HRW)
!+        1997
!+
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| DIM            |   | DIMENSION OF THE FUNCTION
!| ENTREE         |   | INPUT VALUES
!| SORTIE         |   | OUTPUT VALUES
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER DIM,I
      DOUBLE PRECISION ENTREE(dim),SORTIE(dim)
!
      INTRINSIC dabs
!
!-----------------------------------------------------------------------
!     INITIALISES ENTREE
      DO i = 1, dim
        entree(i) = dabs(entree(i))
      ENDDO
      WRITE(lu,*)
!
!-----------------------------------------------------------------------
!
      DO i = 1, dim
!
!-----------------------------------------------------------------------
        IF(i.EQ.1) THEN
          sortie(i)=31.d0*entree(i)+9.d0*entree(i+1)
     &-3.d0*entree(i+2)-5.d0*entree(i+3)+3.d0*entree(i+4)
          sortie(i)=sortie(i)/35.d0
!
!-----------------------------------------------------------------------
        ELSEIF(i.EQ.2) THEN
          sortie(i)=9.d0*entree(i)+13.d0*entree(i+1)
     &+12.d0*entree(i+2)+6.d0*entree(i+3)-5.d0*entree(i+4)
          sortie(i)=sortie(i)/35.d0
!
!-----------------------------------------------------------------------
        ELSEIF(i.EQ.dim-1) THEN
          sortie(i)=-5.d0*entree(dim-4)+6.d0*entree(dim-3)
     &+12.d0*entree(dim-2)+13.d0*entree(dim-1)+9.d0*entree(dim)
          sortie(i)=sortie(i)/35.d0
!
!-----------------------------------------------------------------------
        ELSEIF(i.EQ.dim) THEN
          sortie(i)=3.d0*entree(dim-4)-5.d0*entree(dim-3)
     &-3.d0*entree(dim-2)+9.d0*entree(dim-1)+31.d0*entree(dim)
          sortie(i)=sortie(i)/35.d0
!
!-----------------------------------------------------------------------
        ELSE
          sortie(i)=-3.d0*entree(i-2)+12.d0*entree(i-1)
     &+17.d0*entree(i)+12.d0*entree(i+1)-3.d0*entree(i+2)
          sortie(i)=sortie(i)/35.d0
        ENDIF
      ENDDO
!
!-----------------------------------------------------------------------
      RETURN
      END