From 12e2327660dba4cc66d58bece1b0d5c27d79fd17 Mon Sep 17 00:00:00 2001 From: Tobias Neumann Date: Thu, 28 Mar 2024 16:21:36 +0000 Subject: [PATCH] Remove default 164/ktopanom initial calls + decrease boost of virt pieces relative to real --- src/Integrate/mcfm_vegas_adaptive.f | 9 +++------ 1 file changed, 3 insertions(+), 6 deletions(-) diff --git a/src/Integrate/mcfm_vegas_adaptive.f b/src/Integrate/mcfm_vegas_adaptive.f index 1ecb68d..52b9c1f 100644 --- a/src/Integrate/mcfm_vegas_adaptive.f +++ b/src/Integrate/mcfm_vegas_adaptive.f @@ -544,14 +544,11 @@ c--- ===================================== if (nproc == 164 .or. nproc==169) then iterationStorage(lord,1)%vinfo%damp = [1.5, 0.8] - SETUP(lord,250000) iterationStorage(nloVirt,1)%vinfo%damp = [1.5, 0.8] - SETUP(nloVirt,100000) enableIps(nloReal, 1:2) = .true. iterationStorage(nloReal,1)%vinfo%damp = [1.5, 0.8] - SETUP(nloReal,500000) endif if (nproc == 1281 .or. nproc == 1291 .or. nproc == 1301 .or. @@ -1513,9 +1510,9 @@ c =========================== ! we fake the "real emission" contributions uncertainty ! for the prioritization system here ! to boost the easier to compute contributions - sdMod(nloReal,:) = sdMod(nloReal,:) / 1.5_dp - sdMod(snloAbove,:) = sdMod(snloAbove,:) / 1.5_dp - sdMod(nnloRealAbove,:) = sdMod(nnloRealAbove,:) / 1.5_dp + sdMod(nloReal,:) = sdMod(nloReal,:) / 1.2_dp + sdMod(snloAbove,:) = sdMod(snloAbove,:) / 1.2_dp + sdMod(nnloRealAbove,:) = sdMod(nnloRealAbove,:) / 1.2_dp sdMod(nloResummed,:) = sdMod(nloResummed,:) / 15._dp nextLoc = maxloc(sdMod(:,:))