Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(pdf_kernel) | :: | self | ||||
real(kind=8) | :: | sdev | ||||
integer | :: | N | ||||
real(kind=8) | :: | h |
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