special_functions_gbl::cfp_asyik Interface Reference

Public Member Functions
subroutine	wp_asyik (X, FNU, KODE, FLGIK, RA, ARG, IN, Y)
subroutine	ep_asyik (X, FNU, KODE, FLGIK, RA, ARG, IN, Y)
	Quad prec version of wp_asyik tranlated to F90 standard. More...

Detailed Description

Warning

Quad precision version not implemented yet.

Member Function/Subroutine Documentation

ep_asyik()

subroutine special_functions_gbl::cfp_asyik::ep_asyik	(	real(kind=ep1), intent(in)	X,
		real(kind=ep1), intent(in)	FNU,
		integer, intent(in)	KODE,
		real(kind=ep1), intent(in)	FLGIK,
		real(kind=ep1), intent(out)	RA,
		real(kind=ep1), intent(out)	ARG,
		integer, intent(in)	IN,
		real(kind=ep1), dimension(*), intent(out)	Y
	)

Quad prec version of wp_asyik tranlated to F90 standard.

Here is the call graph for this function:

wp_asyik()

subroutine special_functions_gbl::cfp_asyik::wp_asyik	(	real(kind=wp)	X,
		real(kind=wp)	FNU,
		integer	KODE,
		real(kind=wp)	FLGIK,
		real(kind=wp)	RA,
		real(kind=wp)	ARG,
		integer	IN,
		real(kind=wp), dimension(*)	Y
	)

***SUBSIDIARY
***PURPOSE  Subsidiary to cfp_besi and DBESK
***LIBRARY   SLATEC
***TYPE      REAL(kind=wp) (ASYIK-S, wp_asyik-D)
***AUTHOR  Amos, D. E., (SNLA)
***DESCRIPTION

                    wp_asyik computes Bessel functions I and K
                  for arguments X.GT.0.0 and orders FNU.GE.35
                  on FLGIK = 1 and FLGIK = -1 respectively.

                                    INPUT

      X    - Argument, X.GT.0.00_wp
      FNU  - Order of first Bessel function
      KODE - A parameter to indicate the scaling option
             KODE=1 returns Y(I)=        I/SUB(FNU+I-1)/(X), I=1,IN
                    or      Y(I)=        K/SUB(FNU+I-1)/(X), I=1,IN
                    on FLGIK = 1.00_wp or FLGIK = -1.00_wp
             KODE=2 returns Y(I)=EXP(-X)*I/SUB(FNU+I-1)/(X), I=1,IN
                    or      Y(I)=EXP( X)*K/SUB(FNU+I-1)/(X), I=1,IN
                    on FLGIK = 1.00_wp or FLGIK = -1.00_wp
     FLGIK - Selection parameter for I or K FUNCTION
             FLGIK =  1.00_wp gives the I function
             FLGIK = -1.00_wp gives the K function
        RA - SQRT(1.+Z*Z), Z=X/FNU
       ARG - Argument of the leading exponential
        IN - Number of functions desired, IN=1 or 2

                                    OUTPUT

         Y - A vector whose first IN components contain the sequence

     Abstract  **** A REAL(kind=wp) routine ****
         wp_asyik implements the uniform asymptotic expansion of
         the I and K Bessel functions for FNU.GE.35 and real
         X.GT.0.00_wp. The forms are identical except for a change
         in sign of some of the terms. This change in sign is
         accomplished by means of the FLAG FLGIK = 1 or -1.

***SEE ALSO  cfp_besi, DBESK
***ROUTINES CALLED  F1MACH
***REVISION HISTORY  (YYMMDD)
   750101  DATE WRITTEN
   890531  Changed all specific intrinsics to generic.  (WRB)
   890911  Removed unnecessary intrinsics.  (WRB)
   891214  Prologue converted to Version 4.0 format.  (BAB)
   900328  Added TYPE section.  (WRB)
   910408  Updated the AUTHOR section.  (WRB)

Here is the call graph for this function:

---

The documentation for this interface was generated from the following file:

• /home/jacob/Work/codes/UKRmol/gitlab/ukrmol-in/source/gbtolib/source/special_functions.f90