g1_triweight Function

public pure function g1_triweight(r) result(val)

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: r

Distance from filter center

Return Value real(kind=wp)

Result


Source Code

    pure function g1_triweight(r) result(val)
      ! Triweight filter.
      implicit none
      real(wp), intent(in) :: r                                                !! Distance from filter center
      real(wp)             :: val                                              !! Result

      if ( (abs(r).le. 1.0_wp)) then
        val=(35.0_wp/32.0_wp)*((1-r**2)**3)
      else
        val=0.0_WP
      end if

      return
    end function g1_triweight