22 lines
429 B
Julia
22 lines
429 B
Julia
println("### Test: transpose(U) * U == identity")
|
|
|
|
using LinearAlgebra
|
|
|
|
include("../ho_basis.jl")
|
|
|
|
E_max = 30
|
|
Λ = 0
|
|
|
|
println("No of threads = ", Threads.nthreads())
|
|
|
|
@time "Basis" basis = ho_basis_2B(E_max, Λ)
|
|
|
|
println("Basis size = ", basis.dim)
|
|
|
|
@time "Moshinsky brackets" U = Moshinsky_transform(basis)
|
|
|
|
check = transpose(U) * U - I
|
|
|
|
println("Maximum = ", maximum(abs.(check)))
|
|
println("Norm = ", sum(check .* conj(check)))
|