← Back to branch summary

~maddevelopers/mg5amcnlo/2.6.3_rwgt

« back to all changes in this revision

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

  • Committer: olivier-mattelaer
  • Date: 2017-08-16 21:28:30 UTC
  • mfrom: (272.1.60 2.5.6)
  • Revision ID: olivier-mattelaer-20170816212830-sfl8jfv8pmpm5l2j
pass to 2.6.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/TestCmdMatchBox/MatchBoxOutput/%TEST%SubProcesses%P1_uux_uux%improve_ps.f
780
780
C     Apply the rescaling
781
781
      DO I=1,NEXTERNAL
782
782
        DO J=1,3
783
 
          NEWP(J,I)=NEWP(J,I)*XSCALE
 
783
C         Consider scaling by x**2 for the first particle so that
 
784
C         the algorithm for numerically solving for XSCALE has a
 
785
C          non-vanishing
 
786
C         derivative in the case that all particle are massless.
 
787
          IF (I.EQ.1) THEN
 
788
            NEWP(J,I)=NEWP(J,I)*XSCALE**2
 
789
          ELSE
 
790
            NEWP(J,I)=NEWP(J,I)*XSCALE
 
791
          ENDIF
784
792
        ENDDO
785
793
      ENDDO
786
794
 
964
972
      RES=ZERO
965
973
      BUFF=ZERO
966
974
 
 
975
C     Consider scaling by x**2 for the first particle so that
 
976
C     the algorithm for numerically solving for XSCALE has a
 
977
C      non-vanishing
 
978
C     derivative in the case that all particle are massless.
 
979
 
967
980
      DO I=1,NEXTERNAL
968
981
        IF (I.LE.NINITIAL) THEN
969
982
          FACTOR=-ONE
970
983
        ELSE
971
984
          FACTOR=ONE
972
985
        ENDIF
973
 
        BUFF=MASSES(I)**2+PVECSQ(I)*X**2
 
986
        IF (I.EQ.1) THEN
 
987
          BUFF=MASSES(I)**2+PVECSQ(I)*X**4
 
988
        ELSE
 
989
          BUFF=MASSES(I)**2+PVECSQ(I)*X**2
 
990
        ENDIF
974
991
        IF (BUFF.LT.ZERO) THEN
975
992
          RES=ZERO
976
993
          ERROR = 1
977
994
          GOTO 800
978
995
        ENDIF
979
996
        IF (DERIVATIVE) THEN
980
 
          RES=RES + FACTOR*((X*PVECSQ(I))/SQRT(BUFF))
 
997
          IF (I.EQ.1) THEN
 
998
            RES=RES + FACTOR*((2*X*PVECSQ(I))/SQRT(BUFF))
 
999
          ELSE
 
1000
            RES=RES + FACTOR*((X*PVECSQ(I))/SQRT(BUFF))
 
1001
          ENDIF
981
1002
        ELSE
982
1003
          RES=RES + FACTOR*SQRT(BUFF)
983
1004
        ENDIF