conjugate gradient code will raise error

[nosyndicate]

Conjugate gradient code will complain divides by 0 error if using numpy 1.13.3 in environmental.yaml on Mac OS on anaconda, so we have to use 1.13.1. But 1.13.1 will raise error on linux. yuck!