Skip to content

Commit

Permalink
Revise type stability fixes in Lobatto and order condition functions …
Browse files Browse the repository at this point in the history
…(replace adjoint with transpose).
  • Loading branch information
michakraus committed Nov 28, 2020
1 parent d7538f1 commit a89866a
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 4 deletions.
4 changes: 2 additions & 2 deletions src/order_conditions.jl
Original file line number Diff line number Diff line change
Expand Up @@ -57,7 +57,7 @@ function solve_simplifying_assumption_c(c::AbstractVector{T}) where {T}
M \ r
end

T.(vcat([row(i)' for i in 1:s]...))
vcat([transpose(row(i)) for i in 1:s]...)
end

@doc raw"""
Expand All @@ -77,5 +77,5 @@ function solve_simplifying_assumption_d(b::AbstractVector{T}, c::AbstractVector{
M \ r
end

T.(hcat([row(j) for j in 1:s]...))
hcat([row(j) for j in 1:s]...)
end
4 changes: 2 additions & 2 deletions src/tableaus/lobatto.jl
Original file line number Diff line number Diff line change
Expand Up @@ -89,7 +89,7 @@ function get_lobatto_c_coefficients(s, T=BigFloat)
M \ r
end

T.(hcat(b[1] * ones(T,s), vcat([row(i)' for i in 1:s]...)))
hcat(b[1] * ones(T,s), vcat([transpose(row(i)) for i in 1:s]...))
end

@doc raw"""
Expand All @@ -113,7 +113,7 @@ function get_lobatto_c̄_coefficients(s, T=BigFloat)
M \ r
end

T.(hcat(vcat([row(i)' for i in 1:s]...), zeros(T,s)))
hcat(vcat([transpose(row(i)) for i in 1:s]...), zeros(T,s))
end

get_lobatto_d_coefficients(s, T=BigFloat) = (get_lobatto_c_coefficients(s,T) .+ get_lobatto_c̄_coefficients(s,T)) ./ 2
Expand Down

0 comments on commit a89866a

Please sign in to comment.