cluster_walktrap.hpp

Go to the documentation of this file.

#ifndef SCRAN_CLUSTER_WALKTRAP_HPP
#define SCRAN_CLUSTER_WALKTRAP_HPP
 
#include <vector>
#include <algorithm>
 
#include "raiigraph/raiigraph.hpp"
#include "igraph.h"
 
namespace scran_graph_cluster {
 
struct ClusterWalktrapOptions {
    int steps = 4;
 
    bool report_merges = true;
 
    bool report_modularity = true;
};
 
struct ClusterWalktrapResults {
    int status = 0;
    
    raiigraph::IntegerVector membership;
 
    raiigraph::IntegerMatrix merges;
 
    raiigraph::RealVector modularity;
};
 
inline void cluster_walktrap(const igraph_t* graph, const igraph_vector_t* weights, const ClusterWalktrapOptions& options, ClusterWalktrapResults& output) {
    auto membership = output.membership.get();
    auto modularity = (options.report_modularity ? output.modularity.get() : NULL);
    auto merges = (options.report_merges ? output.merges.get() : NULL);
    output.status = igraph_community_walktrap(graph, weights, options.steps, merges, modularity, membership);
}
 
inline ClusterWalktrapResults cluster_walktrap(const raiigraph::Graph& graph, const std::vector<igraph_real_t>& weights, const ClusterWalktrapOptions& options) {
    // No need to free this, as it's just a view.
    igraph_vector_t weight_view;
    igraph_vector_view(&weight_view, weights.data(), weights.size());
 
    ClusterWalktrapResults output;
    cluster_walktrap(graph.get(), &weight_view, options, output);
    return output;
}
 
}
 
#endif