← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%COLLIER_interface.f

mwrge with 2.7.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%COLLIER_interface.f
171
171
            IF (NCALLS_IN_CACHE(I).NE.-1.AND..NOT. ( NCOLLIERCALLS(I)
172
172
     $       .EQ.NCALLS_IN_CACHE(I).OR. ( CTMODERUN.EQ.
173
173
     $       -1.AND..NOT.COLLIERUSEINTERNALSTABILITYTEST.AND.NROTATIONS
174
 
     $_DP.GT.0.AND.MOD(NCALLS_IN_CACHE(I),2).EQ.0.AND.(NCALLS_IN_CACHE(
175
 
     $I)/2).EQ.NCOLLIERCALLS(I) ) ) ) THEN
 
174
     $_DP.GT.0.AND.MOD(NCALLS_IN_CACHE(I),2)
 
175
     $       .EQ.0.AND.(NCALLS_IN_CACHE(I)/2).EQ.NCOLLIERCALLS(I) ) ) )
 
176
     $        THEN
176
177
              WRITE(*,*) 'WARNING: A consistency check in MadLoop'
177
178
     $         //' failed and, for safety, forced MadLoop to'
178
179
     $         //' reinitialize the global cache of COLLIER. Report'
263
264
C               The expression below is like taking the absolute value
264
265
C                when summing errors linearly
265
266
C               TIRCOEFSERRORS(I,K)=(TIR_COEFS_ERROR(CURR_RANK,ID)/MaxC
266
 
C               oefForRank(CURR_RANK))*DCMPLX( ABS(DBLE(TIRCOEFS(I,K)))
267
 
C               ,ABS(DIMAG(TIRCOEFS(I,K))) )
 
267
C               oefForRank(CURR_RANK))*DCMPLX(
 
268
C                ABS(DBLE(TIRCOEFS(I,K))),ABS(DIMAG(TIRCOEFS(I,K))) )
268
269
C               But empirically, I observed that retaining the
269
270
C                original complex phase leads to slightly more
270
271
C                accurate estimates
486
487
C             The expression below is like taking the absolute value
487
488
C              when summing errors linearly
488
489
C             TIRCOEFSERRORS(I,K)=(TNtenerr(CURR_RANK)/MaxCoefForRank(C
489
 
C             URR_RANK))*DCMPLX( ABS(DBLE(TIRCOEFS(I,K))),ABS(DIMAG(TIR
490
 
C             COEFS(I,K))) )                    
 
490
C             URR_RANK))*DCMPLX(
 
491
C              ABS(DBLE(TIRCOEFS(I,K))),ABS(DIMAG(TIRCOEFS(I,K))) )                     
491
492
C             But empirically, I observed that retaining the original
492
493
C              complex phase leads to slightly more accurate estimates
493
494
              TIRCOEFSERRORS(I,K)=(TNTENERR(CURR_RANK)
595
596
      LOGICAL COLLIERUSEINTERNALSTABILITYTEST_BU
596
597
      INTEGER USERHEL_BU, SQSO_TARGET_BU, COLLIERMODE_BU,CTMODERUN_BU
597
598
      COMMON/ML5_0_COLLIER_CACHE_RELEVANT_PARAMS/USERHEL_BU
598
 
     $ ,SQSO_TARGET_BU,COLLIERMODE_BU,CTMODERUN_BU,COLLIERUSEINTERNALST
599
 
     $ABILITYTEST_BU
 
599
     $ ,SQSO_TARGET_BU,COLLIERMODE_BU,CTMODERUN_BU
 
600
     $ ,COLLIERUSEINTERNALSTABILITYTEST_BU
600
601
 
601
602
C     The common blocks below are to retrieve the necessary
602
603
C      information about
603
 
C     MadLoop running mode and store it in the sCOLLIER_CACHE_RELEVANT_
604
 
C     PARAMS common block.
 
604
C     MadLoop running mode and store it in the
 
605
C      sCOLLIER_CACHE_RELEVANT_PARAMS common block.
605
606
 
606
607
      INCLUDE 'MadLoopParams.inc'
607
608
      INCLUDE 'unique_id.inc'
670
671
        USERHEL_BU         = USERHEL
671
672
        SQSO_TARGET_BU     = SQSO_TARGET
672
673
        COLLIERMODE_BU     = COLLIERMODE
673
 
        COLLIERUSEINTERNALSTABILITYTEST_BU = COLLIERUSEINTERNALSTABILIT
674
 
     $YTEST
 
674
        COLLIERUSEINTERNALSTABILITYTEST_BU =
 
675
     $    COLLIERUSEINTERNALSTABILITYTEST
675
676
        CTMODERUN_BU       = CTMODERUN
676
677
        CALL SWITCHONCACHE_CLL((UNIQUE_ID-1)*4+1)
677
678
        COLLIER_CACHE_ACTIVE = 1
718
719
        USERHEL_BU         = USERHEL
719
720
        SQSO_TARGET_BU     = SQSO_TARGET
720
721
        COLLIERMODE_BU     = COLLIERMODE
721
 
        COLLIERUSEINTERNALSTABILITYTEST_BU = COLLIERUSEINTERNALSTABILIT
722
 
     $YTEST
 
722
        COLLIERUSEINTERNALSTABILITYTEST_BU =
 
723
     $    COLLIERUSEINTERNALSTABILITYTEST
723
724
        CTMODERUN_BU       = CTMODERUN
724
725
        IF (COLLIERGLOBALCACHE.EQ.-1) THEN
725
726
          CALL INITCACHESYSTEM_CLL(N_CACHES*NPROCS,MAXNEXTERNAL)