- #1

- 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: