romg.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file romg.f

C
C
C
  FUNCTION romg(X)
C   ********This gives the eikonal function from a table
C     calculated in the first call
  dimension fr(0:1000)
  DATA i0/0/

        SAVE i0, fr                             !Uzhi

  IF(i0.NE.0) go to 100
  DO 50 i=1,1001
  xr=(i-1)*0.01
  fr(i-1)=omg0(xr)
50  CONTINUE
100 i0=1
  IF(x.GE.10.0) THEN
    romg=0.0
    RETURN
  ENDIF
  ix=int(x*100)
  romg=(fr(ix)*((ix+1)*0.01-x)+fr(ix+1)*(x-ix*0.01))/0.01
  RETURN
  END