From 9a0c3935faf501daa8ceeb9ffa884c3c7548e985 Mon Sep 17 00:00:00 2001 From: vijay gopal chilkuri Date: Tue, 17 Mar 2020 10:44:50 +0100 Subject: [PATCH] Added warning for when the number of extradiagonal terms exceeds the max value. --- examples/tJ_4x4.inp | 2 +- src/extra_diag.irp.f | 11 +++++++++++ 2 files changed, 12 insertions(+), 1 deletion(-) diff --git a/examples/tJ_4x4.inp b/examples/tJ_4x4.inp index da2bc07..a193afb 100644 --- a/examples/tJ_4x4.inp +++ b/examples/tJ_4x4.inp @@ -1,5 +1,5 @@ 16 -11 +33 4 1 0 diff --git a/src/extra_diag.irp.f b/src/extra_diag.irp.f index bd1a2b6..8bbfd04 100644 --- a/src/extra_diag.irp.f +++ b/src/extra_diag.irp.f @@ -96,6 +96,17 @@ tistart=tistart+1 enddo + +! Check for large number of extradiagonal terms + if(detfound + count1 > maxlien)then + write(6,*)" ====================================================" + write(6,*)" *****WARNING*****: The number of extradiagonal terms & + exceedes maximum dimension of :",maxlien, & + " \n The calculation will be unreliable. Decrease the number & + of extradiagonal terms (nnz: second row in input) to solve & + this problem or increase maxlien." + write(6,*)" ====================================================" + endif Touch foundet foundetadr detfound foundadd foundaddh foundetdmat det deth call adrfull()