From 6f10af12d430f4b621943a4c90250551f3f338b9 Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Fri, 12 Apr 2024 16:48:11 -0400
Subject: [PATCH] Moshinsky transform unit tests

---
 test/moshinsky.jl | 21 +++++++++++++++++++++
 1 file changed, 21 insertions(+)
 create mode 100644 test/moshinsky.jl

diff --git a/test/moshinsky.jl b/test/moshinsky.jl
new file mode 100644
index 0000000..a89dd33
--- /dev/null
+++ b/test/moshinsky.jl
@@ -0,0 +1,21 @@
+println("### Test: U' * U == identity")
+
+using LinearAlgebra
+
+include("../ho_basis.jl")
+
+E_max = 15
+Λ = 0
+
+println("No of threads = ", Threads.nthreads())
+
+@time "Basis" Es, n1s, l1s, n2s, l2s = get_2p_basis(E_max)
+
+println("Basis size = ", length(Es))
+
+@time "Moshinsky brackets" U = Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ)
+
+check = conj(U) * U - I
+
+println("Maximum = ", maximum(abs.(check)))
+println("Norm = ", sum(check .* conj(check)))