1 #include <Sequence/SummStats/lHaf.hpp> 9 lHaf(
const SimData &data,
const double l)
15 std::vector<unsigned> dcounts;
16 dcounts.reserve(data.numsites());
17 for (
auto i = data.sbegin(); i < data.send(); ++i)
19 dcounts.push_back(static_cast<unsigned>(
20 std::count(i->second.begin(), i->second.end(),
'1')));
23 std::vector<double> rv;
24 rv.reserve(data.size());
28 = std::find_if(i.cbegin(), i.cend(),
29 [](
const char &ch) {
return ch ==
'1'; });
33 size_t d2 = size_t(j - i.cbegin());
34 score += std::pow(static_cast<double>(dcounts[d2]), l);
35 j = std::find(j + 1, i.cend(),
'1');
The namespace in which this library resides.