const ArrayXd inverseLink::linkInv(const ArrayXd& eta) const {return eta.inverse();}
const ArrayXd inverseLink::muEta( const ArrayXd& eta) const {return -(eta.inverse().square());}
//@{ const ArrayXd inverseLink::linkFun(const ArrayXd& mu) const {return mu.inverse();}