hdiff output

r31324/mc_gbh.F90 2016-10-16 19:30:12.824488431 +0100 r31323/mc_gbh.F90 2016-10-16 19:30:13.068491670 +0100
196:   !DO I=1,N-1196:   !DO I=1,N-1
197:   !   Y = MIN(DEXP(-(VALUES(I)-VALUES(N))/TEMP(1)), 1.0D0)197:   !   Y = MIN(DEXP(-(VALUES(I)-VALUES(N))/TEMP(1)), 1.0D0)
198:   !   X = X + Y198:   !   X = X + Y
199:   !   PSUM(I) = X ! store partial sum199:   !   PSUM(I) = X ! store partial sum
200:   !ENDDO200:   !ENDDO
201:   !PSUM(N) = DBLE(N-1)201:   !PSUM(N) = DBLE(N-1)
202:   !202:   !
203:   ELOWEST=MINVAL(VALUES)203:   ELOWEST=MINVAL(VALUES)
204:   X=0.0D0 ! initialise total sum204:   X=0.0D0 ! initialise total sum
205:   DO I=1,N-1205:   DO I=1,N-1
206:      IF(VALUES(I) > VALUES(N) .AND. VALUES(N) > ELOWEST) THEN206:      IF(VALUES(I) > VALUES(N) .AND. VALUES(N) < ELOWEST) THEN
207:         Y=0.0D0207:         Y=0.0D0
208:      ELSE208:      ELSE
209:         Y=DEXP(-(VALUES(I)-ELOWEST)/TEMP(1))209:         Y=DEXP(-(VALUES(I)-ELOWEST)/TEMP(1))
210:      ENDIF210:      ENDIF
211:      X = X + Y211:      X = X + Y
212:      PSUM(I) = X ! store partial sum212:      PSUM(I) = X ! store partial sum
213:   ENDDO213:   ENDDO
214:   IF(VALUES(N) > ELOWEST) THEN214:   IF(VALUES(N) < ELOWEST) THEN
215:      PSUM(N) = PSUM(N-1) ! Disallow rejection215:      PSUM(N) = PSUM(N-1) ! Disallow rejection
216:   ELSE216:   ELSE
217:      PSUM(N) = DBLE(N-1) ! Allow rejection217:      PSUM(N) = DBLE(N-1) ! Allow rejection
218:   ENDIF218:   ENDIF
219:   !219:   !
220:   X=PSUM(N)*DPRAND()220:   X=PSUM(N)*DPRAND()
221:   !221:   !
222:   I=1222:   I=1
223:   DO WHILE (X > PSUM(I))223:   DO WHILE (X > PSUM(I))
224:      I=I+1224:      I=I+1


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0