Temporary solution for strange problem on Opteron, aparently a bug in gcc
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 4 Nov 2004 16:08:37 +0000 (16:08 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 4 Nov 2004 16:08:37 +0000 (16:08 +0000)
HIJING/hijing1_36/vegas.F

index e6bf3cc..b01f1b2 100644 (file)
@@ -197,7 +197,10 @@ C1114  FORMAT(1X,'**************END NOTICE*************')
       END IF    
       IF(D(I,J).LE.1.0D-18) GO TO 24
       XO=DT(J)/D(I,J)
-      R(I)=((XO-ONE)/XO/DLOG(XO))**ALPH
+C     R(I)=((XO-ONE)/XO/DLOG(XO))**ALPH
+C     The line above doesn't work on Opteron, probably because ALPH is
+C     used in DATA. As a temporary solution we use 1.5D0 directly (PH)
+      R(I)=((XO-ONE)/XO/DLOG(XO))**(1.5D0)
 24    RC=RC+R(I)
       RC=RC/XND
       K=0