/**
* The total exploring function
* @param startingNode
* @return
*/
int* exploreArea(int startingNode) {
reInit();
printf("Calculating odd nodes...\n");
getOddNodes();
int iterations = 400;
printf("First expansion...\n");
initGen();
while (iterations--) firstExpansion();
printf("Improve genetically...\n");
int Z;
for (Z = 0; Z < 5; Z++) improveGenetically();
printf("Create odd paths...\n");
createOddPaths();
printf("Create even paths...\n");
createEvenPaths();
printf("Merge even paths...\n");
mergeEvenPaths();
printf("Exploring...\n");
explore(startingNode);
if (pathReturned == NULL) return NULL;
int len, i, index = 0; for(len = 0; pathReturned[len] != -1; len++);
int* total = malloc(sizeof(int)*len*2);
for (i = 0; i < len-1; i++) {
total[index++] = pathReturned[i]; total[index++] = pathReturned[i+1];
}
total[index++] = -2; total[index++] = -2;
return total;
}
static stList *splitMultipleStubCycle(stList *cycle,
stList *nonZeroWeightAdjacencyEdges, stSortedSet *allAdjacencyEdges,
stList *stubEdges, stList *chainEdges) {
/*
* Takes a simple cycle containing k stub edges and splits into k cycles, each containing 1 stub edge.
*/
/*
* Get sub-components containing only adjacency and chain edges.
*/
stSortedSet *stubAndChainEdgesSet = getSetOfMergedLists(stubEdges,
chainEdges);
stList *adjacencyEdgeMatching =
stList_filterToExclude(cycle, stubAndChainEdgesSet); //Filter out the the non-adjacency edges
//Make it only the chain edges present in the original component
stList *stubFreePaths = getComponents2(adjacencyEdgeMatching, NULL,
chainEdges);
stList_destruct(adjacencyEdgeMatching);
assert(stList_length(stubFreePaths) >= 1);
stList *splitCycles = stList_construct3(0,
(void(*)(void *)) stList_destruct); //The list to return.
if (stList_length(stubFreePaths) > 1) {
/*
* Build the list of adjacency edges acceptable in the merge
*/
stSortedSet *oddNodes = getOddNodes(cycle);
stList *oddToEvenNonZeroWeightAdjacencyEdges =
getOddToEvenAdjacencyEdges(oddNodes,
nonZeroWeightAdjacencyEdges);
stSortedSet *oddToEvenAllAdjacencyEdges = getOddToEvenAdjacencyEdges2(oddNodes, allAdjacencyEdges);
/*
* Merge together the best two components.
*/
stList *l = filterListsToExclude(stubFreePaths, stubAndChainEdgesSet);
doBestMergeOfTwoSimpleCycles(l, oddToEvenNonZeroWeightAdjacencyEdges,
oddToEvenAllAdjacencyEdges); //This is inplace.
stList *l2 = stList_join(l);
stList_destruct(l);
l = getComponents2(l2, stubEdges, chainEdges);
assert(stList_length(l) == 2);
stList_destruct(l2);
/*
* Cleanup
*/
stSortedSet_destruct(oddNodes);
stList_destruct(oddToEvenNonZeroWeightAdjacencyEdges);
stSortedSet_destruct(oddToEvenAllAdjacencyEdges);
/*
* Call procedure recursively.
*/
for (int64_t i = 0; i < stList_length(l); i++) {
/*
* Split into adjacency edges, stub edges and chain edges.
*/
stList *subCycle = stList_get(l, i);
stList *subAdjacencyEdges;
stList *subStubEdges;
stList *subChainEdges;
splitIntoAdjacenciesStubsAndChains(subCycle,
nonZeroWeightAdjacencyEdges, stubEdges, chainEdges,
&subAdjacencyEdges, &subStubEdges, &subChainEdges);
/*
* Call recursively.
*/
l2 = splitMultipleStubCycle(subCycle, subAdjacencyEdges,
allAdjacencyEdges, subStubEdges, subChainEdges);
stList_appendAll(splitCycles, l2);
/*
* Clean up
*/
stList_setDestructor(l2, NULL);
stList_destruct(l2);
stList_destruct(subAdjacencyEdges);
stList_destruct(subStubEdges);
stList_destruct(subChainEdges);
}
stList_destruct(l);
} else {
stList_append(splitCycles, stList_copy(cycle, NULL));
}
stSortedSet_destruct(stubAndChainEdgesSet);
stList_destruct(stubFreePaths);
return splitCycles;
}
