In the algebraic approach to mass action chemical reaction networks there is a very useful factorization of the RHS of the ODE system
\frac{dx}{dt}=Y\cdot A_k \cdot \Phi(x),
where Y is the complex stoichiometric matrix, A_k=-L_k and \Phi(x) is the vector of monomials. In order to define L_k consider the complex graph associated to the CRN as a weighted directed graph with reaction rates as weights and then L_k is its Laplacian matrix.
Y is given by the Catalyst function complexstoichmat(). My question is if there are other functions to compute the other two factors, or at least the possibility to construct the weighted graph.