pdf_sigma_data_2d Subroutine

subroutine pdf_sigma_data_2d(self, data, h)

Arguments

Type IntentOptional Attributes Name
type(pdf_kernel_2d) :: self
real(kind=8), dimension(:,:) :: data
real(kind=8) :: h(2,2)

Calls

proc~~pdf_sigma_data_2d~~CallsGraph proc~pdf_sigma_data_2d pdf_sigma_data_2d get_var get_var proc~pdf_sigma_data_2d->get_var

Source Code

subroutine pdf_sigma_data_2d(self,data,h)
  type(pdf_kernel_2d)    :: self
  real(8),dimension(:,:) :: data
  real(8)                :: h(2,2)
  integer                :: L
  if(.not.self%status)stop "PDF_=SIGMA: PDF not allocated"
  L = size(data,2)
  if(any(shape(data)/=[2,L]))stop "PDF_SET_SIGMA: DATA wrong dimensions"
  h = 0d0
  h(1,1) = get_var(data(1,:))
  h(2,2) = get_var(data(2,:))
  h(1,1) = (1d0/size(data,2))**(2/6d0)*h(1,1) !Silverman's rule of thumb.
  h(2,2) = (1d0/size(data,2))**(2/6d0)*h(2,2) !Silverman's rule of thumb.
end subroutine pdf_sigma_data_2d