From 848ba6bfcefdcd192986dded8d72456eacfb0374 Mon Sep 17 00:00:00 2001 From: Anthony Scemama Date: Sat, 17 Apr 2021 10:36:49 +0200 Subject: [PATCH] Fixed out of bound --- src/csf/sigma_vector.irp.f | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/src/csf/sigma_vector.irp.f b/src/csf/sigma_vector.irp.f index 6703d41b..265d2384 100644 --- a/src/csf/sigma_vector.irp.f +++ b/src/csf/sigma_vector.irp.f @@ -38,7 +38,11 @@ endif ncfg = ncfgpersomo - ncfgprev !detDimperBF = max(1,nint((binom(i,(i+1)/2)))) - dimcsfpercfg = max(1,nint((binom(i-2,(i-2+1)/2)-binom(i-2,((i-2+1)/2)+1)))) + if (i > 2) then + dimcsfpercfg = max(1,nint((binom(i-2,(i-2+1)/2)-binom(i-2,((i-2+1)/2)+1)))) + else + dimcsfpercfg = 1 + endif n_CSF += ncfg * dimcsfpercfg !if(cfg_seniority_index(i+2) == -1) EXIT !if(detDimperBF > maxDetDimPerBF) maxDetDimPerBF = detDimperBF