diff --git a/NEWS.md b/NEWS.md
index 780bf3cc6a..16e457f554 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -7,6 +7,10 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 ## [0.9.10] – 2023-12-27
 
+### Added
+
+* Compatibility with `RecursiveArrayTools` v3.
+
 ### Fixed
 
 * Fixed mixed array index number in-place `parallel_transport_to!` on real `Circle`, on Julia 1.6.
diff --git a/Project.toml b/Project.toml
index b297892cc5..0574f8379e 100644
--- a/Project.toml
+++ b/Project.toml
@@ -59,7 +59,7 @@ Plots = "1"
 Quaternions = "0.5, 0.6, 0.7"
 Random = "1.6"
 RecipesBase = "1.1"
-RecursiveArrayTools = "2"
+RecursiveArrayTools = "2, 3"
 Requires = "0.5, 1"
 SimpleWeightedGraphs = "1.2"
 SpecialFunctions = "0.8, 0.9, 0.10, 1.0, 2"
diff --git a/test/manifolds/product_manifold.jl b/test/manifolds/product_manifold.jl
index f9233a984e..5bbf9cd3cf 100644
--- a/test/manifolds/product_manifold.jl
+++ b/test/manifolds/product_manifold.jl
@@ -350,12 +350,12 @@ using RecursiveArrayTools: ArrayPartition
         Z = -X
         H1 = riemannian_Hessian(N, p, Y, Z, X)
         H2 = ArrayPartition(
-            [riemannian_Hessian(M, p[i, :], Y[i, :], Z[i, :], X[i, :]) for i in 1:2]...,
+            [riemannian_Hessian(M, p.x[i], Y.x[i], Z.x[i], X.x[i]) for i in 1:2]...,
         )
         @test H1 == H2
         V = ArrayPartition([0.2, 0.0, 0.0], [0.0, 0.0, 0.3])
         W1 = Weingarten(N, p, X, V)
-        W2 = ArrayPartition([Weingarten(M, p[i, :], X[i, :], V[i, :]) for i in 1:2]...)
+        W2 = ArrayPartition([Weingarten(M, p.x[i], X.x[i], V.x[i]) for i in 1:2]...)
         @test W1 == W2
     end
     @testset "Manifold volume" begin