- #1
nguyenthanhlam
- 9
- 0
Dear all,
I can not figure out why I got negative infinity in my output, so please help. Here is my calculation for FAMAX:
FAMAG=SQRT(FX**2+FY**2+FZ**2)
FAMAX=MAXVAL(FAMAG(1:NATOM))
and for FPBPMAX:
FMAGPBP(1:NSYN)=SQRT(FXGTASE(1:NSYN)**2+FYGTASE(1:NSYN)**2+FZGTASE(1:NSYN)**2)
FPBPMAX=MAXVAL(FMAGPBP(1:NSYN))
FMAGPBP(1:NSYN)=SQRT(FXTPASE(1:NSYN)**2+FYTPASE(1:NSYN)**2+FZTPASE(1:NSYN)**2)
FPBPMAX=MAX(MAXVAL(FMAGPBP(1:NSYN)),FPBPMAX)
FMAGPBP(1:NSYN)=SQRT(FXEDASE(1:NSYN)**2+FYEDASE(1:NSYN)**2+FZEDASE(1:NSYN)**2)
FPBPMAX=MAX(MAXVAL(FMAGPBP(1:NSYN)),FPBPMAX)
I expected it crash and got either NaN or Infinity for FAMAX and FPBPMAX at some point. But what I did not expect was FAMAX and FPBPMAX got NEGATIVE infinity when it crashed. Does anyone know why?
Thanks,
Lam
I can not figure out why I got negative infinity in my output, so please help. Here is my calculation for FAMAX:
FAMAG=SQRT(FX**2+FY**2+FZ**2)
FAMAX=MAXVAL(FAMAG(1:NATOM))
and for FPBPMAX:
FMAGPBP(1:NSYN)=SQRT(FXGTASE(1:NSYN)**2+FYGTASE(1:NSYN)**2+FZGTASE(1:NSYN)**2)
FPBPMAX=MAXVAL(FMAGPBP(1:NSYN))
FMAGPBP(1:NSYN)=SQRT(FXTPASE(1:NSYN)**2+FYTPASE(1:NSYN)**2+FZTPASE(1:NSYN)**2)
FPBPMAX=MAX(MAXVAL(FMAGPBP(1:NSYN)),FPBPMAX)
FMAGPBP(1:NSYN)=SQRT(FXEDASE(1:NSYN)**2+FYEDASE(1:NSYN)**2+FZEDASE(1:NSYN)**2)
FPBPMAX=MAX(MAXVAL(FMAGPBP(1:NSYN)),FPBPMAX)
I expected it crash and got either NaN or Infinity for FAMAX and FPBPMAX at some point. But what I did not expect was FAMAX and FPBPMAX got NEGATIVE infinity when it crashed. Does anyone know why?
Thanks,
Lam
Last edited: