Update lookup table for binomial coefficients in TDMs.jl

arnab82 · Jan 27, 2024 · 595a740 · 595a740
1 parent 6a570ab
commit 595a740
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/FCI/TDMs.jl b/src/FCI/TDMs.jl
@@ -47,6 +47,9 @@ 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
+                break
+            end
             _binomial[j+1,i+1] = binomial(j,i)
         end
     end