MCPcopy Create free account
hub / github.com/BirolLab/abyss / findRepeats

Function findRepeats

SimpleGraph/SimpleGraph.cpp:237–254  ·  view source on GitHub ↗

Return the set of contigs that appear more than once in a single * solution. */

Source from the content-addressed store, hash-verified

235 * solution.
236 */
237static set<ContigID> findRepeats(ContigID seed,
238 const ContigPaths& solutions)
239{
240 set<ContigID> repeats;
241 for (ContigPaths::const_iterator solIt = solutions.begin();
242 solIt != solutions.end(); ++solIt) {
243 map<ContigID, unsigned> count;
244 count[seed]++;
245 for (ContigPath::const_iterator it = solIt->begin();
246 it != solIt->end(); ++it)
247 count[it->contigIndex()]++;
248 for (map<ContigID, unsigned>::const_iterator
249 it = count.begin(); it != count.end(); ++it)
250 if (it->second > 1)
251 repeats.insert(it->first);
252 }
253 return repeats;
254}
255
256/** The fewest number of pairs in a distance estimate. */
257static unsigned g_minNumPairs = UINT_MAX;

Callers 1

handleEstimateFunction · 0.70

Calls 4

contigIndexMethod · 0.80
beginMethod · 0.45
endMethod · 0.45
insertMethod · 0.45

Tested by

no test coverage detected