X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=HERWIG%2Fherwig6510.f;h=3b38f3402c39d74ba90d162417f81dbf6b49b8e2;hb=e5ef0c3a8d1bac74198c5c84fc0a812890ee5a37;hp=11113180f8f6977f028ec0798932b20de203ef5d;hpb=d08a783265dc8c3ee04e2885de31f9743ba8ae3a;p=u%2Fmrichter%2FAliRoot.git diff --git a/HERWIG/herwig6510.f b/HERWIG/herwig6510.f index 11113180f8f..3b38f3402c3 100644 --- a/HERWIG/herwig6510.f +++ b/HERWIG/herwig6510.f @@ -45908,28 +45908,29 @@ CDECK ID>, HWRGEN. *-- Author : F. James, modified by Mike Seymour *- Split in 3 files by M. Kirsanov. Initial seeds ISEED set in HWUDAT C----------------------------------------------------------------------- - FUNCTION HWRGEN(I) + +C FUNCTION HWRGEN(I) C----------------------------------------------------------------------- C MAIN RANDOM NUMBER GENERATOR C USES METHOD OF l'Ecuyer, (VIA F.JAMES, COMP PHYS COMM 60(1990)329) C----------------------------------------------------------------------- - IMPLICIT NONE - DOUBLE PRECISION HWRGEN - COMMON/HWSEED/ISEED(2) - INTEGER ISEED - INTEGER I,K,IZ -C - K=ISEED(1)/53668 - ISEED(1)=40014*(ISEED(1)-K*53668)-K*12211 - IF (ISEED(1).LT.0) ISEED(1)=ISEED(1)+2147483563 - K=ISEED(2)/52774 - ISEED(2)=40692*(ISEED(2)-K*52774)-K*3791 - IF (ISEED(2).LT.0) ISEED(2)=ISEED(2)+2147483399 - IZ=ISEED(1)-ISEED(2) - IF (IZ.LT.1) IZ=IZ+2147483562 - HWRGEN=DBLE(IZ)*4.656613001013252D-10 +C IMPLICIT NONE +C DOUBLE PRECISION HWRGEN +C COMMON/HWSEED/ISEED(2) +C INTEGER ISEED +C INTEGER I,K,IZ +C +C K=ISEED(1)/53668 +C ISEED(1)=40014*(ISEED(1)-K*53668)-K*12211 +C IF (ISEED(1).LT.0) ISEED(1)=ISEED(1)+2147483563 +C K=ISEED(2)/52774 +C ISEED(2)=40692*(ISEED(2)-K*52774)-K*3791 +C IF (ISEED(2).LT.0) ISEED(2)=ISEED(2)+2147483399 +C IZ=ISEED(1)-ISEED(2) +C IF (IZ.LT.1) IZ=IZ+2147483562 +C HWRGEN=DBLE(IZ)*4.656613001013252D-10 C---> (4.656613001013252D-10 = 1.D0/2147483589) - END +c END CDECK ID>, HWRSET. *CMZ :- -26/04/91 12.42.30 by Federico Carminati *-- Author : F. James, modified by Mike Seymour @@ -60420,6 +60421,7 @@ C---CHECK FOR FATAL ERROR IF (NEGWTS.AND.EVWGT.LT.ZERO) NNEGEV=NNEGEV-1 C---PRINT FIRST MAXPR EVENTS ELSEIF (NEVHEP.LE.MAXPR) THEN + write(6,*) "Calling HWUEPR", NEVHEP, MAXPR CALL HWUEPR END IF END