diff --git a/ho_basis.jl b/ho_basis.jl
index b8d7f22..7641cd5 100644
--- a/ho_basis.jl
+++ b/ho_basis.jl
@@ -181,7 +181,7 @@ function get_jacobi_V2_matrix(V_of_r, basis::ho_basis_2B, μω_global; atol=10^-
 
     V_relative = get_sp_V_matrix(V_relative_elem, basis.n1s, basis.l1s, [basis.n2s, basis.l2s]; dtype=ComplexF64, cache=V_relative_cache) + get_sp_V_matrix(V_relative_elem, basis.n2s, basis.l2s, [basis.n1s, basis.l1s]; dtype=ComplexF64, cache=V_relative_cache)
     U = Moshinsky_transform(basis)
-    V2 =  U' * V_relative * U
+    V2 =  transpose(U) * V_relative * U
 
     return V2
 end
@@ -217,7 +217,7 @@ function get_src_V12_matrix(V_of_r, basis::ho_basis_2B, μω_global; atol=10^-6,
     V_relative = get_sp_V_matrix(V_relative_elem, basis.n1s, basis.l1s, [basis.n2s, basis.l2s]; dtype=ComplexF64, E_max=basis.E_max)
 
     U = Moshinsky_transform(basis)
-    V12 =  U' * V_relative * U
+    V12 =  transpose(U) * V_relative * U
 
     return V12
 end
diff --git a/ho_basis_3body.jl b/ho_basis_3body.jl
index c883e5e..5797bbd 100644
--- a/ho_basis_3body.jl
+++ b/ho_basis_3body.jl
@@ -31,7 +31,7 @@ V_relative_elem(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; μω_gen=μω_global
 @time "V2" V2 = get_sp_V_matrix(V_elem, basis.n2s, basis.l2s, [basis.n1s, basis.l1s]; E_max=E_max)
 @time "V relative" V_relative = get_sp_V_matrix(V_relative_elem, basis.n1s, basis.l1s, [basis.n2s, basis.l2s]; E_max=E_max)
 @time "Moshinsky brackets" U = Moshinsky_transform(basis)
-@time "V12" V12 =  U' * V_relative * U
+@time "V12" V12 =  transpose(U) * V_relative * U
 
 println("Calculating spectrum")
 @time "H" H = T1 + T2 + T_cross + V1 + V2 + V12
diff --git a/test/moshinsky.jl b/test/moshinsky.jl
index 85019dc..73dc0bf 100644
--- a/test/moshinsky.jl
+++ b/test/moshinsky.jl
@@ -1,4 +1,4 @@
-println("### Test: U' * U == identity")
+println("### Test: transpose(U) * U == identity")
 
 using LinearAlgebra
 
@@ -15,7 +15,7 @@ println("Basis size = ", basis.dim)
 
 @time "Moshinsky brackets" U = Moshinsky_transform(basis)
 
-check = U' * U - I
+check = transpose(U) * U - I
 
 println("Maximum = ", maximum(abs.(check)))
 println("Norm = ", sum(check .* conj(check)))