← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%polynomial.f

mwrge with 2.7.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%polynomial.f
78
78
C     
79
79
C     FUNCTIONS
80
80
C     
81
 
      INTEGER MG5_1_ML5SOINDEX_FOR_BORN_AMP, MG5_1_ML5SOINDEX_FOR_LOOP_
82
 
     $AMP, MG5_1_ML5SQSOINDEX
 
81
      INTEGER MG5_1_ML5SOINDEX_FOR_BORN_AMP,
 
82
     $  MG5_1_ML5SOINDEX_FOR_LOOP_AMP, MG5_1_ML5SQSOINDEX
83
83
C     
84
84
C     GLOBAL VARIABLES
85
85
C     
113
113
        CFTOT=CMPLX(CF_N(COLOR_ID,I)/(ONE*ABS(CF_D(COLOR_ID,I))),ZERO
114
114
     $   ,KIND=8)
115
115
        IF(CF_D(COLOR_ID,I).LT.0) CFTOT=CFTOT*IMAG1
116
 
        CONST(MG5_1_ML5SOINDEX_FOR_BORN_AMP(I))=CONST(MG5_1_ML5SOINDEX_
117
 
     $FOR_BORN_AMP(I))+CFTOT*CONJG(AMP(I))
 
116
        CONST(MG5_1_ML5SOINDEX_FOR_BORN_AMP(I))
 
117
     $   =CONST(MG5_1_ML5SOINDEX_FOR_BORN_AMP(I))+CFTOT*CONJG(AMP(I))
118
118
      ENDDO
119
119
 
120
120
      DO I=1,NAMPSO
125
125
            CONST(I)=CONST(I)*GOODHEL(HELCONFIG)
126
126
          ENDIF
127
127
          CALL MG5_1_MERGE_WL(LOOP_WF,RANK,LCUT_SIZE,CONST(I)
128
 
     $     ,LOOPCOEFS(0,MG5_1_ML5SQSOINDEX(I,MG5_1_ML5SOINDEX_FOR_LOOP_
129
 
     $AMP(COLOR_ID)),LOOP_GROUP_NUMBER))
 
128
     $     ,LOOPCOEFS(0,MG5_1_ML5SQSOINDEX(I
 
129
     $     ,MG5_1_ML5SOINDEX_FOR_LOOP_AMP(COLOR_ID)),LOOP_GROUP_NUMBER)
 
130
     $     )
130
131
        ENDIF
131
132
      ENDDO
132
133
 
197
198
C     
198
199
C     FUNCTIONS
199
200
C     
200
 
      INTEGER MG5_1_ML5SOINDEX_FOR_BORN_AMP, MG5_1_ML5SOINDEX_FOR_LOOP_
201
 
     $AMP, MG5_1_ML5SQSOINDEX
 
201
      INTEGER MG5_1_ML5SOINDEX_FOR_BORN_AMP,
 
202
     $  MG5_1_ML5SOINDEX_FOR_LOOP_AMP, MG5_1_ML5SQSOINDEX
202
203
C     
203
204
C     GLOBAL VARIABLES
204
205
C     
232
233
        CFTOT=CMPLX(CF_N(COLOR_ID,I)/(ONE*ABS(CF_D(COLOR_ID,I))),ZERO
233
234
     $   ,KIND=16)
234
235
        IF(CF_D(COLOR_ID,I).LT.0) CFTOT=CFTOT*IMAG1
235
 
        CONST(MG5_1_ML5SOINDEX_FOR_BORN_AMP(I))=CONST(MG5_1_ML5SOINDEX_
236
 
     $FOR_BORN_AMP(I))+CFTOT*CONJG(AMP(I))
 
236
        CONST(MG5_1_ML5SOINDEX_FOR_BORN_AMP(I))
 
237
     $   =CONST(MG5_1_ML5SOINDEX_FOR_BORN_AMP(I))+CFTOT*CONJG(AMP(I))
237
238
      ENDDO
238
239
 
239
240
      DO I=1,NAMPSO
244
245
            CONST(I)=CONST(I)*GOODHEL(HELCONFIG)
245
246
          ENDIF
246
247
          CALL MP_MG5_1_MERGE_WL(LOOP_WF,RANK,LCUT_SIZE,CONST(I)
247
 
     $     ,LOOPCOEFS(0,MG5_1_ML5SQSOINDEX(I,MG5_1_ML5SOINDEX_FOR_LOOP_
248
 
     $AMP(COLOR_ID)),LOOP_GROUP_NUMBER))
 
248
     $     ,LOOPCOEFS(0,MG5_1_ML5SQSOINDEX(I
 
249
     $     ,MG5_1_ML5SOINDEX_FOR_LOOP_AMP(COLOR_ID)),LOOP_GROUP_NUMBER)
 
250
     $     )
249
251
        ENDIF
250
252
      ENDDO
251
253
 
413
415
      ENDDO
414
416
      END
415
417
 
 
418
      SUBROUTINE MG5_1_UPDATE_WL_0_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
 
419
     $ ,OUT)
 
420
      USE MG5_1_POLYNOMIAL_CONSTANTS
 
421
      INTEGER I,J,K
 
422
      COMPLEX*16 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
 
423
      COMPLEX*16 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
 
424
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
 
425
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
 
426
 
 
427
      DO I=1,LCUT_SIZE
 
428
        DO J=1,OUT_SIZE
 
429
          DO K=0,0
 
430
            OUT(J,K,I)=(0.0D0,0.0D0)
 
431
          ENDDO
 
432
          DO K=1,IN_SIZE
 
433
            OUT(J,0,I)=OUT(J,0,I)+A(K,0,I)*B(J,0,K)
 
434
          ENDDO
 
435
        ENDDO
 
436
      ENDDO
 
437
      END
 
438
 
 
439
      SUBROUTINE MP_MG5_1_UPDATE_WL_0_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
 
440
     $ ,OUT)
 
441
      USE MG5_1_POLYNOMIAL_CONSTANTS
 
442
      INTEGER I,J,K
 
443
      COMPLEX*32 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
 
444
      COMPLEX*32 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
 
445
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
 
446
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
 
447
 
 
448
      DO I=1,LCUT_SIZE
 
449
        DO J=1,OUT_SIZE
 
450
          DO K=0,0
 
451
            OUT(J,K,I)=CMPLX(0.0E0_16,0.0E0_16,KIND=16)
 
452
          ENDDO
 
453
          DO K=1,IN_SIZE
 
454
            OUT(J,0,I)=OUT(J,0,I)+A(K,0,I)*B(J,0,K)
 
455
          ENDDO
 
456
        ENDDO
 
457
      ENDDO
 
458
      END
 
459
 
416
460
      SUBROUTINE MG5_1_UPDATE_WL_2_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
417
461
     $ ,OUT)
418
462
      USE MG5_1_POLYNOMIAL_CONSTANTS
533
577
      ENDDO
534
578
      END
535
579
 
536
 
      SUBROUTINE MG5_1_UPDATE_WL_0_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
537
 
     $ ,OUT)
538
 
      USE MG5_1_POLYNOMIAL_CONSTANTS
539
 
      INTEGER I,J,K
540
 
      COMPLEX*16 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
541
 
      COMPLEX*16 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
542
 
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
543
 
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
544
 
 
545
 
      DO I=1,LCUT_SIZE
546
 
        DO J=1,OUT_SIZE
547
 
          DO K=0,0
548
 
            OUT(J,K,I)=(0.0D0,0.0D0)
549
 
          ENDDO
550
 
          DO K=1,IN_SIZE
551
 
            OUT(J,0,I)=OUT(J,0,I)+A(K,0,I)*B(J,0,K)
552
 
          ENDDO
553
 
        ENDDO
554
 
      ENDDO
555
 
      END
556
 
 
557
 
      SUBROUTINE MP_MG5_1_UPDATE_WL_0_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
558
 
     $ ,OUT)
559
 
      USE MG5_1_POLYNOMIAL_CONSTANTS
560
 
      INTEGER I,J,K
561
 
      COMPLEX*32 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
562
 
      COMPLEX*32 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
563
 
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
564
 
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
565
 
 
566
 
      DO I=1,LCUT_SIZE
567
 
        DO J=1,OUT_SIZE
568
 
          DO K=0,0
569
 
            OUT(J,K,I)=CMPLX(0.0E0_16,0.0E0_16,KIND=16)
570
 
          ENDDO
571
 
          DO K=1,IN_SIZE
572
 
            OUT(J,0,I)=OUT(J,0,I)+A(K,0,I)*B(J,0,K)
573
 
          ENDDO
574
 
        ENDDO
575
 
      ENDDO
576
 
      END
577
 
 
578
580
      SUBROUTINE MG5_1_UPDATE_WL_1_1(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
579
581
     $ ,OUT)
580
582
      USE MG5_1_POLYNOMIAL_CONSTANTS
597
599
            OUT(J,4,I)=OUT(J,4,I)+A(K,0,I)*B(J,4,K)+A(K,4,I)*B(J,0,K)
598
600
            OUT(J,5,I)=OUT(J,5,I)+A(K,1,I)*B(J,1,K)
599
601
            OUT(J,6,I)=OUT(J,6,I)+A(K,1,I)*B(J,2,K)+A(K,2,I)*B(J,1,K)
 
602
            OUT(J,8,I)=OUT(J,8,I)+A(K,1,I)*B(J,3,K)+A(K,3,I)*B(J,1,K)
 
603
            OUT(J,11,I)=OUT(J,11,I)+A(K,1,I)*B(J,4,K)+A(K,4,I)*B(J,1,K)
600
604
            OUT(J,7,I)=OUT(J,7,I)+A(K,2,I)*B(J,2,K)
601
 
            OUT(J,8,I)=OUT(J,8,I)+A(K,1,I)*B(J,3,K)+A(K,3,I)*B(J,1,K)
602
605
            OUT(J,9,I)=OUT(J,9,I)+A(K,2,I)*B(J,3,K)+A(K,3,I)*B(J,2,K)
 
606
            OUT(J,12,I)=OUT(J,12,I)+A(K,2,I)*B(J,4,K)+A(K,4,I)*B(J,2,K)
603
607
            OUT(J,10,I)=OUT(J,10,I)+A(K,3,I)*B(J,3,K)
604
 
            OUT(J,11,I)=OUT(J,11,I)+A(K,1,I)*B(J,4,K)+A(K,4,I)*B(J,1,K)
605
 
            OUT(J,12,I)=OUT(J,12,I)+A(K,2,I)*B(J,4,K)+A(K,4,I)*B(J,2,K)
606
608
            OUT(J,13,I)=OUT(J,13,I)+A(K,3,I)*B(J,4,K)+A(K,4,I)*B(J,3,K)
607
609
            OUT(J,14,I)=OUT(J,14,I)+A(K,4,I)*B(J,4,K)
608
610
          ENDDO
632
634
            OUT(J,4,I)=OUT(J,4,I)+A(K,0,I)*B(J,4,K)+A(K,4,I)*B(J,0,K)
633
635
            OUT(J,5,I)=OUT(J,5,I)+A(K,1,I)*B(J,1,K)
634
636
            OUT(J,6,I)=OUT(J,6,I)+A(K,1,I)*B(J,2,K)+A(K,2,I)*B(J,1,K)
 
637
            OUT(J,8,I)=OUT(J,8,I)+A(K,1,I)*B(J,3,K)+A(K,3,I)*B(J,1,K)
 
638
            OUT(J,11,I)=OUT(J,11,I)+A(K,1,I)*B(J,4,K)+A(K,4,I)*B(J,1,K)
635
639
            OUT(J,7,I)=OUT(J,7,I)+A(K,2,I)*B(J,2,K)
636
 
            OUT(J,8,I)=OUT(J,8,I)+A(K,1,I)*B(J,3,K)+A(K,3,I)*B(J,1,K)
637
640
            OUT(J,9,I)=OUT(J,9,I)+A(K,2,I)*B(J,3,K)+A(K,3,I)*B(J,2,K)
 
641
            OUT(J,12,I)=OUT(J,12,I)+A(K,2,I)*B(J,4,K)+A(K,4,I)*B(J,2,K)
638
642
            OUT(J,10,I)=OUT(J,10,I)+A(K,3,I)*B(J,3,K)
639
 
            OUT(J,11,I)=OUT(J,11,I)+A(K,1,I)*B(J,4,K)+A(K,4,I)*B(J,1,K)
640
 
            OUT(J,12,I)=OUT(J,12,I)+A(K,2,I)*B(J,4,K)+A(K,4,I)*B(J,2,K)
641
643
            OUT(J,13,I)=OUT(J,13,I)+A(K,3,I)*B(J,4,K)+A(K,4,I)*B(J,3,K)
642
644
            OUT(J,14,I)=OUT(J,14,I)+A(K,4,I)*B(J,4,K)
643
645
          ENDDO