double WishartModel::logp(const SpdMatrix &W) const{ return dWish(W, sumsq(), nu(), true); }
inline double dWishinv(const Spd &S, const Spd &sumsq, double df, bool logscale){ return dWish(S, sumsq, df, logscale, true); }