dsinint Function

private function dsinint(xvalue) result(fn_val)

Calculate sine integral of xvalue. AUTHOR: Allan MacLeod Dept. of Mathematics and Statistics University of Paisley Scotland (e-mail: macl_ms0@paisley.ac.uk)

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(in) :: xvalue

Return Value real(kind=rk)


Called by

proc~~dsinint~~CalledByGraph proc~dsinint dsinint proc~tem_eval_heaviside_gibbs tem_eval_heaviside_gibbs proc~tem_eval_heaviside_gibbs->proc~dsinint proc~tem_spatial_for_coord tem_spatial_for_coord proc~tem_spatial_for_coord->proc~tem_eval_heaviside_gibbs proc~tem_spatial_scalar_for_index tem_spatial_scalar_for_index proc~tem_spatial_scalar_for_index->proc~tem_spatial_for_coord interface~tem_spatial_for tem_spatial_for interface~tem_spatial_for->proc~tem_spatial_for_coord interface~tem_spatial_for->proc~tem_spatial_scalar_for_index proc~tem_spatial_vector_storeval tem_spatial_vector_storeVal proc~tem_spatial_vector_storeval->interface~tem_spatial_for proc~tem_spacetime_scalar_for_index tem_spacetime_scalar_for_index proc~tem_spacetime_scalar_for_index->interface~tem_spatial_for proc~tem_spacetime_vector_for_index tem_spacetime_vector_for_index proc~tem_spacetime_vector_for_index->interface~tem_spatial_for proc~tem_spacetime_for_coord tem_spacetime_for_coord proc~tem_spacetime_for_coord->interface~tem_spatial_for proc~tem_spacetime_for_treeids tem_spacetime_for_treeIDs proc~tem_spacetime_for_treeids->interface~tem_spatial_for proc~tem_spatial_scalar_storeval tem_spatial_scalar_storeVal proc~tem_spatial_scalar_storeval->interface~tem_spatial_for proc~tem_spacetime_vector_for_coord tem_spacetime_vector_for_coord proc~tem_spacetime_vector_for_coord->interface~tem_spatial_for proc~tem_spacetime_vector_for_treeids tem_spacetime_vector_for_treeIDs proc~tem_spacetime_vector_for_treeids->interface~tem_spatial_for

Contents


Variables

TypeVisibilityAttributesNameInitial
integer, private :: i
integer, private :: indsgn
real(kind=rk), private :: cx
real(kind=rk), private :: fival
real(kind=rk), private :: gival
real(kind=rk), private :: sumden
real(kind=rk), private :: sumnum
real(kind=rk), private :: sx
real(kind=rk), private :: x
real(kind=rk), private :: xhigh
real(kind=rk), private :: xsq
real(kind=rk), private, parameter:: zero =0.0_rk
real(kind=rk), private, parameter:: one =1.0_rk
real(kind=rk), private, parameter:: six =6.0_rk
real(kind=rk), private, parameter:: twelve =12.0_rk
real(kind=rk), private, parameter:: piby2 =1.5707963267948966192_rk
real(kind=rk), private, parameter:: xlow =4.47E-8_rk
real(kind=rk), private, parameter:: xhigh1 =2.32472E8_rk
real(kind=rk), private, parameter:: xhigh2 =9.0072E15_rk
real(kind=rk), private, parameter:: xhigh3 =1.4148475E16_rk
real(kind=rk), private, parameter:: asintn(0:7) =(/1.0_rk, -0.44663998931312457298E-1_rk, 0.11209146443112369449E-2_rk, -0.13276124407928422367E-4_rk, 0.85118014179823463879E-7_rk, -0.29989314303147656479E-9_rk, 0.55401971660186204711E-12_rk, -0.42406353433133212926E-15_rk/)
real(kind=rk), private, parameter:: asintd(0:7) =(/1.0_rk, 0.10891556624243098264E-1_rk, 0.59334456769186835896E-4_rk, 0.21231112954641805908E-6_rk, 0.54747121846510390750E-9_rk, 0.10378561511331814674E-11_rk, 0.13754880327250272679E-14_rk, 0.10223981202236205703E-17_rk/)
real(kind=rk), private, parameter:: afn1(0:7) =(/0.99999999962173909991_rk, 0.36451060338631902917E3_rk, 0.44218548041288440874E5_rk, 0.22467569405961151887E7_rk, 0.49315316723035561922E8_rk, 0.43186795279670283193E9_rk, 0.11847992519956804350E10_rk, 0.45573267593795103181E9_rk/)
real(kind=rk), private, parameter:: afd1(0:7) =(/1.0_rk, 0.36651060273229347594E3_rk, 0.44927569814970692777E5_rk, 0.23285354882204041700E7_rk, 0.53117852017228262911E8_rk, 0.50335310667241870372E9_rk, 0.16575285015623175410E10_rk, 0.11746532837038341076E10_rk/)
real(kind=rk), private, parameter:: agn1(0:8) =(/0.99999999920484901956_rk, 0.51385504875307321394E3_rk, 0.92293483452013810811E5_rk, 0.74071341863359841727E7_rk, 0.28142356162841356551E9_rk, 0.49280890357734623984E10_rk, 0.35524762685554302472E11_rk, 0.79194271662085049376E11_rk, 0.17942522624413898907E11_rk/)
real(kind=rk), private, parameter:: agd1(0:8) =(/1.0_rk, 0.51985504708814870209E3_rk, 0.95292615508125947321E5_rk, 0.79215459679762667578E7_rk, 0.31977567790733781460E9_rk, 0.62273134702439012114E10_rk, 0.54570971054996441467E11_rk, 0.18241750166645704670E12_rk, 0.15407148148861454434E12_rk/)
real(kind=rk), private, parameter:: afn2(0:7) =(/0.19999999999999978257E1_rk, 0.22206119380434958727E4_rk, 0.84749007623988236808E6_rk, 0.13959267954823943232E9_rk, 0.10197205463267975592E11_rk, 0.30229865264524075951E12_rk, 0.27504053804288471142E13_rk, 0.21818989704686874983E13_rk/)
real(kind=rk), private, parameter:: afd2(0:7) =(/1.0_rk, 0.11223059690217167788E4_rk, 0.43685270974851313242E6_rk, 0.74654702140658116258E8_rk, 0.58580034751805687471E10_rk, 0.20157980379272098841E12_rk, 0.26229141857684496445E13_rk, 0.87852907334918467516E13_rk/)
real(kind=rk), private, parameter:: agn2(0:8) =(/0.59999999999999993089E1_rk, 0.96527746044997139158E4_rk, 0.56077626996568834185E7_rk, 0.15022667718927317198E10_rk, 0.19644271064733088465E12_rk, 0.12191368281163225043E14_rk, 0.31924389898645609533E15_rk, 0.25876053010027485934E16_rk, 0.12754978896268878403E16_rk/)
real(kind=rk), private, parameter:: agd2(0:8) =(/1.0_rk, 0.16287957674166143196E4_rk, 0.96636303195787870963E6_rk, 0.26839734750950667021E9_rk, 0.37388510548029219241E11_rk, 0.26028585666152144496E13_rk, 0.85134283716950697226E14_rk, 0.11304079361627952930E16_rk, 0.42519841479489798424E16_rk/)