denom_lj is used in some derived classes, so we need to initialize both denom_lj and denom_lj_inv
this closes #126 for real