indep.c File Reference

#include "mincov_int.h"

Include dependency graph for indep.c:

Go to the source code of this file.

Functions
solution_t *	sm_maximal_independent_set (sm_matrix *A, int *weight)

Function Documentation

sm_maximal_independent_set()

solution_t * sm_maximal_independent_set	(	sm_matrix *	A,
		int *	weight )

Definition at line 44 of file indep.c.

{
    register sm_row *best_row, *prow;
    register sm_element *p;
    int least_weight;
    sm_row *save;
    sm_matrix *B;
    solution_t *indep;
 
    indep = solution_alloc();
    B = build_intersection_matrix(A);
 
    while (B->nrows > 0) {
    /*  Find the row which is disjoint from a maximum number of rows */
    best_row = B->first_row;
    for(prow = B->first_row->next_row; prow != 0; prow = prow->next_row) {
        if (prow->length < best_row->length) {
        best_row = prow;
        }
    }
 
    /* Find which element in this row has least weight */
    if (weight == NIL(int)) {
        least_weight = 1;
    } else {
        prow = sm_get_row(A, best_row->row_num);
        least_weight = weight[prow->first_col->col_num];
        for(p = prow->first_col->next_col; p != 0; p = p->next_col) {
        if (weight[p->col_num] < least_weight) {
            least_weight = weight[p->col_num];
        }
        }
    }
    indep->cost += least_weight;
    (void) sm_row_insert(indep->row, best_row->row_num);
 
    /*  Discard the rows which intersect this row */
    save = sm_row_dup(best_row);
    for(p = save->first_col; p != 0; p = p->next_col) {
        sm_delrow(B, p->col_num);
        sm_delcol(B, p->col_num);
    }
    sm_row_free(save);
    }
 
    sm_free(B);
 
/*
    if (! verify_indep_set(A, indep->row)) {
    fail("sm_maximal_independent_set: row set is not independent");
    }
*/
    return indep;
}

Here is the call graph for this function: