Skip to content

Commit

Permalink
Remove old comments and show [no ci]
Browse files Browse the repository at this point in the history
  • Loading branch information
DanielVandH committed Jul 13, 2024
1 parent 77b2cd5 commit ce92ea0
Show file tree
Hide file tree
Showing 2 changed files with 0 additions and 7 deletions.
5 changes: 0 additions & 5 deletions src/derivatives.jl
Original file line number Diff line number Diff line change
Expand Up @@ -55,13 +55,8 @@ function diff(P::SemiclassicalJacobi{T}; dims=1) where {T}
Dₐ₀ᵪᵃ⁺¹¹ᶜ⁺¹ = diff(Pᵗᵃ⁰ᶜ)
Pᵗᵃ⁺¹¹ᶜ⁺¹ = Dₐ₀ᵪᵃ⁺¹¹ᶜ⁺¹.args[1]
Pᵗᵃ⁺¹⁰ᶜ⁺¹ = SemiclassicalJacobi(P.t, P.a + 1, zero(P.b), P.c + 1, Pᵗᵃ⁰ᶜ)
# Rₐ₊₁₁ᵪ₊₁ᵗᵃ⁺¹⁰ᶜ⁺¹ = ApplyArray(inv, Pᵗᵃ⁺¹¹ᶜ⁺¹ \ Pᵗᵃ⁺¹⁰ᶜ⁺¹) #
# Rₐ₊₁₁ᵪ₊₁ᵗᵃ⁺¹⁰ᶜ⁺¹ = Pᵗᵃ⁺¹⁰ᶜ⁺¹ \ Pᵗᵃ⁺¹¹ᶜ⁺¹
Rₐ₊₁₀ᵪ₊₁ᵗᵃ⁺¹¹ᶜ⁺¹ = Pᵗᵃ⁺¹¹ᶜ⁺¹ \ Pᵗᵃ⁺¹⁰ᶜ⁺¹
# Dₐ₋₁ᵪᵃ⁺¹⁰ᶜ⁺¹ = Rₐ₊₁₁ᵪ₊₁ᵗᵃ⁺¹⁰ᶜ⁺¹ * coefficients(Dₐ₀ᵪᵃ⁺¹¹ᶜ⁺¹) * Rᵦₐ₁ᵪᵗᵃ⁰ᶜ # bidigaonalconjugation data and avoid inv
Dₐ₋₁ᵪᵃ⁺¹⁰ᶜ⁺¹ = BidiagonalConjugation(Rₐ₊₁₀ᵪ₊₁ᵗᵃ⁺¹¹ᶜ⁺¹, coefficients(Dₐ₀ᵪᵃ⁺¹¹ᶜ⁺¹), Rᵦₐ₁ᵪᵗᵃ⁰ᶜ, 'U')
# b2 = Vcat(zero(T), zero(T), Dₐ₋₁ᵪᵃ⁺¹⁰ᶜ⁺¹[band(1)])
# b1 = Vcat(zero(T), Dₐ₋₁ᵪᵃ⁺¹⁰ᶜ⁺¹[band(0)])
b2 = Vcat(zero(T), zero(T), supdiagonaldata(Dₐ₋₁ᵪᵃ⁺¹⁰ᶜ⁺¹))
b1 = Vcat(zero(T), diagonaldata(Dₐ₋₁ᵪᵃ⁺¹⁰ᶜ⁺¹))
data = Hcat(b2, b1)'
Expand Down
2 changes: 0 additions & 2 deletions test/test_neg1b.jl
Original file line number Diff line number Diff line change
Expand Up @@ -77,7 +77,6 @@ end
# Why does this take SO long for Ps[4]? Without them this takes 40 s, but with them it takes 10m!
for P in Ps
P === Ps[4] && continue
@show 1
for (idx, g) in enumerate((x -> exp(x) + sin(x), x -> (1 - x) * cos(x^3), x -> 5.0 + (1 - x)))
f = expand(P, g)
for x in LinRange(0, 1, 100)
Expand Down Expand Up @@ -161,7 +160,6 @@ end
f = expand(P, g)
df = diff(f)
for x in LinRange(0, 1, 100)
@show x
@test df[x] dg(x) atol = 1e-5
end
end
Expand Down

0 comments on commit ce92ea0

Please sign in to comment.