Skip to content

Commit

Permalink
int(t) => 1.0/t
Browse files Browse the repository at this point in the history
  • Loading branch information
MikaelSlevinsky committed Jun 14, 2017
1 parent 3685fcf commit e14d654
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions src/sincosint.jl
Original file line number Diff line number Diff line change
Expand Up @@ -35,7 +35,7 @@ function sinint(x::Float64)
0.13754880327250272679E-14,
0.10223981202236205703E-17)
elseif t 144.0
invt = inv(t)
invt = 1.0/t
return copysign/2, x) - cos(x) *
@horner(invt, 0.99999999962173909991E0,
0.36451060338631902917E3,
Expand Down Expand Up @@ -72,7 +72,7 @@ function sinint(x::Float64)
0.18241750166645704670E12,
0.15407148148861454434E12)
elseif t < Inf
invt = inv(t)
invt = 1.0/t
return copysign/2, x) - cos(x) / x * (1.0 -
@horner(invt, 0.19999999999999978257E1,
0.22206119380434958727E4,
Expand Down Expand Up @@ -153,7 +153,7 @@ function cosint(x::Float64)
0.37718676301688932926E-12,
0.27706844497155995398E-15)
elseif x 12.0
invt = inv(t)
invt = 1.0/t
return sin(x) * @horner(invt, 0.99999999962173909991E0,
0.36451060338631902917E3,
0.44218548041288440874E5,
Expand Down Expand Up @@ -189,7 +189,7 @@ function cosint(x::Float64)
0.18241750166645704670E12,
0.15407148148861454434E12)
elseif x < Inf
invt = inv(t)
invt = 1.0/t
return sin(x)/x * (1.0 - @horner(invt, 0.19999999999999978257E1,
0.22206119380434958727E4,
0.84749007623988236808E6,
Expand Down

0 comments on commit e14d654

Please sign in to comment.