Add a routine to SampleCollection that computes estimate of mode.

[Matthew Parno]

No description provided.

[Matthew Parno]

It’s straightforward to just grab the LogTarget metadata and find the index with the largest value. Given how easy this is, I don’t think we need to include anything in the SampleCollection class to enable this.

‌

More specifically,

‌

Eigen::MatrixXd logTargetDens = samps->GetMeta("LogTarget");

double maxLogDens;
unsigned int maxRow, maxCol;

maxLogDens = logTargetDens.maxCoeff(&maxRow, &maxCol);

‌

[Matthew Parno]

Not necessary given ease of extracting log density information for each sample.