Ah yes, this looks like it - in my example, if I change x_grid = np.linspace(1, 10, 1024) (i.e. so that all values are above the lower value) then the values are finite as expected. OK I think this is a bug, so will make an issue - thanks for helping debug this!
1 Like