pdf_sigma_sdev_1d Subroutine

subroutine pdf_sigma_sdev_1d(self, sdev, N, h)

Arguments

Type IntentOptional Attributes Name
type(pdf_kernel) :: self
real(kind=8) :: sdev
integer :: N
real(kind=8) :: h

Source Code

subroutine pdf_sigma_sdev_1d(self,sdev,N,h)
  type(pdf_kernel) :: self    
  real(8)          :: sdev
  integer          :: N
  real(8)          :: h
  if(.not.self%status)stop "PDF_SET_SIGMA: PDF not allocated"
  h = (4d0/3d0/N)**(1d0/5d0)*sdev !Silverman's rule of thumb.
end subroutine pdf_sigma_sdev_1d