MatrixXd MultivariateFNormalSufficient::compute_PTP() const
{
  LOG( "MVN:   computing PTP" << std::endl);
  VectorXd peps(get_Peps());
  MatrixXd tmp(peps*peps.transpose());
  return tmp;
}
  MatrixXd MultivariateFNormalSufficient::compute_PTP() const
{
  timer_.start(PTP);
  IMP_LOG(TERSE, "MVN:   computing PTP" << std::endl);
  VectorXd peps(get_Peps());
  MatrixXd tmp(peps*peps.transpose());
  timer_.stop(PTP);
  return tmp;
}