Skip to content

Commit

Permalink
fix coordinate allocation
Browse files Browse the repository at this point in the history
  • Loading branch information
mateuszbaran committed Oct 1, 2024
1 parent 1ee0b43 commit cf59553
Showing 1 changed file with 8 additions and 8 deletions.
16 changes: 8 additions & 8 deletions src/bases.jl
Original file line number Diff line number Diff line change
Expand Up @@ -303,15 +303,15 @@ function allocate_result(
f::typeof(get_coordinates),
p,
X,
basis::AbstractBasis,
)
T = allocate_result_type(M, f, (p, X))
basis::AbstractBasis{𝔽},
) where {𝔽}
T = coordinate_eltype(M, p, 𝔽)
return allocate_coordinates(M, p, T, number_of_coordinates(M, basis))
end

@inline function allocate_result_type(
M::AbstractManifold,
f::Union{typeof(get_coordinates),typeof(get_vector)},
f::typeof(get_vector),
args::Tuple{Any,Vararg{Any}},
)
apf = allocation_promotion_function(M, f, args)
Expand Down Expand Up @@ -559,8 +559,8 @@ function _get_coordinates(M::AbstractManifold, p, X, B::DefaultOrthonormalBasis)
return get_coordinates_orthonormal(M, p, X, number_system(B))
end
function get_coordinates_orthonormal(M::AbstractManifold, p, X, N)
Y = allocate_result(M, get_coordinates, p, X, DefaultOrthonormalBasis(N))
return get_coordinates_orthonormal!(M, Y, p, X, N)
c = allocate_result(M, get_coordinates, p, X, DefaultOrthonormalBasis(N))
return get_coordinates_orthonormal!(M, c, p, X, N)
end

function _get_coordinates(M::AbstractManifold, p, X, B::DiagonalizingOrthonormalBasis)
Expand All @@ -572,8 +572,8 @@ function get_coordinates_diagonalizing(
X,
B::DiagonalizingOrthonormalBasis,
)
Y = allocate_result(M, get_coordinates, p, X, B)
return get_coordinates_diagonalizing!(M, Y, p, X, B)
c = allocate_result(M, get_coordinates, p, X, B)
return get_coordinates_diagonalizing!(M, c, p, X, B)
end

function _get_coordinates(M::AbstractManifold, p, X, B::CachedBasis)
Expand Down

0 comments on commit cf59553

Please sign in to comment.