Construct an adjacency matrix for specifying which rates must be updated when a transition occurs, with adjacency determined at the lumped population level.

Constructs an adjacency matrix with rates given in the rows, and reactions on which they each depend given in columns. Thus, when a reaction occurs, all of the non-zero entries in the corresponding column must be updated.

build_rate_adjmat(rates, compartment_codes)

Arguments

rates

intermediate list of rate functions created within the stem_dynamics function

compartment_codes

vector of compartment codes.

Value

adjacency matrix