Skip to content

Commit c82a41c

Browse files
authored
clean up whitespaces (#180)
1 parent cb622ba commit c82a41c

File tree

4 files changed

+28
-28
lines changed

4 files changed

+28
-28
lines changed

src/PaduaTransform.jl

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -123,7 +123,7 @@ Pre-plan a Padua Transform.
123123
"""
124124
function plan_paduatransform!(::Type{T},N::Integer,lex) where T
125125
n=Int(cld(-3+sqrt(1+8N),2))
126-
if N  ((n+1)*(n+2))÷2
126+
if N ((n+1)*(n+2))÷2
127127
error("Padua transforms can only be applied to vectors of length (n+1)*(n+2)/2.")
128128
end
129129
PaduaTransformPlan(Array{T}(undef,n+2,n+1),FFTW.plan_r2r!(Array{T}(undef,n+2,n+1),FFTW.REDFT00),lex)

src/chebyshevtransform.jl

Lines changed: 15 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -364,29 +364,29 @@ plan_chebyshevutransform(x::AbstractArray, dims...; kws...) = plan_chebyshevutra
364364

365365
@inline function _chebu1_prescale!(_, x::AbstractVector{T}) where T
366366
n = length(x)
367-
for k=1:n # sqrt(1-x_j^2) weight
367+
for k=1:n # sqrt(1-x_j^2) weight
368368
x[k] *= sinpi(one(T)/(2n) + (k-one(T))/n)/n
369369
end
370370
x
371371
end
372372

373373
@inline function _chebu1_postscale!(_, x::AbstractVector{T}) where T
374374
n = length(x)
375-
for k=1:n # sqrt(1-x_j^2) weight
375+
for k=1:n # sqrt(1-x_j^2) weight
376376
x[k] /= sinpi(one(T)/(2n) + (k-one(T))/n)/n
377377
end
378378
x
379379
end
380380

381381
function *(P::ChebyshevUTransformPlan{T,1,K,true}, x::AbstractVector{T}) where {T,K}
382-
length(x)  1 && return x
382+
length(x) 1 && return x
383383
_chebu1_prescale!(P.plan.region, x)
384384
P.plan * x
385385
end
386386

387387
function mul!(y::AbstractVector{T}, P::ChebyshevUTransformPlan{T,1,K,false}, x::AbstractVector{T}) where {T,K}
388388
n = length(x)
389-
length(x)  1 && return copyto!(y, x)
389+
length(x) 1 && return copyto!(y, x)
390390
_chebu1_prescale!(P.plan.region, x)
391391
_plan_mul!(y, P.plan, x)
392392
_chebu1_postscale!(P.plan.region, x)
@@ -396,7 +396,7 @@ end
396396
@inline function _chebu2_prescale!(_, x::AbstractVector{T}) where T
397397
n = length(x)
398398
c = one(T)/ (n+1)
399-
for k=1:n # sqrt(1-x_j^2) weight
399+
for k=1:n # sqrt(1-x_j^2) weight
400400
x[k] *= sinpi(k*c)
401401
end
402402
x
@@ -405,22 +405,22 @@ end
405405
@inline function _chebu2_postscale!(_, x::AbstractVector{T}) where T
406406
n = length(x)
407407
c = one(T)/ (n+1)
408-
@inbounds for k=1:n # sqrt(1-x_j^2) weight
408+
@inbounds for k=1:n # sqrt(1-x_j^2) weight
409409
x[k] /= sinpi(k*c)
410410
end
411411
x
412412
end
413413

414414
function *(P::ChebyshevUTransformPlan{T,2,K,true}, x::AbstractVector{T}) where {T,K}
415415
n = length(x)
416-
n  1 && return x
416+
n 1 && return x
417417
_chebu2_prescale!(P.plan.region, x)
418418
lmul!(one(T)/ (n+1), P.plan * x)
419419
end
420420

421421
function mul!(y::AbstractVector{T}, P::ChebyshevUTransformPlan{T,2,K,false}, x::AbstractVector{T}) where {T,K}
422422
n = length(x)
423-
n  1 && return copyto!(y, x)
423+
n 1 && return copyto!(y, x)
424424
_chebu2_prescale!(P.plan.region, x)
425425
_plan_mul!(y, P.plan, x)
426426
_chebu2_postscale!(P.plan.region, x)
@@ -485,14 +485,14 @@ plan_ichebyshevutransform(x::AbstractArray, dims...; kws...) = plan_ichebyshevut
485485

486486
function _ichebyu1_postscale!(_, x::AbstractVector{T}) where T
487487
n = length(x)
488-
@inbounds for k=1:n # sqrt(1-x_j^2) weight
488+
@inbounds for k=1:n # sqrt(1-x_j^2) weight
489489
x[k] /= 2sinpi(one(T)/(2n) + (k-one(T))/n)
490490
end
491491
x
492492
end
493493
function *(P::IChebyshevUTransformPlan{T,1,K,true}, x::AbstractVector{T}) where {T<:fftwNumber,K}
494494
n = length(x)
495-
n  1 && return x
495+
n 1 && return x
496496

497497
x = P.plan * x
498498
_ichebyu1_postscale!(P.plan.region, x)
@@ -501,7 +501,7 @@ end
501501
function mul!(y::AbstractVector{T}, P::IChebyshevUTransformPlan{T,1,K,false}, x::AbstractVector{T}) where {T<:fftwNumber,K}
502502
n = length(x)
503503
length(y) == n || throw(DimensionMismatch("output must match dimension"))
504-
n  1 && return x
504+
n 1 && return x
505505

506506
_plan_mul!(y, P.plan, x)
507507
_ichebyu1_postscale!(P.plan.region, y)
@@ -510,7 +510,7 @@ end
510510
function _ichebu2_rescale!(_, x::AbstractVector{T}) where T
511511
n = length(x)
512512
c = one(T)/ (n+1)
513-
for k=1:n # sqrt(1-x_j^2) weight
513+
for k=1:n # sqrt(1-x_j^2) weight
514514
x[k] /= sinpi(k*c)
515515
end
516516
ldiv!(2, x)
@@ -519,7 +519,7 @@ end
519519

520520
function *(P::IChebyshevUTransformPlan{T,2,K,true}, x::AbstractVector{T}) where {T<:fftwNumber,K}
521521
n = length(x)
522-
n  1 && return x
522+
n 1 && return x
523523

524524
x = P.plan * x
525525
_ichebu2_rescale!(P.plan.region, x)
@@ -528,7 +528,7 @@ end
528528
function mul!(y::AbstractVector{T}, P::IChebyshevUTransformPlan{T,2,K,false}, x::AbstractVector{T}) where {T<:fftwNumber,K}
529529
n = length(x)
530530
length(y) == n || throw(DimensionMismatch("output must match dimension"))
531-
n  1 && return x
531+
n 1 && return x
532532

533533
_plan_mul!(y, P.plan, x)
534534
_ichebu2_rescale!(P.plan.region, y)
@@ -679,4 +679,4 @@ end
679679
*(P::IChebyshevTransformPlan{T,1,Nothing,true,N,R}, x::AbstractVector{T}) where {T,N,R} =
680680
copyto!(x, IChebyshevTransformPlan{T,1,Nothing,false,N,R}() * x)
681681
# *(P::IChebyshevTransformPlan{T,SECONDKIND,false,Nothing}, x::AbstractVector{T}) where T =
682-
# IChebyshevTransformPlan{T,SECONDKIND,true,Nothing}() * copy(x)
682+
# IChebyshevTransformPlan{T,SECONDKIND,true,Nothing}() * copy(x)

test/chebyshevtests.jl

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -191,7 +191,7 @@ using FastTransforms, Test
191191
@test ichebyshevtransform(X,Val(2),2) ichebyshevtransform!(copy(X),Val(2),2) hcat(ichebyshevtransform.([X[k,:] for k=axes(X,1)],Val(2))...)'
192192

193193
@test @inferred(ichebyshevtransform(X)) @inferred(ichebyshevtransform!(copy(X))) ichebyshevtransform(ichebyshevtransform(X,1),2)
194-
@test @inferred(ichebyshevtransform(X,Val(2))) @inferred(ichebyshevtransform!(copy(X),Val(2))) ichebyshevtransform(ichebyshevtransform(X,Val(2),1),Val(2),2)
194+
@test @inferred(ichebyshevtransform(X,Val(2))) @inferred(ichebyshevtransform!(copy(X),Val(2))) ichebyshevtransform(ichebyshevtransform(X,Val(2),1),Val(2),2)
195195

196196
@test ichebyshevtransform(chebyshevtransform(X)) X
197197
@test chebyshevtransform(ichebyshevtransform(X)) X
@@ -208,7 +208,7 @@ using FastTransforms, Test
208208
@test chebyshevtransform([1,2,3], Val(2)) == chebyshevtransform([1.,2,3], Val(2))
209209
@test ichebyshevtransform([1,2,3]) == ichebyshevtransform([1.,2,3])
210210
@test ichebyshevtransform([1,2,3], Val(2)) == ichebyshevtransform([1.,2,3], Val(2))
211-
211+
212212
@test chebyshevutransform([1,2,3]) == chebyshevutransform([1.,2,3])
213213
@test chebyshevutransform([1,2,3], Val(2)) == chebyshevutransform([1.,2,3], Val(2))
214214
@test ichebyshevutransform([1,2,3]) == ichebyshevutransform([1.,2,3])

test/clenshawtests.jl

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -14,9 +14,9 @@ import FastTransforms: clenshaw, clenshaw!, forwardrecurrence!, forwardrecurrenc
1414
@test @inferred(clenshaw(elty[],1)) zero(elty)
1515

1616
x = elty[1,0,0.1]
17-
@test @inferred(clenshaw(c,x)) @inferred(clenshaw!(c,copy(x)))
17+
@test @inferred(clenshaw(c,x)) @inferred(clenshaw!(c,copy(x)))
1818
@inferred(clenshaw!(c,x,similar(x)))
19-
@inferred(clenshaw(cf,x)) @inferred(clenshaw!(cf,copy(x)))
19+
@inferred(clenshaw(cf,x)) @inferred(clenshaw!(cf,copy(x)))
2020
@inferred(clenshaw!(cf,x,similar(x))) elty[6,-2,-1.74]
2121

2222
@testset "Strided" begin
@@ -34,7 +34,7 @@ import FastTransforms: clenshaw, clenshaw!, forwardrecurrence!, forwardrecurrenc
3434
@testset "matrix coefficients" begin
3535
c = [1 2; 3 4; 5 6]
3636
@test clenshaw(c,0.1) [clenshaw(c[:,1],0.1), clenshaw(c[:,2],0.1)]
37-
@test clenshaw(c,[0.1,0.2]) [clenshaw(c[:,1], 0.1) clenshaw(c[:,2], 0.1); clenshaw(c[:,1], 0.2) clenshaw(c[:,2], 0.2)]
37+
@test clenshaw(c,[0.1,0.2]) [clenshaw(c[:,1], 0.1) clenshaw(c[:,2], 0.1); clenshaw(c[:,1], 0.2) clenshaw(c[:,2], 0.2)]
3838
end
3939
end
4040
end
@@ -45,19 +45,19 @@ import FastTransforms: clenshaw, clenshaw!, forwardrecurrence!, forwardrecurrenc
4545
@testset "forwardrecurrence!" begin
4646
@test @inferred(forwardrecurrence(N, A, B, C, 1)) == @inferred(forwardrecurrence!(Vector{Int}(undef,N), A, B, C, 1)) == 1:N
4747
@test forwardrecurrence!(Vector{Int}(undef,N), A, B, C, -1) == (-1) .^ (0:N-1) .* (1:N)
48-
@test forwardrecurrence(N, A, B, C, 0.1) forwardrecurrence!(Vector{Float64}(undef,N), A, B, C, 0.1)
48+
@test forwardrecurrence(N, A, B, C, 0.1) forwardrecurrence!(Vector{Float64}(undef,N), A, B, C, 0.1)
4949
sin.((1:N) .* acos(0.1)) ./ sqrt(1-0.1^2)
5050
end
5151

5252
c = [1,2,3]
53-
@test c'forwardrecurrence(3, A, B, C, 0.1) clenshaw([1,2,3], A, B, C, 0.1)
53+
@test c'forwardrecurrence(3, A, B, C, 0.1) clenshaw([1,2,3], A, B, C, 0.1)
5454
1 + (2sin(2acos(0.1)) + 3sin(3acos(0.1)))/sqrt(1-0.1^2)
5555

5656
@testset "matrix coefficients" begin
5757
c = [1 2; 3 4; 5 6]
5858
@test clenshaw(c,A,B,C,0.1) [clenshaw(c[:,1],A,B,C,0.1), clenshaw(c[:,2],A,B,C,0.1)]
59-
@test clenshaw(c,A,B,C,[0.1,0.2]) [clenshaw(c[:,1], A,B,C,0.1) clenshaw(c[:,2], A,B,C,0.1); clenshaw(c[:,1], A,B,C,0.2) clenshaw(c[:,2], A,B,C,0.2)]
60-
end
59+
@test clenshaw(c,A,B,C,[0.1,0.2]) [clenshaw(c[:,1], A,B,C,0.1) clenshaw(c[:,2], A,B,C,0.1); clenshaw(c[:,1], A,B,C,0.2) clenshaw(c[:,2], A,B,C,0.2)]
60+
end
6161
end
6262

6363
@testset "Chebyshev-as-general" begin
@@ -79,10 +79,10 @@ import FastTransforms: clenshaw, clenshaw!, forwardrecurrence!, forwardrecurrenc
7979
@test @inferred(clenshaw(Float64[], A, B, C, 1)) 0.0
8080

8181
x = [1,0,0.1]
82-
@test @inferred(clenshaw(c, A, B, C, x)) @inferred(clenshaw!(c, A, B, C, copy(x)))
82+
@test @inferred(clenshaw(c, A, B, C, x)) @inferred(clenshaw!(c, A, B, C, copy(x)))
8383
@inferred(clenshaw!(c, A, B, C, x, one.(x), similar(x)))
8484
@inferred(clenshaw!(cf, Af, Bf, Cf, x, one.(x),similar(x)))
85-
@inferred(clenshaw([1.,2,3], A, B, C, x))
85+
@inferred(clenshaw([1.,2,3], A, B, C, x))
8686
@inferred(clenshaw!([1.,2,3], A, B, C, copy(x))) [6,-2,-1.74]
8787
end
8888

@@ -136,4 +136,4 @@ import FastTransforms: clenshaw, clenshaw!, forwardrecurrence!, forwardrecurrenc
136136
c = randn(N)
137137
@test clenshaw(c, A, B, C, 0.1) == clenshaw(c, A, Vector(B), C, 0.1)
138138
end
139-
end
139+
end

0 commit comments

Comments
 (0)