← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%loop_matrix.f

  • Committer: olivier-mattelaer
  • Date: 2017-08-14 09:19:26 UTC
  • mfrom: (272.1.55 2.5.6)
  • Revision ID: olivier-mattelaer-20170814091926-d6hz7b55vz7jprdd
pass to 2.6.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%V0_dxu_wp%loop_matrix.f
353
353
      COMPLEX*32 MPW(20,NWAVEFUNCS)
354
354
      COMMON/MP_W/MPW
355
355
 
356
 
      COMPLEX*16 WL(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE
357
 
     $ ,0:NLOOPWAVEFUNCS)
358
 
      COMPLEX*16 PL(0:3,0:NLOOPWAVEFUNCS)
 
356
      COMPLEX*16 WL(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE,
 
357
     $ -1:NLOOPWAVEFUNCS)
 
358
      COMPLEX*16 PL(0:3,-1:NLOOPWAVEFUNCS)
359
359
      COMMON/WL/WL,PL
360
360
 
361
361
      COMPLEX*16 LOOPCOEFS(0:LOOPMAXCOEFS-1,NSQUAREDSO,NLOOPGROUPS)
583
583
 
584
584
C       SETUP OF THE COMMON STARTING EXTERNAL LOOP WAVEFUNCTION
585
585
C       IT IS ALSO PS POINT INDEPENDENT, SO IT CAN BE DONE HERE.
 
586
C       The index -1 is for the charge-conjugated fermions with
 
587
C        flipped fermion flow.
586
588
        DO I=0,3
 
589
          PL(I,-1)=DCMPLX(0.0D0,0.0D0)
587
590
          PL(I,0)=DCMPLX(0.0D0,0.0D0)
588
591
        ENDDO
589
592
        DO I=1,MAXLWFSIZE
590
593
          DO J=0,LOOPMAXCOEFS-1
591
594
            DO K=1,MAXLWFSIZE
 
595
              WL(I,J,K,-1)=(0.0D0,0.0D0)
592
596
              IF(I.EQ.K.AND.J.EQ.0) THEN
593
597
                WL(I,J,K,0)=(1.0D0,0.0D0)
594
598
              ELSE
603
607
        ENDIF
604
608
      ENDIF
605
609
 
 
610
C     This is the chare conjugate version of the unit 4-currents in
 
611
C      the canonical cartesian basis.
 
612
C     This, for now, is only defined for 4-fermionic currents.
 
613
      WL(1,0,2,-1) = DCMPLX(-1.0D0,0.0D0)
 
614
      WL(2,0,1,-1) = DCMPLX(1.0D0,0.0D0)
 
615
      WL(3,0,4,-1) = DCMPLX(1.0D0,0.0D0)
 
616
      WL(4,0,3,-1) = DCMPLX(-1.0D0,0.0D0)
 
617
 
606
618
C     Make sure that lorentz rotation tests are not used if there is
607
619
C      external loop wavefunction of spin 2 and that one specific
608
620
C      helicity is asked
1927
1939
          ENDIF
1928
1940
        ENDIF
1929
1941
 
 
1942
C       If NaN are present in the evaluation, automatically set the
 
1943
C        accuracy to 1.0d99.
 
1944
        DO I=1,3
 
1945
          DO J=1,MAXSTABILITYLENGTH
 
1946
            IF (ISNAN(FULLLIST(I,K,J))) THEN
 
1947
              ACC(K) = 1.0D99
 
1948
            ENDIF
 
1949
          ENDDO
 
1950
        ENDDO
 
1951
 
1930
1952
      ENDDO
1931
1953
 
1932
1954
      END