int main()
{
Graph *graph = getGraph("input.txt");
if (graph == NULL)
{
printf("Ошибка открытия файла\n");
return 0;
}
int connectedComponents = countConnectedComponents(graph);
printf("Количество компонент связности: %d\n", connectedComponents);
deleteGraph(graph);
}
Locus *extractGraphLoci(Graph * argGraph, ReadSet * reads,
boolean * dubious, ShortLength * lengths,
IDnum * locusCount, boolean scaffolding)
{
Locus *loci;
graph = argGraph;
buildScaffold(graph, reads, dubious, lengths, scaffolding);
puts("Extracting loci from connection graph...");
*locusCount = countConnectedComponents(graph);
velvetLog("Counted %li mRNA loci\n", (long) *locusCount);
loci = extractConnectedComponents(*locusCount);
if (doubleStrandedGraph(graph))
orientLoci(loci, *locusCount);
transitiveReduction();
simplifyLoci(loci, *locusCount);
return loci;
}
