DPM generator : Bug in DPM fortran code? [message #21538] |
Tue, 26 September 2017 15:25 |
Klaus Götzen
Messages: 293 Registered: June 2006 Location: GSI
|
first-grade participant |
From: *gsi.de
|
|
Hi DPM experts,
when running many jobs with background events on the kronos cluster, it happens quite regulary, that the jobs get stuck, apparently in an infinite loop. I think I identified the problematic part of code in dpm_gen.f:
Nrepeat=0.
560 CONTINUE
Nrepeat=Nrepeat+1
if(Nrepeat.ge.100) go to 550
SumMt=0.
SumaPz=0.
do i=1,Nhad
if(AMF(i).le.0.) AMF(i)=-Amass(NREF(i)) ! Uzhi -Amass
SumMt=SumMt+sqrt(AMF(i)**2+PXF(i)**2+PYF(i)**2)
SumaPz=SumaPz+abs(PZF(i))
enddo
if(SumMt.gt.SS) go to 560
In gdb it looks like it sometimes happens, that this loop to 560 runs infinitely, since SumMt > SS all the time. Could somebody take a look (I'm not a FORTRAN guy) what's going wrong there?
Best regards,
Klaus
|
|
|