Skip to content

Commit a050c28

Browse files
heltonmcheltonmc
committed
fix spacing
1 parent 3d2b465 commit a050c28

File tree

12 files changed

+10
-69
lines changed

12 files changed

+10
-69
lines changed

‎src/Bessels.jl

Lines changed: 5 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -2,18 +2,20 @@ module Bessels
22

33
export besselj0
44
export besselj1
5+
56
export bessely0
67
export bessely1
8+
79
export besseli0
810
export besseli0x
911
export besseli1
1012
export besseli1x
13+
1114
export besselk0
1215
export besselk0x
1316
export besselk1
1417
export besselk1x
1518

16-
1719
include("Float32/besseli.jl")
1820
include("Float32/besselj.jl")
1921
include("Float32/besselk.jl")
@@ -32,14 +34,9 @@ include("Float128/besselk.jl")
3234
include("Float128/bessely.jl")
3335
include("Float128/constants.jl")
3436

35-
3637
include("chebyshev.jl")
37-
include("constants.jl")
38-
#include("j0_y0_constants.jl")
39-
#include("j0.jl")
40-
#include("y0.jl")
41-
#include("j1.jl")
42-
38+
include("math_constants.jl")
4339
#include("parse.jl")
40+
#include("hankel.jl")
4441

4542
end

‎src/Float32/besselj.jl

Lines changed: 0 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -25,7 +25,6 @@ function besselj0(x::Float32)
2525
return p
2626
end
2727
end
28-
2928
function besselj1(x::Float32)
3029
x = abs(x)
3130
if iszero(x)

‎src/Float32/bessely.jl

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -25,7 +25,6 @@ function bessely0(x::Float32)
2525
return p
2626
end
2727
end
28-
2928
function bessely1(x::Float32)
3029
if x <= zero(x)
3130
if iszero(x)
@@ -52,4 +51,4 @@ function bessely1(x::Float32)
5251
p = p * sin(xn + x)
5352
return p
5453
end
55-
end
54+
end

‎src/Float32/constants.jl

Lines changed: 0 additions & 22 deletions
Original file line numberDiff line numberDiff line change
@@ -6,13 +6,11 @@ const A_i0(::Type{Float32}) = (
66
-2.37374148058994688156f-2, 4.93052842396707084878f-2, -9.49010970480476444210f-2,
77
1.71620901522208775349f-1, -3.04682672343198398683f-1, 6.76795274409476084995f-1
88
)
9-
109
const B_i0(::Type{Float32}) = (
1110
3.39623202570838634515f-9, 2.26666899049817806459f-8, 2.04891858946906374183f-7,
1211
2.89137052083475648297f-6, 6.88975834691682398426f-5, 3.36911647825569408990f-3,
1312
8.04490411014108831608f-1
1413
)
15-
1614
const A_i1(::Type{Float32}) = (
1715
9.38153738649577178388f-9, -4.44505912879632808065f-8, 2.00329475355213526229f-7,
1816
-8.56872026469545474066f-7, 3.47025130813767847674f-6, -1.32731636560394358279f-5,
@@ -21,93 +19,73 @@ const A_i1(::Type{Float32}) = (
2119
2.47264490306265168283f-2, -5.29459812080949914269f-2, 1.02643658689847095384f-1,
2220
-1.76416518357834055153f-1, 2.52587186443633654823f-1
2321
)
24-
2522
const B_i1(::Type{Float32}) = (
2623
-3.83538038596423702205f-9, -2.63146884688951950684f-8, -2.51223623787020892529f-7,
2724
-3.88256480887769039346f-6, -1.10588938762623716291f-4, -9.76109749136146840777f-3,
2825
7.78576235018280120474f-1
2926
)
30-
31-
3227
const JP_j0(::Type{Float32}) = (
3328
-1.729150680240724f-1, 1.332913422519003f-2, -3.969646342510940f-4,
3429
6.388945720783375f-6, -6.068350350393235f-8
3530
)
36-
3731
const MO_j0(::Type{Float32}) = (
3832
7.978845717621440f-1, -3.355424622293709f-6, -4.969382655296620f-2,
3933
-3.560281861530129f-3, 1.197549369473540f-1, -2.145007480346739f-1,
4034
1.864949361379502f-1, -6.838999669318810f-2
4135
)
42-
4336
const PH_j0(::Type{Float32}) = (
4437
-1.249992184872738f-1, 6.490598792654666f-2, -1.939906941791308f-1,
4538
1.001973420681837f0, -4.974978466280903f0, 1.756221482109099f1,
4639
-3.630592630518434f1, 3.242077816988247f1,
4740
)
48-
49-
5041
const JP32 = (
5142
-3.405537384615824f-2, 1.937383947804541f-3, -4.541343896997497f-5,
5243
6.009061827883699f-7, -4.878788132172128f-9
5344
)
54-
55-
5645
const YP32 = (
5746
4.202369946500099f-2, -2.641785726447862f-3, 6.719543806674249f-5,
5847
-9.496460629917016f-7, 8.061978323326852f-9
5948
)
60-
61-
6249
const MO132 = (
6350
7.978845453073848f-1, 4.976029650847191f-6, 1.493389585089498f-1,
6451
5.435364690523026f-3, -2.102302420403875f-1, 3.138238455499697f-1,
6552
-2.284801500053359f-1,6.913942741265801f-2,
6653
)
67-
6854
const PH132 = (
6955
3.749989509080821f-1, -1.637986776941202f-1, 3.503787691653334f-1,
7056
-1.544842782180211f0, 7.222973196770240f0, -2.485774108720340f1,
7157
5.073465654089319f1, -4.497014141919556f1,
7258
)
73-
7459
const YP_y1(::Type{Float32}) = (
7560
1.707584643733568f-1, -1.584289289821316f-2, 5.344486707214273f-4,
7661
-9.413212653797057f-6, 9.454583683980369f-8
7762
)
78-
7963
const MO_y1(::Type{Float32}) = (
8064
7.978845717621440f-1, -3.355424622293709f-6, -4.969382655296620f-2,
8165
-3.560281861530129f-3, 1.197549369473540f-1, -2.145007480346739f-1,
8266
1.864949361379502f-1, -6.838999669318810f-2
8367
)
84-
8568
const PH_y1(::Type{Float32}) = (
8669
-1.249992184872738f-1, 6.490598792654666f-2, -1.939906941791308f-1,
8770
1.001973420681837f0, -4.974978466280903f0, 1.756221482109099f1,
8871
-3.630592630518434f1, 3.242077816988247f1,
8972
)
90-
9173
const A_k0(::Type{Float32}) = (
9274
1.90451637722020886025f-9, 2.53479107902614945675f-7, 2.28621210311945178607f-5,
9375
1.26461541144692592338f-3, 3.59799365153615016266f-2, 3.44289899924628486886f-1,
9476
-5.35327393233902768720f-1
9577
)
96-
97-
9878
const B_k0(::Type{Float32}) = (
9979
-1.69753450938905987466f-9, 8.57403401741422608519f-9, -4.66048989768794782956f-8,
10080
2.76681363944501510342f-7, -1.83175552271911948767f-6, 1.39498137188764993662f-5,
10181
-1.28495495816278026384f-4, 1.56988388573005337491f-3, -3.14481013119645005427f-2,
10282
2.44030308206595545468f0
10383
)
104-
10584
const A_k1(::Type{Float32}) = (
10685
-2.21338763073472585583f-8, -2.43340614156596823496f-6, -1.73028895751305206302f-4,
10786
-6.97572385963986435018f-3, -1.22611180822657148235f-1, -3.53155960776544875667f-1,
10887
1.52530022733894777053f0
10988
)
110-
11189
const B_k1(::Type{Float32}) = (
11290
2.01504975519703286596f-9, -1.03457624656780970260f-8, 5.74108412545004946722f-8,
11391
-3.50196060308781257119f-7, 2.40648494783721712015f-6, -1.93619797416608296024f-5,

‎src/Float64/besselj.jl

Lines changed: 0 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -55,8 +55,6 @@ function besselj1(x::Float64)
5555
return p * SQ2OPI(T) / sqrt(x)
5656
end
5757
end
58-
59-
6058
function besselj(n::Int, x::Float64)
6159
if n < 0
6260
n = -n

‎src/Float64/besselk.jl

Lines changed: 0 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -58,4 +58,3 @@ function besselk1x(x::Float64)
5858
return chbevl(8.0 / x - 2.0, B_k1(T)) / sqrt(x)
5959
end
6060
end
61-

‎src/Float64/bessely.jl

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -51,7 +51,6 @@ function bessely1(x::Float64)
5151
return p * SQ2OPI(T) / sqrt(x)
5252
end
5353
end
54-
5554
function bessely(n::Int, x)
5655
if n < 0
5756
n = -n
@@ -87,4 +86,4 @@ function bessely(n::Int, x)
8786
end
8887

8988
return sign * an
90-
end
89+
end

‎src/Float64/constants.jl

Lines changed: 0 additions & 24 deletions
Original file line numberDiff line numberDiff line change
@@ -7,7 +7,6 @@ const RQ_j0(::Type{Float64}) = (
77
2.11277520115489217587E12, 1.11855537045356834862E10, 4.84409658339962045305E7,
88
1.73785401676374683123E5, 4.99563147152651017219E2, 1.00000000000000000000E0
99
)
10-
1110
const PP_j0(::Type{Float64}) = (
1211
9.99999999999999997821E-1, 5.30324038235394892183E0, 8.74716500199817011941E0,
1312
5.44725003058768775090E0, 1.23953371646414299388E0, 8.28352392107440799803E-2,
@@ -18,7 +17,6 @@ const PQ_j0(::Type{Float64}) = (
1817
5.47097740330417105182E0, 1.25352743901058953537E0, 8.56288474354474431428E-2,
1918
9.24408810558863637013E-4
2019
)
21-
2220
const YP_y0(::Type{Float64}) = (
2321
-1.84950800436986690637E16, 4.42733268572569800351E16, -3.46628303384729719441E15,
2422
8.75906394395366999549E13, -9.82136065717911466409E11, 5.43526477051876500413E9,
@@ -29,7 +27,6 @@ const YQ_y0(::Type{Float64}) = (
2927
8.64002487103935000337E10, 2.68919633393814121987E8, 6.26107330137134956842E5,
3028
1.04128353664259848412E3, 1.00000000000000000000E0
3129
)
32-
3330
const PP_y0(::Type{Float64}) = (
3431
9.99999999999999997821E-1, 5.30324038235394892183E0, 8.74716500199817011941E0,
3532
5.44725003058768775090E0, 1.23953371646414299388E0, 8.28352392107440799803E-2,
@@ -40,20 +37,16 @@ const PQ_y0(::Type{Float64}) = (
4037
5.47097740330417105182E0, 1.25352743901058953537E0, 8.56288474354474431428E-2,
4138
9.24408810558863637013E-4
4239
)
43-
4440
const QP_y0(::Type{Float64}) = (
4541
-6.05014350600728481186E0, -5.14105326766599330220E1, -1.47077505154951170175E2,
4642
-1.77681167980488050595E2, -9.32060152123768231369E1, -1.95539544257735972385E1,
4743
-1.28252718670509318512E0, -1.13663838898469149931E-2
4844
)
49-
5045
const QQ_y0(::Type{Float64}) = (
5146
2.42005740240291393179E2, 2.06209331660327847417E3, 5.93072701187316984827E3,
5247
7.24046774195652478189E3, 3.88240183605401609683E3, 8.56430025976980587198E2,
5348
6.43178256118178023184E1, 1.00000000000000000000E0
5449
)
55-
56-
5750
const YP_y1(::Type{Float64}) = (
5851
-7.78877196265950026825E17, 2.02439475713594898196E17,
5952
-8.12770255501325109621E15, 1.14509511541823727583E14,
@@ -66,12 +59,10 @@ const YQ_y1(::Type{Float64}) = (
6659
2.35564092943068577943E5, 5.94301592346128195359E2,
6760
1.00000000000000000000E0
6861
)
69-
7062
const RP_j1(::Type{Float64}) = (
7163
3.68295732863852883286E15, -7.27494245221818276015E13,
7264
4.52228297998194034323E11, -8.99971225705559398224E8
7365
)
74-
7566
const RQ_j1(::Type{Float64}) = (
7667
5.32278620332680085395E18, 8.95222336184627338078E16,
7768
7.84369607876235854894E14, 4.74914122079991414898E12,
@@ -85,28 +76,24 @@ const PP_j1(::Type{Float64}) = (
8576
1.12719608129684925192E0, 7.31397056940917570436E-2,
8677
7.62125616208173112003E-4
8778
)
88-
8979
const PQ_j1(::Type{Float64}) = (
9080
9.99999999999999997461E-1, 5.20982848682361821619E0,
9181
8.39985554327604159757E0, 5.07386386128601488557E0,
9282
1.10514232634061696926E0, 6.88455908754495404082E-2,
9383
5.71323128072548699714E-4
9484
)
95-
9685
const QP_j1(::Type{Float64}) = (
9786
2.52070205858023719784E1, 2.11688757100572135698E2,
9887
5.97489612400613639965E2, 7.10856304998926107277E2,
9988
3.66779609360150777800E2, 7.58238284132545283818E1,
10089
4.98213872951233449420E0, 5.10862594750176621635E-2
10190
)
102-
10391
const QQ_j1(::Type{Float64}) = (
10492
3.36093607810698293419E2, 2.82619278517639096600E3,
10593
7.99704160447350683650E3, 9.56231892404756170795E3,
10694
4.98641058337653607651E3, 1.05644886038262816351E3,
10795
7.42373277035675149943E1, 1.00000000000000000000E0
10896
)
109-
11097
const QP_j0(::Type{Float64}) = (
11198
-6.05014350600728481186E0, -5.14105326766599330220E1, -1.47077505154951170175E2,
11299
-1.77681167980488050595E2, -9.32060152123768231369E1, -1.95539544257735972385E1,
@@ -117,7 +104,6 @@ const QQ_j0(::Type{Float64}) = (
117104
7.24046774195652478189E3, 3.88240183605401609683E3, 8.56430025976980587198E2,
118105
6.43178256118178023184E1, 1.00000000000000000000E0
119106
)
120-
121107
const A_i0(::Type{Float64}) = (
122108
-4.41534164647933937950E-18, 3.33079451882223809783E-17, -2.43127984654795469359E-16,
123109
1.71539128555513303061E-15, -1.16853328779934516808E-14, 7.67618549860493561688E-14,
@@ -130,7 +116,6 @@ const A_i0(::Type{Float64}) = (
130116
-2.37374148058994688156E-2, 4.93052842396707084878E-2, -9.49010970480476444210E-2,
131117
1.71620901522208775349E-1, -3.04682672343198398683E-1, 6.76795274409476084995E-1
132118
)
133-
134119
const B_i0(::Type{Float64}) = (
135120
-7.23318048787475395456E-18, -4.83050448594418207126E-18, 4.46562142029675999901E-17,
136121
3.46122286769746109310E-17, -2.82762398051658348494E-16, -3.42548561967721913462E-16,
@@ -142,8 +127,6 @@ const B_i0(::Type{Float64}) = (
142127
2.89137052083475648297E-6, 6.88975834691682398426E-5, 3.36911647825569408990E-3,
143128
8.04490411014108831608E-1
144129
)
145-
146-
147130
const A_i1(::Type{Float64}) = (
148131
2.77791411276104639959E-18, -2.11142121435816608115E-17, 1.55363195773620046921E-16,
149132
-1.10559694773538630805E-15, 7.60068429473540693410E-15, -5.04218550472791168711E-14,
@@ -156,7 +139,6 @@ const A_i1(::Type{Float64}) = (
156139
2.47264490306265168283E-2, -5.29459812080949914269E-2, 1.02643658689847095384E-1,
157140
-1.76416518357834055153E-1, 2.52587186443633654823E-1
158141
)
159-
160142
const B_i1(::Type{Float64}) = (
161143
7.51729631084210481353E-18, 4.41434832307170791151E-18, -4.65030536848935832153E-17,
162144
-3.20952592199342395980E-17, 2.96262899764595013876E-16, 3.30820231092092828324E-16,
@@ -168,15 +150,12 @@ const B_i1(::Type{Float64}) = (
168150
-3.88256480887769039346E-6, -1.10588938762623716291E-4, -9.76109749136146840777E-3,
169151
7.78576235018280120474E-1
170152
)
171-
172-
173153
const A_k0(::Type{Float64}) = (
174154
1.37446543561352307156E-16, 4.25981614279661018399E-14, 1.03496952576338420167E-11,
175155
1.90451637722020886025E-9, 2.53479107902614945675E-7, 2.28621210311945178607E-5,
176156
1.26461541144692592338E-3, 3.59799365153615016266E-2, 3.44289899924628486886E-1,
177157
-5.35327393233902768720E-1
178158
)
179-
180159
const B_k0(::Type{Float64}) = (
181160
5.30043377268626276149E-18, -1.64758043015242134646E-17, 5.21039150503902756861E-17,
182161
-1.67823109680541210385E-16, 5.51205597852431940784E-16, -1.84859337734377901440E-15,
@@ -188,15 +167,12 @@ const B_k0(::Type{Float64}) = (
188167
-1.28495495816278026384E-4, 1.56988388573005337491E-3, -3.14481013119645005427E-2,
189168
2.44030308206595545468E0
190169
)
191-
192-
193170
const A_k1(::Type{Float64}) = (
194171
-7.02386347938628759343E-18, -2.42744985051936593393E-15, -6.66690169419932900609E-13,
195172
-1.41148839263352776110E-10, -2.21338763073472585583E-8, -2.43340614156596823496E-6,
196173
-1.73028895751305206302E-4, -6.97572385963986435018E-3, -1.22611180822657148235E-1,
197174
-3.53155960776544875667E-1, 1.52530022733894777053E0
198175
)
199-
200176
const B_k1(::Type{Float64}) = (
201177
-5.75674448366501715755E-18, 1.79405087314755922667E-17, -5.68946255844285935196E-17,
202178
1.83809354436663880070E-16, -6.05704724837331885336E-16, 2.03870316562433424052E-15,

‎src/constants.jl renamed to ‎src/math_constants.jl

File renamed without changes.

‎test/besseli_test.jl

Lines changed: 1 addition & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -40,7 +40,6 @@ i1_32 = besseli1.(Float32.(x))
4040
@test i1_64[1] isa Float64
4141
@test i1_32[1] isa Float32
4242

43-
4443
# test against SpecialFunctions.jl
4544
@test i1_64 ≈ i1_SpecialFunctions
4645
@test i1_32 ≈ i1_SpecialFunctions
@@ -55,7 +54,6 @@ i1x_32 = besseli1x.(Float32.(x))
5554
@test i1x_64[1] isa Float64
5655
@test i1x_32[1] isa Float32
5756

58-
5957
# test against SpecialFunctions.jl
6058
@test i1x_64 ≈ i1x_SpecialFunctions
61-
@test i1x_32 ≈ i1x_SpecialFunctions
59+
@test i1x_32 ≈ i1x_SpecialFunctions

0 commit comments

Comments
 (0)