From d987b022906ad42526a0e061c67b2631c7d640c0 Mon Sep 17 00:00:00 2001
From: ysyapa <ysyapa@ncsu.edu>
Date: Tue, 16 Apr 2024 19:32:40 +0000
Subject: [PATCH] COSMO tested

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

diff --git a/test/cosmo.jl b/test/cosmo.jl
new file mode 100644
index 0000000..8f89c44
--- /dev/null
+++ b/test/cosmo.jl
@@ -0,0 +1,19 @@
+# depends on a fork of alvolya/cosmo
+
+println("COSMO test against Brody et al.")
+
+get_bracket(n, l, N, L, n1, l1, n2, l2, lambda, d) = @ccall "../cosmo/bin/shared.so"._Z24MoshinskyBracketsWrapperiiiiiiiiid(n::Cint,l::Cint,N::Cint,L::Cint,n1::Cint,l1::Cint,n2::Cint,l2::Cint,lambda::Cint,d::Cdouble)::Cdouble
+
+# Table 1 from Buck et al.
+n1 = [0,0,0,0,0,0,2,2,2,2];
+l1 = [0,1,1,2,2,2,2,2,2,2];
+n2 = [0,0,0,0,0,0,1,1,1,1];
+l2 = [0,3,5,2,4,5,3,3,4,4];
+N = [0,0,0,0,1,0,0,1,0,3];
+L = [0,2,1,1,3,5,3,0,2,2];
+n = [0,1,0,0,0,0,1,2,4,0];
+l = [0,0,5,3,1,2,6,5,2,4];
+Λ = [0,2,6,4,3,4,4,5,2,4];
+
+bracs = get_bracket.(n, l, N, L, n1, l1, n2, l2, Λ, 1.0)
+display(bracs)
\ No newline at end of file