[BART-user] safelog10

[aj-foster]

The following is an archived message from the BART-user mailing list, which has now closed.

From: sa***@kn*** (sally.blumenthal)
Date: Mon, 15 Jun 2015 17:10:44 +0000
Subject: [BART-user] safelog10

In makeatm.radpress(), it uses interp1d with the input of np.log10(press). If press == 0, it will break. I know, I know, zero pressure- what gives? I don't know but using the Venot code outputs there is press==0 so may I suggest this which Joe and I made encountering a similar problem in the CEA project:

import numpy as np

lowerlimit = <eps?>

def safelog10(x):

y = np.log10(x)

y[np.where(~np.isfinite(y))] = np.log10(lowerlimit)

return y

• Sally

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://physics.ucf.edu/pipermail/bart-user/attachments/20150615/32c49bf3/attachment.html