ksampRSA.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	ksamprsa (costheta)

Function/Subroutine Documentation

ksamprsa()

subroutine ksamprsa

(

real*8

costheta

)

Definition at line 8 of file ksampRSA.f.

References d0, and rndc().

Referenced by ccohs().

      implicit none
!      sample cos from  (1+cos^2)dcos
      real*8  costheta  ! output  cos ofsampled angle  
      real*8  u1, u2, u3

      call rndc(u1)
      if(u1 .lt. 0.75d0) then
         call rndc(u2)
         costheta = 2.0*u2 -1.
      else
!         second term. take max of 3 u's.
         call rndc(u1)
         call rndc(u2)
         call rndc(u3)
         costheta = max(u1, u2, u3)
         call rndc(u1)
         if(u1 .lt. 0.5) then
            costheta = -costheta
         endif
      endif

Here is the call graph for this function:

Here is the caller graph for this function: