* * $Id: c300gt.F,v 1.1.1.1 1996/02/15 17:48:39 mclareni Exp $ * * $Log: c300gt.F,v $ * Revision 1.1.1.1 1996/02/15 17:48:39 mclareni * Kernlib * * #include "kernnumt/pilot.h" #if defined(CERNLIB_NEVER) SUBROUTINE C300GT(I,ISIGN,X,EREF,ECREF,F2REF) IMPLICIT DOUBLE PRECISION (A-H,O-Z) DIMENSION XTAB(100),YTAB(100) DATA NTAB/100/ C C IF ISIGN.EQ.1 SETS X=XTAB(I). C IF ISIGN.EQ.2 SETS X=-XTAB(I). C SETS EREF=ERF(X), ECREF=ERFC(X), F2REF=FREQ(SQRT(2)*X), WHERE THE C XTAB(I) AND THE FUNCTION VALUES ARE TAKEN FROM THE TABLE OF DATA C STATEMENTS (PREVIOUSLY CONSTRUCTED BY PROGRAM ECWRITE). C C TABLE. DATA XTAB( 1)/.917534D+01/, YTAB( 1)/ .1676371373263653D-37/ DATA XTAB( 2)/.443708D+00/, YTAB( 2)/ .4696670689209019D+00/ DATA XTAB( 3)/.325435D+00/, YTAB( 3)/ .3546521550410488D+00/ DATA XTAB( 4)/.953785D+00/, YTAB( 4)/ .1773833253302947D+00/ DATA XTAB( 5)/.241720D-01/, YTAB( 5)/ .2726986999298563D-01/ DATA XTAB( 6)/.170710D+00/, YTAB( 6)/ .1907706944951034D+00/ DATA XTAB( 7)/.496750D+00/, YTAB( 7)/ .5176391984092610D+00/ DATA XTAB( 8)/.113611D+00/, YTAB( 8)/ .1276468520075838D+00/ DATA XTAB( 9)/.170011D+00/, YTAB( 9)/ .1900045197721693D+00/ DATA XTAB( 10)/.118468D+01/, YTAB( 10)/ .9385763479811183D-01/ DATA XTAB( 11)/.400029D+00/, YTAB( 11)/ .4284202394208670D+00/ DATA XTAB( 12)/.463354D+00/, YTAB( 12)/ .4877134774720710D+00/ DATA XTAB( 13)/.784774D+01/, YTAB( 13)/ .1277377188008478D-27/ DATA XTAB( 14)/.409844D+00/, YTAB( 14)/ .4378202902930709D+00/ DATA XTAB( 15)/.257202D+00/, YTAB( 15)/ .2839467513613266D+00/ DATA XTAB( 16)/.584129D+01/, YTAB( 16)/ .1446611273515512D-15/ DATA XTAB( 17)/.165450D+00/, YTAB( 17)/ .1850007630301463D+00/ DATA XTAB( 18)/.405904D+00/, YTAB( 18)/ .4340558435695968D+00/ DATA XTAB( 19)/.678301D+01/, YTAB( 19)/ .8587216485742127D-21/ DATA XTAB( 20)/.131452D+00/, YTAB( 20)/ .1474777597992756D+00/ DATA XTAB( 21)/.399609D+00/, YTAB( 21)/ .4280163336112963D+00/ DATA XTAB( 22)/.719236D-01/, YTAB( 22)/ .8101736677809329D-01/ DATA XTAB( 23)/.801624D-01/, YTAB( 23)/ .9026020320336464D-01/ DATA XTAB( 24)/.316115D+00/, YTAB( 24)/ .3451640158076886D+00/ DATA XTAB( 25)/.110325D+01/, YTAB( 25)/ .1187052757341097D+00/ DATA XTAB( 26)/.183969D+00/, YTAB( 26)/ .2052684780282969D+00/ DATA XTAB( 27)/.200078D+00/, YTAB( 27)/ .2227871504045401D+00/ DATA XTAB( 28)/.597662D+01/, YTAB( 28)/ .2858225799027289D-16/ DATA XTAB( 29)/.228525D+00/, YTAB( 29)/ .2534434621564421D+00/ DATA XTAB( 30)/.231465D+00/, YTAB( 30)/ .2565899708120994D+00/ DATA XTAB( 31)/.168515D+01/, YTAB( 31)/ .1716463961002001D-01/ DATA XTAB( 32)/.256622D+00/, YTAB( 32)/ .2833340935618121D+00/ DATA XTAB( 33)/.471115D+00/, YTAB( 33)/ .4947532969229083D+00/ DATA XTAB( 34)/.164010D+01/, YTAB( 34)/ .2037055912411078D-01/ DATA XTAB( 35)/.276759D+00/, YTAB( 35)/ .3044956963291484D+00/ DATA XTAB( 36)/.233049D+00/, YTAB( 36)/ .2582834617301759D+00/ DATA XTAB( 37)/.482360D+01/, YTAB( 37)/ .9003172523514829D-11/ DATA XTAB( 38)/.298440D+00/, YTAB( 38)/ .3270172406081429D+00/ DATA XTAB( 39)/.135617D+00/, YTAB( 39)/ .1520943928223444D+00/ DATA XTAB( 40)/.133854D+01/, YTAB( 40)/ .5836034014253724D-01/ DATA XTAB( 41)/.334745D+00/, YTAB( 41)/ .3640728256113021D+00/ DATA XTAB( 42)/.438535D+00/, YTAB( 42)/ .4648621244359200D+00/ DATA XTAB( 43)/.880946D+01/, YTAB( 43)/ .1257794187408745D-34/ DATA XTAB( 44)/.189212D+00/, YTAB( 44)/ .2109821274364948D+00/ DATA XTAB( 45)/.937724D-01/, YTAB( 45)/ .1055014982462454D+00/ DATA XTAB( 46)/.123848D+00/, YTAB( 46)/ .1390362812367687D+00/ DATA XTAB( 47)/.250532D+00/, YTAB( 47)/ .2768902426361637D+00/ DATA XTAB( 48)/.435657D+00/, YTAB( 48)/ .4621794224069148D+00/ DATA XTAB( 49)/.211924D+00/, YTAB( 49)/ .2355984220670371D+00/ DATA XTAB( 50)/.257939D+00/, YTAB( 50)/ .2847249855776238D+00/ DATA XTAB( 51)/.468113D+00/, YTAB( 51)/ .4920363107479127D+00/ DATA XTAB( 52)/.684849D+01/, YTAB( 52)/ .3484321315524656D-21/ DATA XTAB( 53)/.428405D-01/, YTAB( 53)/ .4831077084644222D-01/ DATA XTAB( 54)/.349104D-02/, YTAB( 54)/ .3939200804670809D-02/ DATA XTAB( 55)/.529556D+01/, YTAB( 55)/ .6937244552912752D-13/ DATA XTAB( 56)/.242891D+00/, YTAB( 56)/ .2687774684053429D+00/ DATA XTAB( 57)/.361056D+00/, YTAB( 57)/ .3903760322920853D+00/ DATA XTAB( 58)/.341639D+01/, YTAB( 58)/ .1355047657224369D-05/ DATA XTAB( 59)/.193337D+00/, YTAB( 59)/ .2154694715212153D+00/ DATA XTAB( 60)/.490651D+00/, YTAB( 60)/ .5122458560929465D+00/ DATA XTAB( 61)/.781475D+01/, YTAB( 61)/ .2150415374192762D-27/ DATA XTAB( 62)/.418332D+00/, YTAB( 62)/ .4458887518137371D+00/ DATA XTAB( 63)/.365366D-02/, YTAB( 63)/ .4122695482703374D-02/ DATA XTAB( 64)/.816031D+00/, YTAB( 64)/ .2484829061485597D+00/ DATA XTAB( 65)/.308077D+00/, YTAB( 65)/ .3369359227204185D+00/ DATA XTAB( 66)/.295853D+00/, YTAB( 66)/ .3243448200464001D+00/ DATA XTAB( 67)/.196144D+00/, YTAB( 67)/ .2185189606135350D+00/ DATA XTAB( 68)/.214419D+00/, YTAB( 68)/ .2382886566807954D+00/ DATA XTAB( 69)/.206053D+00/, YTAB( 69)/ .2292568392047439D+00/ DATA XTAB( 70)/.487353D+01/, YTAB( 70)/ .5493133987860958D-11/ DATA XTAB( 71)/.577190D-01/, YTAB( 71)/ .6505666408085011D-01/ DATA XTAB( 72)/.124564D+00/, YTAB( 72)/ .1398318324566679D+00/ DATA XTAB( 73)/.835124D+01/, YTAB( 73)/ .3447652187669018D-31/ DATA XTAB( 74)/.270668D+00/, YTAB( 74)/ .2981188555376397D+00/ DATA XTAB( 75)/.149457D+00/, YTAB( 75)/ .1673968448379378D+00/ DATA XTAB( 76)/.465863D+01/, YTAB( 76)/ .4448978249281597D-10/ DATA XTAB( 77)/.338135D+00/, YTAB( 77)/ .3674886523173324D+00/ DATA XTAB( 78)/.352951D+00/, YTAB( 78)/ .3823249396434660D+00/ DATA XTAB( 79)/.161540D+01/, YTAB( 79)/ .2234095399787408D-01/ DATA XTAB( 80)/.202682D+00/, YTAB( 80)/ .2256086726378163D+00/ DATA XTAB( 81)/.457367D+00/, YTAB( 81)/ .4822480672698303D+00/ DATA XTAB( 82)/.101656D+01/, YTAB( 82)/ .1505382298224469D+00/ DATA XTAB( 83)/.406001D+00/, YTAB( 83)/ .4341486666930002D+00/ DATA XTAB( 84)/.862220D-01/, YTAB( 84)/ .9705055036924218D-01/ DATA XTAB( 85)/.146177D+01/, YTAB( 85)/ .3871013595458409D-01/ DATA XTAB( 86)/.354186D+00/, YTAB( 86)/ .3835547284199477D+00/ DATA XTAB( 87)/.451218D+00/, YTAB( 87)/ .4766035302472280D+00/ DATA XTAB( 88)/.344042D+00/, YTAB( 88)/ .3734219296483067D+00/ DATA XTAB( 89)/.315927D+00/, YTAB( 89)/ .3449720427714967D+00/ DATA XTAB( 90)/.100722D+00/, YTAB( 90)/ .1132694411057261D+00/ DATA XTAB( 91)/.550044D+00/, YTAB( 91)/ .4366399457502890D+00/ #if !defined(CERNLIB_NUMD38) DATA XTAB( 92)/.113232D+02/, YTAB( 92)/ .1029847777887544D-56/ DATA XTAB( 93)/.108756D+02/, YTAB( 93)/ .2215026448145411D-52/ DATA XTAB( 94)/.100892D+02/, YTAB( 94)/ .3449654485465253D-45/ DATA XTAB( 95)/.113914D+02/, YTAB( 95)/ .2174722536864140D-57/ DATA XTAB( 96)/.104325D+02/, YTAB( 96)/ .2908914569100026D-48/ DATA XTAB( 97)/.102246D+02/, YTAB( 97)/ .2175298337170668D-46/ DATA XTAB( 98)/.104371D+02/, YTAB( 98)/ .2641489538998184D-48/ DATA XTAB( 99)/.107497D+02/, YTAB( 99)/ .3410472643765532D-51/ #endif DATA XTAB( 100)/.400000D+02/, YTAB( 100)/0. / C C START. IF(I.GT.NTAB) GO TO 100 X=XTAB(I) Y=YTAB(I) IF(ISIGN.EQ.2) X=-X IF(ABS(X).GT.0.5D0) GO TO 1 C C (ABS(X).LE.0.5). IF(ISIGN.EQ.2) Y=-Y EREF=Y ECREF=1D0-Y F2REF=0.5D0+0.5D0*Y RETURN C C (ABS(X).GT.0.5). 1 IF(ISIGN.EQ.2) GO TO 2 C (X.GE.0). EREF=1D0-Y ECREF=Y F2REF=1D0-0.5D0*Y RETURN C (X.LT.0). 2 EREF=Y-1D0 ECREF=2D0-Y F2REF=0.5D0*Y RETURN C C ERROR MESSAGE. 100 WRITE(*,2000) I,NTAB RETURN C 2000 FORMAT( // 21H ***** C300GT ... I =, I6, * 17H EXCEEDS NTAB =, I6 / ) END #if defined(CERNLIB_NUMTIME) #endif #endif