Merge remote-tracking branch 'origin/master' into mbaran/cholesky-bases

.JuliaFormatter.toml

@@ -8,4 +8,4 @@ whitespace_in_kwargs = false
 remove_extra_newlines = true
 annotate_untyped_fields_with_any = false
 conditional_to_if = false
-ignore = ["tutorials"]
+ignore = ["tutorials", ".git"]

NEWS.md

@@ -5,13 +5,17 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
-## [0.10.12] - unreleased
+## [0.10.12] - 2025-01-10
 
 ### Added
 
 * Orthonormal bases for `CholeskySpace` and `LogCholesky` metric for `SymmetricPositiveDefinite`.
 * `rand` for `CholeskySpace`.
 
+### Changed
+
+* Improved performance of selected `get_vector` and `get_coordinates` methods for complex `Euclidean` manifold.
+
 ## [0.10.11] - 2025-01-02
 
 ### Added

src/manifolds/Euclidean.jl

@@ -251,41 +251,39 @@ function get_coordinates_induced_basis!(
     X,
     ::InducedBasis{ℝ,TangentSpaceType,<:RetractionAtlas},
 )
-    S = representation_size(M)
-    PS = prod(S)
-    copyto!(c, reshape(X, PS))
+    copyto!(c, vec(X))
     return c
 end
 
-function get_coordinates_orthonormal!(M::Euclidean{<:Any,ℂ}, c, ::Any, X, ::RealNumbers)
-    S = representation_size(M)
-    PS = prod(S)
-    c .= [reshape(real.(X), PS)..., reshape(imag(X), PS)...]
+function get_coordinates_orthonormal!(::Euclidean{<:Any,ℂ}, c, ::Any, X, ::RealNumbers)
+    Xvec = vec(X)
+    d = div(length(c), 2)
+    view(c, 1:d) .= real.(Xvec)
+    view(c, (d + 1):(2d)) .= imag.(Xvec)
     return c
 end
 
 function get_coordinates_diagonalizing!(
-    M::Euclidean{<:Any,ℂ},
+    ::Euclidean{<:Any,ℂ},
     c,
     ::Any,
     X,
     ::DiagonalizingOrthonormalBasis{ℝ},
 )
-    S = representation_size(M)
-    PS = prod(S)
-    c .= [reshape(real.(X), PS)..., reshape(imag(X), PS)...]
+    Xvec = vec(X)
+    d = div(length(c), 2)
+    view(c, 1:d) .= real.(Xvec)
+    view(c, (d + 1):(2d)) .= imag.(Xvec)
     return c
 end
 function get_coordinates_diagonalizing!(
-    M::Euclidean{<:Any,𝔽},
+    ::Euclidean{<:Any,𝔽},
     c,
     p,
     X,
     ::DiagonalizingOrthonormalBasis{𝔽},
 ) where {𝔽}
-    S = representation_size(M)
-    PS = prod(S)
-    copyto!(c, reshape(X, PS))
+    copyto!(c, vec(X))
     return c
 end
 
@@ -377,7 +375,7 @@ end
 function get_vector_orthonormal!(M::Euclidean{<:Any,ℂ}, Y, ::Any, c, ::RealNumbers)
     S = representation_size(M)
     N = div(length(c), 2)
-    copyto!(Y, reshape(c[1:N] + im * c[(N + 1):end], S))
+    copyto!(Y, reshape(c[1:N] .+ im .* c[(N + 1):end], S))
     return Y
 end
 function get_vector_diagonalizing!(
@@ -389,7 +387,7 @@ function get_vector_diagonalizing!(
 )
     S = representation_size(M)
     N = div(length(c), 2)
-    copyto!(Y, reshape(c[1:N] + im * c[(N + 1):end], S))
+    copyto!(Y, reshape(c[1:N] .+ im .* c[(N + 1):end], S))
     return Y
 end