Update lookup table for binomial coefficients in TDMs.jl

nmayhall-vt · Jan 27, 2024 · 0e9ecd6 · 0e9ecd6
1 parent 595a740
commit 0e9ecd6
Showing 1 changed file with 22 additions and 1 deletion.
diff --git a/src/FCI/TDMs.jl b/src/FCI/TDMs.jl
@@ -47,7 +47,7 @@ function compute_creation(no::Integer, bra_na, bra_nb, ket_na, ket_nb, bra_v::Ma
     _binomial = Array{Int,2}(undef,no+1,no+1)
     for i in 0:no
         for j in i:no
-            if binomial(j,i)>1E12
+            if binomial(j,i)>1E15
                 break
             end
             _binomial[j+1,i+1] = binomial(j,i)
@@ -161,6 +161,9 @@ function compute_AA(no::Integer, bra_na, bra_nb, ket_na, ket_nb, bra_v::Matrix,
     _binomial = Array{Int,2}(undef,no+1,no+1)
     for i in 0:no
         for j in i:no
+            if binomial(j,i)>1E15
+                break
+            end
             _binomial[j+1,i+1] = binomial(j,i)
         end
     end
@@ -267,6 +270,9 @@ function compute_Aa(no::Integer, bra_na, bra_nb, ket_na, ket_nb, bra_v::Matrix,
     _binomial = Array{Int,2}(undef,no+1,no+1)
     for i in 0:no
         for j in i:no
+            if binomial(j,i)>1E15
+                break
+            end
             _binomial[j+1,i+1] = binomial(j,i)
         end
     end
@@ -372,6 +378,9 @@ function compute_Ab(no::Integer, bra_na, bra_nb, ket_na, ket_nb, bra_v::Matrix,
     _binomial = Array{Int,2}(undef,no+1,no+1)
     for i in 0:no
         for j in i:no
+            if binomial(j,i)>1E15
+                break
+            end
             _binomial[j+1,i+1] = binomial(j,i)
         end
     end
@@ -486,6 +495,9 @@ function compute_AB(no::Integer, bra_na, bra_nb, ket_na, ket_nb, bra_v::Matrix,
     _binomial = Array{Int,2}(undef,no+1,no+1)
     for i in 0:no
         for j in i:no
+            if binomial(j,i)>1E15
+                break
+            end
             _binomial[j+1,i+1] = binomial(j,i)
         end
     end
@@ -598,6 +610,9 @@ function compute_AAa(no::Integer, bra_na, bra_nb, ket_na, ket_nb, bra_v::Matrix,
     _binomial = Array{Int,2}(undef,no+1,no+1)
     for i in 0:no
         for j in i:no
+            if binomial(j,i)>1E15
+                break
+            end
             _binomial[j+1,i+1] = binomial(j,i)
         end
     end
@@ -699,6 +714,9 @@ function compute_ABa(no::Integer, bra_na, bra_nb, ket_na, ket_nb, bra_v::Matrix,
     _binomial = Array{Int,2}(undef,no+1,no+1)
     for i in 0:no
         for j in i:no
+            if binomial(j,i)>1E15
+                break
+            end
             _binomial[j+1,i+1] = binomial(j,i)
         end
     end
@@ -811,6 +829,9 @@ function compute_ABb(no::Integer, bra_na, bra_nb, ket_na, ket_nb, bra_v::Matrix,
     _binomial = Array{Int,2}(undef,no+1,no+1)
     for i in 0:no
         for j in i:no
+            if binomial(j,i)>1E15
+                break
+            end
             _binomial[j+1,i+1] = binomial(j,i)
         end
     end