ksampAF.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	ksampaf0 (x, y, n, coef, nc, yi, total, coef2)
subroutine	ksampaf0_b (x, y, n, coef, nc)
subroutine	ksampaf (x, yi, n, coef2, nc, xs)

Function/Subroutine Documentation

ksampaf()

subroutine ksampaf	(	real*8, dimension(n)	x,
		real*8, dimension(n)	yi,
		integer	n,
		real*8, dimension(nc, 3)	coef2,
		integer	nc,
		real*8	xs
	)

Definition at line 120 of file ksampAF.f.

References kcsplintp(), and rndc().

Referenced by csampaf().

      implicit none
      integer n
      real*8 x(n), yi(n)
      integer nc
      real*8 coef2(nc, 3)
      real*8 xs  !output.  sampled x
      real*8 u

      call rndc(u)
      call kcsplintp(yi, x, n, coef2, nc, u,  xs)

Here is the call graph for this function:

Here is the caller graph for this function:

ksampaf0()

subroutine ksampaf0	(	real*8, dimension(n)	x,
		real*8, dimension(n)	y,
		integer	n,
		real*8, dimension(nc, 3)	coef,
		integer	nc,
		real*8, dimension(n)	yi,
		real*8	total,
		real*8, dimension(nc, 3)	coef2
	)

Definition at line 55 of file ksampAF.f.

References kcsplcoef(), and kcsplinteg().

Referenced by csampaf0().

      implicit none
      integer n  !  input number of data points
      real*8 x(n), y(n) ! input data points (x,y)
!                          zero is assumed at x<x(1) or x>x(n)
      integer nc  ! input. must be >=n-1
      real*8 coef(nc, 3) ! output. to keep the spline coefficents
      real*8 yi(n) ! output. yi(i)= Normalized integral from x=x(1) to x(i)
!                   (i = 1, n)      (yia(1)=0. yi(n)= 1.0) 
      real*8  total  ! integral valeu of the function
!
      real*8 coef2(nc, 3)  !ouptut. coef2 to be used for (yi, x)
!                   spline interploation      
!



      real*8 dx, a, b, sum
      integer i
      real*8 x1, x2
 
      if( nc .lt. n-1 ) then
         write(0,*) ' nc< n-1 in ksamplAF0'
         stop 1234
      endif

      call kcsplcoef(x, y, n, coef, nc)
      a = x(1)
      b = x(n)
      call kcsplinteg(x, y, n, coef, nc, a, b, total)

      yi(1) = 0.
      do i = 2, n
         b = x(i)
         call kcsplinteg(x, y, n, coef, nc, a, b, sum)
         yi(i) = sum/total
      enddo
!          to assure yi(n) = 1.
      yi(n) = 1.
      call kcsplcoef(yi, x, n, coef2, nc)

Here is the call graph for this function:

Here is the caller graph for this function:

ksampaf0_b()

subroutine ksampaf0_b	(	real*8, dimension(n)	x,
		real*8, dimension(n)	y,
		integer	n,
		real*8, dimension(nc, 3)	coef,
		integer	nc
	)

Definition at line 98 of file ksampAF.f.

References kcsplcoef().

      implicit none
!              this is only for getting interplation
      integer n  !  input number of data points
      real*8 x(n), y(n) ! input data points (x,y)
!                          zero is assumed at x<x(1) or x>x(n)
      integer nc  ! input. must be >=n-1
      real*8 coef(nc, 3) ! output. to keep the spline coefficents


      real*8 dx, a, b
      integer i
      real*8 x1, x2
 
      if( nc .lt. n-1 ) then
         write(0,*) ' nc< n-1 in ksamplAF0'
         stop 1234
      endif

      call kcsplcoef(x, y, n, coef, nc)

Here is the call graph for this function: