Return the set of contigs that appear more than once in a single * solution. */
| 235 | * solution. |
| 236 | */ |
| 237 | static 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. */ |
| 257 | static unsigned g_minNumPairs = UINT_MAX; |
no test coverage detected