← 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

  • Committer: olivier-mattelaer
  • Date: 2020-08-21 09:16:56 UTC
  • mfrom: (284.2.24 python3)
  • Revision ID: olivier-mattelaer-20200821091656-iizux2mj94tkssuo
pass to 2.8.0 (and move to python3 branch)

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
415
415
      ENDDO
416
416
      END
417
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
 
418
460
      SUBROUTINE MG5_1_UPDATE_WL_2_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
419
461
     $ ,OUT)
420
462
      USE MG5_1_POLYNOMIAL_CONSTANTS
535
577
      ENDDO
536
578
      END
537
579
 
538
 
      SUBROUTINE MG5_1_UPDATE_WL_0_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
539
 
     $ ,OUT)
540
 
      USE MG5_1_POLYNOMIAL_CONSTANTS
541
 
      INTEGER I,J,K
542
 
      COMPLEX*16 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
543
 
      COMPLEX*16 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
544
 
      COMPLEX*16 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
545
 
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
546
 
 
547
 
      DO I=1,LCUT_SIZE
548
 
        DO J=1,OUT_SIZE
549
 
          DO K=0,0
550
 
            OUT(J,K,I)=(0.0D0,0.0D0)
551
 
          ENDDO
552
 
          DO K=1,IN_SIZE
553
 
            OUT(J,0,I)=OUT(J,0,I)+A(K,0,I)*B(J,0,K)
554
 
          ENDDO
555
 
        ENDDO
556
 
      ENDDO
557
 
      END
558
 
 
559
 
      SUBROUTINE MP_MG5_1_UPDATE_WL_0_0(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
560
 
     $ ,OUT)
561
 
      USE MG5_1_POLYNOMIAL_CONSTANTS
562
 
      INTEGER I,J,K
563
 
      COMPLEX*32 A(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
564
 
      COMPLEX*32 B(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
565
 
      COMPLEX*32 OUT(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE)
566
 
      INTEGER LCUT_SIZE,IN_SIZE,OUT_SIZE
567
 
 
568
 
      DO I=1,LCUT_SIZE
569
 
        DO J=1,OUT_SIZE
570
 
          DO K=0,0
571
 
            OUT(J,K,I)=CMPLX(0.0E0_16,0.0E0_16,KIND=16)
572
 
          ENDDO
573
 
          DO K=1,IN_SIZE
574
 
            OUT(J,0,I)=OUT(J,0,I)+A(K,0,I)*B(J,0,K)
575
 
          ENDDO
576
 
        ENDDO
577
 
      ENDDO
578
 
      END
579
 
 
580
580
      SUBROUTINE MG5_1_UPDATE_WL_1_1(A,LCUT_SIZE,B,IN_SIZE,OUT_SIZE
581
581
     $ ,OUT)
582
582
      USE MG5_1_POLYNOMIAL_CONSTANTS
599
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)
600
600
            OUT(J,5,I)=OUT(J,5,I)+A(K,1,I)*B(J,1,K)
601
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)
602
604
            OUT(J,7,I)=OUT(J,7,I)+A(K,2,I)*B(J,2,K)
603
 
            OUT(J,8,I)=OUT(J,8,I)+A(K,1,I)*B(J,3,K)+A(K,3,I)*B(J,1,K)
604
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)
605
607
            OUT(J,10,I)=OUT(J,10,I)+A(K,3,I)*B(J,3,K)
606
 
            OUT(J,11,I)=OUT(J,11,I)+A(K,1,I)*B(J,4,K)+A(K,4,I)*B(J,1,K)
607
 
            OUT(J,12,I)=OUT(J,12,I)+A(K,2,I)*B(J,4,K)+A(K,4,I)*B(J,2,K)
608
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)
609
609
            OUT(J,14,I)=OUT(J,14,I)+A(K,4,I)*B(J,4,K)
610
610
          ENDDO
634
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)
635
635
            OUT(J,5,I)=OUT(J,5,I)+A(K,1,I)*B(J,1,K)
636
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)
637
639
            OUT(J,7,I)=OUT(J,7,I)+A(K,2,I)*B(J,2,K)
638
 
            OUT(J,8,I)=OUT(J,8,I)+A(K,1,I)*B(J,3,K)+A(K,3,I)*B(J,1,K)
639
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)
640
642
            OUT(J,10,I)=OUT(J,10,I)+A(K,3,I)*B(J,3,K)
641
 
            OUT(J,11,I)=OUT(J,11,I)+A(K,1,I)*B(J,4,K)+A(K,4,I)*B(J,1,K)
642
 
            OUT(J,12,I)=OUT(J,12,I)+A(K,2,I)*B(J,4,K)+A(K,4,I)*B(J,2,K)
643
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)
644
644
            OUT(J,14,I)=OUT(J,14,I)+A(K,4,I)*B(J,4,K)
645
645
          ENDDO